A method and device for identifying an intra-platform beach zone of an ancient carbonate platform

Abstract

The application provides a method and device for identifying an intra-platform shoal zone of an ancient carbonate platform, and the corresponding method comprises the following steps: establishing a thickness-valley number-cloudy ash ratio relationship model according to the formation thickness, the valley number and the cloudy ash ratio determined by first logging data in a target work area; determining the cloudy ash ratio corresponding to the seismic data according to the formation thickness and the valley number determined by the seismic data in the target work area based on the thickness-valley number-cloudy ash ratio relationship model, and then generating a cloudy ash ratio planar distribution map of the target work area; determining a cloudy ash ratio variation interface according to the cloudy ash ratio planar distribution map; determining the seismic facies characteristics corresponding to the cloudy ash ratio variation interface according to the seismic data; and identifying the ancient carbonate platform intra-platform shoal zone according to the seismic facies characteristics. The application uses seismic data to quickly identify the shoal zone, improves the identification accuracy and efficiency, and is significantly faster than manual seismic facies interpretation. The application has a significant effect on studying the distribution rule of the intra-platform shoal zone and finding a scale reservoir.

Description

Technical Field

[0001] This application belongs to the field of oil and gas exploration technology, particularly the field of oil and gas exploration data processing technology, and specifically relates to a method and device for identifying hilly and shoal zones within ancient carbonate platforms. Background Technology

[0002] Mound-shoal bodies are complexes of carbonate biomass mounds and grain shoals, and have attracted widespread attention in the industry due to their frequent use as oil and gas reservoirs. Mound-shoal bodies are one of the main target types for carbonate oil and gas exploration and have become a very important exploration area.

[0003] The lithology of ancient carbonate platforms is dominated by dolomite. Seismic reflections of the hillocks and shoals within these platforms exhibit multiple interpretations, making identification challenging. These hillocks and shoals typically occur in clusters and belts, referred to in the industry as hillock-shoal zones, which are closely related to sedimentary environments such as paleogeography, water depth, and hydrodynamics. Therefore, identifying these hillock-shoal development zones can significantly improve the accuracy of identifying individual hillocks and shoals, and also pinpoint favorable areas for the development of large-scale reservoirs, thus laying the foundation for the search for large oil and gas fields. Taking the carbonate rock sector of the Tarim Basin as an example, large oil and gas fields such as the Tahe field controlled by weathering crust and the Shunbei field controlled by faults have been discovered. However, large-scale breakthroughs in facies-controlled areas have yet to be seen. The main factor restricting this exploration progress is the lack of identification of facies-controlled large-scale reservoirs.

[0004] Existing technologies have yielded significant work on the identification and internal characterization of hill-shoal bodies, but research on their development patterns and the prediction of their distribution and scale is relatively limited. Identifying hill-shoal zones is a crucial issue in oil and gas exploration. As an important reservoir type, hill-shoal bodies have been extensively studied in terms of their external morphology and internal structure. For example, previous researchers have effectively highlighted the envelope features of hill-shoal bodies and identified their internal layered sedimentary characteristics using spectral inversion high-resolution processing techniques. However, this method primarily relies on seismic data for the identification of individual hill-shoal bodies, leading to ambiguity and uncertainty in its interpretation. Summary of the Invention

[0005] One object of the present invention is to provide a method that can accurately predict the planar distribution characteristics of hill and shoal development zones and improve the accuracy of hill and shoal prediction in areas with low exploration levels.

[0006] Another object of the present invention is to provide a device for identifying the hilly and shoal zones within ancient carbonate platform areas. A further object of the present invention is to provide an electronic device comprising a memory and a processor, wherein the memory stores a computer program, and the processor, when executing the computer program, implements the steps of the aforementioned method for identifying the hilly and shoal zones within ancient carbonate platform areas. A further object of the present invention is to provide a readable medium storing a computer program thereon, wherein the computer program, when executed by a processor, implements the steps of the aforementioned method for identifying the hilly and shoal zones within ancient carbonate platform areas.

[0007] To address the technical problems in the background section of this application, the present invention provides the following technical solutions:

[0008] In a first aspect, the present invention provides a method for identifying the hilly shoal zone within an ancient carbonate platform, comprising:

[0009] Based on the formation thickness, number of valleys, and cloud-to-gray ratio determined by the first well logging data within the target work area, a model relating thickness, number of valleys, and cloud-to-gray ratio is established.

[0010] Based on the thickness-valley number-cloud-gray ratio relationship model, the cloud-gray ratio corresponding to the seismic data is determined according to the stratum thickness and valley number determined by the seismic data in the target work area, and then a cloud-gray ratio planar distribution map of the target work area is generated.

[0011] The cloud-to-gray ratio variation interface is determined based on the cloud-to-gray ratio planar distribution map.

[0012] Based on the earthquake data, determine the seismic facies characteristics corresponding to the cloud-to-gray ratio change interface;

[0013] The ancient carbonate platform inner hilly and shoal zone was identified based on the seismic facies characteristics.

[0014] In some embodiments of the present invention, determining the seismic facies characteristics corresponding to the cloud-to-gray ratio change interface based on the seismic data includes:

[0015] Based on the seismic data corresponding to the cloud-to-gray ratio variation interface, the seismic facies characteristics of the windward side and the leeward side of the hilly-shoal zone within the ancient carbonate platform are determined.

[0016] In some embodiments of the present invention, the seismic phase characteristics of the windward surface include a change in the seismic reflection phase axis from disordered, thickened reflections to parallel reflections;

[0017] The seismic phase characteristics of the leeward side include the transformation of the seismic reflection phase axis from parallel reflection to complex waves and multiple axes.

[0018] In some embodiments of the present invention, the speed at which the chaotic, thickened reflection transforms into parallel reflection is greater than the speed at which the parallel reflection transforms into complex waves or multi-axis reflection.

[0019] In some embodiments of the present invention, the first logging data is formation element logging data; determining the cloud-to-ash ratio based on the first logging data includes:

[0020] The corresponding captured gamma spectrum is determined based on the well logging data of the formation elements.

[0021] The calcium and magnesium content was determined based on the captured gamma spectrum.

[0022] The dolomite content and limestone content are determined based on the calcium content and the magnesium content;

[0023] The dolomite-lime ratio is determined based on the dolomite content and the limestone content.

[0024] In some embodiments of the present invention, before determining the cloud-to-gray ratio change interface based on the cloud-to-gray ratio planar distribution map, the method further includes:

[0025] The cloud-to-gray ratio plane distribution map is corrected based on the second logging data and well logging data within the target work area.

[0026] Secondly, the present invention provides an identification device for the hilly shoal zone within an ancient carbonate platform, the device comprising:

[0027] The three-parameter model building module is used to establish a thickness-valley number-cloud-ash ratio relationship model based on the formation thickness, valley number, and cloud-ash ratio determined by the first well logging data in the target work area.

[0028] The planar distribution map generation module is used to determine the cloud-gray ratio corresponding to the seismic data based on the thickness-valley number-cloud-gray ratio relationship model, according to the stratum thickness and valley number determined by the seismic data in the target work area, and then generate a planar distribution map of the cloud-gray ratio of the target work area.

