Method and device for evaluating geological storage of carbon dioxide based on euler number

Abstract

The application provides a carbon dioxide geological storage evaluation method and device based on Euler characteristic number, and the method comprises the following steps: obtaining three-dimensional pore grid data of a carbon dioxide geological storage reservoir; inputting the three-dimensional pore grid data into a pre-trained topological classification branch model to obtain the Euler characteristic number of pores; determining the pore topological category based on the Euler characteristic number and generating a topological normalization coefficient; obtaining local interface parameters according to the pore topological category; wherein the local interface parameters comprise pore deficiency curvature and the number of unit volume three-phase intersection lines; obtaining a dynamic wetting angle based on the topological normalization coefficient and the local interface parameters; and obtaining the carbon dioxide geological storage suitability evaluation result of the target carbon dioxide geological storage reservoir according to the dynamic wetting angle and the pore topological category. Through the technical scheme, the fine degree and the result accuracy of the carbon dioxide geological storage reservoir evaluation can be improved.

Description

Technical Field

[0001] This application relates to the field of carbon dioxide geological storage technology, and in particular to a method, apparatus, equipment and storage medium for evaluating carbon dioxide geological storage based on Euler characteristic number. Background Technology

[0002] Reservoir wettability is a key indicator for regulating the behavior of the carbon dioxide-brine-rock three-phase interface, the efficiency of carbon dioxide retention and plugging, and the safety of storage. Pore topology is the core factor determining wettability characteristics. Existing technologies for evaluating the wettability of carbon dioxide storage reservoirs and analyzing their suitability for storage only perform simple coupled calculations based on local interface and topological parameters such as wetting angle, defect curvature, and the number of three-phase intersections per unit volume. The topological representation is relatively one-sided, using only local parameters to describe the pore structure. It cannot distinguish between pores with different topological natures, such as simply connected, connected tubular, and network-like multi-connected pores, resulting in poor versatility of the evaluation models. Summary of the Invention

[0003] This application aims to at least partially address one of the technical problems in the related art.

[0004] In a first aspect, this application proposes a method for evaluating carbon dioxide geological storage based on Euler characteristic numbers. The method includes: acquiring three-dimensional pore mesh data of a carbon dioxide geological storage layer; inputting the three-dimensional pore mesh data into a pre-trained topology classification branch model to obtain the Euler characteristic numbers of the pores; determining the pore topology category based on the Euler characteristic numbers and generating topology normalization coefficients; obtaining local interface parameters according to the pore topology category; wherein the local interface parameters include pore deficiency curvature and the number of three-phase intersections per unit volume; obtaining a dynamic wetting angle based on the topology normalization coefficients and the local interface parameters; and obtaining a carbon dioxide geological storage suitability evaluation result for the target carbon dioxide geological storage layer according to the dynamic wetting angle and the pore topology category.

[0005] In one implementation, determining the pore topology category based on the Euler characteristic number and generating a topology normalization coefficient includes: establishing a mapping relationship between different Euler characteristic number numerical ranges and different preset pore topology types; obtaining the pore topology category corresponding to the Euler characteristic number based on the mapping relationship; and obtaining the corresponding topology normalization coefficient based on the topology benchmark threshold and reservoir geological background parameters corresponding to the pore topology category.

[0006] In one implementation, obtaining local interface parameters based on the pore topology category includes: establishing a three-dimensional pore mesh model of the carbon dioxide geological sealing reservoir; based on the three-dimensional pore mesh model, and combining the pore topology category to partition and locate the pore wall, pore boundary and fluid contact interface region, so as to divide the carbon dioxide geological sealing reservoir into multiple calculation partitions; performing surface geometry operations on each calculation partition to obtain the corresponding pore deficit curvature; and based on the pore deficit curvature, obtaining the number of gas-liquid-solid three-phase intersection lines per unit reservoir volume.

