Microseismic-based CCUS / CCS geologic body integrity evaluation method and device and electronic equipment

Abstract

The embodiment of the invention provides a micro-seismic-based CCUS / CCS geologic body integrity evaluation method and device and electronic equipment. The microseism-based CCUS / CCS geologic body integrity evaluation method comprises the following steps: responding to a plurality of seismic signal detectors arranged in a preset range of a target CCUS / CCS geologic body, detecting data information of microseism signals from the target CCUS / CCS geologic body, and according to the data information of the microseism signals, evaluating the integrity of the microseism-based CCUS / CCS geologic body according to the data information of the microseism signals; obtaining a micro-seismic event attribute corresponding to the micro-seismic signal and pore pressure of a stratum generating the micro-seismic signal; obtaining a geologic body integrity evaluation safety threshold value of the target CCUS / CCS geologic body according to the microseismic event attribute and the pore pressure; and evaluating the geologic body integrity of the target CCUS / CCS geologic body according to the geologic body integrity evaluation threshold value.

Description

Technical Field

[0001] The present application relates to the technical field of petroleum geophysical exploration, and particularly relates to a method, device and electronic device for evaluating the integrity of CCUS / CCS geological bodies based on microseismicity. Background Art

[0002] The problem of global warming is becoming increasingly serious, and the international community's demand for greenhouse gas emissions reduction is becoming increasingly urgent. Carbon dioxide (CO2) emissions are one of the main targets of greenhouse gas emissions reduction. During the process of oil extraction, a large amount of carbon dioxide emissions are generated, which not only have a negative impact on the environment, but also restrict the sustainable development of oilfields. Carbon dioxide can reach the supercritical state under certain temperature and pressure conditions. Carbon dioxide in the supercritical state can be compressed like a gas, which helps the carbon dioxide to diffuse and migrate in the reservoir, and has a large density like a liquid, which is conducive to storing more carbon dioxide in the same space, realizing geological storage and oil displacement of carbon dioxide. Therefore, the carbon capture, utilization and storage technology of carbon dioxide can be used to reduce carbon dioxide emissions, and at the same time, carbon dioxide can be converted into valuable resources for utilization. The carbon capture, utilization and storage technology of carbon dioxide includes CCS (Carbon Capture and Storage), CCUS (Carbon Capture, Utilization and Storage), etc.

[0003] In the related art, during the process of geological storage, migration and oil displacement of carbon dioxide underground, there are problems such as the gas drive front, the injection sweep range, the dominant injection direction, fracture zone crossflow, unclear understanding of the injection rate reduction and the injection sweep range, which directly affect the effectiveness of the regulation and control of injection-production well patterns and injection-production parameters; the activation of caprock and faults will lead to the failure of the integrity of geological bodies.

[0004] Therefore, how to monitor the integrity of the geological body storing carbon dioxide using CCUS / CCS technology in real time and evaluate the risk of carbon dioxide leakage is a technical problem to be solved urgently. Summary of the Invention

[0005] The embodiments of the present application provide a method, device and electronic device for evaluating the integrity of CCUS / CCS geological bodies based on microseismicity, which are used to monitor the integrity of the geological body storing carbon dioxide using CCUS / CCS technology in real time and evaluate the risk of carbon dioxide leakage.

[0006] One embodiment of the present application provides a method for evaluating the geological body integrity of CCUS / CCS based on microseismicity. The method includes: in response to a plurality of seismic signal detectors arranged within a preset range of a target CCUS / CCS geological body, detecting data information of microseismic signals from the target CCUS / CCS geological body, and based on the data information of the microseismic signals, obtaining the microseismic event attributes corresponding to the microseismic signals and the pore pressure of the formation where the microseismic signals are generated; obtaining a safety valve threshold for evaluating the geological body integrity of the target CCUS / CCS geological body according to the microseismic event attributes and the pore pressure; and evaluating the geological body integrity of the target CCUS / CCS geological body according to the safety valve threshold for evaluating the geological body integrity.

[0007] In some embodiments, the microseismic event attributes include the microseismic event magnitude; the obtaining of the microseismic event attributes corresponding to the microseismic signals according to the data information of the microseismic signals includes: calculating the value of the microseismic event magnitude M using the following formula w :

[0008]

[0009] where v p is the longitudinal wave propagation velocity of the microseismic signal; Ω0 is the zero-frequency limit value of the source spectrum; F is the average propagation coefficient of the seismic wave; d is the average value of the distances from the source of the microseismic signal to each seismic signal detector, which is obtained using a multi-wave joint location algorithm based on the data information of the microseismic signals.