[0029] The change interface determination module is used to determine the cloud-to-gray ratio change interface based on the cloud-to-gray ratio planar distribution map.

[0030] The seismic facies feature determination module is used to determine the seismic facies features corresponding to the cloud-to-gray ratio change interface based on the seismic data.

[0031] The hill and shoal zone identification module is used to identify the hill and shoal zone within the ancient carbonate platform based on the seismic facies characteristics.

[0032] In some embodiments of the present invention, the seismic facies feature determination module includes:

[0033] The seismic facies feature determination unit is used to determine the seismic facies features of the windward side and the leeward side of the hilly and shoal zone within the ancient carbonate platform based on the seismic data corresponding to the cloud-to-gray ratio change interface.

[0034] In some embodiments of the present invention, the seismic phase characteristics of the windward surface include a change in the seismic reflection phase axis from disordered, thickened reflections to parallel reflections;

[0035] The seismic phase characteristics of the leeward side include the transformation of the seismic reflection phase axis from parallel reflection to complex waves and multiple axes.

[0036] In some embodiments of the present invention, the speed at which the chaotic, thickened reflection transforms into parallel reflection is greater than the speed at which the parallel reflection transforms into complex waves or multi-axis reflection.

[0037] In some embodiments of the present invention, the first logging data is formation element logging data; the three-parameter model establishment module includes:

[0038] The energy spectrum determination unit is used to determine the corresponding captured gamma energy spectrum based on the formation element logging data.

[0039] A calcium and magnesium content determination unit is used to determine the calcium and magnesium content based on the captured gamma spectrum.

[0040] A rock content determination unit is used to determine the dolomite content and limestone content based on the calcium content and the magnesium content;

[0041] The cloud-to-ash ratio determination unit is used to determine the cloud-to-ash ratio based on the dolomite content and the limestone content.

[0042] In some embodiments of the present invention, an identification device for an inner hilly shoal zone of an ancient carbonate platform further includes:

[0043] The planar distribution map correction module is used to correct the cloud-to-gray ratio planar distribution map based on the second logging data and well logging data within the target work area.

[0044] Thirdly, the present invention provides a computer program product, including a computer program / instruction, which, when executed by a processor, implements the steps of a method for identifying an inner hilly shoal zone of an ancient carbonate platform.

[0045] Fourthly, the present invention provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the steps of a method for identifying the hilly and shoal belts within an ancient carbonate platform.

[0046] Fifthly, the present invention provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of a method for identifying the hilly and shoal zones within an ancient carbonate platform.

[0047] As described above, embodiments of the present invention provide a method and apparatus for identifying the hilly and shoal zones within an ancient carbonate platform. The method for identifying the hilly and shoal zones within an ancient carbonate platform includes: First, establishing a thickness-trough number-cloud-ash ratio relationship model based on the stratigraphic thickness, number of valleys, and cloud-ash ratio determined by the first well logging data within the target work area; Next, based on the thickness-trough number-cloud-ash ratio relationship model, determining the cloud-ash ratio corresponding to the seismic data based on the stratigraphic thickness and number of valleys determined by the seismic data within the target work area, and then generating a cloud-ash ratio planar distribution map of the target work area; Determining the cloud-ash ratio variation interface based on the cloud-ash ratio planar distribution map; Determining the seismic facies characteristics corresponding to the cloud-ash ratio variation interface based on the seismic data; Finally, identifying the hilly and shoal zones within an ancient carbonate platform based on the seismic facies characteristics.

[0048] The corresponding identification device for the hill-shoal zone within the ancient carbonate platform includes: a three-parameter model establishment module, used to establish a thickness-valley number-cloud-ash ratio relationship model based on the stratigraphic thickness, valley number, and cloud-ash ratio determined by the first well logging data within the target area; a planar distribution map generation module, used to determine the cloud-ash ratio corresponding to the seismic data based on the thickness-valley number-cloud-ash ratio relationship model, the stratigraphic thickness and valley number determined by the seismic data within the target area, and then generate a planar distribution map of the cloud-ash ratio for the target area; a change interface determination module, used to determine the cloud-ash ratio change interface based on the cloud-ash ratio planar distribution map; a seismic facies characteristic determination module, used to determine the seismic facies characteristics corresponding to the cloud-ash ratio change interface based on the seismic data; and a hill-shoal zone identification module, used to identify the hill-shoal zone within the ancient carbonate platform based on the seismic facies characteristics.

[0049] In summary, this invention improves the accuracy of hill and shoal prediction in areas with low exploration levels by statistically analyzing the cloud-to-gray ratio, formation thickness, and number of extracted valleys of the target layer, quantitatively establishing the relationship between the three, and predicting the planar distribution characteristics of hill and shoal development zones based on this relationship. Attached Figure Description

[0050] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0051] Figure 1This is a flowchart illustrating a method for identifying an inner hilly shoal zone of an ancient carbonate platform according to an embodiment of the present invention.

[0052] Figure 2 This is a flowchart illustrating step 100 of a method for identifying an inner hilly shoal zone of an ancient carbonate platform according to an embodiment of the present invention.

[0053] Figure 3 This is a flowchart illustrating a method for identifying an inner hilly shoal zone of an ancient carbonate platform according to a specific embodiment of the present invention.

[0054] Figure 4 This is a schematic diagram illustrating the cloud-to-ash ratio, thickness, and number of troughs of a typical drilling formation in the study area, as described in a specific embodiment of the present invention.

[0055] Figure 5 This is a schematic diagram illustrating the relationship between cloud dust ratio and thickness / valley number in a typical well in the Qiulitag Formation of the Tarim Basin, according to a specific embodiment of the present invention.

[0056] Figure 6 Thickness map of the Upper Cambrian Qiulitag Formation in the northern Tarim region, as shown in a specific embodiment of the present invention;

[0057] Figure 7 This is a wave valley number attribute map of the Upper Cambrian Qiulitag Formation in the northern Tarim region, as described in a specific embodiment of the present invention.

[0058] Figure 8 This is a cloud-gray ratio distribution map of the Upper Cambrian transition in the northern Tarim region, as shown in a specific embodiment of the present invention.

[0059] Figure 9 This is a cloud-to-gray ratio distribution map of the Upper Cambrian period in the northern Tarim region, as shown in a specific embodiment of the present invention.

[0060] Figure 10 This is a schematic diagram of the seismic facies analysis results of a typical Upper Cambrian section in the northern Tarim region according to a specific embodiment of the present invention;

[0061] Figure 11 This is a distribution map of the Upper Cambrian hilly and shoal zone in the northern Tarim region, as described in a specific embodiment of the present invention.

[0062] Figure 12 A block diagram of an identification device for an inner hilly shoal zone of an ancient carbonate platform, as described in an embodiment of the present invention.

[0063] Figure 13 This is a schematic diagram of the structure of an electronic device in an embodiment of the present invention. Detailed Implementation

