Karst ancient landform detail drawing compilation method, apparatus and device, and computer program

Abstract

The invention relates to a karst ancient landform detail drawing compilation method and device, equipment, a storage medium and a computer program. The method comprises the following steps: determining ancient landform categories contained in a region according to a geological background; according to the unconformity development characteristics, the upper stratum thickness and the lower stratum thickness of the karst unconformity surface are determined, and the thickness ratio of the upper stratum to the lower stratum at all drilling positions in the area and the ancient landform type at each drilling position are determined; determining a thickness ratio interval corresponding to each ancient landform type according to a thickness ratio between the upper stratum and the lower stratum at all drilling positions; performing tracking interpretation on the karst unconformity surface according to the seismic data so as to determine a weathering crust of a seismic interpretation layer and a ratio distribution diagram between the upper stratum thickness and the lower stratum thickness of the weathering crust; and determining a karst ancient landform detail drawing according to the ratio distribution diagram of the calibrated drilling position and the thickness ratio interval. The ancient landform is finely depicted in a well-seismic combination mode, operation is convenient, and efficiency is high.

Description

Technical Field

[0001] This disclosure relates to the field of geological exploration technology, and in particular to a method, apparatus, equipment, storage medium, and computer program for compiling detailed maps of karst paleomorphology. Background Technology

[0002] Paleomorphology refers to the geomorphic features of a stratum in a sedimentary area before deposition. Scholars' research on paleomorphology within sedimentary basins is still in the qualitative to semi-quantitative stage, mainly focusing on morphological reconstruction and unit division of paleomorphology during or before deposition.

[0003] Under current technology, paleogeographic restoration methods mainly include: impression method, residual thickness method, sequence stratigraphy method, two-interface method, and comprehensive geological analysis method.

[0004] For example, Zhang et al. (2020) used a combination of well and seismic methods to compile residual thickness maps of the Maokou Formation and thickness maps of the overlying Longtan Formation. Then, based on the residual thickness method and the impression method, they used a comprehensive geological analysis method to characterize the karst paleomorphology of the Maokou Formation in the Sichuan Basin. Zhou et al. (2015), taking the Tahe 2 area as an example, compiled a thickness distribution map of the overlying strata using a combination of well and seismic methods and applied the impression method to restore the early Hercynian karst paleomorphology of the Ordovician strata in the area and its surroundings. Pu et al. (2021), based on the impression method, supplemented by a combination of well and seismic methods, and using 3D modeling, comprehensively applied these three methods to restore the pre-Jurassic paleomorphological features of northern Shaanxi. Feng Lei et al. (2012), addressing the uneven distribution of wells in the western Liaohe River tidal flat region, proposed using seismic data for inter-well paleogeographic trend compensation based on stratigraphic thickness compaction correction of well data. This organically integrated well and seismic data to reconstruct the pre-depositional paleogeography of the Sha-1 section in the western Liaohe River tidal flat region. Liu Yongtao et al. (2019), taking the Tianhuan Depression in the Ordos Basin as an example, proposed a new method for reconstructing the paleogeography of the critical hydrocarbon generation period using double-interface seismic layer flattening technology.

[0005] Currently, with the increasing refinement of exploration requirements and the concealment of exploration targets, paleogeomorphological reconstruction has expanded from macroscopic absolute paleogeomorphology to fine-grained relative paleogeomorphology. Recent studies have found that traditional methods for reconstructing fine-grained relative paleogeomorphology have many limitations. For example, when marker beds are easily penetrated using the impression method, it is difficult to account for compaction when reconstructing pre-depositional paleogeomorphology, and it fails to consider tectonic changes when reconstructing depositional paleogeomorphology, resulting in significant errors. The residual thickness method does not consider the influence of pre-depositional paleotopography and erosion differences, making it less applicable to large-scale reconstructions. Sequence stratigraphy is influenced by complex factors, requires consideration of numerous factors, and involves a large workload for basic mapping. The two-interface method is suitable for smaller study areas, and similar to the geological synthesis method, it requires a comprehensive analysis of the sedimentary and tectonic background of the reconstructed object to determine the paleogeomorphological trend.

[0006] Therefore, in order to overcome the shortcomings of current commonly used methods and explore new methods suitable for fine relative paleogeographic restoration, it is particularly important to be able to perform comprehensive, accurate and rapid relative paleogeographic restoration. Summary of the Invention

[0007] This disclosure provides a method, apparatus, equipment, storage medium, and computer program for compiling detailed maps of karst paleogeography.

[0008] To address the problems of current paleogeographic restoration methods, which require consideration of compaction, tectonic activity, pre-depositional paleotopography, and differential erosion, and suffer from large errors, heavy workload, and unsuitability for large study areas, this disclosure proposes a method for compiling detailed karst paleogeographic maps. This method, based on a ratio-based well-seismic combination approach, identifies karst paleogeographic unit types and compiles paleogeographic overview maps. This method uses thickness ratios as its core and employs a well-seismic combination approach. It determines paleogeographic categories through geological understanding, then extends seismic tracking to the entire area. This well-seismic combination achieves a detailed depiction of relative paleogeography. It is not only convenient and efficient to operate, but also unrestricted by traditional paleogeographic methods such as compaction correction and differential erosion, and the compiled paleogeographic maps are accurate and reliable.

[0009] Firstly, this disclosure provides a method for compiling detailed maps of karst paleogeography, including:

[0010] Based on the geological background of the region where the karst paleolandforms are located, determine the paleolandform categories included in the region;

[0011] Based on the unconformity development characteristics of the area where the karst paleolandform is located, the thickness of the upper and lower strata of the karst unconformity surface is determined.

[0012] Based on the thickness of the upper stratum and the thickness of the lower stratum, determine the thickness ratio between the upper and lower stratum at all drilling locations within the region, and determine the paleogeographic type at each drilling location based on the thickness ratio between the upper and lower stratum at each drilling location and the paleogeographic category.

[0013] Based on the thickness ratio between the upper and lower strata at all drilling locations, the thickness ratio range corresponding to each paleogeographic type is determined.

[0014] Based on the seismic data of the region, the karst unconformity is traced and interpreted to determine the weathering crust of the seismic interpretation layer, and the distribution map of the ratio between the thickness of the upper strata and the thickness of the lower strata of the weathering crust.

[0015] All well locations are marked on the ratio distribution map, and based on the ratio distribution map of the marked well locations and the thickness ratio range, a detailed karst paleogeographic map of the region is determined, thereby completing the compilation of the detailed karst paleogeographic map.

[0016] In some embodiments, it also includes:

[0017] Based on the seismic data and the drilling calibration data of the region, the extent of the gully distribution in the region is traced and interpreted to determine the gully distribution map;

[0018] The karst paleogeographic detailed map is updated based on the trench distribution map to obtain the updated karst paleogeographic detailed map.