[0007] In one implementation, obtaining the dynamic wetting angle based on the topology normalization coefficient and the local interface parameters includes: constructing a fitting model based on the topology normalization coefficient and the topology wetting coupling relationship; substituting the local interface parameters into the fitting model to calculate the dynamic wetting angle.

[0008] In one implementation, the formula for calculating the topology normalization coefficient is:

[0009] ; In the formula, These are the topological normalization coefficients. It is an Euler characteristic.

[0010] In one alternative implementation, the formula for calculating the porosity defect curvature is: ; In the formula: For the pore defect curvature, The ideal equivalent radius of curvature of the pores. The normal height of the pore wall. The tangential length of the pore wall; the formula for calculating the number of three-phase intersection lines per unit volume is: ; In the formula, The number of three-phase intersection lines per unit volume; The total length of the three-phase intersection line within the pore. This represents the pore volume.

[0011] Secondly, this application proposes a carbon dioxide geological storage assessment device based on Euler characteristic number, the device comprising: an acquisition module for acquiring three-dimensional pore grid data of carbon dioxide geological storage layer; The first processing module is used to input the three-dimensional pore mesh data into a pre-trained topology classification branch model to obtain the Euler characteristic number of the pores; the second processing module is used to determine the pore topology category based on the Euler characteristic number and generate topology normalization coefficients; the third processing module is used to obtain local interface parameters according to the pore topology category; wherein, the local interface parameters include pore defect curvature and the number of three-phase intersections per unit volume; the fourth processing module is used to obtain the dynamic wetting angle based on the topology normalization coefficients and the local interface parameters; the fifth processing module is used to obtain the carbon dioxide geological storage suitability assessment result of the target carbon dioxide geological storage layer according to the dynamic wetting angle and the pore topology category.

[0012] In one implementation, the second processing module can be used to: establish a mapping relationship between different Euler characteristic number numerical ranges and different preset pore topology types; obtain the pore topology category corresponding to the Euler characteristic number based on the mapping relationship; and obtain the corresponding topology normalization coefficient based on the topology benchmark threshold and reservoir geological background parameters corresponding to the pore topology category.

[0013] In one implementation, the third processing module can be used to: establish a three-dimensional pore mesh model of the carbon dioxide geological sealing reservoir; based on the three-dimensional pore mesh model, and combined with the pore topology partitioning, locate the pore walls, pore boundaries and fluid contact interface areas to divide the carbon dioxide geological sealing reservoir into multiple computational partitions; perform surface geometry operations on each computational partition to obtain the corresponding pore deficit curvature; and based on the pore deficit curvature, obtain the number of gas-liquid-solid three-phase intersection lines per unit reservoir volume.

[0014] In one implementation, the fourth processing module can be used to: construct a fitting model based on the topology normalization coefficient and the topology wetting coupling relationship; substitute the local interface parameters into the fitting model to calculate the dynamic wetting angle.

[0015] In one implementation, the formula for calculating the topology normalization coefficient is:

[0016] In the formula, These are the topological normalization coefficients. It is an Euler characteristic.

[0017] In one alternative implementation, the formula for calculating the porosity defect curvature is:

[0018] In the formula: For the pore defect curvature, The ideal equivalent radius of curvature of the pores. The normal height of the pore wall. The tangential length of the pore wall; the formula for calculating the number of three-phase intersection lines per unit volume is:

[0019] In the formula, The number of three-phase intersection lines per unit volume; The total length of the three-phase intersection line within the pore. This represents the pore volume.

[0020] Thirdly, this application provides an electronic device comprising: at least one processor; and a memory communicatively connected to the at least one processor; wherein the memory stores instructions executable by the at least one processor, the instructions being executed by the at least one processor to enable the at least one processor to perform the method as described in the first aspect.

[0021] Fourthly, this application proposes a storage medium storing instructions that, when executed on an electronic device, cause the electronic device to perform the method described in the first aspect.