[0064] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0065] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0066] It should be noted that the terms "comprising" and "having," and any variations thereof, in the specification, claims, and accompanying drawings of this application are intended to cover a non-exclusive inclusion. For example, a process, method, system, product, or apparatus that includes a series of steps or units is not necessarily limited to those explicitly listed, but may include other steps or units not explicitly listed or inherent to these processes, methods, products, or apparatuses. Without conflict, the embodiments and features in the embodiments of this application can be combined with each other. This application will now be described in detail with reference to the accompanying drawings and embodiments.

[0067] The invention patent with publication number CN112130205A discloses a method for predicting the distribution of carbonate hills and shoals based on seismic data. The method includes: obtaining a geological model and rock velocity and density by performing profile measurements on outcrops of carbonate hills and shoals; obtaining a velocity model based on the geological model and the rock velocity and density; using the velocity model to perform forward modeling simulation of the seismic reflection characteristics of carbonate hills and shoals at different frequencies; determining quantitative frequency prediction parameters based on the seismic reflection characteristics of carbonate hills and shoals at different frequencies and the geological model; establishing a stratigraphic framework for carbonate hills and shoals based on three-dimensional seismic data and drilling data; obtaining a data volume for predicting carbonate hills and shoals based on the quantitative frequency prediction parameters; and identifying carbonate hills and shoals based on the data volume.

[0068] The invention patent with publication number CN116009094A discloses a method for characterizing ultra-deep carbonate rock mounds and shoals based on frequency division properties, including the following steps:

[0069] 1: By loading drilling and logging data into the seismic interpretation area, combined with core and thin section data and field outcrop data into the seismic interpretation area, single-well sequence stratigraphy is carried out to establish a sequence stratigraphic framework;

[0070] 2. After calculating and synthesizing seismic records from drilling logging data using sonic waves and density curves, fine well-seismic calibration is performed using stratigraphic data to obtain time-depth relationships. After establishing a regional seismic framework profile using time-depth relationships, the three-dimensional seismic horizons of the entire region are precisely traced and interpreted.

[0071] 3. Based on drilling and logging data, field outcrop data, core data, and the established sequence stratigraphic framework, determine the scale, petrophysical parameters, and longitudinal and lateral distribution patterns of the hill and shoal bodies in the study area, and establish a geological prediction model for the hill and shoal bodies.

[0072] 4. Based on the established geological prediction model of the hilly area, conduct seismic forward modeling and obtain the seismic forward modeling results;

[0073] 5. The seismic forward modeling results are processed by seismic frequency division to determine the seismic reflection characteristics of hills and shoals of different sizes in different frequency bands. Based on the seismic reflection characteristics, the seismic response model of the hills and shoals in the study area is established. Based on the seismic response model, the optimal identification frequency corresponding to the hills and shoals of different sizes is matched.

[0074] 6. Obtain the dominant frequency and effective bandwidth of the seismic data in the study area, perform seismic frequency division processing on the three-dimensional post-stack seismic data, and obtain multiple single-band data volumes based on the optimal identification frequency obtained in step 5;

[0075] 7. Based on the three-dimensional seismic horizon of the whole region, attribute extraction and algorithm optimization are performed on each single-frequency band data volume at different frequencies to obtain the optimized frequency-division attributes and sequence stratigraphic thickness;

[0076] 8. After matching and statistically correlating the data on hill and shoal bodies encountered in actual drilling and the data related to the development of hill and shoal bodies with the selected frequency division attributes and sequence formation thickness, establish the conversion relationship between frequency division attributes and hill and shoal bodies;

[0077] 9. Remove outliers from the frequency division attributes. Based on the conversion relationship between the frequency division attributes and the hills and beaches established in step 8, integrate the frequency division attributes to characterize the planar distribution of the hills and beaches in the study area.

[0078] Understandably, earthquake prediction methods mainly include earthquake attribute analysis, earthquake facies analysis, and earthquake inversion. These two methods are mainly based on geophysical theories such as the kinematics, dynamics, and geometry of seismic waves, and extract information from earthquake data to make predictions. However, due to the lack of constraints from sedimentary models, the prediction results often lack geological regularity in many regions.

[0079] To address these issues, some geological-seismic methods have emerged in recent years, such as phase-controlled seismic inversion and geological model-based inversion. However, these methods primarily use geological laws as constraints for earthquake prediction to limit or correct the prediction results, without incorporating geological parameters as input conditions. Consequently, they still suffer from poor prediction accuracy and a lack of geological regularity.

[0080] In summary, existing methods mainly rely on seismic data to identify individual hills and shoals, and their interpretation is subject to multiple solutions and uncertainties.

[0081] The acquisition, storage, use, and processing of data in this application comply with relevant laws and regulations.

[0082] Example 1:

[0083] For the reasons stated above, embodiments of the present invention provide a specific implementation method for identifying the hilly shoal zone within an ancient carbonate platform, see [link to relevant documentation]. Figure 1 Specifically, it includes the following:

[0084] Step 100: Based on the formation thickness, number of troughs, and cloud-to-ash ratio determined by the first well logging data within the target work area, establish a model relating thickness, number of troughs, and cloud-to-ash ratio.

[0085] Step 200: Based on the thickness-valley number-cloud-gray ratio relationship model, determine the cloud-gray ratio corresponding to the seismic data according to the stratum thickness and valley number determined by the seismic data in the target work area, and then generate a cloud-gray ratio planar distribution map of the target work area.

[0086] Step 300: Determine the cloud-to-gray ratio variation interface based on the cloud-to-gray ratio planar distribution map;

[0087] Step 400: Determine the seismic facies characteristics corresponding to the cloud-to-gray ratio change interface based on the seismic data;

[0088] Step 500: Identify the hilly and shoal zone within the ancient carbonate platform based on the seismic facies characteristics.

[0089] It is understandable that dolomite is generally more developed in more confined and enclosed carbonate platforms, while limestone is more developed in open environments. Based on the needs of oil and gas exploration research, this invention proposes a method for identifying hill-shoal zones within dolomite platforms. By leveraging the barrier effect of these zones on seawater, which causes differences in dolomite development, a quantitative relationship is established between the dolomite-to-lime ratio, the number of troughs, and stratigraphic thickness. Through the extraction of trough number attributes, combined with stratigraphic thickness, the boundary between dolomite and limestone-bearing dolomite within the platform is identified, thereby determining the development zone of hill-shoal bodies, locating large-scale reservoirs, and providing targets for oil and gas exploration in ancient carbonate platforms. Based on this, embodiments of the present invention provide a method for identifying the hilly and shoal zone within an ancient carbonate platform, comprising: establishing a cloud-to-grain ratio determination model based on the thickness of the target layer, the number of troughs, and the cloud-to-grain ratio of at least one well in the target work area; then, generating a cloud-to-grain ratio planar distribution map of the target work area based on the thickness of the target layer, the number of troughs, and the cloud-to-grain ratio determination model; and finally, identifying the hilly and shoal zone within the ancient carbonate platform of the target work area based on the cloud-to-grain ratio planar distribution map.