[0019] In some embodiments, determining the thickness of the upper and lower strata of the karst unconformity surface based on the unconformity development characteristics of the region where the karst paleogeography is located includes:

[0020] Based on the unconformity development characteristics of the region, the upper and lower marker layers of the karst unconformity surface are determined;

[0021] The thickness between the karst unconformity surface and the upper marker layer is taken as the thickness of the upper stratum, and the thickness between the karst unconformity surface and the lower marker layer is taken as the thickness of the lower stratum.

[0022] In some embodiments, the strata between the karst unconformity surface and the upper marker layer are the upper strata, and the strata between the karst unconformity surface and the lower marker layer are the lower strata; the ratio between the thickness of the upper strata and the thickness of the lower strata is used as the thickness ratio.

[0023] In some embodiments, updating the karst paleogeographic map based on the trench distribution map includes:

[0024] The trench distribution map and the detailed karst paleogeographic map are overlaid to update the detailed karst paleogeographic map.

[0025] In some embodiments, the step of tracking and interpreting the karst unconformity surface based on seismic data of the region to determine the weathering crust of the seismically interpreted layer and a distribution map of the ratio between the thickness of the upper and lower strata of the weathering crust includes:

[0026] Based on the seismic data of the region, the karst unconformity, the upper marker layer of the karst unconformity, and the lower marker layer of the karst unconformity are traced and interpreted to determine the weathering crust, the upper marker layer, and the lower marker layer of the seismic interpretation layer.

[0027] The thickness of the strata between the weathering crust and the upper marker layer of the seismic interpretation plane is taken as the upper stratum thickness of the seismic interpretation plane, and the thickness of the strata between the weathering crust and the lower marker layer of the seismic interpretation plane is taken as the lower stratum thickness of the seismic interpretation plane.

[0028] Based on the thickness of the upper strata and the thickness of the lower strata of the seismic interpretation plane, a distribution map of the ratio between the thickness of the upper strata and the thickness of the lower strata of the weathering crust is determined.

[0029] Secondly, this disclosure provides a device for compiling detailed karst paleogeographic maps, including:

[0030] The paleomorphology category determination module is used to determine the paleomorphology category of the region based on the geological background of the area where the karst paleomorphology is located.

[0031] The stratigraphic thickness determination module is used to determine the upper and lower marker layers of the karst unconformity surface based on the unconformity development characteristics of the area where the karst paleogeography is located, and to determine the upper stratigraphic thickness between the karst unconformity surface and the upper marker layer, as well as the lower stratigraphic thickness between the karst unconformity surface and the lower marker layer.

[0032] The paleogeographic type determination module is used to determine the thickness ratio between the upper and lower strata at all drilling locations within the region based on the thickness of the upper strata and the thickness of the lower strata, and to determine the paleogeographic type at each drilling location based on the thickness ratio between the upper and lower strata at each drilling location and the paleogeographic category.

[0033] The thickness ratio range determination module is used to determine the thickness ratio range corresponding to each paleogeographic type based on the thickness ratio between the upper and lower strata at all drilling locations.

[0034] The ratio distribution map determination module is used to trace and interpret the karst unconformity, the upper marker layer, and the lower marker layer based on the seismic data of the region, so as to determine the weathering crust, the second upper marker layer, and the second lower marker layer of the seismic interpretation plane, and to determine the thickness of the upper stratum between the weathering crust and the second upper marker layer and the thickness of the lower stratum between the weathering crust and the second lower marker layer, and to determine the ratio distribution map between the thickness of the upper stratum and the thickness of the lower stratum based on the thickness of the upper stratum and the thickness of the lower stratum;

[0035] The paleogeographic mapping module is used to mark all well locations on the ratio distribution map and determine the karst paleogeographic map of the region based on the ratio distribution map after the well locations have been marked and the thickness ratio range.

[0036] Thirdly, this disclosure provides a computer device including a memory, a processor, and a computer program stored in the memory, wherein the processor executes the computer program to implement the steps of the method described in the foregoing aspects.

[0037] Fourthly, this disclosure provides a computer-readable storage medium having a computer program stored thereon that, when executed by a processor, implements the steps of the methods described in the above aspects.

[0038] Fifthly, this disclosure provides a computer program product, including a computer program / instructions that, when executed by a processor, implement the steps of the methods described in the foregoing aspects.

[0039] This application circumvents the problems of traditional karst paleogeographic reconstruction by using the stratigraphic thickness ratio method. It proposes a method for identifying karst paleogeographic unit types and compiling paleogeographic maps based on a combination of well and seismic analysis using the thickness ratio method. This method uses thickness ratios as the core, employing a combined well and seismic approach. It determines paleogeographic categories through geological understanding, and then extends this method to the entire region through seismic tracking. This well-seismic approach achieves a detailed characterization of the relative paleogeography. This method has been applied to the reconstruction of pre-Carboniferous paleogeographic features in two blocks, resulting in detailed and reliable paleogeographic distribution maps. In practical work, it has guided well placement and achieved good results.

[0040] This disclosure provides a method, apparatus, equipment, storage medium, and computer program for compiling detailed karst paleogeographic maps. Addressing the problems of current paleogeographic reconstruction methods, which require consideration of compaction, tectonic activity, pre-depositional paleotopography, and erosion differences, and suffer from significant errors, large workloads, and unsuitability for large study areas, this disclosure proposes a method for compiling detailed karst paleogeographic maps. This method uses a ratio-based well-seismic combination approach to identify karst paleogeographic unit types and compile paleogeographic overview maps. This method uses thickness ratios as its core and employs a well-seismic combination approach. It determines paleogeographic categories through geological understanding, then extends seismic tracking to the entire area. This well-seismic combination achieves a detailed depiction of relative paleogeography. It is not only convenient and efficient to operate, but also unrestricted by traditional paleogeographic methods such as compaction correction and differential erosion, and the compiled paleogeographic maps are accurate and reliable. Attached Figure Description

[0041] The present disclosure will be described in more detail below based on embodiments and with reference to the accompanying drawings:

[0042] Figure 1 A flowchart illustrating a method for compiling detailed karst paleogeographic maps according to an embodiment of this disclosure;

[0043] Figure 2 A schematic diagram of key elements for compiling detailed karst paleogeographic maps based on the ratio method, provided in this embodiment of the disclosure;

[0044] Figure 3 This disclosure provides an embodiment of karst paleogeographic features and typical drilling diagrams for a certain block A.

[0045] Figure 4A distribution map of the K0 ratio of upper and lower strata thicknesses provided in this embodiment of the present disclosure;

[0046] Figure 5 This is a trench layout diagram of a certain block A provided in an embodiment of the present disclosure;

[0047] Figure 6 A detailed karst paleogeographic map of a certain block A, compiled based on the ratio method and well-seismic analysis, is provided as an embodiment of this disclosure.