[0022] Fifthly, this application proposes a program product comprising at least one of a program and instructions, wherein when the program and instructions are executed by an electronic device, they implement the steps of the method described in the first aspect.

[0023] The method, apparatus, equipment, and storage medium for assessing carbon dioxide geological storage based on Euler characteristic numbers provided in this application can acquire three-dimensional pore mesh data of reservoirs, calculate pore Euler characteristic numbers using a pre-trained model, classify pore topology categories based on Euler characteristic numbers, and construct topology normalization coefficients. Then, by combining topological partitioning, local interface parameters such as pore deficiency curvature and the number of three-phase intersections per unit volume are extracted. Dynamic wetting angles are inverted through coupling topology coefficients and interface parameters, and finally, a comprehensive assessment of reservoir CO2 storage suitability is completed by combining the wetting angle and pore topology characteristics. This can improve the refinement and accuracy of carbon dioxide geological storage reservoir assessment.

[0024] Additional aspects and advantages of this application will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of this application. Attached Figure Description

[0025] The above and / or additional aspects and advantages of this application will become apparent and readily understood from the following description of the embodiments taken in conjunction with the accompanying drawings, wherein: Figure 1 This is a flowchart illustrating a method for evaluating carbon dioxide geological sequestration based on Euler characteristic numbers, provided in an embodiment of this application. Figure 2 This is a schematic diagram of the structure of a carbon dioxide geological sequestration assessment device based on Euler characteristic numbers provided in an embodiment of this application; Figure 3 This is a schematic diagram of the structure of the electronic device provided in the embodiments of this application. Detailed Implementation

[0026] The embodiments of this application are described in detail below. Examples of these embodiments are shown in the accompanying drawings, wherein the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The embodiments described below with reference to the accompanying drawings are exemplary and intended to explain this application, and should not be construed as limiting this application.

[0027] The following describes, with reference to the accompanying drawings, a method and apparatus for evaluating carbon dioxide geological sequestration based on Euler characteristic numbers, according to embodiments of this application.

[0028] Figure 1 This is a schematic flowchart illustrating a carbon dioxide geological sequestration assessment method based on Euler characteristic numbers provided in an embodiment of this application. Figure 1 As shown, the method may include, but is not limited to, the following steps: S1. Obtain three-dimensional pore mesh data of the carbon dioxide geological sealing reservoir.

[0029] For example, numerical modeling or geological 3D mesh generation techniques can be used to establish 3D pore mesh data of carbon dioxide geological storage layers.

[0030] S2. Input the three-dimensional pore mesh data into the pre-trained topology classification branch model to obtain the Euler characteristic number of the pores.

[0031] For example, after preprocessing the three-dimensional porous mesh data such as normalization, noise reduction, and mesh clipping, it is input into a pre-trained deep learning model for topological classification to obtain the Euler characteristic number topological index of the three-dimensional porous medium.

[0032] As an example, the formula for calculating the Euler characteristic can be expressed as:

[0033] In the formula, Here is the Euler characteristic number for porosity. The number of vertices in the porous mesh. This represents the number of edges in the pore mesh. The number of pore mesh surfaces. This represents the number of pore mesh cells.

[0034] S3. Determine the pore topology category based on the Euler characteristic number and generate the topology normalization coefficient.

[0035] For example, for closed curved face gaps, the Euler characteristic can be directly related to the topological genus, and the aforementioned formula for calculating the Euler characteristic can be simplified to:

[0036] In the formula, It is a topological genus. Thus, the correlation between the Euler characteristic and the porosity topological category can be simplified using this formula.

[0037] At this point, the logic for determining the pore topology category based on the Euler characteristic number is as follows: When Euler characteristic > 0 There are no connected pores, the pore structure is regular, and the pore topology is determined to be a simple topological pore with quasi-spherical single connectivity; When the Euler characteristic is 0 The pores are single-channel connected and tubular or annular, and the topology of the pores is determined to be a medium-topological pore with cylindrical / annular connectivity. When the Euler characteristic is less than 0 The pore structure is extremely complex, with multiple channels interwoven, and the pore topology is determined to be a complex topological pore with multiple interconnections in a mesh.