[0090] Specifically, this invention first establishes a thickness-valley number-cloud-ash ratio relationship model based on the formation thickness, valley number, and cloud-ash ratio determined by the first well logging data within the target work area; then, based on the thickness-valley number-cloud-ash ratio relationship model, the cloud-ash ratio corresponding to the seismic data is determined according to the formation thickness and valley number determined by the seismic data within the target work area, thereby generating a cloud-ash ratio planar distribution map of the target work area; the cloud-ash ratio variation interface is determined according to the cloud-ash ratio planar distribution map; the seismic facies characteristics corresponding to the cloud-ash ratio variation interface are determined according to the seismic data; finally, the hilly and shoal zones within the ancient carbonate platform are identified according to the seismic facies characteristics.

[0091] Example 2:

[0092] For step 100, the thickness and cloud cover ratio of the target layer can be determined using data from at least one typical well. Seismic data analysis and well-seismic calibration further clarify parameters such as the number of troughs in the target layer. Furthermore, it is understood that variations in the number of troughs are closely related to sedimentary facies changes; therefore, by extracting and analyzing the seismic trough number attributes, the distribution of hill-shoal bodies can be effectively characterized.

[0093] Regarding step 200, since seismic data is continuous on a plane, the thickness of the target layer and the number of valleys on the plane can be obtained from the seismic data, and a cloud-to-gray ratio plane distribution map can be generated from this.

[0094] When implementing step 300, a seismic facies study is conducted near the cloud-to-gray ratio change interface to clarify the seismic reflection boundaries on the windward and leeward sides. The area between the two reflection boundaries is the development range on the hill and beach development zone plane.

[0095] For step 400, near the identified cloud-to-grain ratio transition interface, seismic facies identification using the phase surface method (used to describe the behavior of seismic waves propagating within the Earth, including changes in velocity, direction, and propagation mode) is performed. Well-seismic calibration clarifies the reflection characteristics of strata with different cloud-to-grain ratios during seismic events.

[0096] Example 3:

[0097] In some embodiments of the present invention, determining the seismic facies characteristics corresponding to the cloud-to-gray ratio change interface based on the seismic data includes:

[0098] Based on the seismic data corresponding to the cloud-to-gray ratio variation interface, the seismic facies characteristics of the windward side and the leeward side of the hilly-shoal zone within the ancient carbonate platform are determined.

[0099] Specifically, firstly, the cloud-to-grain ratio (VGR) variation interface is determined on the cloud-to-grain ratio planar distribution map. Then, near the identified VGR transition interface, seismic facies identification using the phase surface method (used to describe the behavior of seismic waves propagating within the Earth, including changes in velocity, direction, and propagation mode) is conducted. Well-seismic calibration clarifies the reflection characteristics of strata with different VGRs in seismic events.

[0100] Interpreting and tracing seismic anomalies in hill-like bodies reveals that they often exhibit low, flat-topped, convex features with chaotic internal reflections. Combining these characteristics with the seismic reflection features on the windward and leeward sides, the distribution range of seismic anomalies in hill-like bodies is characterized, thereby determining the extent of the hill-like zone within the platform.

[0101] In some embodiments of the present invention, the seismic phase characteristics of the windward surface include a change in the seismic reflection phase axis from disordered, thickened reflections to parallel reflections;

[0102] The seismic phase characteristics of the leeward side include the transformation of the seismic reflection phase axis from parallel reflection to complex waves and multiple axes.

[0103] The seismic facies characteristics of the cloud-to-grain ratio variation interface are specifically manifested as low-frequency characteristics of dolomite and high-frequency characteristics of grain-bearing dolomite. Leeward interface characteristics: Seismic reflection phase axes show a gradual transition from parallel reflections to complex, multi-axis characteristics. Windward interface characteristics: Seismic reflection phase axes show a rapid transition from chaotic, thickened, mound-like reflections to parallel reflections.

[0104] Example 4:

[0105] In some embodiments of the present invention, the speed at which the chaotic, thickened reflection transforms into parallel reflection is greater than the speed at which the parallel reflection transforms into complex waves or multi-axis reflection.

[0106] In some embodiments of the present invention, the first logging data is formation element logging data; see [link to relevant documentation]. Figure 2Step 100, determining the cloud-to-ash ratio based on the first well logging data, includes:

[0107] Step 101: Determine the corresponding captured gamma spectrum based on the formation element logging data;

[0108] Formation elemental logging refers to the use of logging techniques to assess the content and distribution of elements in underground rocks. By measuring the formation during drilling using logging instruments, information about the various elements contained in the formation rocks can be obtained.

[0109] Formation element logging techniques mainly include radioactive logging, neutron logging, and gamma-ray logging. These methods can be used to identify parameters such as the type, content, and porosity of formation rocks. By obtaining information on the elemental content of underground rocks through formation element logging, it is possible to help determine the properties, lithology, porosity, and other parameters of the formation, providing important geological data and basis for exploration and development work.

[0110] Step 102: Determine the calcium and magnesium content based on the captured gamma spectrum;

[0111] Step 103: Determine the dolomite content and limestone content based on the calcium and magnesium content;

[0112] Step 104: Determine the dolomite-lime ratio based on the dolomite content and limestone content.

[0113] In steps 101 to 104, formation elemental logging compares the results of stripping and analyzing the inelastic scattering and captured gamma-ray spectra recorded by the instrument with experimental standard spectra to obtain the elemental composition of the formation. It then uses oxide closure models and cluster analysis to determine the types and contents of minerals in the formation, ultimately evaluating the formation. This type of logging method can provide raw data on the contents of formation elements such as silicon, aluminum, calcium, iron, magnesium, and gadolinium. Its comprehensive interpretation results can provide a formation lithological profile and further obtain the dolomite-to-limestone ratio (the ratio of dolomite to limestone).

[0114] In some embodiments of the present invention, prior to step 300, a method for identifying an inner hilly shoal zone of an ancient carbonate platform further includes:

[0115] The cloud-to-gray ratio plane distribution map is corrected based on the second logging data and well logging data within the target work area.

[0116] The cloud-to-grain ratio distribution map is corrected based on multiple actual drilled wells. Specifically, the cloud-to-grain ratio of the target layer is calculated using the logging and well logging data of each actual drilled well. Then, the cloud-to-grain ratio at the drilling location is numerically corrected, so that the cloud-to-grain ratio distribution map is closer to the actual geological conditions, thus clarifying the interface where the plane cloud-to-grain ratio changes.

[0117] Well logging data primarily records information extracted during the drilling process. It includes data obtained from drilling mud (i.e., the fluid circulating during drilling):

[0118] Rock cuttings description: During drilling, rock cuttings (fragments of rock that are broken off from the drill bit) are brought to the surface and analyzed to determine the type of formation that the drilling passed through.

[0119] Gas detection: Gases released from drilling mud are collected and analyzed to identify any hydrocarbon gases that may be present.

[0120] Temperature and pressure data: Record various physical parameters during the drilling process, such as bottom hole temperature and pressure.

[0121] Drilling parameters, including drilling speed, drill pipe rotation speed, and pump speed, are important data for evaluating drilling efficiency and optimizing drilling operations.