[0010] In some embodiments, the microseismic event attributes further include the microseismic event energy; the obtaining of the microseismic event attributes corresponding to the microseismic signals according to the data information of the microseismic signals includes:

[0011] calculating the value of the microseismic event energy E using the following formula nergy :

[0012]

[0013] where ρ is the density of the medium of the formation where the microseismic signal is generated; t is the average value of the first arrival times of the microseismic signal picked up by each seismic signal detector; W rms,i is the root mean square amplitude of the microseismic signal detected by the i-th seismic signal detector; K is the number of seismic signal detectors that detect the microseismic signal.

[0014] In some embodiments, the microseismic event attribute further includes the signal-to-noise ratio of the microseismic event; obtaining the microseismic event attribute corresponding to the microseismic signal according to the data information of the microseismic signal includes:

[0015] Using the following formula to calculate the signal-to-noise ratio R SN of the microseismic event:

[0016]

[0017] where RMS signal is the root mean square amplitude of the effective signal detected by the seismic signal detector; RMS noise is the root mean square amplitude of the noise detected by the seismic signal detector; wherein, RMS signal is obtained using the following formula:

[0018]

[0019] where N is the number of seismic signal detectors detecting the microseismic signal; A i is the amplitude value of the microseismic signal detected by the i-th seismic signal detector.

[0020] In some embodiments, obtaining the pore pressure of the formation that generates the microseismic signal according to the data information of the microseismic signal includes: using the following formula to calculate the pore pressure P p :

[0021]

[0022] where P o is the overburden pressure of the target CCUS / CCS geological body, calculated by density logging; h is the vertical depth of the formation that generates the microseismic signal, obtained using the multi-wave joint positioning algorithm according to the data information of the microseismic signal; P h is the hydrostatic pressure of the formation that generates the microseismic signal; DT n is the acoustic wave travel time of mudstone under normal pressure; DT is the measured acoustic wave travel time of the formation by logging; b is the Eaton index; ρ is the medium density of the formation space where the microseismic signal is generated.

[0023] In some embodiments, obtaining the geological body integrity evaluation safety valve threshold of the target CCUS / CCS geological body according to the microseismic event attributes and the pore pressure includes: respectively normalizing the microseismic event magnitude, the microseismic event energy, the microseismic event signal-to-noise ratio, and the pore pressure; according to the normalized microseismic event magnitude, the microseismic event energy, the microseismic event signal-to-noise ratio, and the pore pressure, using the weighted analysis function shown below, calculating to obtain the geological body integrity evaluation safety valve threshold S ccus / ccs :

[0024] S ccus / ccs = w1M w1 + w2E nergy2 + w3R SN3 + w4P p4

[0025] where w1, w2, w3, and w4 are the weight coefficients of the microseismic event magnitude, the microseismic event energy, the microseismic event signal-to-noise ratio, and the pore pressure respectively; M w1 、E nergy2 、R SN3 、P p4 are the values after normalization of the microseismic event magnitude, the values after normalization of the microseismic event energy, the values after normalization of the microseismic event signal-to-noise ratio, and the values after normalization of the pore pressure respectively.

[0026] In some embodiments, evaluating the geological body integrity of the target CCUS / CCS geological body according to the geological body integrity evaluation safety valve threshold includes: in response to the geological body integrity evaluation safety valve threshold being within the first preset threshold range, determining that the carbon dioxide leakage risk of the target CCUS / CCS geological body is a high risk level; in response to the geological body integrity evaluation safety valve threshold being within the second preset threshold range, determining that the carbon dioxide leakage risk of the target CCUS / CCS geological body is a medium risk level; in response to the geological body integrity evaluation safety valve threshold being within the third preset threshold range, determining that the carbon dioxide leakage risk of the target CCUS / CCS geological body is a low risk level.

[0027] One embodiment of the present application provides a device for evaluating the integrity of a CCUS / CCS geological body based on microseismicity. The device includes: a first acquisition module configured to, in response to a plurality of seismic signal detectors disposed within a preset range of a target CCUS / CCS geological body, detect data information of microseismic signals from the target CCUS / CCS geological body, and obtain microseismic event attributes corresponding to the microseismic signals and the pore pressure of the formation that generates the microseismic signals according to the data information of the microseismic signals; a second acquisition module configured to obtain a safety valve threshold for evaluating the integrity of the target CCUS / CCS geological body according to the microseismic event attributes and the pore pressure; and an evaluation module configured to evaluate the integrity of the target CCUS / CCS geological body according to the safety valve threshold for evaluating the integrity of the geological body.