[0038] In one implementation, the pore topology category is determined based on the Euler characteristic number, and a topology normalization coefficient is generated. This includes: establishing a mapping relationship between different Euler characteristic number numerical ranges and different preset pore topology types; obtaining the pore topology category corresponding to the Euler characteristic number based on the mapping relationship; and obtaining the corresponding topology normalization coefficient based on the topology benchmark threshold and reservoir geological background parameters corresponding to the pore topology category.

[0039] For example, a one-to-one mapping relationship is established between different Euler characteristic number ranges and various preset pore topology types. This mapping relationship is used to match and determine the range to which each pore Euler characteristic number belongs in the aforementioned steps, thereby determining the pore topology category of the target reservoir. Then, the benchmark thresholds corresponding to various pore topology types are combined with the actual geological background parameters of the reservoir for correction to obtain the topology normalization coefficient.

[0040] In one alternative implementation, the topology normalization coefficient is calculated as follows: ; In the formula, These are the topological normalization coefficients. It is an Euler characteristic.

[0041] It should be noted that this application can achieve accurate classification of the pore topology by introducing Euler characteristic as a global topological invariant, solve the problem of one-sided local parameter representation, and distinguish the differences in wetting behavior of pores with different topological structures.

[0042] S4. Obtain local interface parameters based on the pore topology category.

[0043] In the embodiments of this application, the local interface parameters include the pore defect curvature and the number of three-phase intersections per unit volume.

[0044] In one implementation, local interface parameters are obtained based on the pore topology category, including: establishing a three-dimensional pore mesh model of the carbon dioxide geological sealing reservoir; based on the three-dimensional pore mesh model, and combined with the pore topology category, locating the pore wall, pore boundary and fluid contact interface region to divide the carbon dioxide geological sealing reservoir into multiple computational partitions; performing surface geometry operations on each computational partition to obtain the corresponding pore deficit curvature; and based on the pore deficit curvature, obtaining the number of gas-liquid-solid three-phase intersection lines per unit reservoir volume.

[0045] For example, a three-dimensional pore grid model of a carbon dioxide geological reservoir is constructed. Based on the pre-defined pore topology categories, the pore walls, pore boundaries, and fluid contact interfaces are located, thereby dividing the reservoir into multiple zones. Subsequently, surface geometry operations are performed on each zone to obtain the corresponding pore deficit curvature. By combining the inherent correlation laws of differential topology geometry of porous media, a quantitative mapping relationship between pore deficit curvature and the number of gas-liquid-solid three-phase intersection lines is established. The curvature characteristics of all zones are integrated and normalized and homogenized by combining the volume parameters of the reservoir grid unit. The number of gas-liquid-solid three-phase intersection lines per unit reservoir volume is then statistically obtained.

[0046] In one alternative implementation, the formula for calculating the porosity defect curvature is:

[0047] In the formula: For the pore defect curvature, The ideal equivalent radius of curvature of the pores. The normal height of the pore wall. The tangential length of the pore wall; The formula for calculating the number of three-phase intersection lines per unit volume is:

[0048] In the formula, The number of three-phase intersection lines per unit volume; The total length of the three-phase intersection line within the pore. This represents the pore volume.

[0049] S5. Obtain the dynamic wetting angle based on the topology normalization coefficient and local interface parameters.

[0050] In one implementation, the dynamic wetting angle is obtained based on the topology normalization coefficient and local interface parameters, including: constructing a fitting model based on the topology normalization coefficient and the topology wetting coupling relationship; substituting the local interface parameters into the fitting model to calculate the dynamic wetting angle.