[0048] Figure 7 This disclosure provides an embodiment of karst paleogeographic features and typical well diagrams for a certain Block B.

[0049] Figure 8 A K0 diagram showing the thickness ratio of upper and lower strata in a certain block B, provided in an embodiment of this disclosure;

[0050] Figure 9 A distribution map of karst paleogullies in a certain block B provided in this embodiment of the disclosure;

[0051] Figure 10 This invention provides a detailed karst paleogeographic map compiled based on the ratio method and well-seismic analysis, which is an embodiment of the present disclosure.

[0052] In the accompanying drawings, the same parts are referred to by the same reference numerals, and the drawings are not drawn to scale. Detailed Implementation

[0053] To enable those skilled in the art to better understand the technical solutions of this disclosure, and to fully understand and implement the process of how this disclosure applies technical means to solve technical problems and achieve corresponding technical effects, the technical solutions in the embodiments of this disclosure will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this disclosure, not all embodiments. The embodiments of this disclosure and the various features within them can be combined with each other without conflict, and the resulting technical solutions are all within the protection scope of this disclosure. All other embodiments obtained by those skilled in the art based on the embodiments of this disclosure without creative effort should fall within the protection scope of this disclosure.

[0054] It should be noted that the terms "first," "second," etc., in the specification, claims, and accompanying drawings of this disclosure are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that the embodiments of this disclosure described herein can be implemented in orders other than those illustrated or described herein. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover non-exclusive inclusion; for example, a process, method, system, product, or apparatus that comprises a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, methods, products, or apparatus.

[0055] It should be noted that the steps shown in the flowchart in the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions, and although a logical order is shown in the flowchart, in some cases the steps shown or described may be executed in a different order than that shown here.

[0056] As the background technology indicates, paleogeography refers to the geomorphic features of a stratum in a sedimentary area before deposition. Scholars' research on paleogeography within sedimentary basins is still in the qualitative to semi-quantitative stage, mainly focusing on morphological reconstruction and unit division of paleogeography during or before deposition.

[0057] Under current technology, paleogeographic restoration methods mainly include: impression method, residual thickness method, sequence stratigraphy method, two-interface method, and comprehensive geological analysis method.

[0058] For example, Zhang et al. (2020) used a combination of well and seismic methods to compile residual thickness maps of the Maokou Formation and thickness maps of the overlying Longtan Formation. Then, based on the residual thickness method and the impression method, they used a comprehensive geological analysis method to characterize the karst paleomorphology of the Maokou Formation in the Sichuan Basin. Zhou et al. (2015), taking the Tahe 2 area as an example, compiled a thickness distribution map of the overlying strata using a combination of well and seismic methods and applied the impression method to restore the early Hercynian karst paleomorphology of the Ordovician strata in the area and its surroundings. Pu et al. (2021), based on the impression method, supplemented by a combination of well and seismic methods, and using 3D modeling, comprehensively applied these three methods to restore the pre-Jurassic paleomorphological features of northern Shaanxi. Feng Lei et al. (2012), addressing the uneven distribution of wells in the western Liaohe River tidal flat region, proposed using seismic data for inter-well paleogeographic trend compensation based on stratigraphic thickness compaction correction of well data. This organically integrated well and seismic data to reconstruct the pre-depositional paleogeography of the Sha-1 section in the western Liaohe River tidal flat region. Liu Yongtao et al. (2019), taking the Tianhuan Depression in the Ordos Basin as an example, proposed a new method for reconstructing the paleogeography of the critical hydrocarbon generation period using double-interface seismic layer flattening technology.

[0059] Currently, with the increasing refinement of exploration requirements and the concealment of exploration targets, paleogeomorphological reconstruction has expanded from macroscopic absolute paleogeomorphology to fine-grained relative paleogeomorphology. Recent studies have found that traditional methods for reconstructing fine-grained relative paleogeomorphology have many limitations. For example, when marker beds are easily penetrated using the impression method, it is difficult to account for compaction when reconstructing pre-depositional paleogeomorphology, and it fails to consider tectonic changes when reconstructing depositional paleogeomorphology, resulting in significant errors. The residual thickness method does not consider the influence of pre-depositional paleotopography and erosion differences, making it less applicable to large-scale reconstructions. Sequence stratigraphy is influenced by complex factors, requires consideration of numerous factors, and involves a large workload for basic mapping. The two-interface method is suitable for smaller study areas, and similar to the geological synthesis method, it requires a comprehensive analysis of the sedimentary and tectonic background of the reconstructed object to determine the paleogeomorphological trend.

[0060] Therefore, to overcome the shortcomings of current commonly used methods and explore new methods suitable for detailed relative paleogeomorphic reconstruction, it is particularly important to be able to perform comprehensive, accurate, and rapid relative paleogeomorphic reconstruction. This application discloses a method, apparatus, equipment, storage medium, and computer program for compiling detailed karst paleogeomorphic maps. Addressing the problems of current paleogeomorphic reconstruction methods needing to consider the effects of compaction, tectonic activity, pre-depositional paleotopography, and erosion differences, and suffering from large errors, heavy workload, and unsuitability for large study areas, this disclosure proposes a method for compiling detailed karst paleogeomorphic maps. This method, based on the ratio method and combining well and seismic analysis, identifies karst paleogeomorphic unit types and compiles paleogeomorphic general maps. This method uses thickness ratios as the core and employs a well-seismic combination approach. It determines paleogeomorphic categories through geological understanding, and then extends seismic tracking to the entire area. This well-seismic combination achieves a detailed depiction of relative paleogeomorphic features. It is not only convenient and efficient to operate, but also not limited by traditional paleogeomorphic methods such as compaction correction and differential erosion, and the compiled paleogeomorphic maps are accurate and reliable.

[0061] Example 1

[0062] This embodiment discloses a method for compiling detailed maps of karst paleogeography. Figure 1 This is a flowchart illustrating a method for compiling detailed karst paleogeographic maps according to an embodiment of this disclosure. Figure 1 As shown, the method for compiling detailed karst paleogeographic maps disclosed in this embodiment includes:

[0063] Step 110: Determine the paleolandform category of the area based on the geological background of the karst paleolandform.

[0064] It is understandable that the geological background of the area where the karst paleolandform is located is based on known, pre-defined data.

[0065] Step 120: Based on the unconformity development characteristics of the area where the karst paleolandform is located, determine the thickness of the upper and lower strata of the karst unconformity surface.

[0066] It is understandable that the unconformity development characteristics of the karst paleolandform area are based on known, pre-defined data.

[0067] In some embodiments, determining the thickness of the upper and lower strata of the karst unconformity surface based on the unconformity development characteristics of the region where the karst paleogeography is located includes:

[0068] Based on the unconformity development characteristics of the region, the upper and lower marker layers of the karst unconformity surface are determined;

[0069] The thickness between the karst unconformity surface and the upper marker layer is taken as the thickness of the upper stratum, and the thickness between the karst unconformity surface and the lower marker layer is taken as the thickness of the lower stratum.