[0122] Well logging data is collected using various logging tools, which are typically lowered into the well after drilling to measure physical, chemical, and other properties to assess the formation and surrounding environment. The main types include:

[0123] Resistivity logging: Identifying the type of fluid (oil, water, gas) in a formation by measuring its resistivity.

[0124] Sonic logging: Measuring the sound velocity in a rock formation to infer the rock's density and porosity.

[0125] Density logging: assessing rock porosity and fluid saturation by measuring formation density.

[0126] Nuclear magnetic resonance logging: assesses the pore structure and fluid type of rocks, and is particularly suitable for the analysis of complex reservoirs.

[0127] Imaging logging: Provides detailed images of the formation, helping to identify fractures, formation dip, and rock texture.

[0128] As described above, well logging data and well logging data are closer to the actual underground geological data. Therefore, well logging data and well logging data can be used to correct the cloud-to-gray ratio plane distribution map.

[0129] As described above, this invention provides a method for identifying hill-shoal zones within ancient carbonate platforms. The method includes: First, statistically analyzing parameters such as cloud-to-grain ratio, thickness, and number of troughs in the target strata of typical wells in the study area, and performing correlation analysis to clarify the relationship between these three parameters. Next, predicting the cloud-to-grain ratio distribution based on the stratigraphic thickness and number of troughs in the 3D seismic survey area, and then correcting the cloud-to-grain ratio distribution map based on existing well data. This clarifies the transition interfaces between lithologies with different cloud-to-grain ratios on the plane. Near the identified cloud-to-grain ratio transition interfaces, seismic facies identification using the facies surface method is conducted. Tracking seismic anomalies such as hill-shaped bodies, and combining windward and leeward features, the distribution range is characterized, thereby determining the extent of the hill-shoal zone within the platform. This invention quantitatively establishes the relationship between the cloud-to-grain ratio of the target strata, stratigraphic thickness, and the extracted trough number attribute, thereby predicting the planar distribution characteristics of hill-shoal development zones and improving the accuracy of hill-shoal prediction in areas with low exploration levels.

[0130] Understandably, identifying mound-shoal bodies within large dolomite formations is challenging. The platform margins and intraplatform mound-shoal zones are the primary locations for the development of large-scale carbonate reservoirs. Platform margins are relatively easy to identify due to their distinct geomorphological differences, while intraplatform mound-shoal zones are more difficult to identify.

[0131] This invention comprehensively utilizes drilling and seismic data, establishing an effective method for predicting cloud-to-ash ratio distribution by analyzing the interrelationships between cloud-to-ash ratio elemental logging, formation thickness, and wave trough number attributes. Based on the principle that hill-shoal zones separate platforms and cause differences in cloud formation, seismic facies identification of hill-shoal zones is conducted near the boundary of abrupt changes in cloud-to-ash ratio. Furthermore, the planar distribution of hill-shoal zones is clarified based on the differences in sedimentary structure-seismic reflection on the windward and leeward sides. Particularly in low-exploration areas with few or no wells, seismic data is used to rapidly identify hill-shoal zones, significantly improving accuracy and efficiency compared to artificial seismic facies interpretation. This has significant scientific and practical implications for studying the distribution patterns of hill-shoal zones within platforms and locating large-scale reservoirs. This invention has strong practical applications and can be used in research fields such as sedimentary facies, lithofacies paleogeography, large-scale reservoir prediction, and oil and gas exploration, facilitating its widespread application.

[0132] Example 5:

[0133] To further illustrate the scheme, this invention also provides a specific implementation method for identifying the hilly-shoal zone within an inner platform of a paleocrustate platform, taking the Upper Cambrian Qiulitag Formation in the Tarim Basin as an example. See [link to relevant documentation]. Figure 3 Specifically, it includes the following:

[0134] The Upper Cambrian Qiulitag Formation in the Tarim Basin developed a (semi-)restricted platform sedimentary environment. During the depositional period, the sea level fluctuated at high altitudes, resulting in thick-layered dolomite. The Upper Cambrian Qiulitag Formation is mainly composed of grayish-white to dark gray thick-layered to massive dolomite, often containing chert bands and nodules, interbedded with stromatolites and algal dolomite, and including oolitic, sandy, algal, and bamboo-leaf-shaped clastic dolomite. Based on field outcrops and the latest drilling, logging, and seismic data, the Qiulitag Formation is divided into two third-order sequences.

[0135] S1: Statistical study of the target layer thickness, cloud-to-ash ratio, and number of troughs in the drilling area.

[0136] Specifically, the thickness and cloud cover ratio of the target layer are determined using typical drilling data. Seismic data analysis and well-seismic calibration clarify parameters such as the number of wave troughs in the target layer.

[0137] The thickness and cloud-to-ash ratio of the Qiulitag Formation in ten wells were statistically analyzed (Table 1). Based on the sequence division of single wells, well-seismic calibration was carried out, followed by seismic sequence interpretation and comparison in different facies zones, and the number of wave troughs corresponding to the Qiulitag Formation was counted.

[0138] Table 1. Parameters such as drilling thickness, cloud-to-ash ratio, and trough number in the Qiulitag Formation.

[0139]

[0140] The X1-X4 wells drilled through the Qiulitag Formation, which consists of extremely thick layered dolomite, interbedded with medium to thick layered argillaceous dolomite, and developed in a restricted platform environment. For example, the X2 (Zhonghan 2) well drilled through the Qiulitag Formation, which is composed of extremely thick layered dark gray dolomite, interbedded with similarly colored medium to thick layered argillaceous dolomite. The thickness is 691 meters, the dolomite-to-gray ratio is 22.89, and the number of wave troughs is 3. Seismically, it exhibits subparallel, medium-to-strong amplitude, low-frequency characteristics. Figure 4 ).

[0141] The X5-X10 wells drilled through the Qiulitag Formation, which contains dolomite and argillaceous dolomite, with increasing argillaceous content upwards, indicating a restricted platform environment. For example, the X8 (Yu Shen 1) well penetrated the Qiulitag Formation, which contains argillaceous dolomite with increasing argillaceous content upwards, reaching a thickness of 687.00 meters, a dolomite-to-lime ratio of 2.18, and 7 wave troughs. Seismically, it exhibits strong amplitude, high frequency, and continuous reflection characteristics. Statistical analysis of the thickness, dolomite-to-lime ratio, and wave trough number of these Qiulitag Formation parameters (Table 1) can provide data support for finding correlations.

[0142] S2: Quantitatively establish the correlation between cloud-gray ratio, number of troughs, and thickness.

[0143] Correlation analysis was conducted on the cloud-to-ash ratio, thickness, and trough number of typical wells in the work area to determine the fitting formula. The interrelationship among the three was explored by statistically analyzing the characteristics of cloud-to-ash ratio, thickness, and trough number of typical wells.