[0051] For example, based on the topological normalization coefficient, a fitting calculation model is constructed in combination with the topological wetting coupling mechanism. Then, the extracted local interface parameters such as the pore defect curvature and the number of three-phase intersections per unit volume are uniformly substituted into the fitting model. Through parameter coupling calculation and model solving, the dynamic wetting angle value corresponding to the reservoir pores is calculated.

[0052] The formula for calculating the dynamic wetting angle is:

[0053] In the formula, For dynamic wetting angle, The function is a generalized function of the number of three-phase intersection lines per unit volume, which adaptively switches according to the pore topology type (linear / logarithmic / power function). All are global average coupling coefficients. This is a global correction item. The topological normalization coefficient can be calculated using the following formula:

[0054] S6. Based on the dynamic wetting angle and pore topology category, obtain the suitability assessment results for carbon dioxide geological storage of the target carbon dioxide geological storage layer.

[0055] In one implementation, obtaining the suitability assessment result for carbon dioxide geological storage of a target carbon dioxide geological storage reservoir based on the dynamic wetting angle and pore topology category may include the following steps: determining the reservoir hydrophilicity level based on the dynamic wetting angle; and obtaining the suitability assessment result for carbon dioxide geological storage based on the hydrophilicity level.

[0056] Understandably, the stronger the hydrophilicity, the better the CO2 capillary retention and blocking ability, and the higher the corresponding storage suitability level; conversely, the lower the suitability.

[0057] In one implementation, the hydrophilicity level includes, from high to low, strong hydrophilicity, weak hydrophilicity, and hydrophobicity. The determination process for the reservoir hydrophilicity level based on the dynamic wetting angle includes: when the dynamic wetting angle < a first threshold, the hydrophilicity level is determined to be strong hydrophilic; when the first threshold ≤ the dynamic wetting angle ≤ a second threshold, the hydrophilicity level is determined to be weak hydrophilic; when the dynamic wetting angle > the second threshold, the hydrophilicity level is determined to be hydrophobic; when the target carbon dioxide geological sealing reservoir has complex topological pores and exhibits strong hydrophilicity, the wettability is determined. The evaluation result is excellent; when the target carbon dioxide geological seal reservoir has medium topological porosity and is strongly or weakly hydrophilic, or when the target carbon dioxide geological seal reservoir has complex topological porosity and is weakly hydrophilic, the wettability evaluation result is good to above average; when the target carbon dioxide geological seal reservoir has simple topological porosity and is strongly hydrophilic, or when the target carbon dioxide geological seal reservoir has medium topological porosity and is weakly hydrophilic, the wettability evaluation result is good to below average; in the remaining cases, the wettability evaluation result is poor.

[0058] By implementing the embodiments of this application, three-dimensional pore mesh data of the reservoir can be obtained. The Euler characteristic number of the pores can be calculated using a pre-trained model, and pore topology categories can be classified based on the Euler characteristic number to construct topology normalization coefficients. Then, local interface parameters such as pore deficiency curvature and the number of three-phase intersections per unit volume can be extracted by combining topology partitioning. Dynamic wetting angles can be inverted through coupling topology coefficients and interface parameters. Finally, a comprehensive assessment of the reservoir's CO2 storage suitability can be completed by combining the wetting angle and pore topology characteristics. This can improve the refinement and accuracy of CO2 geological storage reservoir assessments.

[0059] Please see Figure 2 , Figure 2 This is a schematic diagram of a carbon dioxide geological sequestration assessment device based on Euler characteristic numbers provided in an embodiment of this application. Figure 2 As shown, the device 200 includes: an acquisition module 201 for acquiring three-dimensional pore mesh data of a carbon dioxide geological storage layer; a first processing module 202 for inputting the three-dimensional pore mesh data into a pre-trained topology classification branch model to acquire the Euler characteristic number of the pores; a second processing module 203 for determining the pore topology category based on the Euler characteristic number and generating topology normalization coefficients; a third processing module 204 for acquiring local interface parameters according to the pore topology category, wherein the local interface parameters include pore defect curvature and the number of three-phase intersections per unit volume; a fourth processing module 205 for acquiring the dynamic wetting angle based on the topology normalization coefficients and local interface parameters; and a fifth processing module 206 for acquiring the carbon dioxide geological storage suitability assessment result of the target carbon dioxide geological storage layer according to the dynamic wetting angle and the pore topology category.