[0028] An embodiment of the present application provides an electronic device. The electronic device includes a memory and a processor. The memory stores a computer program, and the processor executes the method described above when running the program.

[0029] An embodiment of the present application provides a storage medium for storing a computer-readable program. When the computer-readable program is run, it executes the method described above.

[0030] The above technical solutions provided by the embodiments of the present application have at least the following advantages compared with the prior art:

[0031] In the embodiments provided by the present application, according to the data information of the microseismic signals, microseismic event attributes corresponding to the microseismic signals and the pore pressure of the formation that generates the microseismic signals are obtained; according to the microseismic event attributes and the pore pressure, a safety valve threshold for evaluating the integrity of the target CCUS / CCS geological body is obtained; and according to the safety valve threshold for evaluating the integrity of the geological body, the integrity of the target CCUS / CCS geological body is evaluated. Thus, the integrity of the geological body for storing carbon dioxide using the CCUS / CCS technology can be monitored in real time, and the risk of carbon dioxide leakage can be evaluated. BRIEF DESCRIPTION OF THE DRAWINGS

[0032] The present application will be further described in an exemplary embodiment, and these exemplary embodiments will be described in detail through the drawings. These embodiments are not restrictive. In these embodiments, the same numbers represent the same structures, where:

[0033] Figure 1 is an exemplary flowchart of a method for evaluating the integrity of a CCUS / CCS geological body based on microseismicity according to some embodiments of the present application;

[0034] Figure 2 is an exemplary diagram of a microseismic signal according to some embodiments of the present application;

[0035] Figure 3 is an exemplary schematic diagram of the forward microseismic event location result shown in some embodiments of the present application;

[0036] Figure 4 is an exemplary schematic diagram of the CCUS / CCS geological body integrity evaluation effect obtained according to some embodiments of the present application;

[0037] Figure 5 is an exemplary schematic diagram of the CCUS / CCS geological body integrity evaluation device based on microseismic shown in some embodiments of the present application;

[0038] Figure 6 is an exemplary structural schematic diagram of an electronic device shown in some embodiments of the present application. Detailed implementation manners

[0039] To more clearly illustrate the technical solutions of the embodiments of the present application, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings in the following description are only some examples or embodiments of the present application. For those of ordinary skill in the art, without creative efforts, the present application can also be applied to other similar scenarios based on these drawings. Unless obvious from the language context or otherwise stated, the same reference numerals in the figures represent the same structure or operation.

[0040] It should be understood that the "system", "device", "unit" and / or "module" used herein is a method for distinguishing different components, elements, parts, portions or assemblies at different levels. However, if other words can achieve the same purpose, the words can be replaced by other expressions.

[0041] As shown in the present application and the claims, unless the context clearly indicates an exceptional situation, words such as "a", "an", "one" and / or "the" are not specifically singular and may also include the plural. Generally speaking, the terms "comprising" and "including" only indicate the inclusion of the clearly identified steps and elements, and these steps and elements do not constitute an exclusive list. The method or device may also include other steps or elements.

[0042] Flowcharts are used in the present application to illustrate the operations performed by the systems according to the embodiments of the present application. It should be understood that the operations before or after may not necessarily be executed precisely in sequence. On the contrary, the steps can be processed in reverse order or simultaneously. At the same time, other operations can also be added to these processes, or one or several steps can be removed from these processes.

[0043] For ease of understanding, the technical solutions of the present application are introduced below in conjunction with the accompanying drawings and embodiments.

[0044] Figure 1 is an exemplary flowchart of a microseismic-based CCUS / CCS geological body integrity evaluation method according to some embodiments of the present application. As Figure 1 shown, the microseismic-based CCUS / CCS geological body integrity evaluation method includes the following steps:

[0045] Step S110, in response to a plurality of seismic signal detectors arranged within a preset range of the target CCUS / CCS geological body, detecting data information of microseismic signals from the target CCUS / CCS geological body, and obtaining microseismic event attributes corresponding to the microseismic signals and the pore pressure of the formation where the microseismic signals are generated according to the data information of the microseismic signals.

[0046] The target CCUS / CCS geological body is a geological body that needs to monitor its operation safety and stores carbon dioxide using CCUS / CCS technology.