[0144] Specifically, the relationship between the cloud-to-ash ratio, thickness, and number of troughs in typical wells was established by statistically analyzing these characteristics: Y (cloud-to-ash ratio) = 0.1556 x (thickness / number of troughs) - 15.072; R 2 =0.9064 ( Figure 5 The cloud-to-grain ratio and the stratum thickness / valley number have a good correlation. Therefore, the cloud-to-grain ratio parameter can be calculated using stratum thickness and valley number obtained from 3D seismic data.

[0145] S3: Extract the valley number and thickness attributes, and convert them into a cloud-to-gray ratio planar distribution map.

[0146] Based on the target layer thickness and trough number in the three-dimensional work area, the cloud-to-gray ratio is calculated using a fitting formula, and the planar distribution of the cloud-to-gray ratio is determined, and a planar distribution map of the cloud-to-gray ratio is drawn.

[0147] First, wellbore calibration was performed on the above wells. The top and bottom interfaces of the Qiulitag Formation were traced using 3D seismic data, and the thickness of the Upper Cambrian Qiulitag Formation in the North Tarim region was analyzed. Figure 6 When the amount of erosion is small, the thickness can reflect the ancient landform to some extent.

[0148] Next, a wave trough number attribute analysis was performed on the three-dimensional region of the Upper Cambrian Qiulitag Formation in the northern Tarim Basin. Figure 7 The variation in the number of seismic troughs is closely related to the variation in sedimentary facies. By extracting and analyzing the seismic wave trough number attributes, the distribution of hill and shoal bodies can be effectively characterized.

[0149] Finally, the analysis results of the wave trough number attributes and thickness of the Upper Cambrian Qiulitag Formation in the North Tarim region were converted into a cloud-gray ratio planar distribution map. Figure 8 ).

[0150] S4: Cloud-to-ash ratio distribution map corrected using cloud-to-ash ratio from actual well drilling.

[0151] The cloud-to-grain ratio distribution map was corrected based on data from multiple drilled wells. Specifically, the cloud-to-grain ratio of the target layer was calculated using logging and well logging data from each drilled well. Then, the cloud-to-grain ratio at the drilling location was numerically corrected using dual-arc software, making the cloud-to-grain ratio distribution map more closely reflect the actual geological conditions. The interface at which the planar cloud-to-grain ratio changes was also clearly defined.

[0152] The obtained cloud-to-ash ratio distribution map was corrected using the cloud-to-ash ratio values ​​measured in the drilled wells. Figure 9 This is done to make the cloud-to-gray ratio distribution map more closely reflect the actual situation.

[0153] S5: Define the interface for varying cloud-to-gray ratio.

[0154] Specifically, seismic facies identification using the phase surface method was conducted near the identified cloud-to-grain ratio transition interface. Well-seismic calibration clarified the seismic reflection characteristics of strata with different cloud-to-grain ratios, specifically revealing low-frequency characteristics in dolomite and high-frequency characteristics in ash-bearing dolomite. Leeward interface characteristics: Seismic reflection phase axes showed a gradual transition from parallel reflections to complex, multi-axis characteristics. Windward interface characteristics: Seismic reflection phase axes showed a rapid transition from chaotic, thickened reflections to parallel reflections.

[0155] Specifically, seismic facies identification is conducted near the cloud-to-gray ratio transition boundary. It is understandable that seismic facies reflect sedimentary facies distribution. It should be noted that this must be considered in conjunction with the geological model of the hill-shoal zone, where multiple phases of hill-shoal bodies develop. The windward slope is steeper, dominated by wave-resistant skeletal hills, and the seismic reflection phase axis shows a rapid transition from chaotic, thickened, hill-like reflections to parallel reflections. The leeward slope is gentler, with developed back-hill shoals, and the seismic reflection phase axis shows a gradual transition from parallel reflections to complex, multi-axis characteristics. Figure 10 The area west of this hilly area is a dolomitic plateau, and the area east of it is a grey dolomitic plateau.

[0156] S6: Conduct seismic facies studies near the cloud-to-gray ratio change interface to clarify the seismic reflection boundaries on the windward and leeward sides.

[0157] Artificial interpretation and tracking of seismic anomalies such as those found in hilly areas are employed. These hilly areas often exhibit low, flat-topped, convex features and chaotic internal reflection characteristics. Combining the aforementioned seismic reflection characteristics of the windward and leeward sides, the distribution range of seismic anomalies in hilly areas is characterized, thereby determining the extent of the hilly area within the platform.

[0158] S7: The area between the boundaries is the extent of the hill and beach development zone on the plane.

[0159] Based on the identified boundaries of the hill and beach zones, draw a distribution map of the hill and beach zones. Figure 11 ).

[0160] As described above, the specific embodiments of the present invention provide a method for identifying the hilly and shoal zones within an ancient carbonate platform, comprising: establishing a thickness-trough number-cloud-ash ratio relationship model based on the stratigraphic thickness, number of troughs, and cloud-ash ratio determined by first well logging data within the target work area; determining the cloud-ash ratio corresponding to the seismic data based on the thickness-trough number-cloud-ash ratio relationship model, according to the stratigraphic thickness and number of troughs determined by seismic data within the target work area, and thereby generating a cloud-ash ratio planar distribution map of the target work area; determining the cloud-ash ratio variation interface based on the cloud-ash ratio planar distribution map; determining the seismic facies characteristics corresponding to the cloud-ash ratio variation interface based on the seismic data; and identifying the hilly and shoal zones within the ancient carbonate platform based on the seismic facies characteristics.

[0161] This invention is based on the principle that the confinement of seawater in the platform's hilly shoal zone promotes the differentiation of dolomite and pyrolithic dolomite within the platform. Generally, open seawater environments are conducive to limestone deposition, while confined and enclosed environments, due to the influence of Mg... 2+ Increased ion concentration favors the development of dolomite. By utilizing the correlation between parameters such as cloud-to-ash ratio and wave trough number, the boundaries of intra-station facies zones are clearly defined, and seismic facies identification is carried out near these boundaries. Generally, the windward side exhibits accretion and a steeper characteristic, while the leeward side exhibits migration and a gentler characteristic of hill-shoal bodies. The location where complex waves begin is the nearshore boundary of the hill-shoal zone, and the location where the accretion of the hill-shoal body disappears is the nearshore boundary of the hill-shoal zone, thus determining the distribution location and extent of the hill-shoal zone. This invention effectively solves the problems of inconsistent standards, multiple interpretations, and irregular distribution in previous hill-shoal seismic facies identification methods.

[0162] Example 6:

[0163] Based on the same inventive concept, this application also provides an identification device for the inner hilly and shoal zones of ancient carbonate platform, which can be used to implement the method described in the above embodiments, as shown in the following embodiments. Since the principle of the identification device for the inner hilly and shoal zones of ancient carbonate platform is similar to the identification method for the inner hilly and shoal zones of ancient carbonate platform, the implementation of the identification device for the inner hilly and shoal zones of ancient carbonate platform can refer to the implementation of the identification method for the inner hilly and shoal zones of ancient carbonate platform, and repeated details will not be elaborated further. As used below, the term "unit" or "module" can refer to a combination of software and / or hardware that performs a predetermined function. Although the system described in the following embodiments is preferably implemented in software, hardware implementation, or a combination of software and hardware, is also possible and contemplated.