[0070] In some embodiments, the strata between the karst unconformity surface and the upper marker layer are the upper strata, and the strata between the karst unconformity surface and the lower marker layer are the lower strata; the ratio between the thickness of the upper strata and the thickness of the lower strata is used as the thickness ratio.

[0071] Step 130: Based on the thickness of the upper stratum and the thickness of the lower stratum, determine the thickness ratio between the upper stratum and the lower stratum at all drilling locations within the region, and determine the paleogeographic type at each drilling location based on the thickness ratio between the upper stratum and the lower stratum at each drilling location and the paleogeographic type.

[0072] Step 140: Determine the thickness ratio range corresponding to each paleogeographic type based on the thickness ratio between the upper and lower strata at all drilling locations.

[0073] Step 150: Based on the seismic data of the area, trace and interpret the karst unconformity to determine the weathering crust of the seismic interpretation layer, and the distribution map of the ratio between the thickness of the upper stratum and the thickness of the lower stratum of the weathering crust.

[0074] In some embodiments, the step of tracking and interpreting the karst unconformity surface based on seismic data of the region to determine the weathering crust of the seismically interpreted layer and a distribution map of the ratio between the thickness of the upper and lower strata of the weathering crust includes:

[0075] Based on the seismic data of the region, the karst unconformity, the upper marker layer of the karst unconformity, and the lower marker layer of the karst unconformity are traced and interpreted to determine the weathering crust, the upper marker layer, and the lower marker layer of the seismic interpretation layer.

[0076] The thickness of the strata between the weathering crust and the upper marker layer of the seismic interpretation plane is taken as the upper stratum thickness of the seismic interpretation plane, and the thickness of the strata between the weathering crust and the lower marker layer of the seismic interpretation plane is taken as the lower stratum thickness of the seismic interpretation plane.

[0077] Based on the thickness of the upper strata and the thickness of the lower strata of the seismic interpretation plane, a distribution map of the ratio between the thickness of the upper strata and the thickness of the lower strata of the weathering crust is determined.

[0078] Step 160: Mark all drilling locations on the ratio distribution map, and determine the detailed karst paleogeographic map of the area based on the ratio distribution map of the marked drilling locations and the thickness ratio range, thereby completing the compilation of the detailed karst paleogeographic map.

[0079] In some embodiments, it also includes:

[0080] Step 170: Based on the seismic data and the drilling calibration data of the region, trace and interpret the range of the valley distribution in the region to determine the trench distribution map;

[0081] Step 180: Update the detailed karst paleogeographic map based on the trench distribution map to obtain the updated detailed karst paleogeographic map.

[0082] It is understandable that the drilling calibration data of the karst paleogeographic area is known, pre-defined data.

[0083] In some embodiments, updating the karst paleogeographic map based on the trench distribution map includes:

[0084] The trench distribution map and the detailed karst paleogeographic map are overlaid to update the detailed karst paleogeographic map.

[0085] Example 2

[0086] Based on the above embodiments, this embodiment further explains and illustrates the method for compiling detailed karst paleogeographic maps disclosed in this application.

[0087] Step 110: Determine the paleolandform category of the area based on the geological background of the karst paleolandform.

[0088] It is understood that the paleolandform categories contained in the area include multiple specific paleolandform types.

[0089] Optionally, based on the regional geological background, the types of paleogeomorphic units can be determined, and the stratigraphic characteristics of different paleogeomorphic unit types can be clarified. For example, a certain block of karst paleogeomorphic units may include residual hills, terraces, plateaus, and gullies. The thickness characteristics of the upper and lower strata for different paleogeomorphic units can be clarified. Residual hills are characterized by thinner upper layers and thicker lower layers; terraces are characterized by thinner upper and middle layers or thinner middle layers and thicker lower layers; and plateaus are characterized by thicker upper layers, thinner lower layers, thicker lower layers, and thicker lower layers. Key elements of detailed karst paleogeomorphic maps can be referenced. Figure 2 .

[0090] It should be noted that the characteristics of different paleolandform types are related to their definitions. For example, residual hills are defined as having a thin upper layer and a thicker lower layer; similarly, other paleolandform types are the same. The paleolandform types determined here are based on the geological background of different regions and research needs. Therefore, the main purpose of this action is classification, not to divide them into several categories. Of course, after classification, we will know how many categories there are, but this knowledge is only known after the classification is completed, not initially with the intention of dividing them into several categories.

[0091] Step 120: Based on the unconformity development characteristics of the area where the karst paleolandform is located, determine the thickness of the upper and lower strata of the karst unconformity surface.

[0092] Optionally, based on the characteristics of unconformity development in the region, the upper marker layer C1 and the lower marker layer C3 of the karst unconformity surface are determined respectively; then, C1 and C3 are matched with the unconformity surface C2 to obtain the upper stratum thickness X and the lower stratum thickness Y respectively, where X = C2 - C1 and Y = C3 - C2.

[0093] In some embodiments, determining the thickness of the upper and lower strata of the karst unconformity surface based on the unconformity development characteristics of the region where the karst paleogeography is located includes:

[0094] Based on the unconformity development characteristics of the region, the upper and lower marker layers of the karst unconformity surface are determined;

[0095] The thickness between the karst unconformity surface and the upper marker layer is taken as the thickness of the upper stratum, and the thickness between the karst unconformity surface and the lower marker layer is taken as the thickness of the lower stratum.

[0096] In some embodiments, the strata between the karst unconformity surface and the upper marker layer are the upper strata, and the strata between the karst unconformity surface and the lower marker layer are the lower strata; the ratio between the thickness of the upper strata and the thickness of the lower strata is used as the thickness ratio.

[0097] Step 130: Based on the thickness of the upper stratum and the thickness of the lower stratum, determine the thickness ratio between the upper stratum and the lower stratum at all drilling locations within the region, and determine the paleogeographic type at each drilling location based on the thickness ratio between the upper stratum and the lower stratum at each drilling location and the paleogeographic type.

[0098] Optionally, based on all wells in the area, determine the formation thicknesses K1, K2, K3, ~Kn above and below all wells, where Kn = Xn / Yn.

[0099] Furthermore, based on the thickness ratio between the upper and lower strata and the paleogeographic category at each drilling location, the specific paleogeographic type at each drilling location is determined.

[0100] Step 140: Determine the thickness ratio range corresponding to each paleogeographic type based on the thickness ratio between the upper and lower strata at all drilling locations.

[0101] Optionally, based on the thickness ratio between the upper and lower strata at the drilling location determined in step 130 and the paleogeographic type at the drilling location, the distribution range of the thickness ratio between the upper and lower strata for different paleogeographic types is statistically determined. For example, the ratio range for residual hills in a certain block is Kmin-A1, the ratio range for terraces is A1-A2, the ratio range for step plateaus is A2-A3, and the ratio range for trenches is A3-Kmax, where the values ​​of A1, A2, and A3 are between Kmin and Kmax.