[0047] In a specific implementation process, a plurality of seismic signal detectors can be arranged within the preset range of the target CCUS / CCS geological body to monitor microseismic signals from the target CCUS / CCS geological body (as Figure 2 shown).

[0048] In a specific implementation process, the data information of microseismic signals can be collected using seismic signal detectors. The data information of microseismic signals includes but is not limited to: the amplitude value of the microseismic signal picked up by each seismic signal detector, the first arrival time of the microseismic signal picked up by each seismic signal detector, etc.

[0049] The microseismic event attribute is a characteristic that can reflect the intensity of the microseismic signal. The microseismic event attribute can include but is not limited to the microseismic event magnitude, microseismic event energy, microseismic event signal-to-noise ratio, etc.

[0050] When there are pores and fractures in the rock mass, the change of pore pressure will cause the stress state in the rock mass to change, thus generating microseismic signals. Therefore, the pore pressure of the formation where the microseismic signals are generated can also reflect the intensity of the microseismic signals.

[0051] In some embodiments, the microseismic event attribute includes the microseismic event magnitude, and the following formula can be used to calculate the numerical value of the microseismic event magnitude M w of:

[0052]

[0053] where, v p$v_p$ is the longitudinal wave propagation velocity of the microseismic signal, which can be calculated from the first arrival data picked up by the geophones; $\Omega_0$ is the zero-frequency limit value of the source spectrum; $F$ is the average propagation coefficient of the seismic wave. The average propagation coefficient $F$ of the longitudinal wave is 0.52, and the average propagation coefficient $F$ of the shear wave is 0.63; $d$ is the average value of the distance from the source of the microseismic signal to each geophone.

[0054] In the specific implementation process, according to the data information of the microseismic signal, the distance from the source of the microseismic signal to each geophone can be obtained by using the multi-wave joint location algorithm.

[0055] The multi-wave joint location algorithm can make full use of clear and effective signals to accurately locate microseismic events. Its principle is to perform polarization rotation on the three-component signals, determine the azimuth of the microseismic event through the vibration direction of the P-wave, and then use the P-S wave time difference and velocity model to determine the spatial position of the microseismic event.

[0056] Let the source position coordinates of the $k$-th microseismic signal be $Q$ k (x qk , y qk , z qk ), and the coordinates of the $i$-th geophone be $P$ i (x pi , y pi , z pi ), then the distance between the source and the observation point is:

[0057]

[0058] where $d$ ki represents the distance between the source and the $i$-th geophone, in meters; $x$ qk , y qk , z qk respectively represent the coordinates of the source position of the $k$-th microseismic signal on the x, y, and z axes in three-dimensional space, in meters; $x$ pi , y pi , z pi respectively represent the coordinates of the $i$-th geophone on the x, y, and z axes in three-dimensional space, in meters.

[0059] Given the longitudinal and transverse wave propagation velocities $v$ p and $v$ s , then the P-S wave arrival time difference $\Delta T$ i recorded at point $p$ ki can be expressed by the following formula, that is

[0060]

[0061]

[0062] Among them, T Spick is the first arrival travel time of the S-wave (shear wave) of the microseismic signal, and T Ppick is the first arrival travel time of the P-wave (compressional wave) of the microseismic event.

[0063] Then:

[0064]

[0065] Among them, ΔT ki represents that the source coordinates recorded at point p i are Q i (x qk , y qk , z qk ), and the shear-wave and compressional-wave travel time difference of the microseismic signal.

[0066] When the number of geophones receiving the seismic signal is greater than 3, the above formula can be solved to obtain the microseismic positioning result as shown in Figure 3 .

[0067] After obtaining the distance from the source of the microseismic signal to each geophone using the above method, the average value of all distances can be calculated to obtain the average value d of the distance from the source of the microseismic signal to each geophone in formula (1).

[0068] In some embodiments, the microseismic event attribute includes the microseismic event energy, and the following formula can be used to calculate the value of the microseismic event energy E nergy :

[0069]

[0070] Among them, ρ is the density of the medium of the formation where the microseismic signal is generated; t is the average value of the first arrival time of the microseismic signal picked up by each geophone; W rms,i is the root mean square amplitude of the microseismic signal detected by the i-th geophone; K is the number of geophones detecting the microseismic signal.

[0071] In some embodiments, the microseismic event attribute includes the microseismic event signal-to-noise ratio, and the following formula can be used to calculate the value of the microseismic event signal-to-noise ratio R SN :