[0164] The present invention provides a specific implementation of a device for identifying the inner hilly and shoal zones of ancient carbonate platforms, capable of realizing a method for identifying these zones. See [link to specific implementation details]. Figure 12 An identification device for an inner hilly shoal zone of an ancient carbonate platform, comprising:

[0165] The three-parameter model establishment module 10 is used to establish a thickness-valley number-cloud-ash ratio relationship model based on the formation thickness, valley number and cloud-ash ratio determined by the first well logging data in the target work area.

[0166] The planar distribution map generation module 20 is used to determine the cloud-gray ratio corresponding to the seismic data based on the thickness-valley number-cloud-gray ratio relationship model, according to the stratum thickness and valley number determined by the seismic data in the target work area, and then generate a cloud-gray ratio planar distribution map of the target work area.

[0167] The change interface determination module 30 is used to determine the cloud-to-gray ratio change interface based on the cloud-to-gray ratio plane distribution map.

[0168] The seismic facies feature determination module 40 is used to determine the seismic facies features corresponding to the cloud-to-gray ratio change interface based on the seismic data.

[0169] The hill and shoal zone identification module 50 is used to identify the hill and shoal zone within the ancient carbonate platform based on the seismic facies characteristics.

[0170] In some embodiments of the present invention, the seismic facies feature determination module includes:

[0171] The seismic facies feature determination unit is used to determine the seismic facies features of the windward side and the leeward side of the hilly and shoal zone within the ancient carbonate platform based on the seismic data corresponding to the cloud-to-gray ratio change interface.

[0172] In some embodiments of the present invention, the seismic phase characteristics of the windward surface include a change in the seismic reflection phase axis from disordered, thickened reflections to parallel reflections;

[0173] The seismic phase characteristics of the leeward side include the transformation of the seismic reflection phase axis from parallel reflection to complex waves and multiple axes.

[0174] In some embodiments of the present invention, the speed at which the chaotic, thickened reflection transforms into parallel reflection is greater than the speed at which the parallel reflection transforms into complex waves or multi-axis reflection.

[0175] In some embodiments of the present invention, the first logging data is formation element logging data; the three-parameter model establishment module includes:

[0176] The energy spectrum determination unit is used to determine the corresponding captured gamma energy spectrum based on the formation element logging data.

[0177] A calcium and magnesium content determination unit is used to determine the calcium and magnesium content based on the captured gamma spectrum.

[0178] A rock content determination unit is used to determine the dolomite content and limestone content based on the calcium content and the magnesium content;

[0179] The cloud-to-ash ratio determination unit is used to determine the cloud-to-ash ratio based on the dolomite content and the limestone content.

[0180] In some embodiments of the present invention, an identification device for an inner hilly shoal zone of an ancient carbonate platform further includes:

[0181] The planar distribution map correction module is used to correct the cloud-to-gray ratio planar distribution map based on the second logging data and well logging data within the target work area.

[0182] As described above, this invention provides an identification device for hill-shoal zones within an ancient carbonate platform, comprising: a three-parameter model establishment module, used to establish a thickness-valley number-cloud-ash ratio relationship model based on the stratigraphic thickness, valley number, and cloud-ash ratio determined by first well logging data within the target work area; a planar distribution map generation module, used to determine the cloud-ash ratio corresponding to the seismic data based on the thickness-valley number-cloud-ash ratio relationship model, the stratigraphic thickness, and valley number determined by seismic data within the target work area, and thereby generate a planar distribution map of the cloud-ash ratio of the target work area; a change interface determination module, used to determine the cloud-ash ratio change interface based on the cloud-ash ratio planar distribution map; a seismic facies characteristic determination module, used to determine the seismic facies characteristics corresponding to the cloud-ash ratio change interface based on the seismic data; and a hill-shoal zone identification module, used to identify the hill-shoal zone within the ancient carbonate platform based on the seismic facies characteristics.

[0183] In summary, this invention improves the accuracy of hill and shoal prediction in areas with low exploration levels by statistically analyzing the cloud-to-gray ratio, formation thickness, and number of extracted valleys of the target layer, quantitatively establishing the relationship between the three, and predicting the planar distribution characteristics of hill and shoal development zones based on this relationship.

[0184] This application also provides a specific implementation of an electronic device capable of performing all steps in the method for identifying the inner hilly shoal zone of the ancient carbonate platform described in the above embodiments. See [link to implementation details]. Figure 13 The electronic devices specifically include the following:

[0185] Processor 1201, memory 1202, communications interface 1203, and bus 1204;

[0186] The processor 1201, memory 1202, and communication interface 1203 communicate with each other via bus 1204; the communication interface 1203 is used to realize information transmission between server-side devices and client-side devices and other related devices.

[0187] The processor 1201 is used to call the computer program in the memory 1202. When the processor executes the computer program, it implements all the steps in the method for identifying the hilly and shoal zones of the ancient carbonate platform in the above embodiment. For example, when the processor executes the computer program, it implements the following steps:

[0188] Based on the formation thickness, number of valleys, and cloud-to-gray ratio determined by the first well logging data within the target work area, a model relating thickness, number of valleys, and cloud-to-gray ratio is established.

[0189] Based on the thickness-valley number-cloud-gray ratio relationship model, the cloud-gray ratio corresponding to the seismic data is determined according to the stratum thickness and valley number determined by the seismic data in the target work area, and then a cloud-gray ratio planar distribution map of the target work area is generated.

[0190] The cloud-to-gray ratio variation interface is determined based on the cloud-to-gray ratio planar distribution map.

[0191] Based on the earthquake data, determine the seismic facies characteristics corresponding to the cloud-to-gray ratio change interface;

[0192] The ancient carbonate platform inner hilly and shoal zone was identified based on the seismic facies characteristics.

[0193] Embodiments of this application also provide a computer-readable storage medium capable of implementing all steps of the method for identifying the inner hilly and shoal zones of ancient carbonate platforms in the above embodiments. The computer-readable storage medium stores a computer program that, when executed by a processor, implements all steps of the method for identifying the inner hilly and shoal zones of ancient carbonate platforms in the above embodiments. For example, when the processor executes the computer program, it implements the following steps:

[0194] Based on the formation thickness, number of valleys, and cloud-to-gray ratio determined by the first well logging data within the target work area, a model relating thickness, number of valleys, and cloud-to-gray ratio is established.

[0195] Based on the thickness-valley number-cloud-gray ratio relationship model, the cloud-gray ratio corresponding to the seismic data is determined according to the stratum thickness and valley number determined by the seismic data in the target work area, and then a cloud-gray ratio planar distribution map of the target work area is generated.

[0196] The cloud-to-gray ratio variation interface is determined based on the cloud-to-gray ratio planar distribution map.

[0197] Based on the earthquake data, determine the seismic facies characteristics corresponding to the cloud-to-gray ratio change interface;

[0198] The ancient carbonate platform inner hilly and shoal zone was identified based on the seismic facies characteristics.