[0102] Step 150: Based on the seismic data of the area, trace and interpret the karst unconformity to determine the weathering crust of the seismic interpretation layer, and the distribution map of the ratio between the thickness of the upper stratum and the thickness of the lower stratum of the weathering crust.

[0103] In some embodiments, the step of tracking and interpreting the karst unconformity surface based on seismic data of the region to determine the weathering crust of the seismically interpreted layer and a distribution map of the ratio between the thickness of the upper and lower strata of the weathering crust includes:

[0104] Based on the seismic data of the region, the karst unconformity, the upper marker layer of the karst unconformity, and the lower marker layer of the karst unconformity are traced and interpreted to determine the weathering crust, the upper marker layer, and the lower marker layer of the seismic interpretation layer.

[0105] The thickness of the strata between the weathering crust and the upper marker layer of the seismic interpretation plane is taken as the upper stratum thickness of the seismic interpretation plane, and the thickness of the strata between the weathering crust and the lower marker layer of the seismic interpretation plane is taken as the lower stratum thickness of the seismic interpretation plane.

[0106] Based on the thickness of the upper strata and the thickness of the lower strata of the seismic interpretation plane, a distribution map of the ratio between the thickness of the upper strata and the thickness of the lower strata of the weathering crust is determined.

[0107] Optionally, by tracing and interpreting the upper marker layer C1, lower marker layer C3, and unconformity surface C2 determined in the aforementioned steps on the seismic data of the karst paleogeographic region, the upper marker layer C10, lower marker layer C30, and weathering crust C20 of the seismic interpretation layer can be obtained.

[0108] Furthermore, based on the upper marker layer C10, lower marker layer C30, and weathering crust C20 of the seismic interpretation layer, the upper stratum thickness map X0, the lower stratum thickness map Y0, and the ratio distribution map K0 of the upper and lower stratum thicknesses of the weathering crust are calculated again. Among them, the upper and lower strata are extracted from the seismic data by the method of tracking interpretation C1 and C3; in a single well, the X and Y values ​​are obtained according to step (1), and in the seismic data, the same layer is used to obtain X0 and Y0. X0 = C20 - C10 (here it should be C20 and C10, where C20 is the weathering crust of the seismic tracking interpretation and C10 is the upper stratum of the seismic tracking interpretation), similarly Y0 = C30 - C20, where C30 is the lower marker layer of the seismic interpretation, and K0 = X0 / Y0.

[0109] Step 160: Mark all drilling locations on the ratio distribution map, and determine the detailed karst paleogeographic map of the area based on the ratio distribution map of the marked drilling locations and the thickness ratio range, thereby completing the compilation of the detailed karst paleogeographic map.

[0110] In some embodiments, steps 170 and 180 are further included; wherein,

[0111] Step 170: Based on the seismic data and the drilling calibration data of the region, trace and interpret the range of the valley distribution in the region to determine the trench distribution map;

[0112] Step 180: Update the detailed karst paleogeographic map based on the trench distribution map to obtain the updated detailed karst paleogeographic map.

[0113] It should be noted that trenches exhibit a significant response to earthquakes and can be interpreted independently even without well calibration data. However, the interpretation is more accurate when well calibration data is available. Of course, in some cases, well calibration data can be disregarded.

[0114] In some embodiments, updating the karst paleogeographic map based on the trench distribution map includes:

[0115] The trench distribution map and the detailed karst paleogeographic map are overlaid to update the detailed karst paleogeographic map.

[0116] Example 3

[0117] Based on the above embodiments, this embodiment provides a specific example.

[0118] In this embodiment, taking the restoration of karst paleogeography in a certain block of a basin as an example, the top of Tai 1 section of Taiyuan Formation and the bottom of Ma 5 5 are selected as the upper and lower marker layers, respectively.

[0119] Based on the regional geological background, the main paleolandform types in this area are identified as residual hills, terraces, plateaus, and valleys. Among them, residual hills are characterized by being thinner at the top and thicker at the bottom; terraces are characterized by being thinner in the upper and middle sections or thinner in the middle and thicker in the lower middle; and plateaus are characterized by being thicker in the upper middle, thinner in the lower middle, and thicker in the lower middle to the bottom. For details, please refer to [reference needed]. Figure 3 .

[0120] Based on data from 41 wells in the area, the distance X between the top of the Tai 1 section of the Taiyuan Formation and the weathered crust, and the distance Y between the bottom of the Ma 5 5 section and the weathered crust were calculated. The distances X1, X2, X3...X41, Y1, Y2, Y3,...Y41 of the 41 wells were obtained, and the K value of each well was calculated.

[0121] Where K1=X1 / Y1, K2=X2 / Y2, K3=X3 / Y3…………, K41=X41 / Y41, thus obtaining the upper and lower formation thickness ratios K1, K2, K3…………K41 for 41 wells.

[0122] Based on the paleogeographic types identified in the second step of the study, the paleogeographic unit types encountered in each well were determined. According to the paleogeographic types, the distribution range of K values ​​for each paleogeographic type was statistically analyzed. The paleogeographic types and their K value distribution ranges for a certain block A are shown in Table 1.

[0123] Table 1. Statistical table of distribution range of upper and lower stratum thickness ratios for different karst paleogeographic types in Block A.

[0124]

[0125]

[0126] Conduct seismic horizon delineation and horizon tracking interpretation work, in accordance with Figure 3 As shown in the paleogeographic reconstruction features of Block A, the three strata—the top of Tai 1, the weathering crust, and the bottom of Ma 5 5—are interpreted separately. The upper stratum thickness map X0 and the lower stratum thickness map Y0 are calculated within the seismic data body, and the ratio of the two thickness maps is calculated to obtain the upper and lower stratum thickness ratio distribution map K0, where K0 = X0 / Y0. Figure 4 As shown.

[0127] On the K0 ratio distribution map above, the contour lines of different ratios exhibit a regular and gradual distribution characteristic. Combining the paleogeographic types already obtained from the well sites at various locations, and following the steps above to obtain the distribution range of different paleogeographic ratios, a detailed paleogeographic map without gullies can be compiled.

[0128] Based on seismic data, combined with well calibration, characteristics such as the distribution range of valleys are traced and interpreted to obtain valley distribution maps, such as... Figure 5 As shown.

[0129] By overlaying the above detailed paleogeographic map without gullies with the gully distribution map, the final detailed paleogeographic map is obtained, as shown below. Figure 6 As shown.

[0130] Example 4

[0131] Based on the above embodiments, this embodiment provides a specific example.

[0132] Taking the restoration of karst paleogeography in Block B of a basin as an example, based on the regional geological background, the bottom boundary of He Yiduan and the bottom of Ma Sanduan were selected as the upper and lower marker layers, respectively.