[0072]

[0073] Among them, RMS signal is the root mean square amplitude of the effective signal detected by the geophone; RMS noise is the root mean square amplitude of the noise detected by the geophone;

[0074] Among them, the RMS can be obtained using the following formula signal :

[0075]

[0076] where N is the number of seismic signal detectors that detect microseismic signals; A i is the amplitude value of the microseismic signal detected by the i-th seismic signal detector.

[0077] In some embodiments, the pore pressure P can be calculated using the following formula p :

[0078]

[0079] where P o is the overburden pressure of the target CCUS / CCS geological body, in megapascals (MPa); obtained by density logging; h is the vertical depth of the formation where the microseismic signal is generated, in meters, obtained using the multi-wave joint location algorithm based on the data information of the microseismic signal; P h is the hydrostatic pressure of the formation where the microseismic signal is generated, in megapascals (MPa), P h = 0.0098ρh; DT n is the acoustic wave transit time of shale under normal pressure (i.e., the acoustic wave transit time of the normal trend line), in microseconds per inch (μs / ft); DT is the measured acoustic wave transit time of the formation by logging, in microseconds per inch (μs / ft); b is the Eaton index (usually taken as the constant 3); ρ is the medium density of the formation space where the microseismic signal is generated, in kilograms per cubic meter.

[0080] Step S120, obtain the geological body integrity evaluation safety threshold of the target CCUS / CCS geological body according to the microseismic event attributes and pore pressure.

[0081] In the specific implementation process, in order to facilitate the intuitive evaluation of the geological body integrity of the target CCUS / CCS geological body, it is necessary to perform normalization processing on the microseismic event attributes and pore pressure obtained in step S110 respectively. The normalization processing formula is as follows

[0082]

[0083] where x norm is the normalized data, x norm ∈[0 1]; x is the original data of the microseismic event attributes and pore pressure; x max and x min are the maximum and minimum values in the original data of the microseismic event attributes and pore pressure respectively.

[0084] In some embodiments, based on the normalized microseismic event attributes and pore pressure, the weighted analysis function shown below can be used to calculate the geological body integrity evaluation safety valve threshold S ccus / ccs :

[0085] S ccus / ccs = w1M w1 + w2E nergy2 + w3R SN3 + w4P p4 (11)

[0086] where w1, w2, w3, and w4 are the weight coefficients of the microseismic event magnitude, microseismic event energy, microseismic event signal-to-noise ratio, and pore pressure respectively, {w1w2w3w4} ∈ [0 1]; M w1 、E nergy2 、R SN3 、P p4 are the normalized values of the microseismic event magnitude, the normalized value of the microseismic event energy, the normalized value of the microseismic event signal-to-noise ratio, and the normalized value of the pore pressure respectively, {M w1 E nergy2 R SN3 P p4} ∈ [0 1].

[0087] Step S130, evaluate the geological body integrity of the target CCUS / CCS geological body according to the geological body integrity evaluation safety valve threshold.

[0088] Merely as an example, as Figure 4 shown, when the carbon dioxide in the target CCUS / CCS geological body is above the first interface, the geological body integrity evaluation safety valve threshold range is the first preset threshold range: [0.8 1]; when the carbon dioxide in the target CCUS / CCS geological body is between the first interface and the second interface, the geological body integrity evaluation safety valve threshold range is the second preset threshold range: [0.6 0.8]; when the carbon dioxide in the target CCUS / CCS geological body is below the second interface, the corresponding geological body integrity evaluation safety valve threshold range is the third preset threshold range: [0 0.6].

[0089] In some embodiments, when the geological body integrity evaluation safety valve threshold is within the first preset threshold range, it is determined that the carbon dioxide leakage risk of the target CCUS / CCS geological body is at a high risk level. In this case, the target CCUS / CCS geological body is in a high-risk state of carbon dioxide leakage, with a relatively high failure risk, and measures need to be taken in a timely manner to avoid carbon dioxide leakage.

[0090] When the safety valve threshold value for geological body integrity evaluation is within the range of the second preset threshold value, the carbon dioxide leakage risk of the target CCUS / CCS geological body is determined to be at the medium risk level. In this case, the target CCUS / CCS geological body is in a medium-risk state of carbon dioxide leakage. At this time, the gas injection parameters should be adjusted to avoid the transformation of the carbon dioxide leakage risk of the target CCUS / CCS geological body to the high-risk level.