[0060] In one implementation, the second processing module 203 can be used to: establish a mapping relationship between different Euler characteristic number numerical ranges and different preset pore topology types; obtain the pore topology category corresponding to the Euler characteristic number based on the mapping relationship; and obtain the corresponding topology normalization coefficient based on the topology benchmark threshold and reservoir geological background parameters corresponding to the pore topology category.

[0061] In one implementation, the third processing module 204 can be used to: establish a three-dimensional pore mesh model of the carbon dioxide geological sealing reservoir; based on the three-dimensional pore mesh model, and combined with pore topology partitioning, locate the pore wall, pore boundary and fluid contact interface regions to divide the carbon dioxide geological sealing reservoir into multiple calculation partitions; perform surface geometry operations on each calculation partition to obtain the corresponding pore deficit curvature; and based on the pore deficit curvature, obtain the number of gas-liquid-solid three-phase intersection lines per unit reservoir volume.

[0062] In one implementation, the fourth processing module 205 can be used to: construct a fitting model based on the topological normalization coefficient and the topological wetting coupling relationship; substitute the local interface parameters into the fitting model to calculate the dynamic wetting angle.

[0063] In one implementation, the topological normalization coefficient is calculated as follows:

[0064] In the formula, These are the topological normalization coefficients. It is an Euler characteristic.

[0065] In one alternative implementation, the formula for calculating the porosity defect curvature is:

[0066] In the formula: For the pore defect curvature, The ideal equivalent radius of curvature of the pores. The normal height of the pore wall. The tangential length of the pore wall; the formula for calculating the number of three-phase intersection lines per unit volume is:

[0067] In the formula, The number of three-phase intersection lines per unit volume; The total length of the three-phase intersection line within the pore. This represents the pore volume.

[0068] It should be noted that the foregoing explanation of the embodiment of the carbon dioxide geological storage assessment method based on Euler characteristic number also applies to the carbon dioxide geological storage assessment device based on Euler characteristic number in this embodiment, and will not be repeated here.

[0069] To implement the above embodiments, this application also proposes an electronic device. Please see [link to relevant documentation]. Figure 3 , Figure 3 This is a schematic diagram of the structure of the electronic device provided in an embodiment of this application. For example... Figure 3 As shown, the electronic device 300 includes: a processor 301 and a memory 302 communicatively connected to the processor 301; the memory 302 stores computer execution instructions; the processor 301 executes the computer execution instructions stored in the memory to implement the method provided in the foregoing embodiments.

[0070] To implement the above embodiments, this application also proposes a storage medium storing instructions that, when executed on an electronic device, cause the electronic device to perform the methods provided in the foregoing embodiments.

[0071] To implement the above embodiments, this application also proposes a program product, including at least one of a program and instructions, wherein when the program and instructions are executed by an electronic device, they implement the steps of the method provided in the foregoing embodiments.

[0072] It should be noted that the acquisition, transmission, storage, use, and processing of data in this application comply with the relevant provisions of national laws and regulations and do not violate public order and good morals.

[0073] It should be noted that the information (including but not limited to user device information, user personal information, etc.), data (including but not limited to data used for analysis, data stored, data displayed, etc.) and signals involved in this application are all authorized by the user or fully authorized by all parties, and the collection, use and processing of related data must comply with the relevant laws, regulations and standards of the relevant countries and regions.

