Lithologic trap evaluation method, apparatus and device, medium and computer program
By obtaining distribution maps of seismic facies, reservoirs, and sealing layers and combining them with wave impedance relationships, the accuracy problem of lithological trap evaluation was solved, and the exploration success rate was improved.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- CHINA PETROLEUM & CHEMICAL CORP
- Filing Date
- 2024-11-25
- Publication Date
- 2026-05-26
AI Technical Summary
The lack of quantitative evaluation methods in existing technologies leads to poor accuracy in lithological trap evaluation results, which in turn affects the exploration success rate.
In the case of single-well sedimentary facies, by obtaining the seismic facies planar distribution map, reservoir porosity planar distribution map, and sealing layer porosity planar distribution map, and combining the wave impedance relationship, quantitative prediction of seismic facies, reservoir, and sealing layer can be achieved, thereby determining the evaluation results of lithological traps.
It improves the accuracy of lithological trap evaluation, increases the exploration success rate, and provides a quantitative basis for lithological trap evaluation.
Smart Images
Figure CN122087294A_ABST
Abstract
Description
Technical Field
[0001] This disclosure relates to the field of oil and gas exploration technology, and in particular to a method, apparatus, equipment, medium and computer program for evaluating lithological traps. Background Technology
[0002] Trap evaluation is a crucial task in petroleum exploration and a primary basis for exploration planning. The content of trap evaluation mainly includes geological assessment of the trap, resource calculation, and, where conditions permit, economic evaluation and exploration risk analysis. With the gradual improvement of the trap evaluation system, the order of trap evaluation steps is as follows: trap identification, trap accumulation condition evaluation, trap resource calculation, trap economic evaluation, and comprehensive trap evaluation.
[0003] Due to the unique characteristics of hydrocarbon accumulation conditions in lithological traps, special attention must be paid to evaluating these unique conditions in addition to standard evaluation criteria. Current techniques typically rely on qualitative analysis based on seismic attributes, lacking quantitative evaluation methods and approaches. This results in poor accuracy in lithological trap evaluations, leading to a low success rate in lithological trap exploration. Summary of the Invention
[0004] This disclosure provides a method, apparatus, equipment, medium, and computer program for evaluating lithological traps to overcome the problems existing in related technologies.
[0005] Firstly, this disclosure provides a method for evaluating lithological traps, including:
[0006] When the sedimentary facies type is a single-well sedimentary facies, a seismic facies plane distribution map corresponding to the target area is obtained. The seismic facies plane distribution map is used to characterize the distribution of multiple seismic facies in the target area.
[0007] Based on the first correspondence and the seismic phase plane distribution map, a reservoir porosity plane distribution map is obtained. The first correspondence is the correspondence between reservoir porosity and wave impedance. The reservoir porosity plane distribution map is used to characterize the distribution of the reservoir in each seismic phase.
[0008] Based on the second correspondence and the seismic phase plane distribution map, a sealing layer porosity plane distribution map is obtained. The second correspondence is the correspondence between sealing layer porosity and wave impedance. The sealing layer porosity plane distribution map is used to characterize the distribution of the sealing layer in each seismic phase.
[0009] Based on the seismic facies plane distribution, the reservoir porosity plane distribution map, and the sealing layer porosity plane distribution map, the evaluation results of the lithological traps are determined.
[0010] In some embodiments, obtaining the seismic facies plane distribution map corresponding to the target area includes:
[0011] Based on the seismic data of the target area, a well-through reservoir seismic facies analysis was performed to obtain the first seismic reflection characteristics;
[0012] Evolution simulations are performed based on reservoir seismic forward modeling and relevant reservoir information to obtain second seismic reflection characteristics. The relevant reservoir information includes the reservoir's sedimentary structure, lithology, physical properties, thickness, and rock physical parameters.
[0013] By comparing the first seismic reflection feature and the second seismic reflection feature, a target seismic reflection feature is determined. The target seismic reflection feature is used to characterize the reservoir identification result of each reservoir in multiple reservoirs in the target area.
[0014] Based on the target earthquake reflection characteristics, a target earthquake attribute is selected from multiple preset earthquake attributes. The correlation between the target earthquake attribute and the target earthquake reflection characteristics is greater than the correlation between the target earthquake attribute and other preset earthquake attributes among the multiple preset earthquake attributes.
[0015] Based on the first seismic reflection characteristics and the target seismic attributes, the seismic phase plane distribution map is determined.
[0016] In some embodiments, obtaining the reservoir porosity plane distribution map based on the first correspondence and the seismic facies plane distribution map includes:
[0017] Based on the logging curves, the reservoir porosity of each reservoir in the target area is determined;
[0018] Based on the first correspondence, the reservoir porosity and porosity threshold of each reservoir, a first wave impedance threshold is determined. The porosity threshold is used to characterize the minimum reservoir porosity that the reservoir needs to meet, and the first wave impedance is used to characterize the wave impedance condition that the reservoir needs to meet.
[0019] Based on the first wave impedance threshold, the porosity threshold, and the seismic phase plane distribution map, the reservoir porosity plane distribution map is obtained.
[0020] In some embodiments, obtaining the porosity plane distribution map of the sealing layer based on the second correspondence and the seismic facies plane distribution map includes:
[0021] Based on the well logging curves, the porosity of each packer in the target area is determined;
[0022] Based on the second correspondence, the porosity and porosity threshold of each sealing layer, a second wave impedance threshold is determined. The porosity threshold is used to characterize the maximum sealing layer porosity that the sealing layer needs to meet, and the second wave impedance threshold is used to characterize the wave impedance condition that the sealing layer needs to meet.
[0023] Based on the second wave impedance threshold, the porosity threshold, and the seismic phase plane distribution map, the porosity plane distribution map of the sealing layer is obtained.
[0024] In some embodiments, the method further includes:
[0025] Based on the third correspondence, the breakthrough pressure corresponding to the porosity of each sealing layer is determined. The third correspondence is the correspondence between the porosity of the sealing layer and the breakthrough pressure.
[0026] The determination of the second wave impedance threshold based on the second correspondence, the porosity of each sealing layer, and the porosity threshold of each sealing layer includes:
[0027] Based on the second correspondence, the porosity of each sealing layer, the porosity threshold, the fourth correspondence, the breakthrough pressure and the breakthrough pressure threshold of each sealing layer, the second wave impedance threshold is determined. The fourth correspondence is the correspondence between breakthrough pressure and wave impedance. The breakthrough pressure threshold is used to characterize the minimum breakthrough pressure that the sealing layer needs to meet.
[0028] In some embodiments, determining the evaluation results of lithological traps based on the seismic facies plane distribution, the reservoir porosity plane distribution map, and the sealing layer porosity plane distribution map includes:
[0029] Based on the preset boundary setting information of lithological traps, the seismic facies plane distribution, the reservoir porosity plane distribution map and the sealing layer porosity plane distribution map are overlaid to obtain a lithological trap plane distribution map. The lithological trap plane distribution map is used to characterize the distribution of multiple lithological traps in the target area.
[0030] The evaluation result of each lithological trap is determined based on the ratio between the area of the target seismic facies and the target reservoir in each lithological trap and the area of the lithological trap.
[0031] Secondly, this disclosure provides an evaluation device for lithological traps, comprising:
[0032] The first distribution map acquisition module is configured to acquire a seismic facies plane distribution map corresponding to the target area when the sedimentary facies type is a single-well sedimentary facies. The seismic facies plane distribution map is used to characterize the distribution of multiple seismic facies in the target area.
[0033] The second distribution map acquisition module is configured to acquire a reservoir porosity plane distribution map based on a first correspondence and the seismic phase plane distribution map. The first correspondence is the correspondence between reservoir porosity and wave impedance. The reservoir porosity plane distribution map is used to characterize the distribution of the reservoir in each seismic phase.
[0034] The third distribution map acquisition module is configured to acquire a sealing layer porosity plane distribution map based on the second correspondence and the seismic phase plane distribution map. The second correspondence is the correspondence between sealing layer porosity and wave impedance. The sealing layer porosity plane distribution map is used to characterize the distribution of the sealing layer in each seismic phase.
[0035] The evaluation module is configured to determine the evaluation results of lithological traps based on the seismic facies plane distribution, the reservoir porosity plane distribution map, and the sealing layer porosity plane distribution map.
[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 lithological trap evaluation method described in the first aspect above.
[0037] Fourthly, this disclosure provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the lithological trap evaluation method described in the first aspect.
[0038] Fifthly, this disclosure provides a computer program product, including a computer program / instructions, which, when executed by a processor, implements the steps of the lithological trap evaluation method described in the first aspect.
[0039] This disclosure provides a method, apparatus, equipment, medium, and computer program for evaluating lithological traps. Under single-well sedimentary facies conditions, it acquires seismic facies planar distribution maps, reservoir porosity planar distribution maps, and sealing layer porosity planar distribution maps, enabling prediction of the distribution of seismic facies, reservoir, and sealing layer. Furthermore, when predicting the distribution of reservoir and sealing layer, it considers the wave impedance corresponding to porosity, achieving quantitative prediction. Therefore, based on the distribution of seismic facies, reservoir, and sealing layer, it obtains evaluation results for lithological traps, improving the accuracy of lithological trap evaluation results and thus increasing the success rate of lithological trap exploration. Attached Figure Description
[0040] The present disclosure will be described in more detail below based on embodiments and with reference to the accompanying drawings:
[0041] Figure 1This is a flowchart illustrating a method for evaluating lithological traps provided in an embodiment of this disclosure.
[0042] Figure 2 This is a schematic diagram of a seismic phase plane distribution map provided in an embodiment of the present disclosure.
[0043] Figure 3 A schematic diagram illustrating the relationship between porosity and wave impedance provided in an embodiment of this disclosure;
[0044] Figure 4 This is a schematic diagram of the planar distribution of reservoir porosity and sealing layer porosity provided in an embodiment of this disclosure.
[0045] Figure 5 This is a schematic diagram of a planar distribution map of lithological traps provided in an embodiment of this disclosure.
[0046] Figure 6 This is a schematic diagram of a system for evaluating lithological traps, provided as an embodiment of the present disclosure.
[0047] Figure 7 This is a block diagram of an evaluation device for lithological traps provided in an embodiment of this disclosure.
[0048] 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
[0049] 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.
[0050] 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.
[0051] 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.
[0052] Due to the unique formation conditions of lithological traps, the evaluation of these traps requires special attention beyond standard assessment criteria. Previous evaluations primarily focused on clastic strata; however, recent discoveries of numerous lithological traps in marine carbonate rocks highlight their importance. Carbonate traps are often associated with high-energy sedimentary facies zones and source-transfer faults, exhibiting characteristics of "lithological body-controlled reservoirs and localized structural enrichment." Exploration success rates are high in structurally high areas, but vary significantly on structural slopes and lower areas. This is mainly because reservoirs are often controlled by multiple factors, including sedimentary facies, paleokarst, and fractures, resulting in strong heterogeneity and making it difficult to predict the sealing conditions of lithological traps in structurally lower areas. Current technologies typically rely on qualitative analysis based on seismic attributes, lacking quantitative evaluation methods and approaches, leading to poor accuracy in lithological trap evaluations and consequently low exploration success rates. Therefore, a more targeted evaluation method is needed to accurately identify and evaluate lithological traps.
[0053] The lithological trap evaluation method provided in this disclosure, under the condition of single-well sedimentary facies, obtains seismic facies planar distribution maps, reservoir porosity planar distribution maps, and sealing layer porosity planar distribution maps, respectively, to predict the distribution of seismic facies, reservoir, and sealing layer. Furthermore, when predicting the distribution of reservoir and sealing layer, the wave impedance corresponding to porosity is considered, achieving quantitative prediction. Thus, based on the distribution of seismic facies, reservoir, and sealing layer, the evaluation results of lithological traps are obtained, improving the accuracy of lithological trap evaluation results and thereby increasing the success rate of lithological trap exploration.
[0054] The method provided in this disclosure is executed by a computer device, which may be a mobile phone, tablet computer, laptop computer, desktop computer, or other similar device. This disclosure does not limit the specific type of computer device.
[0055] Example 1
[0056] Figure 1 This is a flowchart illustrating a method for evaluating lithological traps provided in an embodiment of this disclosure, executed by a computer device. See [link / reference]. Figure 1 The method includes the following steps:
[0057] Step S101: When the sedimentary facies type is single-well sedimentary facies, obtain the seismic facies plane distribution map corresponding to the target area. The seismic facies plane distribution map is used to characterize the distribution of multiple seismic facies in the target area.
[0058] Step S102: Based on the first correspondence and the seismic phase plane distribution map, obtain the reservoir porosity plane distribution map. The first correspondence is the correspondence between reservoir porosity and wave impedance. The reservoir porosity plane distribution map is used to characterize the distribution of the reservoir in each seismic phase.
[0059] Step S103: Based on the second correspondence and the seismic phase plane distribution map, obtain the sealing layer porosity plane distribution map. The second correspondence is the correspondence between the sealing layer porosity and wave impedance. The sealing layer porosity plane distribution map is used to characterize the distribution of the sealing layer in each seismic phase.
[0060] Step 104: Based on the seismic facies plane distribution map, reservoir porosity plane distribution map, and sealing layer porosity plane distribution map, determine the evaluation results of the lithological traps.
[0061] In this embodiment of the disclosure, the current sedimentary facies type is first determined to be a single-well sedimentary facies through core observation description and well logging sedimentary facies analysis. Then, based on the single-well sedimentary facies, the lithological traps are evaluated.
[0062] In step S101, a seismic facies plane distribution map corresponding to the target area is obtained. This seismic facies plane distribution map characterizes the distribution of multiple seismic facies in the target area, which is the area to be explored. The seismic facies include a first seismic facies favorable to reservoir development and a second seismic facies unfavorable to reservoir development.
[0063] In some embodiments, the process of obtaining a seismic facies plane distribution map includes the following steps:
[0064] Step 1-1: Perform well-through reservoir seismic facies analysis based on seismic data of the target area to obtain the first seismic reflection characteristics.
[0065] Among them, seismic data is used to characterize the actual geological conditions of the target area, and the first seismic reflection feature is used to characterize the actual reservoir identification results of each reservoir in multiple reservoirs in the target area.
[0066] Optionally, the seismic facies parameter characteristics of the reservoir obtained through the well can be used to select the first seismic reflection feature that is favorable to the reservoir. These seismic facies parameters include reflection structure, continuity, external geometry, amplitude, frequency, and layer velocity.
[0067] Steps 1-2 involve performing evolution simulations based on the reservoir seismic forward model and relevant reservoir information to obtain the second seismic reflection characteristics.
[0068] Among them, reservoir-related information includes the reservoir's sedimentary structure, lithology, physical properties, thickness, and rock physical parameters. The second seismic reflection feature is used to characterize the reservoir identification results of the evolution simulation of each reservoir in multiple reservoirs in the target area.
[0069] Optionally, a reservoir seismic forward model can be established. Based on relevant reservoir information, an evolutionary simulation can be performed using the non-homogeneous medium wave equation and staggered grid finite difference calculation to obtain the second seismic reflection characteristics. During the evolutionary simulation, new second seismic reflection characteristics can also be obtained by changing factors such as reservoir structure morphology, spatial development pattern, and development location (top, middle, and bottom of carbonate rocks).
[0070] Steps 1-3: Compare the first and second earthquake reflection characteristics to determine the target earthquake reflection characteristics.
[0071] Since the first and second seismic reflection features respectively characterize the actual and simulated reservoir identification results, a target seismic reflection feature is obtained by comparing the first and second seismic reflection features and comprehensively considering the actual and simulated reservoir identification results. This target seismic reflection feature is used to characterize the reservoir identification result of each reservoir in multiple reservoirs in the target area. This target seismic reflection feature is more accurate than a single first seismic reflection feature or a single second seismic reflection feature. Therefore, a more comprehensive and accurate reservoir identification marker can be established based on this target seismic reflection feature, thereby achieving more accurate reservoir identification.
[0072] Steps 1-4: Select the target seismic attribute from multiple preset seismic attributes based on the target seismic reflection characteristics.
[0073] Among them, multiple preset seismic attributes are seismic attributes that can reflect the reservoir. The correlation between the target seismic attribute and the target seismic reflection characteristics is greater than the correlation between the other preset seismic attributes (excluding the target seismic attribute) and the target seismic reflection characteristics. In other words, the target seismic attribute is the seismic attribute that can better reflect the reservoir among the multiple preset seismic attributes.
[0074] Optionally, based on the first seismic reflection characteristics, the second seismic reflection characteristics, and the target seismic reflection characteristics, multiple preset seismic attributes are initially selected to obtain multiple candidate seismic attributes. Then, seismic attributes of well-side or interconnected well profiles are calculated, and the overall anomaly distribution pattern of the seismic attributes is analyzed using methods such as cross-plotting. Based on the overall anomaly distribution pattern of the seismic attributes, candidate seismic attributes that have a clear correspondence with reservoir characteristics are preprocessed to form a seismic attribute set. Then, cluster analysis methods such as pattern recognition, cluster analysis, and neural networks are used to further analyze the correlation between the seismic attributes in the seismic attribute set and the target seismic reflection characteristics, and the target seismic attribute is selected. For example, the selected target seismic attribute can be frequency-division reflection intensity and wave trough amplitude.
[0075] Steps 1-5: Based on the first earthquake reflection characteristics and the target earthquake attributes, determine the earthquake phase plane distribution map.
[0076] Since the first seismic reflection feature reflects the actual reservoir identification result for each reservoir, a more accurate seismic facies plane distribution map is determined based on the first seismic reflection feature and combined with the target seismic attributes related to reservoir height. In one example, the seismic facies plane distribution map is as follows: Figure 2 As shown, based on the sedimentary facies division of a single well and the waveform characteristics analysis of the well-through seismic profile, the inter-hill reservoir in the target area is determined to have the characteristics of "micro-amplitude hillock-like, mid-to-low frequency, discontinuous reflection," while the inter-hill sealing layer has the characteristics of "parallel, low frequency, strong peak reflection." This yields the first seismic reflection characteristic. Then, seismic forward modeling analysis is performed to obtain the second seismic reflection characteristic. Based on the first and second seismic reflection characteristics, two target seismic attributes—frequency-based reflection intensity and trough amplitude—are selected. Finally, seismic facies distribution prediction is performed, resulting in the following... Figure 2 The diagram shows the seismic phase plane distribution.
[0077] In step S102, based on the obtained seismic facies plane distribution map, the reservoir is quantitatively predicted.
[0078] In some embodiments, the process of obtaining a reservoir porosity planar distribution map includes the following steps:
[0079] Step 2-1: Based on the logging curves, determine the reservoir porosity of each reservoir in the target area.
[0080] To evaluate the rock properties of the reservoir development zone, the reservoir porosity of each reservoir in the target area is calculated based on the logging curves. Optionally, if core sampling is available, the reservoir porosity calculated from the logging is corrected using the core-measured porosity.
[0081] Step 2-2: Based on the first correspondence, the reservoir porosity and porosity threshold of each reservoir, determine the first wave impedance threshold.
[0082] The porosity threshold is used to characterize the minimum reservoir porosity that the reservoir must meet, and the first wave impedance is used to characterize the wave impedance condition that the reservoir must meet. For example, the porosity threshold is 2%, or other smaller values.
[0083] Optionally, using the single-well porosity calculation results (reservoir porosity for each reservoir) and pseudo-acoustic logging curves, a high-quality reservoir section of a single lithology is selected, and then a first correspondence is established. Based on the first correspondence, regions with reservoir porosity > 2% are selected to obtain the first wave impedance threshold. In one example, taking wave impedance as the P-wave impedance, a schematic diagram of the correspondence between porosity and wave impedance is shown below. Figure 3 As shown, Figure 3 This includes the relationship between reservoir porosity and wave impedance, as well as the relationship between sealing layer porosity and wave impedance, from... Figure 3 As can be seen, the region with reservoir porosity > 2% represents the correlation between reservoir porosity and wave impedance (i.e., Figure 3 (The darker spots in the middle).
[0084] Steps 2-3: Based on the first wave impedance threshold, porosity threshold, and seismic facies plane distribution map, obtain the reservoir porosity plane distribution map.
[0085] By overlaying the seismic facies plane distribution map onto the reservoir impedance and porosity inversion results, a reservoir porosity plane distribution map is obtained. In one example, the reservoir porosity and sealing layer porosity plane distribution maps are as follows: Figure 4 As shown, that is Figure 4 The schematic diagram shown is a superimposed diagram of the reservoir porosity planar distribution and the sealing layer porosity planar distribution.
[0086] When evaluating lithological traps, the predicted results (reservoir identification results) of reservoirs within seismic facies favorable to reservoir development on the reservoir porosity plane distribution map can truly reflect the geological conditions. In areas outside the seismic facies favorable to reservoir development (seismic facies unfavorable to reservoir development), false images may occur due to lithological changes, and can only be used as a reference when evaluating lithological traps.
[0087] In step S103, based on the obtained seismic phase plane distribution map, the sealing layer is quantitatively predicted.
[0088] In some embodiments, the process of obtaining a planar distribution map of the porosity of the sealing layer includes the following steps:
[0089] Step 3-1: Based on the logging curves, determine the porosity of each packer in the target area.
[0090] To evaluate the rock properties of the sealing layer development section, the sealing layer porosity of each sealing layer in the target area is calculated based on the well logging curves.
[0091] Step 3-2: Based on the second correspondence, the porosity of each sealing layer and the porosity threshold, determine the second wave impedance threshold.
[0092] Among them, the porosity threshold is used to characterize the maximum porosity that the sealing layer must meet, and the second wave impedance threshold is used to characterize the wave impedance condition that the sealing layer must meet.
[0093] Optionally, a second correspondence is established using the single-well porosity calculation results (porosity of each compartment) and the pseudo-acoustic logging curves. Based on this second correspondence, regions with reservoir porosity <2% are selected to obtain the second wave impedance threshold. In one example, taking wave impedance as the P-wave impedance, a schematic diagram of the correspondence between porosity and wave impedance is shown below. Figure 3 As shown, Figure 3 This includes the relationship between reservoir porosity and wave impedance, as well as the relationship between sealing layer porosity and wave impedance, from... Figure 3 As can be seen, the region with reservoir porosity <2% represents the correlation between the porosity and wave impedance of the sealing layer (i.e., Figure 3 (Lighter-colored dots in the middle).
[0094] Optionally, when coring is available, the breakthrough pressure corresponding to the porosity of each packer is determined based on the third correspondence, which is the relationship between packer porosity and breakthrough pressure. That is, the breakthrough pressure of each packer is determined based on the logging curves and the third correspondence. Then, based on the second correspondence, the packer porosity and porosity threshold of each packer, the fourth correspondence, the breakthrough pressure of each packer, and the breakthrough pressure threshold, the second wave impedance threshold is determined. The fourth correspondence is the relationship between breakthrough pressure and wave impedance, and the breakthrough pressure threshold characterizes the minimum breakthrough pressure that the packer must meet. For example, the breakthrough pressure threshold is 5 MPa, or another value.
[0095] Step 3-3: Based on the second wave impedance threshold, porosity threshold, and seismic phase plane distribution map, obtain the porosity plane distribution map of the sealing layer.
[0096] By overlaying the seismic facies plane distribution map onto the results of wave impedance and porosity (and breakthrough pressure) inversion, a porosity plane distribution map of the sealing layer is obtained. In one example, the reservoir porosity and sealing layer porosity plane distribution maps are as follows: Figure 4 As shown, that is Figure 4 The schematic diagram shown is a superimposed diagram of the reservoir porosity planar distribution and the sealing layer porosity planar distribution.
[0097] When evaluating lithological traps, the predicted results of the sealing layer (non-reservoir) porosity distribution map in seismic facies unfavorable to reservoir development (sealing layer identification results) can truly reflect the geological conditions. In seismic facies favorable to reservoir development, the results may be misleading due to changes in lithology and physical properties, and can only be used as a reference when evaluating lithological traps.
[0098] In step S104, based on the obtained seismic facies plane distribution map, reservoir porosity plane distribution map, and sealing layer porosity plane distribution map, the lithological traps are evaluated to obtain the evaluation results of the lithological traps.
[0099] In some embodiments, the process for evaluating lithological traps includes the following steps:
[0100] Step 4-1: Based on the preset boundary setting information of lithological traps, the seismic facies plane distribution map, the reservoir porosity plane distribution map, and the sealing layer porosity plane distribution map are overlaid to obtain the lithological trap plane distribution map.
[0101] The preset boundary setting information refers to setting the boundary of each lithological trap to the reservoir porosity as the porosity threshold, or the reservoir porosity as the porosity threshold and the breakthrough pressure as the breakthrough pressure threshold. Based on the preset boundary setting information, the range of the lithological trap can be delineated in the planar distribution map obtained by overlaying the seismic facies plane distribution map, the reservoir porosity plane distribution map, and the sealing layer porosity plane distribution map, thus obtaining the lithological trap sealing plane distribution map. In one example, the lithological trap plane distribution map is as follows: Figure 5 As shown.
[0102] Among them, the lithological trap planar distribution map is used to characterize the distribution of multiple lithological traps in the target area.
[0103] Step 4-2: Determine the evaluation result for each lithological trap based on the ratio between the area of the target seismic facies and the target reservoir and the area of the lithological trap.
[0104] The larger the ratio, the more valuable the lithological trap is to explore; the smaller the ratio, the less valuable the lithological trap is to explore. Therefore, based on this ratio, we can determine whether each lithological trap is worth exploring. In other words, the evaluation result of each lithological trap indicates whether the lithological trap is worth exploring.
[0105] The lithological trap evaluation method provided in this disclosure, under the condition of single-well sedimentary facies, obtains seismic facies planar distribution maps, reservoir porosity planar distribution maps, and sealing layer porosity planar distribution maps, respectively, to predict the distribution of seismic facies, reservoir, and sealing layer. In predicting the distribution of reservoir and sealing layer, the wave impedance corresponding to porosity is considered, realizing quantitative prediction. Therefore, based on the distribution of seismic facies, reservoir, and sealing layer, the evaluation results of lithological traps are obtained, improving the accuracy of the evaluation results of lithological traps, making the evaluation results of lithological traps more reliable, and thus improving the exploration success rate of lithological traps.
[0106] Furthermore, the method provided in this disclosure combines waveform attributes and physical property inversion to identify the sealing layer of lithological traps and quantitatively evaluate their sealing performance, providing a technical basis for the exploration of lithological traps. Compared with traditional qualitative evaluation methods, this disclosure directly provides quantitative prediction results of the sealing performance of low-impedance sealing layers under phase-controlled constraints, thus having higher reliability and good prospects for widespread application.
[0107] Example 2
[0108] Based on the above embodiments, see Figure 6 The diagram shown is a schematic of a system for evaluating lithological traps. The system for evaluating lithological traps includes a well-through reservoir seismic facies analysis subsystem, a seismic attribute optimization subsystem, a seismic facies distribution prediction subsystem, a single-well reservoir evaluation subsystem, a reservoir physical property quantitative prediction subsystem, a single-well isolation layer evaluation subsystem, an isolation layer physical property quantitative prediction subsystem, and a lithological trap evaluation and optimization subsystem.
[0109] Among them, the well-through reservoir seismic facies analysis subsystem is used to determine the first seismic reflection characteristics, the seismic attribute optimization subsystem is used to determine the target seismic attributes, the seismic facies distribution prediction subsystem is used to determine the seismic facies plane distribution map, the single-well reservoir evaluation subsystem is used to determine the reservoir porosity, the reservoir physical property quantitative prediction subsystem is used to determine the reservoir porosity plane distribution map, the single-well sealing layer evaluation subsystem is used to determine the sealing layer porosity, the sealing layer physical property quantitative prediction subsystem is used to determine the sealing layer porosity screen distribution map, and the lithological trap evaluation and optimization subsystem is used to determine the evaluation results of the lithological trap.
[0110] Example 3
[0111] Based on the above embodiments, this embodiment provides an application example.
[0112] The evaluation process for lithological traps is as follows:
[0113] (1) Seismic facies analysis of reservoirs through wells
[0114] The sedimentary facies of individual wells were determined through core observation and description and well logging analysis. Based on the sedimentary facies of individual wells, seismic sedimentology studies were conducted using seismic data to analyze the seismic facies parameter characteristics of the reservoirs through the wells, and to summarize the seismic reflection characteristics (first seismic reflection characteristics) of favorable reservoirs and their geological significance.
[0115] (2) Seismic attribute optimization
[0116] Based on (1), a reservoir seismic forward model is established. Based on the reservoir seismic forward model, according to the sedimentary structure, lithology, physical properties and thickness of the target layer in the study area, the rock physical parameters are set. The non-homogeneous medium wave equation and the staggered grid finite difference calculation are used to carry out numerical simulation to obtain the second seismic reflection characteristics. By changing the reservoir structure morphology, spatial development law and development location (top, middle and bottom of carbonate rocks) and other factors, the geophysical response characteristics (target seismic reflection characteristics) of the target reservoir are studied, so as to establish a more comprehensive and accurate reservoir identification mark.
[0117] Based on geophysical response characteristics, and combining the first and second seismic reflection characteristics, initial seismic attributes are selected. Seismic attributes are calculated for well-side or interconnected well profiles, and cross-plotting is used. Based on the overall anomaly distribution pattern of seismic attributes, necessary preprocessing is performed on attributes that have a clear correspondence with reservoir characteristics to form a seismic attribute set. Using clustering analysis methods such as pattern recognition, cluster analysis, and neural networks, the correspondence between seismic attributes and reservoir characteristics is further analyzed, and the seismic attribute set is optimized to obtain the target seismic attributes.
[0118] (3) Seismic phase distribution prediction
[0119] Based on the results of (1) and (2), the seismic facies plane distribution of the study area was divided, and a seismic facies plane distribution map was obtained.
[0120] (4) Single-well reservoir evaluation
[0121] Based on the sedimentary facies of a single well, the rock properties of the reservoir development section are evaluated, and the reservoir porosity is calculated using well logging curves. If core sampling is available, the reservoir porosity calculated by well logging is corrected using the measured porosity from the core.
[0122] (5) Quantitative prediction of reservoir properties
[0123] Using the single-well porosity calculation results (reservoir porosity) and pseudo-acoustic logging curves, high-quality reservoir sections with a single lithology are selected, and the mathematical relationship between reservoir porosity and wave impedance is established. Areas with reservoir porosity > 2% are preferred, and the threshold value of wave impedance (first wave impedance threshold) is obtained.
[0124] By overlaying the seismic facies plane distribution map onto the results of wave impedance and reservoir porosity inversion, a reservoir porosity plane distribution map is obtained. Only the reservoir physical property prediction results within the favorable reservoir development facies zone can truly reflect the geological conditions. In areas outside the favorable facies zone, lithological changes may produce false images, which can only be used as a reference for trap evaluation.
[0125] (6) Evaluation of single-well sealing layer
[0126] Based on the sedimentary facies of a single well, the rock properties of the non-reservoir development section are evaluated, and the porosity of the sealing layer is calculated using logging curves. If core sampling is available, the breakthrough pressure measured in the core is used to establish the relationship between the breakthrough pressure and the porosity of the sealing layer. Then, the breakthrough pressure of the non-reservoir section sealing layer is calculated using logging curves.
[0127] (7) Quantitative prediction of sealing layer properties
[0128] Using the calculation results of the packing porosity (and breakthrough pressure) and the pseudo-acoustic logging curves, a mathematical relationship chart between packing porosity (and breakthrough pressure) and wave impedance is established. The region with packing porosity <2% (and breakthrough pressure >5MPa) is selected to determine the threshold value of packing wave impedance (second wave impedance threshold).
[0129] By overlaying the seismic facies plane distribution map onto the results of wave impedance and sealing layer porosity (and breakthrough pressure) inversion, a sealing layer porosity plane distribution map is obtained. Only the prediction results in non-reservoir facies zones can truly reflect the geological conditions. In favorable reservoir facies zones, the difference in lithology and physical properties may produce false images, which can only be used as a reference for trap evaluation.
[0130] (8) Evaluation and selection of lithological traps
[0131] By overlaying the seismic facies plane distribution map, the reservoir porosity plane distribution map, and the sealing layer property distribution map together, and setting the lithological trap boundary to a reservoir porosity of 2% (or a breakthrough pressure of 5 MPa), the lithological trap range is delineated, and a lithological trap plane distribution map is obtained.
[0132] The lithological traps are numbered according to their area size. Based on the ratio of the area occupied by the dominant facies and high-porosity reservoirs in the trap (based on the ratio between the area of the target seismic facies and the target reservoir in each lithological trap and the area of the lithological trap), the lithological traps are ranked and selected to be more favorable for exploration.
[0133] In one example, based on the sedimentary facies division of a single well and the waveform characteristics analysis of the seismic profile, the reservoir between the hills and shoals was identified as having "micro-amplitude hill-shaped, mid-to-low frequency, discontinuous reflections," while the sealing layer between the hills and shoals exhibited "parallel, low-frequency, strong peak reflections." Combined with seismic forward modeling analysis, two attributes—frequency-based reflection intensity and trough amplitude—were selected for seismic facies distribution prediction, resulting in a seismic facies plane distribution map. The core and physical property characteristics of reservoirs and sealing layers from multiple drilled wells in the adjacent area were comprehensively analyzed. Under the constraint of measured data, the porosity of the reservoir and sealing layers was recalculated using well logging curves. A statistical relationship between porosity and pseudo-acoustic impedance was established (e.g.,...). Figure 3 As shown), the evaluation wave impedances of the reservoir and the sealing layer were obtained respectively. Then, under phase control constraints, wave impedance and porosity inversion were performed to obtain the planar distribution maps of reservoir porosity and sealing layer porosity (as shown). Figure 4 (As shown); the seismic facies plane distribution map, reservoir porosity plane distribution map, and sealing layer porosity plane distribution map are shown. Figure 3 By overlaying planar distribution maps, a planar distribution map of lithological traps can be obtained (e.g., ...). Figure 5 As shown in the figure, by utilizing the degree of reservoir development and the degree of lithological trap confirmation within the trap, the traps were optimized and ranked, the most favorable lithological traps were determined, and risk wells were deployed. The drilling results showed that the gas layer thickness was 60m, which confirmed that the Dengying Formation in the lower part of the structure has the geological conditions for lithological reservoir formation, and also proved that the quantitative evaluation method of lithological trap sealing under phase control constraints is effective.
[0134] Example 4
[0135] Figure 7 A block diagram of an evaluation device for lithological traps provided in this disclosure, configured on a computer device, is shown below. Figure 7 The device includes:
[0136] The first distribution map acquisition module 701 is configured to acquire the seismic facies plane distribution map corresponding to the target area when the sedimentary facies type is single-well sedimentary facies. The seismic facies plane distribution map is used to characterize the distribution of multiple seismic facies in the target area.
[0137] The second distribution map acquisition module 702 is configured to acquire a reservoir porosity plane distribution map based on a first correspondence and a seismic phase plane distribution map. The first correspondence is the correspondence between reservoir porosity and wave impedance. The reservoir porosity plane distribution map is used to characterize the distribution of the reservoir in each seismic phase.
[0138] The third distribution map acquisition module 703 is configured to acquire a sealing layer porosity plane distribution map based on the second correspondence and the seismic phase plane distribution map. The second correspondence is the correspondence between sealing layer porosity and wave impedance. The sealing layer porosity plane distribution map is used to characterize the distribution of the sealing layer in each seismic phase.
[0139] Evaluation module 704 is configured to determine the evaluation results of lithological traps based on seismic facies plane distribution, reservoir porosity plane distribution map, and sealing layer porosity plane distribution map.
[0140] In some embodiments, the first distribution map acquisition module 701 is configured to:
[0141] Based on seismic data of the target area, well-through reservoir seismic facies analysis was performed to obtain the first seismic reflection characteristics;
[0142] Evolution simulations were performed based on reservoir seismic forward modeling and relevant reservoir information to obtain second seismic reflection characteristics. Reservoir information included sedimentary structure, lithology, physical properties, thickness, and rock physical parameters.
[0143] By comparing the first and second seismic reflection features, the target seismic reflection features are determined. The target seismic reflection features are used to characterize the reservoir identification results of each reservoir in multiple reservoirs in the target area.
[0144] Based on the target earthquake reflection characteristics, the target earthquake attribute is selected from multiple preset earthquake attributes. The correlation between the target earthquake attribute and the target earthquake reflection characteristics is greater than the correlation between the target earthquake attribute and other preset earthquake attributes among the multiple preset earthquake attributes.
[0145] Based on the first earthquake reflection characteristics and the target earthquake attributes, the earthquake phase plane distribution map is determined.
[0146] In some embodiments, the second distribution map acquisition module 702 is configured to:
[0147] Based on well logging curves, determine the reservoir porosity of each reservoir in the target area;
[0148] Based on the first correspondence, the reservoir porosity and porosity threshold of each reservoir, the first wave impedance threshold is determined. The porosity threshold is used to characterize the minimum reservoir porosity that the reservoir needs to meet, and the first wave impedance is used to characterize the wave impedance condition that the reservoir needs to meet.
[0149] Based on the first wave impedance threshold, porosity threshold, and seismic facies plane distribution map, the reservoir porosity plane distribution map is obtained.
[0150] In some embodiments, the third distribution map acquisition module 703 is configured to:
[0151] Based on well logging curves, determine the porosity of each packer layer in the target area;
[0152] Based on the second correspondence, the porosity of each sealing layer and the porosity threshold, the second wave impedance threshold is determined. The porosity threshold is used to characterize the maximum sealing layer porosity that the sealing layer needs to meet, and the second wave impedance threshold is used to characterize the wave impedance condition that the sealing layer needs to meet.
[0153] Based on the second wave impedance threshold, porosity threshold, and seismic phase plane distribution map, the porosity plane distribution map of the sealing layer is obtained.
[0154] In some embodiments, the third distribution map acquisition module 703 is configured to:
[0155] Based on the third correspondence, the breakthrough pressure corresponding to the porosity of each sealing layer is determined. The third correspondence is the relationship between the porosity of the sealing layer and the breakthrough pressure.
[0156] Based on the second correspondence, the porosity of each sealing layer, the porosity threshold, the fourth correspondence, the breakthrough pressure of each sealing layer, and the breakthrough pressure threshold, the second wave impedance threshold is determined. The fourth correspondence is the correspondence between breakthrough pressure and wave impedance. The breakthrough pressure threshold is used to characterize the minimum breakthrough pressure that the sealing layer needs to meet.
[0157] In some embodiments, the evaluation module 704 is configured to:
[0158] Based on the preset boundary setting information of lithological traps, the seismic facies plane distribution map, reservoir porosity plane distribution map and the sealing layer porosity plane distribution map are overlaid to obtain the lithological trap plane distribution map. The lithological trap plane distribution map is used to characterize the distribution of multiple lithological traps in the target area.
[0159] The evaluation result for each lithological trap is determined based on the ratio between the area of the target seismic facies and the target reservoir and the area of the lithological trap.
[0160] Example 5
[0161] Based on the above embodiments, this embodiment provides a computer device, including a memory, a processor, and a computer program stored in the memory. The processor executes the computer program to implement the steps of the lithological trap evaluation method described in the above embodiments.
[0162] In some embodiments of this example, a computer-readable storage medium is provided, on which a computer program is stored, characterized in that, when the computer program is executed by a processor, it implements the steps of the lithological trap evaluation method described in the above embodiments.
[0163] In some embodiments of this example, a computer program product is provided, including a computer program / instructions, characterized in that, when the computer program is executed by a processor, it implements the steps of the lithological trap evaluation method described in the above embodiments.
[0164] 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.
[0165] 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.).
[0166] 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.
[0167] 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.).
[0168] The processor can communicate with external devices via the I / O bus through wired or wireless networks.
[0169] 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.
[0170] 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.
[0171] 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.
[0172] 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 evaluating lithological traps, characterized in that, include: When the sedimentary facies type is a single-well sedimentary facies, a seismic facies plane distribution map corresponding to the target area is obtained. The seismic facies plane distribution map is used to characterize the distribution of multiple seismic facies in the target area. Based on the first correspondence and the seismic phase plane distribution map, a reservoir porosity plane distribution map is obtained. The first correspondence is the correspondence between reservoir porosity and wave impedance. The reservoir porosity plane distribution map is used to characterize the distribution of the reservoir in each seismic phase. Based on the second correspondence and the seismic phase plane distribution map, a sealing layer porosity plane distribution map is obtained. The second correspondence is the correspondence between sealing layer porosity and wave impedance. The sealing layer porosity plane distribution map is used to characterize the distribution of the sealing layer in each seismic phase. Based on the seismic facies plane distribution map, the reservoir porosity plane distribution map, and the sealing layer porosity plane distribution map, the evaluation results of the lithological traps are determined.
2. The method according to claim 1, characterized in that, The acquisition of the seismic facies plane distribution map corresponding to the target area includes: Based on the seismic data of the target area, a well-through reservoir seismic facies analysis was performed to obtain the first seismic reflection characteristics; Evolution simulations are performed based on reservoir seismic forward modeling and relevant reservoir information to obtain second seismic reflection characteristics. The relevant reservoir information includes the reservoir's sedimentary structure, lithology, physical properties, thickness, and rock physical parameters. By comparing the first seismic reflection feature and the second seismic reflection feature, a target seismic reflection feature is determined. The target seismic reflection feature is used to characterize the reservoir identification result of each reservoir in multiple reservoirs in the target area. Based on the target earthquake reflection characteristics, a target earthquake attribute is selected from multiple preset earthquake attributes. The correlation between the target earthquake attribute and the target earthquake reflection characteristics is greater than the correlation between the target earthquake attribute and other preset earthquake attributes among the multiple preset earthquake attributes. Based on the first seismic reflection characteristics and the target seismic attributes, the seismic phase plane distribution map is determined.
3. The method according to claim 1, characterized in that, The step of obtaining the reservoir porosity plane distribution map based on the first correspondence and the seismic facies plane distribution map includes: Based on the logging curves, the reservoir porosity of each reservoir in the target area is determined; Based on the first correspondence, the reservoir porosity and porosity threshold of each reservoir, a first wave impedance threshold is determined. The porosity threshold is used to characterize the minimum reservoir porosity that the reservoir needs to meet, and the first wave impedance is used to characterize the wave impedance condition that the reservoir needs to meet. Based on the first wave impedance threshold, the porosity threshold, and the seismic phase plane distribution map, the reservoir porosity plane distribution map is obtained.
4. The method according to claim 1, characterized in that, The step of obtaining the porosity plane distribution map of the sealing layer based on the second correspondence and the seismic facies plane distribution map includes: Based on the well logging curves, the porosity of each packer in the target area is determined; Based on the second correspondence, the porosity and porosity threshold of each sealing layer, a second wave impedance threshold is determined. The porosity threshold is used to characterize the maximum sealing layer porosity that the sealing layer needs to meet, and the second wave impedance threshold is used to characterize the wave impedance condition that the sealing layer needs to meet. Based on the second wave impedance threshold, the porosity threshold, and the seismic phase plane distribution map, the porosity plane distribution map of the sealing layer is obtained.
5. The method according to claim 4, characterized in that, The method further includes: Based on the third correspondence, the breakthrough pressure corresponding to the porosity of each sealing layer is determined. The third correspondence is the correspondence between the porosity of the sealing layer and the breakthrough pressure. The determination of the second wave impedance threshold based on the second correspondence, the porosity of each sealing layer, and the porosity threshold of each sealing layer includes: Based on the second correspondence, the porosity of each sealing layer, the porosity threshold, the fourth correspondence, the breakthrough pressure and the breakthrough pressure threshold of each sealing layer, the second wave impedance threshold is determined. The fourth correspondence is the correspondence between breakthrough pressure and wave impedance. The breakthrough pressure threshold is used to characterize the minimum breakthrough pressure that the sealing layer needs to meet.
6. The method according to claim 1, characterized in that, The evaluation results for determining lithological traps based on the seismic facies plane distribution, the reservoir porosity plane distribution map, and the sealing layer porosity plane distribution map include: Based on the preset boundary setting information of lithological traps, the seismic facies plane distribution, the reservoir porosity plane distribution map and the sealing layer porosity plane distribution map are overlaid to obtain a lithological trap plane distribution map. The lithological trap plane distribution map is used to characterize the distribution of multiple lithological traps in the target area. The evaluation result of each lithological trap is determined based on the ratio between the area of the target seismic facies and the target reservoir in each lithological trap and the area of the lithological trap.
7. An evaluation device for lithological traps, characterized in that, include: The first distribution map acquisition module is configured to acquire a seismic facies plane distribution map corresponding to the target area when the sedimentary facies type is a single-well sedimentary facies. The seismic facies plane distribution map is used to characterize the distribution of multiple seismic facies in the target area. The second distribution map acquisition module is configured to acquire a reservoir porosity plane distribution map based on a first correspondence and the seismic phase plane distribution map. The first correspondence is the correspondence between reservoir porosity and wave impedance. The reservoir porosity plane distribution map is used to characterize the distribution of the reservoir in each seismic phase. The third distribution map acquisition module is configured to acquire a sealing layer porosity plane distribution map based on the second correspondence and the seismic phase plane distribution map. The second correspondence is the correspondence between sealing layer porosity and wave impedance. The sealing layer porosity plane distribution map is used to characterize the distribution of the sealing layer in each seismic phase. The evaluation module is configured to determine the evaluation results of lithological traps based on the seismic facies plane distribution, the reservoir porosity plane distribution map, and the sealing layer porosity plane distribution 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 for evaluating lithological traps 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 evaluation method for lithological traps as described in any one of claims 1 to 6.
10. A computer program product, comprising a computer program, characterized in that, When executed by a processor, the computer program implements the steps of the evaluation method for lithological traps as described in any one of claims 1 to 6.