[0199] The various embodiments in this specification are described in a progressive manner. Similar or identical parts between embodiments can be referred to interchangeably. Each embodiment focuses on its differences from other embodiments. In particular, hardware + program embodiments are relatively simple in description because they are fundamentally similar to method embodiments; relevant parts can be referred to the descriptions in the method embodiments.

[0200] The foregoing has described specific embodiments of this specification. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recited in the claims may be performed in a different order than that shown in the embodiments and may still achieve the desired result. Furthermore, the processes depicted in the drawings do not necessarily require the specific or sequential order shown to achieve the desired result. In some embodiments, multitasking and parallel processing are possible or may be advantageous.

[0201] While this application provides method operation steps as shown in the embodiments or flowcharts, more or fewer operation steps may be included based on conventional or non-inventive labor. The order of steps listed in the embodiments is merely one possible execution order among many and does not represent the only execution order. In actual device or client product execution, the method can be executed sequentially as shown in the embodiments or drawings, or in parallel (e.g., in a parallel processor or multi-threaded processing environment).

[0202] For ease of description, the above devices are described in terms of function, divided into various modules. Of course, in implementing the embodiments of this specification, the functions of each module can be implemented in one or more software and / or hardware components, or a module that performs the same function can be implemented by a combination of multiple sub-modules or sub-units. The device embodiments described above are merely illustrative. For example, the division of units is only a logical functional division; in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Furthermore, the coupling or direct coupling or communication connection shown or discussed may be through some interfaces, or indirect coupling or communication connection between devices or units, and may be electrical, mechanical, or other forms.

[0203] Those skilled in the art will also know that, besides implementing the controller using purely computer-readable program code, the same functions can be achieved by logically programming the method steps, making the controller function as logic gates, switches, application-specific integrated circuits (ASICs), programmable logic controllers (PLCs), and embedded microcontrollers. Therefore, such a controller can be considered a hardware component, and the devices within it used to implement various functions can also be considered structures within that hardware component. Alternatively, the devices used to implement various functions can be considered as both software modules implementing the method and structures within a hardware component.

[0204] In a typical configuration, a computing device includes one or more processors (CPU), input / output interfaces, network interfaces, and memory.

[0205] Memory may include non-persistent storage in computer-readable media, such as random access memory (RAM) and / or non-volatile memory, such as read-only memory (ROM) or flash RAM. Memory is an example of computer-readable media.

[0206] The various embodiments in this specification are described in a progressive manner. Similar or identical parts between embodiments can be referred to mutually. Each embodiment focuses on describing the differences from other embodiments. In particular, system embodiments are basically similar to method embodiments, so the description is relatively simple; relevant parts can be referred to the descriptions in the method embodiments. In the description of this specification, the terms "one embodiment," "some embodiments," "example," "specific example," or "some examples," etc., refer to specific features, structures, materials, or characteristics described in connection with that embodiment or example, which are included in at least one embodiment or example of the embodiments in this specification. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described can be combined in any suitable manner in one or more embodiments or examples. Moreover, without contradiction, those skilled in the art can combine and integrate the different embodiments or examples described in this specification and the features of different embodiments or examples.

[0207] The above description is merely an embodiment of the present specification and is not intended to limit the embodiments of the present specification. For those skilled in the art, various modifications and variations can be made to the embodiments of the present specification. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principle of the embodiments of the present specification should be included within the scope of the claims of the embodiments of the present specification.

Claims

1. A method for identifying hilly and shoal zones within ancient carbonate platform terrain, characterized in that, include: Based on the formation thickness, number of valleys, and cloud-to-gray ratio determined by the first well logging data within the target work area, a model relating thickness, number of valleys, and cloud-to-gray ratio is established. Based on the thickness-valley number-cloud-gray ratio relationship model, the cloud-gray ratio corresponding to the seismic data is determined according to the stratum thickness and valley number determined by the seismic data in the target work area, and then a cloud-gray ratio planar distribution map of the target work area is generated. The cloud-to-gray ratio variation interface is determined based on the cloud-to-gray ratio planar distribution map. Based on the earthquake data, determine the seismic facies characteristics corresponding to the cloud-to-gray ratio change interface; The ancient carbonate platform inner hill and shoal zone was identified based on the seismic facies characteristics. Based on the seismic data, the seismic facies characteristics corresponding to the cloud-to-gray ratio change interface are determined, including: Based on the seismic data corresponding to the cloud-to-gray ratio variation interface, the seismic facies characteristics of the windward side and the leeward side of the hilly and shoal zone within the ancient carbonate platform are determined. The seismic phase characteristics of the windward side include the transformation of the seismic reflection phase axis from disordered, thickened, mound-like reflections to parallel reflections; The seismic phase characteristics of the leeward side include the transformation of the seismic reflection in-phase axis from parallel reflection to complex wave and multi-axis reflection; The first logging data is formation element logging data; determining the cloud-to-ash ratio based on the first logging data includes: The corresponding captured gamma spectrum is determined based on the well logging data of the formation elements. The calcium and magnesium content was determined based on the captured gamma spectrum. The dolomite content and limestone content are determined based on the calcium content and the magnesium content; The dolomite-lime ratio is determined based on the dolomite content and the limestone content.

2. The method for identifying the hilly and shoal zone within an ancient carbonate platform according to claim 1, characterized in that, The speed at which the chaotic, thickened reflection transforms into parallel reflection is greater than the speed at which parallel reflection transforms into complex waves or multi-axis reflection.

3. The method for identifying the hilly and shoal zone within an ancient carbonate platform according to any one of claims 1 to 2, characterized in that, Before determining the cloud-to-gray ratio variation interface based on the cloud-to-gray ratio planar distribution map, the method further includes: The cloud-to-gray ratio plane distribution map is corrected based on the second logging data and well logging data within the target work area.

4. An identification device for the ancient carbonate platform inner hilly and shoal zone used to implement the identification method for the ancient carbonate platform inner hilly and shoal zone according to any one of claims 1-3, characterized in that, include: The three-parameter model building module is used to establish a thickness-valley number-cloud-ash ratio relationship model based on the formation thickness, valley number, and cloud-ash ratio determined by the first well logging data in the target work area. The planar distribution map generation module is used to determine the cloud-gray ratio corresponding to the seismic data based on the thickness-valley number-cloud-gray ratio relationship model, according to the stratum thickness and valley number determined by the seismic data in the target work area, and then generate a planar distribution map of the cloud-gray ratio of the target work area. The change interface determination module is used to determine the cloud-to-gray ratio change interface based on the cloud-to-gray ratio planar distribution map. The seismic facies feature determination module is used to determine the seismic facies features corresponding to the cloud-to-gray ratio change interface based on the seismic data. The hill and shoal zone identification module is used to identify the hill and shoal zone within the ancient carbonate platform based on the seismic facies characteristics.

5. 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 for identifying the hilly and shoal zone within the ancient carbonate platform as described in any one of claims 1 to 3.

6. 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 for identifying the hilly and shoal zones within the ancient carbonate platform as described in any one of claims 1 to 3.

7. 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 for identifying the hilly and shoal zones within the ancient carbonate platform as described in any one of claims 1 to 3.

Images