[0091] When the safety valve threshold value for geological body integrity evaluation is within the range of the third preset threshold value, the carbon dioxide leakage risk of the target CCUS / CCS geological body is determined to be at the low risk level. In this case, the target CCUS / CCS geological body is in a low-risk state of carbon dioxide leakage, and carbon dioxide can be injected normally according to the current gas injection parameters.

[0092] In the embodiments provided in the present application, according to the microseismic event attributes and pore pressure, the safety valve threshold value for geological body integrity evaluation of the target CCUS / CCS geological body is obtained; according to the safety valve threshold value for geological body integrity evaluation, the geological body integrity of the target CCUS / CCS geological body is evaluated. Thus, the real-time monitoring of the geological body integrity of CCUS / CCS can be effectively carried out, the risk assessment of carbon dioxide leakage of the target CCUS / CCS geological body can be carried out, the safe operation of the geological body of the target CCUS / CCS geological body can be guaranteed, and a solid technical support can be provided for carbon dioxide storage and utilization.

[0093] Figure 5 It is an exemplary schematic diagram of a microseismic-based CCUS / CCS geological body integrity evaluation device shown in some embodiments of the present application.

[0094] As Figure 5 shown, the microseismic-based CCUS / CCS geological body integrity evaluation device includes: a first acquisition module 510, a second acquisition module 520, and an evaluation module 530.

[0095] The first acquisition module 510 is configured to, in response to a plurality of seismic signal detectors disposed within a preset range of the target CCUS / CCS geological body, detect data information of microseismic signals from the target CCUS / CCS geological body, and obtain microseismic event attributes corresponding to the microseismic signals and pore pressure of the formation where the microseismic signals are generated according to the data information of the microseismic signals.

[0096] The second acquisition module 520 is configured to obtain the safety valve threshold value for geological body integrity evaluation of the target CCUS / CCS geological body according to the microseismic event attributes and the pore pressure.

[0097] An evaluation module 530 is configured to evaluate the geological body integrity of the target CCUS / CCS geological body according to the safety valve threshold for evaluating the geological body integrity.

[0098] In the embodiments of the above-described microseismic-based CCUS / CCS geological body integrity evaluation device, the specific processing of each module and the resulting technical effects can be respectively referred to the relevant descriptions in the corresponding method embodiments, which will not be elaborated herein.

[0099] Figure 6 It is an exemplary structural schematic diagram of an electronic device according to some embodiments of the present application.

[0100] As Figure 6 shown, the electronic device includes: at least one processor 601, at least one communication interface 602, at least one memory 603, and at least one communication bus 604; optionally, the communication interface 602 may be an interface of a communication module, such as an interface of a GSM module; the processor 601 may be a processor CPU, or a specific integrated circuit ASIC (Application Specific Integrated Circuit), or one or more integrated circuits configured to implement the embodiments of the present invention. The memory 603 may include a high-speed RAM memory, and may also include a non-volatile memory, such as at least one disk memory. Among them, the memory 603 stores a program, and the processor 601 calls the program stored in the memory 603 to execute some or all of the above method embodiments.

[0101] The present application relates to a storage medium for storing a computer-readable program, which when run, executes some or all of the above method embodiments.

[0102] Optionally, the storage medium may be a non-temporary computer-readable storage medium. For example, the non-temporary computer-readable storage medium may be a ROM, a random access memory (RAM), a CD-ROM, a magnetic tape, a floppy disk, and an optical data storage device, etc.

[0103] Based on the same inventive concept, the embodiments of the present application also provide a computer program product, including a computer program, which when executed by a processor, implements some or all of the above method embodiments.

[0104] The basic concepts have been described above. Obviously, for those skilled in the art, the above detailed disclosure is only an example and does not constitute a limitation to this application. Although not explicitly stated here, those skilled in the art may make various modifications, improvements, and corrections to this application. Such modifications, improvements, and corrections are proposed in this application, so such modifications, improvements, and corrections still fall within the spirit and scope of the exemplary embodiments of this application.

[0105] At the same time, this application uses specific terms to describe the embodiments of this application. For example, "one embodiment", "an embodiment", and / or "some embodiments" mean a certain feature, structure, or characteristic related to at least one embodiment of this application. Therefore, it should be emphasized and noted that the "one embodiment" or "an embodiment" or "an alternative embodiment" mentioned twice or more at different positions in this application does not necessarily refer to the same embodiment. In addition, certain features, structures, or characteristics in one or more embodiments of this application can be appropriately combined.