[0074] It is worth noting that in the embodiments of this application, certain software, components, models and other existing solutions in the industry may be mentioned. These should be regarded as exemplary and are only intended to illustrate the feasibility of implementing the technical solution of this application. However, they do not mean that the applicant has used or necessarily used the solution.

[0075] In the description of this application, unless otherwise stated, " / " means "or", for example, A / B can mean A or B; "and / or" in this document is merely a description of the relationship between related objects, indicating that there can be three relationships, for example, A and / or B can mean: A exists alone, A and B exist simultaneously, and B exists alone.

[0076] In the foregoing descriptions of the embodiments, 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 this application. 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 may 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, as well as the features of different embodiments or examples.

[0077] Furthermore, the terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one of that feature. In the description of this application, "multiple" means at least two, such as two, three, etc., unless otherwise explicitly specified.

[0078] Any process or method description in the flowchart or otherwise herein can be understood as representing a module, segment, or portion of code comprising one or more executable instructions for implementing custom logic functions or processes, and the scope of the preferred embodiments of this application includes additional implementations in which functions may be performed not in the order shown or discussed, including substantially simultaneously or in reverse order depending on the functions involved, as should be understood by those skilled in the art to which embodiments of this application pertain.

[0079] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as a sequenced list of executable instructions for implementing logical functions, and can be embodied in any computer-readable medium for use by, or in conjunction with, an instruction execution system, apparatus, or device (such as a computer-based system, a processor-included system, or other system that can fetch and execute instructions from, an instruction execution system, apparatus, or device). For the purposes of this specification, "computer-readable medium" can be any means that can contain, store, communicate, propagate, or transmit programs for use by, or in conjunction with, an instruction execution system, apparatus, or device. More specific examples (a non-exhaustive list) of computer-readable media include: an electrical connection having one or more wires (electronic device), a portable computer disk drive (magnetic device), random access memory (RAM), read-only memory (ROM), erasable and editable read-only memory (EPROM or flash memory), fiber optic devices, and portable optical disc read-only memory (CDROM). Alternatively, the computer-readable medium may be paper or other suitable media on which the program can be printed, since the program can be obtained electronically, for example, by optically scanning the paper or other medium, followed by editing, interpreting, or otherwise processing as necessary, and then stored in a computer memory.

[0080] It should be understood that various parts of this application can be implemented using hardware, software, firmware, or a combination thereof. In the above embodiments, multiple steps or methods can be implemented using software or firmware stored in memory and executed by a suitable instruction execution system. For example, if implemented in hardware as in another embodiment, it can be implemented using any one or a combination of the following techniques known in the art: discrete logic circuits having logic gates for implementing logical functions on data signals, application-specific integrated circuits (ASICs) having suitable combinational logic gates, programmable gate arrays (PGAs), field-programmable gate arrays (FPGAs), etc.

[0081] Those skilled in the art will understand that all or part of the steps of the methods in the above embodiments can be implemented by a program instructing related hardware. The program can be stored in a computer-readable storage medium, and when executed, the program includes one or a combination of the steps of the method embodiments.

[0082] Furthermore, the functional units in the various embodiments of this application can be integrated into a processing module, or each unit can exist physically separately, or two or more units can be integrated into a module. The integrated module can be implemented in hardware or as a software functional module. If the integrated module is implemented as a software functional module and sold or used as an independent product, it can also be stored in a computer-readable storage medium.

[0083] The storage medium mentioned above can be a read-only memory, a disk, or an optical disk, etc. Although embodiments of this application have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting this application. Those skilled in the art can make changes, modifications, substitutions, and variations to the above embodiments within the scope of this application.

Claims