[0133] Based on geological characteristics, the paleogeographic types in this area are mainly identified as four types: residual hills, terraces, plateaus, and valleys. Among them, residual hills are characterized by being thinner at the top and thicker at the bottom; terraces are characterized by being thinner in the upper and middle sections or thinner in the middle and thicker in the lower middle; and plateaus are characterized by being thicker in the upper middle, thinner in the lower middle, thicker in the lower middle, and thicker at the bottom. (Reference) Figure 7 .

[0134] Based on data from 17 wells in the area, the distances H1 and H2 between the bottom boundary of the first section of the He and the bottom of the third section of the Ma and the weathered crust were calculated, and the K value for each well was calculated.

[0135] Where K1 = X1 / Y1, K2 = X2 / Y2, K3 = X3 / Y3…………, K17 = H17 / Y17, thus obtaining the upper and lower formation thickness ratios K1, K2, K3…………K17 for the 17 wells.

[0136] Based on the geological background of the area, the paleogeomorphic unit type encountered by each well was determined. According to the paleogeomorphic type, the distribution range of K value for each paleogeomorphic type was statistically analyzed. The paleogeomorphic type and its K value distribution range of a certain block B are shown in Table 2.

[0137] Table 2. Statistical table of distribution range of upper and lower stratum thickness ratios for different karst paleogeographic types in Block B.

[0138]

[0139]

[0140] Seismic stratigraphic identification and interpretation were conducted. Based on the paleogeomorphological reconstruction elements of Block B, the stratigraphic levels of the bottom boundary of the He-1 member, the weathering crust, and the bottom of the Ma-3 member were interpreted. The thickness maps X0 (upper strata) and Y0 (lower strata) were calculated. The ratio of these two thickness maps was calculated to obtain the distribution map K0 of the upper and lower strata thickness ratio, where K0 = X0 / Y0. Figure 8 As shown.

[0141] On the K0 map above, the contour lines of different ratios exhibit a regular, gradual distribution. Combining this with the paleogeographic types already obtained from the well sites at various locations, and following the steps described above, the distribution range of different paleogeographic ratios is obtained, and a detailed paleogeographic map without grooves is compiled.

[0142] Based on seismic data, combined with well calibration, characteristics such as the distribution range of valleys are traced and interpreted to obtain valley distribution maps, which are then used as references. Figure 9 .

[0143] Furthermore, by overlaying the above detailed paleogeographic map without gullies with the gully distribution map, a final detailed paleogeographic map is obtained, for reference. Figure 10 .

[0144] Example 5

[0145] Based on the above embodiments, this embodiment provides a device for determining the cleanliness of a wellbore. This device embodiment can be used to execute the method embodiments of this application. For details not disclosed in this device embodiment, please refer to the method embodiments of this application. The device disclosed in this embodiment includes:

[0146] The paleolandform category determination module is used to determine the paleolandform category contained in the region based on the geological background of the area where the karst paleolandform is located;

[0147] The stratigraphic thickness determination module is used to determine the upper and lower marker layers of the karst unconformity surface based on the unconformity development characteristics of the area where the karst paleogeography is located, and to determine the upper stratigraphic thickness between the karst unconformity surface and the upper marker layer, as well as the lower stratigraphic thickness between the karst unconformity surface and the lower marker layer.

[0148] The paleogeographic type determination module is used to determine the thickness ratio between the upper and lower strata at all drilling locations within the region based on the thickness of the upper strata and the thickness of the lower strata, and to determine the paleogeographic type at each drilling location based on the thickness ratio between the upper and lower strata at each drilling location and the paleogeographic category.

[0149] The thickness ratio range determination module is used to determine the thickness ratio range corresponding to each paleogeographic type based on the thickness ratio between the upper and lower strata at all drilling locations.

[0150] The ratio distribution map determination module is used to trace and interpret the karst unconformity, the upper marker layer, and the lower marker layer based on the seismic data of the region, so as to determine the weathering crust, the second upper marker layer, and the second lower marker layer of the seismic interpretation plane, and to determine the thickness of the upper stratum between the weathering crust and the second upper marker layer and the thickness of the lower stratum between the weathering crust and the second lower marker layer, and to determine the ratio distribution map between the thickness of the upper stratum and the thickness of the lower stratum based on the thickness of the upper stratum and the thickness of the lower stratum;

[0151] The paleogeographic mapping module is used to mark all well locations on the ratio distribution map and determine the karst paleogeographic map of the region based on the ratio distribution map after the well locations have been marked and the thickness ratio range.

[0152] In some embodiments, it also includes:

[0153] The trench distribution map determination module is used to trace and interpret the trench distribution range of the region based on the seismic data and the drilling calibration data of the region, so as to determine the trench distribution map;

[0154] The paleogeographic map update module is used to update the karst paleogeographic map based on the trench distribution map to obtain the updated karst paleogeographic map.

[0155] In some embodiments, the paleogeographic map update module is used to overlay the trench distribution map with the karst paleogeographic map to update the karst paleogeographic map.

[0156] In some embodiments, the formation thickness determination module includes a marker layer determination unit and a formation thickness determination unit; wherein...

[0157] The marker layer determination unit is used to determine the upper and lower marker layers of the karst unconformity surface based on the unconformity development characteristics of the region.

[0158] The stratum thickness determination unit is used to take the thickness between the karst unconformity surface and the upper marker layer as the upper stratum thickness, and the thickness between the karst unconformity surface and the lower marker layer as the lower stratum thickness.

[0159] In some embodiments, the strata between the karst unconformity surface and the upper marker layer are the upper strata, and the strata between the karst unconformity surface and the lower marker layer are the lower strata; the ratio between the thickness of the upper strata and the thickness of the lower strata is used as the thickness ratio.

[0160] In some embodiments, the ratio distribution map determination module includes: a marker layer determination unit, a formation thickness determination unit, and a ratio distribution map determination unit; wherein...

[0161] The marker layer determination unit is used to trace and interpret the karst unconformity, the upper marker layer of the karst unconformity, and the lower marker layer of the karst unconformity based on the seismic data of the region, so as to determine the weathering crust of the seismic interpretation layer, the upper marker layer of the seismic interpretation layer, and the lower marker layer of the seismic interpretation layer.

[0162] The stratigraphic thickness determination unit is used to take the stratigraphic thickness between the weathering crust and the upper marker layer of the seismic interpretation plane as the upper stratigraphic thickness of the seismic interpretation plane, and to take the stratigraphic thickness between the weathering crust and the lower marker layer of the seismic interpretation plane as the lower stratigraphic thickness of the seismic interpretation plane.

[0163] The ratio distribution map determination unit is used to determine the ratio distribution map between the upper and lower strata thicknesses of the weathering crust based on the upper strata thickness and the lower strata thickness of the seismic interpretation plane.