[0106] In addition, unless explicitly stated in the claims, the order of the processing elements and sequences, the use of numbers and letters, or the use of other names in this application are not used to limit the order of the processes and methods of this application. Although some currently considered useful invention embodiments have been discussed through various examples in the above disclosure, it should be understood that such details only serve the purpose of illustration. The appended claims are not limited to the disclosed embodiments. On the contrary, the claims are intended to cover all modifications and equivalent combinations that conform to the essence and scope of the embodiments of this application. For example, although the system components described above can be implemented by hardware devices, they can also be implemented only through software solutions, such as installing the described system on existing servers or mobile devices.

[0107] Similarly, it should be noted that, in order to simplify the expression of the disclosure of this application and thus help the understanding of one or more invention embodiments, in the previous description of the embodiments of this application, sometimes multiple features are merged into one embodiment, drawing, or description thereof. However, this disclosure method does not mean that the features required by the object of this application are more than those mentioned in the claims. In fact, the features of the embodiments are fewer than all the features of the individual embodiments disclosed above.

[0108] In some embodiments, numbers are used to describe components and the quantity of attributes. It should be understood that such numbers used in the description of embodiments are, in some examples, modified by the modifiers "about", "approximately" or "substantially". Unless otherwise specified, "about", "approximately" or "substantially" indicate that the said numbers are allowed to have a variation of ±20%. Accordingly, in some embodiments, the numerical parameters used in the specification and claims are approximate values, and such approximate values may vary according to the characteristics required by individual embodiments. In some embodiments, the numerical parameters should consider the specified significant digits and adopt the method of retaining the general number of digits. Although the numerical ranges and parameters used in some embodiments of the present application to confirm the breadth of their scope are approximate values, in specific embodiments, such numerical settings are as precise as possible within the feasible range.

[0109] For each patent, patent application, patent application publication and other materials cited in the present application, such as articles, books, specifications, publications, documents, etc., their entire contents are hereby incorporated into the present application by reference. Except for the application history documents that are inconsistent with or conflict with the content of the present application, and also except for the documents that limit the broadest scope of the claims of the present application (currently or subsequently appended to the present application). It should be noted that if there are inconsistencies or conflicts between the descriptions, definitions, and / or uses of terms in the attached materials of the present application and the content described in the present application, the descriptions, definitions, and / or uses of terms in the present application shall prevail.

[0110] Finally, it should be understood that the embodiments described in the present application are only used to illustrate the principles of the embodiments of the present application. Other variations may also fall within the scope of the present application. Therefore, by way of example and not limitation, alternative configurations of the embodiments of the present application may be considered to be consistent with the teachings of the present application. Accordingly, the embodiments of the present application are not limited to the embodiments explicitly introduced and described in the present application.

Claims

1. A method for evaluating the integrity of CCUS / CCS geological bodies based on microseismicity, characterized in that, The method includes: In response to a plurality of seismic signal detectors arranged within a preset range of a target CCUS / CCS geological body, detecting data information of microseismic signals from the target CCUS / CCS geological body, and obtaining microseismic event attributes corresponding to the microseismic signals and the pore pressure of the formation where the microseismic signals are generated according to the data information of the microseismic signals; Obtaining a geological body integrity evaluation safety valve threshold of the target CCUS / CCS geological body according to the microseismic event attributes and the pore pressure; Evaluating the geological body integrity of the target CCUS / CCS geological body according to the geological body integrity evaluation safety valve threshold.

2. The method according to claim 1, wherein The microseismic event attributes include microseismic event magnitude; The obtaining of the microseismic event attributes corresponding to the microseismic signals according to the data information of the microseismic signals includes: The magnitude M of the microseismic event is calculated using the following formula w for the value of where v p is the longitudinal wave propagation velocity of the microseismic signal; Ω0 is the zero-frequency limit value of the source spectrum; F is the average propagation coefficient of the seismic wave; d is the average value of the distances from the source of the microseismic signal to each seismic signal detector, which is obtained by using a multi-wave joint location algorithm according to the data information of the microseismic signal.

3. The method according to claim 2, wherein The microseismic event attributes further include microseismic event energy; The obtaining of the microseismic event attributes corresponding to the microseismic signals according to the data information of the microseismic signals includes: The energy E of the microseismic event is calculated using the following formula nergy for the numerical value of where ρ is the density of the medium of the formation that generates the microseismic signal; t is the average of the first arrival times of each microseismic signal picked up by the seismic signal detectors; W rms,i is the root mean square amplitude of the microseismic signal detected by the i-th seismic signal detector; K is the number of seismic signal detectors that detect the microseismic signal.