1. A method for assessing carbon dioxide geological sequestration based on Euler characteristic numbers, characterized in that, include: Obtain three-dimensional pore mesh data of carbon dioxide geological storage layers; The three-dimensional pore mesh data is input into a pre-trained topological classification branch model to obtain the Euler characteristic number of the pores; The pore topology category is determined based on the Euler characteristic number, and topology normalization coefficients are generated. Based on the pore topology category, local interface parameters are obtained; wherein, the local interface parameters include pore defect curvature and the number of three-phase intersections per unit volume; The dynamic wetting angle is obtained based on the topology normalization coefficient and the local interface parameters; Based on the dynamic wetting angle and the pore topology category, the suitability assessment results for carbon dioxide geological storage of the target carbon dioxide geological storage layer are obtained.

2. The method according to claim 1, characterized in that, The step of determining the pore topology category based on the Euler characteristic number and generating topology normalization coefficients includes: Establish a mapping relationship between different ranges of Euler characteristic numbers and different preset pore topology types; Based on the mapping relationship, the porosity topology category corresponding to the Euler characteristic number is obtained; Based on the topological benchmark threshold corresponding to the aforementioned pore topology category and the reservoir geological background parameters, the corresponding topological normalization coefficient is obtained.

3. The method according to claim 1, characterized in that, The step of obtaining local interface parameters based on the pore topology category includes: Establish a three-dimensional pore mesh model of the carbon dioxide geological sealing storage layer; Based on the three-dimensional pore mesh model, and combined with the pore topology partitioning, the pore walls, pore boundaries and fluid contact interface areas are located to divide the carbon dioxide geological sealing storage layer into multiple computational partitions. Perform surface geometry operations on each computational partition to obtain the corresponding porosity deficiency curvature; Based on the pore defect curvature, the number of gas-liquid-solid three-phase intersection lines per unit reservoir volume is obtained.

4. The method according to claim 1, characterized in that, The process of obtaining the dynamic wetting angle based on the topology normalization coefficient and the local interface parameters includes: Based on the aforementioned topology normalization coefficients, a fitting model is constructed by combining the topology wetting coupling relationship. The dynamic wetting angle is calculated by substituting the local interface parameters into the fitting model.

5. The method according to claim 1, characterized in that, The formula for calculating the topology normalization coefficient is: In the formula, These are the topological normalization coefficients. It is an Euler characteristic.

6. The method according to claim 3, characterized in that, The formula for calculating the pore defect curvature is: In the formula: For the pore defect curvature, The ideal equivalent radius of curvature of the pores. The normal height of the pore wall. The tangential length of the pore wall; The formula for calculating the number of three-phase intersection lines per unit volume is: In the formula, The number of three-phase intersection lines per unit volume; The total length of the three-phase intersection line within the pore. This represents the pore volume.

7. A carbon dioxide geological sequestration assessment device based on Euler characteristic, characterized in that, include: The acquisition module is used to acquire three-dimensional pore mesh data of carbon dioxide geological sealing reservoirs; The first processing module is used to input the three-dimensional pore mesh data into a pre-trained topological classification branch model to obtain the Euler characteristic number of the pores; The second processing module is used to determine the pore topology category based on the Euler characteristic number and generate topology normalization coefficients; The third processing module is used to obtain local interface parameters based on the pore topology category; wherein, the local interface parameters include pore defect curvature and the number of three-phase intersections per unit volume. The fourth processing module is used to obtain the dynamic wetting angle based on the topology normalization coefficient and the local interface parameters; The fifth processing module is used to obtain the suitability assessment results for carbon dioxide geological storage of the target carbon dioxide geological storage layer based on the dynamic wetting angle and the pore topology category.

8. An electronic device, characterized in that, include: A processor, and a memory communicatively connected to the processor; The memory stores computer-executed instructions; The processor executes computer execution instructions stored in the memory to implement the method as described in any one of claims 1 to 6.

9. A storage medium storing instructions, characterized in that, When the instructions are executed on an electronic device, the electronic device causes the electronic device to perform the method of any one of claims 1 to 6.

10. A program product comprising at least one of a program and instructions, characterized in that, When at least one of the program or instructions is executed by an electronic device, it implements the steps of the method according to any one of claims 1 to 6.

Images