[0164] Those skilled in the art will understand that the modules or steps described above can be implemented using general-purpose computing devices, either centralized on a single computing device or distributed across a network of multiple computing devices. Optionally, they can be implemented using computer-executable program code, thereby allowing them to be stored in a storage device for execution by a computing device. Furthermore, in some cases, the steps shown or described can be performed in a different order than presented herein, or they can be fabricated as separate integrated circuit modules, or multiple modules or steps can be fabricated as a single integrated circuit module.

[0165] Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the specific working process of each module in the wellbore cleanliness determination device can be referred to the corresponding process in the aforementioned method embodiment, and will not be repeated here.

[0166] Example 6

[0167] Based on the above embodiments, this embodiment provides a computer-readable storage medium. This computer-readable storage medium stores a computer program, which, when executed by a processor, can implement the method steps as described in the foregoing method embodiments:

[0168] Based on the geological background of the region where the karst paleolandforms are located, determine the paleolandform categories included in the region;

[0169] Based on the unconformity development characteristics of the area where the karst paleolandform is located, the thickness of the upper and lower strata of the karst unconformity surface is determined.

[0170] Based on the thickness of the upper stratum and the thickness of the lower stratum, determine the thickness ratio between the upper and lower stratum at all drilling locations within the region, and determine the paleogeographic type at each drilling location based on the thickness ratio between the upper and lower stratum at each drilling location and the paleogeographic category.

[0171] Based on the thickness ratio between the upper and lower strata at all drilling locations, the thickness ratio range corresponding to each paleogeographic type is determined.

[0172] Based on the seismic data of the region, the karst unconformity is traced and interpreted to determine the weathering crust of the seismic interpretation layer, and the distribution map of the ratio between the thickness of the upper strata and the thickness of the lower strata of the weathering crust.

[0173] All well locations are marked on the ratio distribution map, and based on the ratio distribution map of the marked well locations and the thickness ratio range, a detailed karst paleogeographic map of the region is determined, thereby completing the compilation of the detailed karst paleogeographic map.

[0174] Computer-readable storage media may individually include computer programs, data files, data structures, etc., or combinations thereof. The computer-readable storage media or computer program may be specifically designed and understood by those skilled in the art of computer software, or the computer-readable storage media may be known and available to those skilled in the art of computer software. Examples of computer-readable storage media include: magnetic media, such as hard disks, floppy disks, and magnetic tapes; optical media, such as CD-ROMs and DVDs; magneto-optical media, such as optical discs; and hardware devices specifically configured to store and execute computer programs, such as read-only memory (ROM), random access memory (RAM), flash memory; or servers, application stores, etc. Examples of computer programs include machine code (e.g., code generated by a compiler) and files containing high-level code that can be executed by a computer using an interpreter. The described hardware devices may be configured to function as one or more software modules to perform the operations and methods described above, and vice versa. Furthermore, computer-readable storage media may be distributed across networked computer systems, allowing for the decentralized storage and execution of program code or computer programs.

[0175] Example 7

[0176] Based on the above embodiments, this embodiment provides a computer program product. The computer program product includes a computer program or instructions, which, when executed by a processor, implement all or part of the steps of the methods described in the foregoing method embodiments:

[0177] Based on the geological background of the region where the karst paleolandforms are located, determine the paleolandform categories included in the region;

[0178] Based on the unconformity development characteristics of the area where the karst paleolandform is located, the thickness of the upper and lower strata of the karst unconformity surface is determined.

[0179] Based on the thickness of the upper stratum and the thickness of the lower stratum, determine the thickness ratio between the upper and lower stratum at all drilling locations within the region, and determine the paleogeographic type at each drilling location based on the thickness ratio between the upper and lower stratum at each drilling location and the paleogeographic category.

[0180] Based on the thickness ratio between the upper and lower strata at all drilling locations, the thickness ratio range corresponding to each paleogeographic type is determined.

[0181] Based on the seismic data of the region, the karst unconformity is traced and interpreted to determine the weathering crust of the seismic interpretation layer, and the distribution map of the ratio between the thickness of the upper strata and the thickness of the lower strata of the weathering crust.

[0182] All well locations are marked on the ratio distribution map, and based on the ratio distribution map of the marked well locations and the thickness ratio range, a detailed karst paleogeographic map of the region is determined, thereby completing the compilation of the detailed karst paleogeographic map.

[0183] Furthermore, the computer program product may include one or more computer-executable components configured to perform the embodiments when the program is run; the computer program product may also include a computer program tangibly contained on a readable medium thereof, the computer program containing program code for performing any of the methods in the embodiments of this disclosure. In such embodiments, the computer program may be downloaded and installed from a network via a communication component, and / or installed from a removable medium.

[0184] Example 8

[0185] Based on the above embodiments, this embodiment provides a computer device, including a memory, a processor, and a computer program stored in the memory, wherein the processor executes the computer program to implement the steps of the method described in the above embodiments.

[0186] In some embodiments of this example, a computer-readable storage medium is provided, on which a computer program is stored, characterized in that the computer program, when executed by a processor, implements the steps of the method described in the above embodiments.

[0187] In some embodiments of this example, a computer program product is provided, including a computer program / instructions, characterized in that the computer program, when executed by a processor, implements the steps of the method described in the above embodiments.

[0188] The processor may include, but is not limited to, one or more processors or microprocessors. Each processor may be implemented as an Application Specific Integrated Circuit (ASIC), Digital Signal Processor (DSP), Digital Signal Processing Device (DSPD), Programmable Logic Device (PLD), Field Programmable Gate Array (FPGA), controller, microcontroller, microprocessor, or other electronic component, for executing the methods described in the above embodiments.

[0189] Computer-readable storage media can be implemented by any type of volatile or non-volatile storage device or a combination thereof. Computer-readable storage media may include, but are not limited to, random access memory (RAM), read-only memory (ROM), flash memory, EPROM memory, EEPROM memory, registers, and computer storage media (e.g., hard disks, floppy disks, solid-state drives, removable disks, CD-ROMs, DVD-ROMs, Blu-ray discs, etc.).

[0190] Computer-readable storage media may also store at least one computer-executable program / instruction, such as computer-readable instructions. Computer-readable storage media include, but are not limited to, volatile memory and / or non-volatile memory. Volatile memory may include, for example, random access memory (RAM) and / or cache memory. Computer-readable storage media may include, for example, read-only memory (ROM), hard disk, flash memory, etc. For example, a non-transitory computer-readable storage medium may be connected to a computing device such as a computer, and then, when the computing device executes the computer-readable instructions stored on the computer-readable storage medium, the various methods described above can be performed.

[0191] In addition, the computer device may include (but is not limited to) a data bus, an input / output (I / O) bus, a display, and input / output devices (e.g., keyboard, mouse, speakers, etc.).

[0192] The processor can communicate with external devices via the I / O bus through wired or wireless networks.