4. The method according to claim 3, wherein The microseismic event attributes further include microseismic event signal-to-noise ratio; The obtaining of the microseismic event attributes corresponding to the microseismic signals according to the data information of the microseismic signals includes: The signal-to-noise ratio R of the microseismic event is calculated using the following formula SN to obtain the value of Among them, RMS signal is the root mean square amplitude of the effective signal detected by the seismic signal detector; RMS noise is the root mean square amplitude of the noise detected by the seismic signal detector; Among them, RMS is obtained using the following formula signal :[[]]END]] Where N is the number of geophones that detect microseismic signals; A i is the amplitude value of the microseismic signal detected by the i-th geophone.

5. The method according to claim 4, wherein The obtaining of the pore pressure of the formation where the microseismic signals are generated according to the data information of the microseismic signals includes: The pore pressure P is calculated using the following formula p :[[-END]] Among them, P o is the overburden pressure of the target CCUS / CCS geological body, calculated from density logging; h is the vertical depth of the formation where the microseismic signal is generated, obtained by using the multi-wave joint location algorithm based on the data information of the microseismic signal; P h is the hydrostatic pressure of the formation where the microseismic signal is generated; DT n is the acoustic travel time difference of shale under normal pressure; DT is the measured acoustic travel time difference of the formation by logging; b is the Eaton index; ρ is the medium density of the formation space where the microseismic signal is generated.

6. The method according to claim 5, wherein The obtaining of the geological body integrity evaluation safety valve threshold of the target CCUS / CCS geological body according to the microseismic event attributes and the pore pressure includes: Performing normalization processing on the microseismic event magnitude, the microseismic event energy, the microseismic event signal-to-noise ratio, and the pore pressure respectively; According to the magnitude of the microseismic event, the energy of the microseismic event, the signal-to-noise ratio of the microseismic event, and the pore pressure after normalization, the safety valve threshold S for evaluating the integrity of the geological body is calculated using the weighted analysis function shown below ccus / ccs : S ccus / cvs = w1M w1 + w2E nergy2 + w3R SN3 + w4P p4 Among them, w1, w2, w3, and w4 are the weight coefficients of the microseismic event magnitude, the microseismic event energy, the microseismic event signal-to-noise ratio, and the pore pressure, respectively; M w1 , E nergy2 , R SN3 , P p4 are the normalized values of the microseismic event magnitude, the normalized value of the microseismic event energy, the normalized value of the microseismic event signal-to-noise ratio, and the normalized value of the pore pressure, respectively.

7. The method according to claim 6, characterized in that The evaluating of the geological body integrity of the target CCUS / CCS geological body according to the geological body integrity evaluation safety valve threshold includes: In response to the geological body integrity evaluation safety valve threshold being within a first preset threshold range, determining that the carbon dioxide leakage risk of the target CCUS / CCS geological body is a high risk level; In response to the geological body integrity evaluation safety valve threshold being within a second preset threshold range, determining that the carbon dioxide leakage risk of the target CCUS / CCS geological body is a medium risk level; In response to the geological body integrity evaluation safety valve threshold being within a third preset threshold range, determining that the carbon dioxide leakage risk of the target CCUS / CCS geological body is a low risk level.

8. An integrity evaluation device for CCUS / CCS geological bodies based on microseismicity, characterized in that, The device includes: A first acquisition module, configured to, in response to a plurality of seismic signal detectors arranged within a preset range of a target CCUS / CCS geological body, detect data information of microseismic signals from the target CCUS / CCS geological body, and obtain microseismic event attributes corresponding to the microseismic signals and the pore pressure of the formation where the microseismic signals are generated according to the data information of the microseismic signals; A second acquisition module, configured to obtain a geological body integrity evaluation safety valve threshold of the target CCUS / CCS geological body according to the microseismic event attributes and the pore pressure; An evaluation module, configured to evaluate the geological body integrity of the target CCUS / CCS geological body according to the safety valve threshold for evaluating the geological body integrity.

9. An electronic device, comprising a memory and a processor, wherein the memory stores a computer program, and when the processor runs the program, it executes the method according to any one of claims 1 to 7.

10. A storage medium, configured to store a computer-readable program, and when the computer-readable program is run, it executes the method according to any one of claims 1 to 7.

Images