[0193] In one embodiment, the at least one computer-executable instruction may also be compiled into or comprise a software product / computer program product, wherein one or more computer-executable instructions are executed by a processor to perform the steps of the various functions and / or methods in the embodiments described herein.

[0194] In the embodiments provided in this disclosure, it should be understood that the disclosed apparatus and methods can also be implemented in other ways. The apparatus embodiments described above are merely illustrative; for example, the flowcharts and block diagrams in the accompanying drawings illustrate the architecture, functionality, and operation of possible implementations of apparatus, methods, and computer program products according to various embodiments of this disclosure. In this regard, each block in a flowchart or block diagram may represent a module, segment, or portion of code containing one or more executable instructions for implementing a specified logical function. It should also be noted that in some alternative implementations, the functions marked in the blocks may occur in a different order than those marked in the drawings. For example, two consecutive blocks may actually be executed substantially in parallel, and they may sometimes be executed in reverse order, depending on the functions involved. It should also be noted that each block in a block diagram and / or flowchart, and combinations of blocks in block diagrams and / or flowcharts, can be implemented using a dedicated hardware-based system that performs the specified function or action, or using a combination of dedicated hardware and computer instructions.

[0195] It should be noted that, in this disclosure, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Without further limitation, an element limited by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes that element.

[0196] While the embodiments disclosed herein are as described above, the foregoing content is merely for the purpose of facilitating understanding of this disclosure and is not intended to limit this disclosure. Any person skilled in the art to which this disclosure pertains may make any modifications and changes in form and detail of the implementation without departing from the spirit and scope of this disclosure; however, the scope of patent protection of this disclosure shall still be determined by the scope defined in the appended claims.

Claims

1. A method for compiling detailed maps of karst paleogeography, characterized in that, include: Based on the geological background of the region where the karst paleolandforms are located, determine the paleolandform categories included in the region; Based on the unconformity development characteristics of the area where the karst paleolandform is located, the thickness of the upper and lower strata of the karst unconformity surface is determined. Based on the thickness of the upper stratum and the thickness of the lower stratum, determine the thickness ratio between the upper and lower stratum at all drilling locations within the region, and determine the paleogeographic type at each drilling location based on the thickness ratio between the upper and lower stratum at each drilling location and the paleogeographic category. Based on the thickness ratio between the upper and lower strata at all drilling locations, the thickness ratio range corresponding to each paleogeographic type is determined. Based on the seismic data of the region, the karst unconformity is traced and interpreted to determine the weathering crust of the seismic interpretation layer, and the distribution map of the ratio between the thickness of the upper strata and the thickness of the lower strata of the weathering crust. All well locations are marked on the ratio distribution map, and based on the ratio distribution map of the marked well locations and the thickness ratio range, a detailed karst paleogeographic map of the region is determined, thereby completing the compilation of the detailed karst paleogeographic map.

2. The method according to claim 1, characterized in that, Also includes: Based on the seismic data and the drilling calibration data of the region, the extent of the gully distribution in the region is traced and interpreted to determine the gully distribution map; The karst paleogeographic detailed map is updated based on the trench distribution map to obtain the updated karst paleogeographic detailed map.

3. The method according to claim 1, characterized in that, The determination of the thickness of the upper and lower strata of the karst unconformity surface based on the unconformity development characteristics of the region where the karst paleogeography is located includes: Based on the unconformity development characteristics of the region, the upper and lower marker layers of the karst unconformity surface are determined; The thickness between the karst unconformity surface and the upper marker layer is taken as the thickness of the upper stratum, and the thickness between the karst unconformity surface and the lower marker layer is taken as the thickness of the lower stratum.

4. The method according to claim 3, characterized in that, The strata between the karst unconformity surface and the upper marker layer are the upper strata, and the strata between the karst unconformity surface and the lower marker layer are the lower strata; the ratio between the thickness of the upper strata and the thickness of the lower strata is used as the thickness ratio.

5. The method according to claim 2, characterized in that, The step of updating the detailed karst paleogeographic map based on the trench distribution map includes: The trench distribution map and the detailed karst paleogeographic map are overlaid to update the detailed karst paleogeographic map.

6. The method according to claim 3, characterized in that, The step of tracking and interpreting the karst unconformity surface based on seismic data of the region to determine the weathering crust of the seismically interpreted layer, and the distribution map of the ratio between the thickness of the upper and lower strata of the weathering crust, includes: Based on the seismic data of the region, the karst unconformity, the upper marker layer of the karst unconformity, and the lower marker layer of the karst unconformity are traced and interpreted to determine the weathering crust, the upper marker layer, and the lower marker layer of the seismic interpretation layer. The thickness of the strata between the weathering crust and the upper marker layer of the seismic interpretation plane is taken as the upper stratum thickness of the seismic interpretation plane, and the thickness of the strata between the weathering crust and the lower marker layer of the seismic interpretation plane is taken as the lower stratum thickness of the seismic interpretation plane. Based on the thickness of the upper strata and the thickness of the lower strata of the seismic interpretation plane, a distribution map of the ratio between the thickness of the upper strata and the thickness of the lower strata of the weathering crust is determined.

7. A device for compiling detailed maps of karst paleogeography, characterized in that, include: The paleomorphology category determination module is used to determine the paleomorphology category of the region based on the geological background of the area where the karst paleomorphology is located. The stratigraphic thickness determination module is used to determine the thickness of the upper and lower stratigraphic layers of the karst unconformity surface based on the unconformity development characteristics of the area where the karst paleogeography is located. The paleogeographic type determination module is used to determine the thickness ratio between the upper and lower strata at all drilling locations within the region based on the thickness of the upper strata and the thickness of the lower strata, and to determine the paleogeographic type at each drilling location based on the thickness ratio between the upper and lower strata at each drilling location and the paleogeographic category. The thickness ratio range determination module is used to determine the thickness ratio range corresponding to each paleogeographic type based on the thickness ratio between the upper and lower strata at all drilling locations. The ratio distribution map determination module is used to trace and interpret the karst unconformity surface based on the seismic data of the region, so as to determine the weathering crust of the seismic interpretation layer and the ratio distribution map between the thickness of the upper stratum and the thickness of the lower stratum of the weathering crust. The paleogeographic map determination module is used to mark all drilling locations on the ratio distribution map, and determine the karst paleogeographic map of the area based on the ratio distribution map of the marked drilling locations and the thickness ratio range, thereby completing the compilation of the karst paleogeographic map.

8. A computer device, comprising a memory, a processor, and a computer program stored in the memory, characterized in that, The processor executes the computer program to implement the steps of the method according to any one of claims 1 to 6.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that, When executed by a processor, the computer program implements the steps of the method according to any one of claims 1 to 6.

10. A computer program product comprising a computer program / instructions, characterized in that, When executed by a processor, the computer program implements the steps of the method according to any one of claims 1 to 6.

Images