CN113987851B - Fluid conductance characteristic analysis method, system and equipment for carbon dioxide migration channel - Google Patents

Fluid conductance characteristic analysis method, system and equipment for carbon dioxide migration channel Download PDF

Info

Publication number
CN113987851B
CN113987851B CN202111620381.XA CN202111620381A CN113987851B CN 113987851 B CN113987851 B CN 113987851B CN 202111620381 A CN202111620381 A CN 202111620381A CN 113987851 B CN113987851 B CN 113987851B
Authority
CN
China
Prior art keywords
carbon dioxide
migration channel
crack
channel
dioxide migration
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202111620381.XA
Other languages
Chinese (zh)
Other versions
CN113987851A (en
Inventor
付博烨
张旺
单小彩
汪文洋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Geology and Geophysics of CAS
Original Assignee
Institute of Geology and Geophysics of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Geology and Geophysics of CAS filed Critical Institute of Geology and Geophysics of CAS
Priority to CN202111620381.XA priority Critical patent/CN113987851B/en
Publication of CN113987851A publication Critical patent/CN113987851A/en
Application granted granted Critical
Publication of CN113987851B publication Critical patent/CN113987851B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/08—Probabilistic or stochastic CAD
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Physical Or Chemical Processes And Apparatus (AREA)

Abstract

The invention belongs to the field of rock fluid migration, and particularly relates to a method, a system and equipment for analyzing fluid conductance characteristics of a carbon dioxide migration channel, aiming at solving the problem that the prior art can only finish the evaluation of the permeability of an underground structure and cannot describe the change rule of carbon dioxide migration capacity along with time. The invention comprises the following steps: acquiring the surface morphology, the length and the aperture of the crack; carrying out parametric characterization on the crack surface to form a power spectrum function and generating different random crack surfaces; forming a closed crack space as a carbon dioxide migration channel; obtaining a carbon dioxide flow rate change curve along with time, migration channel permeability and elastic modulus based on numerical simulation; and generating an analysis gauge plate for inverting the carbon dioxide migration channel form and the carbon dioxide migration capacity based on the elastic characteristics of the carbon dioxide migration channel. The invention can not only complete the evaluation of the permeability of the underground structure. The change rule of the carbon dioxide transport capacity along with time can also be described.

Description

Fluid conductance characteristic analysis method, system and equipment for carbon dioxide migration channel
Technical Field
The invention belongs to the field of rock fluid migration, and particularly relates to a method, a system and equipment for analyzing fluid conductance characteristics of a carbon dioxide migration channel.
Background
The underground sealing of carbon dioxide is one of important means for realizing carbon neutralization, but the precondition for realizing the underground sealing of carbon dioxide is that the carbon dioxide transportation and conduction capability of carbon dioxide migration channels such as underground cracks and the like is comprehensively known. Currently, the detection of other characteristics of the subsurface medium based on its seismic wave properties is one of the common approaches.
However, the existing quantitative method for judging the carbon dioxide transportation and conduction capacity of the underground medium by using the seismic wave attributes is still incomplete, most methods can only complete the evaluation of the permeability of the underground structure, the real carbon dioxide transportation process is strongly influenced by time, and the change rule of the carbon dioxide transportation capacity along with the time cannot be described in the prior art.
Disclosure of Invention
In order to solve the above problems in the prior art, that is, the prior art can only complete the evaluation of the permeability of the underground structure, but cannot describe the change rule of the carbon dioxide transport capacity along with the time, the invention provides a fluid transport characteristic analysis method of a carbon dioxide transport channel, which comprises the following steps:
step S10, acquiring morphological parameter characteristics forming a carbon dioxide migration channel based on crack outcrop, core observation data and logging data; the morphological parameter characteristics comprise fracture surface morphology, fracture length and fracture aperture;
step S20, based on the morphological parameter characteristics, performing statistical analysis on the crack surface morphology, the crack length and the crack aperture to obtain the carbon dioxide migration channel characteristics; the carbon dioxide migration channel characteristics comprise an affine function for representing the height level spatial distribution of the surface of the crack, a probability density function for representing that the crack has a rough surface and morphological parameters for representing the surface morphology of the crack;
step S30, constructing a power spectrum function of the carbon dioxide migration channel form based on the carbon dioxide migration channel characteristics; wherein the morphological parameters comprise self-affine coefficients, mean and standard deviation of height distribution;
step S40, setting different self-affine coefficients, mean values and standard deviations of height distribution for the power spectrum function of the carbon dioxide migration channel form, and generating different random crack surfaces through Fourier transformation;
step S50, taking the highest plane of the surface of the crack as the top interface of the crack to form a closed crack space as a carbon dioxide migration channel, and obtaining the carbon dioxide migration channel with different self-affine coefficients, average values of height distribution and standard deviations;
step S60, acquiring the elastic modulus and the carbon dioxide transportation capacity of the carbon dioxide migration channel through numerical simulation by using a local cubic law, a Navigneaux equation and a Hooke law; the carbon dioxide transportation capacity comprises a curve of the permeability of a carbon dioxide migration channel and the change of carbon dioxide flow along with time;
and step S70, generating an analysis gauge plate for inverting the form of the carbon dioxide migration channel and the carbon dioxide migration capacity based on the elastic characteristics of the carbon dioxide migration channel based on the relation among the elastic modulus, the form parameters and the carbon dioxide transportation capacity of the carbon dioxide migration channel.
In some preferred embodiments, the carbon dioxide migration channel morphology has a power spectrum function expressed as:
Figure DEST_PATH_IMAGE001
wherein the content of the first and second substances,
Figure 7917DEST_PATH_IMAGE002
as a function of the power spectrum of the carbon dioxide migration channel morphology,
Figure DEST_PATH_IMAGE003
the wave number of the surface of the carbon dioxide migration channel in the wave number domain,
Figure 668093DEST_PATH_IMAGE004
is the self-affine coefficient of the carbon dioxide migration channel form,
Figure DEST_PATH_IMAGE005
respectively representing the mean and standard deviation of the height distribution of the carbon dioxide migration channel morphology,
Figure 539097DEST_PATH_IMAGE006
is a probability density function satisfied by the carbon dioxide migration channel shape,
Figure DEST_PATH_IMAGE007
to characterize the height level of the fracture surfaceA spatially distributed self-affine function.
In some preferred embodiments, the different random fracture surfaces are generated in step S40 by fourier transformation by:
Figure 660506DEST_PATH_IMAGE008
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE009
are data points
Figure 185028DEST_PATH_IMAGE010
The FFT represents the fourier transform.
In some preferred embodiments, the curve of the carbon dioxide flow rate with time is obtained by a random walk around particle tracing method;
step S601, when the total number is
Figure DEST_PATH_IMAGE011
When the carbon dioxide fluid unit passes through the carbon dioxide transport passage, the second
Figure 611461DEST_PATH_IMAGE012
The time required for each cell was:
Figure 535424DEST_PATH_IMAGE013
Figure 589968DEST_PATH_IMAGE014
Figure 488654DEST_PATH_IMAGE015
wherein the content of the first and second substances,
Figure 136804DEST_PATH_IMAGE016
is as follows
Figure 208665DEST_PATH_IMAGE012
A unit arrives at a node
Figure 507928DEST_PATH_IMAGE017
The time required for the operation of the apparatus,
Figure 843094DEST_PATH_IMAGE018
representing slave nodes
Figure 837595DEST_PATH_IMAGE017
Arriving node
Figure 854093DEST_PATH_IMAGE019
At any one point in time, the time of the day,
Figure 883229DEST_PATH_IMAGE020
is as follows
Figure 42159DEST_PATH_IMAGE012
A volume of carbon dioxide comprising a carbon dioxide fluid unit,
Figure 992797DEST_PATH_IMAGE021
as slave nodes in unit time
Figure 406461DEST_PATH_IMAGE022
Any point arrives at node
Figure 165470DEST_PATH_IMAGE017
The volume of carbon dioxide in the gas mixture,
Figure 639176DEST_PATH_IMAGE023
for the mesh width of the plane subdivision of the carbon dioxide migration channel,
Figure 732903DEST_PATH_IMAGE024
for nodes in carbon dioxide channels
Figure 153520DEST_PATH_IMAGE017
And node
Figure 891669DEST_PATH_IMAGE025
The geometric mean of the height cube,Lrepresentative node
Figure 411643DEST_PATH_IMAGE017
And node
Figure 133612DEST_PATH_IMAGE025
The distance between the two or more of the two or more,
Figure 13712DEST_PATH_IMAGE026
and
Figure DEST_PATH_IMAGE027
is a node
Figure 871947DEST_PATH_IMAGE017
And node
Figure 562822DEST_PATH_IMAGE025
The pore pressure of (a) is set,
Figure 240928DEST_PATH_IMAGE028
is the viscosity coefficient of carbon dioxide;
step S602, obtaining a time-varying curve of the carbon dioxide flow rate of the carbon dioxide transport passage by using the time as a variable through the method of step S601.
In some preferred embodiments, the permeability of the carbon dioxide migration channel is expressed as:
Figure DEST_PATH_IMAGE029
Figure 659140DEST_PATH_IMAGE030
Figure DEST_PATH_IMAGE031
wherein the content of the first and second substances,
Figure 43985DEST_PATH_IMAGE032
is the permeability of the carbon dioxide migration channel,
Figure DEST_PATH_IMAGE033
representing the pressure difference across the channel driving carbon dioxide transport,
Figure 764816DEST_PATH_IMAGE034
and
Figure DEST_PATH_IMAGE035
representing the length and width of the carbon dioxide migration channel.
In some preferred embodiments, the surface of the carbon dioxide migration channel is sectioned into a fine cylindrical system, and the elastic modulus of the carbon dioxide migration channel is calculated using hooke's law, which is expressed as:
Figure 120099DEST_PATH_IMAGE036
Figure DEST_PATH_IMAGE037
wherein the content of the first and second substances,
Figure 827155DEST_PATH_IMAGE038
is the elastic modulus of the carbon dioxide migration channel,
Figure DEST_PATH_IMAGE039
in order to realize the maximum length in the cylindrical system when the plane of the carbon dioxide migration channel is split,
Figure 315774DEST_PATH_IMAGE040
is the first in a cylinder system
Figure 676348DEST_PATH_IMAGE017
Length shrinkage of cylinder，
Figure 329046DEST_PATH_IMAGE041
The total number of cylinders in the cylinder system.
In some preferred embodiments, step S70 includes:
step S71, generating a chart form based on the elastic modulus, the form parameters and the carbon dioxide transmission capacity of the carbon dioxide migration channel;
and step S72, reversely deducing the morphological parameters of the carbon dioxide migration channel and the carbon dioxide transportation capacity of the carbon dioxide migration channel based on the elastic modulus of the carbon dioxide migration channel according to the data trend in the chart form.
In another aspect of the present invention, an analysis system for fluid conductance characteristics of a carbon dioxide transport channel is provided, the analysis system comprising the following modules:
the data acquisition module is configured to acquire morphological parameter characteristics forming a carbon dioxide migration channel based on crack outcrop, core observation data and logging data; the morphological parameter characteristics comprise fracture surface morphology, fracture length and fracture aperture;
the carbon dioxide migration channel feature extraction module is configured to perform statistical analysis on the crack surface morphology, the crack length and the crack aperture based on the morphological parameter features to obtain carbon dioxide migration channel features; the carbon dioxide migration channel characteristics comprise an affine function for representing the height level spatial distribution of the surface of the crack, a probability density function for representing that the crack has a rough surface and morphological parameters for representing the surface morphology of the crack;
the power spectrum function construction module is configured to construct a power spectrum function of a carbon dioxide migration channel form based on the carbon dioxide migration channel characteristics; wherein the morphological parameters comprise self-affine coefficients, mean and standard deviation of height distribution;
the crack surface generation module is configured to set different self-affine coefficients, mean values and standard deviations of height distribution for the power spectrum function of the carbon dioxide migration channel form, and generate different random crack surfaces through Fourier transformation;
the migration channel building module is configured to form a closed crack space by taking the highest plane of the crack surface as a top interface of the crack, and obtain the carbon dioxide migration channel with different self-affine coefficients, mean values of height distribution and standard deviation as a carbon dioxide migration channel;
the migration channel parameter calculation module is configured to obtain the elastic modulus and the carbon dioxide transmission capacity of the carbon dioxide migration channel through numerical simulation by utilizing a local cubic law, a Navigneaux equation and a Hooke law; the carbon dioxide transportation capacity comprises a curve of the permeability of a carbon dioxide migration channel and the change of carbon dioxide flow along with time;
and the migration channel analysis gauge plate generation module is configured to generate an analysis gauge plate for inverting the carbon dioxide migration channel form and the carbon dioxide migration capacity based on the elastic characteristics of the carbon dioxide migration channel based on the relation among the elastic modulus, the form parameters and the carbon dioxide transportation capacity of the carbon dioxide migration channel.
In a third aspect of the present invention, an electronic device is provided, including:
at least one processor; and
a memory communicatively coupled to at least one of the processors; wherein the content of the first and second substances,
the memory stores instructions executable by the processor for execution by the processor to implement the fluid conductance characterization method of the carbon dioxide migration channel described above.
In a fourth aspect of the present invention, a computer-readable storage medium is provided, which stores computer instructions for being executed by the computer to implement the fluid conductance analyzing method for the carbon dioxide transport channel.
The invention has the beneficial effects that:
the invention relates to a fluid transport characteristic analysis method of a carbon dioxide migration channel, which establishes a perfect evaluation method of carbon dioxide transport capacity based on underground medium seismic waves (elastic modulus), and evaluates the permeability of an underground transport system and the change of the carbon dioxide transport capacity along with time through the elastic modulus attribute of an underground medium.
Drawings
Other features, objects and advantages of the present application will become more apparent upon reading of the following detailed description of non-limiting embodiments thereof, made with reference to the accompanying drawings in which:
FIG. 1 is a schematic flow chart of a method for analyzing the fluid conducting characteristics of a carbon dioxide transport channel according to the present invention;
FIG. 2 is a graph of random fracture surfaces with self-affine coefficients of 2 and 0.5 for one embodiment of the method for analyzing fluid conductance characteristics of a carbon dioxide migration channel of the present invention;
FIG. 3 is a calculation of the fluid conductance characteristics analysis method of the carbon dioxide transport channel according to an embodiment of the present invention
Figure 105372DEST_PATH_IMAGE042
A schematic diagram of the upper left corner of the temporal planar subdivision method;
FIG. 4 is a graph of the carbon dioxide flow rate of the carbon dioxide migration channel with time varying according to different self-affine coefficients in an embodiment of the method for analyzing the fluid conductance characteristics of the carbon dioxide migration channel according to the present invention;
FIG. 5 is a graph showing the permeability variation curves of carbon dioxide migration channels with different self-affine coefficients, height distribution mean values and standard deviations obtained by using the local cubic law according to an embodiment of the method for analyzing the fluid conductance characteristics of the carbon dioxide migration channels of the present invention;
FIG. 6 is a diagram illustrating an example of a method for analyzing fluid conductance characteristics of a carbon dioxide migration channel according to the present invention, which utilizes Hooke's law to obtain an elastic modulus variation curve of the carbon dioxide migration channel with different self-affine coefficients based on numerical simulation;
FIG. 7 is a schematic diagram of a fracture surface being sectioned into a cylinder system according to an embodiment of the method for analyzing fluid transport characteristics of a carbon dioxide transport channel according to the present invention;
fig. 8 is a rock physical quantity plate including the elastic property of the migration channel and the carbon dioxide transportation capability finally obtained by the fluid transportation characteristic analysis method of the carbon dioxide migration channel according to an embodiment of the invention.
Detailed Description
The present application will be described in further detail with reference to the following drawings and examples. It is to be understood that the specific embodiments described herein are merely illustrative of the relevant invention and not restrictive of the invention. It should be noted that, for convenience of description, only the portions related to the related invention are shown in the drawings.
It should be noted that the embodiments and features of the embodiments in the present application may be combined with each other without conflict. The present application will be described in detail below with reference to the embodiments with reference to the attached drawings.
The invention provides a fluid transportation characteristic analysis method of a carbon dioxide migration channel, further establishes a perfect evaluation method of CO2 (carbon dioxide) transportation capability based on underground medium seismic waves (elastic modulus), evaluates the permeability of an underground transportation system and the change of CO2 transportation capability along with time through the elastic modulus attribute of an underground medium, and provides a theoretical basis for CO2 sequestration site selection.
One application scenario of the invention is as follows: in the evaluation of the fluid migration capacity of the fracture with the known shape, the shape of the underground fracture can be obtained by a seismic exploration method, but after the shape of the underground fracture is obtained, the fluid conductivity of the fracture can not be evaluated by an effective means in the existing method.
The invention discloses a method for analyzing fluid conducting characteristics of a carbon dioxide migration channel, which comprises the following steps:
step S10, acquiring morphological parameter characteristics forming a carbon dioxide migration channel based on crack outcrop, core observation data and logging data; the morphological parameter characteristics comprise fracture surface morphology, fracture length and fracture aperture;
step S20, based on the morphological parameter characteristics, performing statistical analysis on the crack surface morphology, the crack length and the crack aperture to obtain the carbon dioxide migration channel characteristics; the carbon dioxide migration channel characteristics comprise an affine function for representing the height level spatial distribution of the surface of the crack, a probability density function for representing that the crack has a rough surface and morphological parameters for representing the surface morphology of the crack;
step S30, constructing a power spectrum function of the carbon dioxide migration channel form based on the carbon dioxide migration channel characteristics; wherein the morphological parameters comprise self-affine coefficients, mean and standard deviation of height distribution;
step S40, setting different self-affine coefficients, mean values and standard deviations of height distribution for the power spectrum function of the carbon dioxide migration channel form, and generating different random crack surfaces through Fourier transformation;
step S50, taking the highest plane of the surface of the crack as the top interface of the crack to form a closed crack space as a carbon dioxide migration channel, and obtaining the carbon dioxide migration channel with different self-affine coefficients, average values of height distribution and standard deviations;
step S60, acquiring the elastic modulus and the carbon dioxide transportation capacity of the carbon dioxide migration channel through numerical simulation by using a local cubic law, a Navigneaux equation and a Hooke law; the carbon dioxide transportation capacity comprises a curve of the permeability of a carbon dioxide migration channel and the change of carbon dioxide flow along with time;
and step S70, generating an analysis gauge plate for inverting the form of the carbon dioxide migration channel and the carbon dioxide migration capacity based on the elastic characteristics of the carbon dioxide migration channel based on the relation among the elastic modulus, the form parameters and the carbon dioxide transportation capacity of the carbon dioxide migration channel.
In order to more clearly explain the method for analyzing the fluid conducting characteristics of the carbon dioxide transport channel according to the present invention, the following describes the steps in the embodiment of the present invention in detail with reference to fig. 1.
The method for analyzing the fluid conducting characteristics of the carbon dioxide transport passage according to the first embodiment of the present invention includes steps S10-S70, and the steps are described in detail as follows:
step S10, acquiring morphological parameter characteristics forming a carbon dioxide migration channel based on crack outcrop, core observation data and logging data; the morphological parameter characteristics comprise fracture surface morphology, fracture length and fracture aperture.
Step S20, based on the morphological parameter characteristics, performing statistical analysis on the crack surface morphology, the crack length and the crack aperture to obtain the carbon dioxide migration channel characteristics; the carbon dioxide migration channel characteristics comprise an affine function for representing the height level spatial distribution of the surface of the crack, a probability density function for representing that the crack has a rough surface and morphological parameters for representing the surface morphology of the crack.
Step S30, constructing a power spectrum function of the carbon dioxide migration channel form based on the carbon dioxide migration channel characteristics; wherein the morphological parameters include self-affine coefficients, mean and standard deviation of height distribution.
The power spectrum function of the carbon dioxide migration channel form is shown as the formula (1):
Figure 58285DEST_PATH_IMAGE043
wherein the content of the first and second substances,
Figure 511132DEST_PATH_IMAGE002
as a function of the power spectrum of the carbon dioxide migration channel morphology,
Figure 854388DEST_PATH_IMAGE003
the wave number of the surface of the carbon dioxide migration channel in the wave number domain,
Figure 293460DEST_PATH_IMAGE004
is the self-affine coefficient of the carbon dioxide migration channel form,
Figure 976245DEST_PATH_IMAGE005
respectively representing the mean and standard deviation of the height distribution of the carbon dioxide migration channel morphology,
Figure 475360DEST_PATH_IMAGE006
is a probability density function satisfied by the carbon dioxide migration channel shape,
Figure 305912DEST_PATH_IMAGE007
is an affine function characterizing the height level spatial distribution of the fracture surface.
For carbon dioxide migration channels that are typically formed by fractures,
Figure 673309DEST_PATH_IMAGE044
is a probability density function of a lognormal distribution.
Step S40, setting different self-affine coefficients, mean values and standard deviations of height distribution for the power spectrum function of the carbon dioxide migration channel form, and generating different random crack surfaces through Fourier transformation, as shown in formula (2):
Figure 600813DEST_PATH_IMAGE045
wherein the content of the first and second substances,
Figure 146195DEST_PATH_IMAGE009
are data points
Figure 198465DEST_PATH_IMAGE010
The FFT represents the fourier transform. The height of the top interface of the carbon dioxide migration channel is
Figure 244918DEST_PATH_IMAGE009
Maximum value of (2).
As shown in fig. 2, random fracture surfaces with self-affine coefficients of 2 and 0.5 are shown in the left diagram of fig. 2, and random fracture surfaces with self-affine coefficients of 0.5 are shown in the right diagram of fig. 2.
And S50, taking the plane at the highest position of the surface of the crack as the top interface of the crack to form a closed crack space as a carbon dioxide migration channel, and obtaining the carbon dioxide migration channel with different self-affine coefficients, average values of height distribution and standard deviations.
In fig. 2, the space between the upper and lower surfaces is the transport path for carbon dioxide.
Step S60, acquiring the elastic modulus and the carbon dioxide transportation capacity of the carbon dioxide migration channel through numerical simulation by using a local cubic law, a Navigneaux equation and a Hooke law; the carbon dioxide conducting capacity comprises a curve of the permeability of the carbon dioxide migration channel and the change of the carbon dioxide flow with time.
A curve of the carbon dioxide flow changing along with the time is obtained by adopting a random walk around particle tracing method;
step S601, when the total number is
Figure 718992DEST_PATH_IMAGE011
When the carbon dioxide fluid unit passes through the carbon dioxide transport passage, the second
Figure 825488DEST_PATH_IMAGE012
The time required for each unit is as shown in formula (3) to formula (5):
Figure 37158DEST_PATH_IMAGE046
Figure 824985DEST_PATH_IMAGE047
Figure 195924DEST_PATH_IMAGE048
wherein the content of the first and second substances,
Figure 597955DEST_PATH_IMAGE016
is as follows
Figure 687134DEST_PATH_IMAGE012
A unit arrives at a node
Figure 685177DEST_PATH_IMAGE017
The time required for the operation of the apparatus,
Figure 176201DEST_PATH_IMAGE049
representing slave nodes
Figure 483554DEST_PATH_IMAGE017
Arriving node
Figure 997712DEST_PATH_IMAGE022
At any one point in time, the time of the day,
Figure 924080DEST_PATH_IMAGE050
is as follows
Figure 410556DEST_PATH_IMAGE012
A volume of carbon dioxide comprising a carbon dioxide fluid unit,
Figure 764177DEST_PATH_IMAGE051
as slave nodes in unit time
Figure 765631DEST_PATH_IMAGE022
Any point arrives at node
Figure 620324DEST_PATH_IMAGE017
The volume of carbon dioxide in the gas mixture,
Figure 85940DEST_PATH_IMAGE023
for the mesh width of the plane subdivision of the carbon dioxide migration channel,
Figure 485829DEST_PATH_IMAGE024
for nodes in carbon dioxide channels
Figure 974579DEST_PATH_IMAGE017
And node
Figure 508328DEST_PATH_IMAGE025
The geometric mean of the height cube,Lrepresentative node
Figure 690435DEST_PATH_IMAGE017
And node
Figure 651438DEST_PATH_IMAGE025
The distance between the two or more of the two or more,
Figure 299588DEST_PATH_IMAGE026
and
Figure 574712DEST_PATH_IMAGE052
is a node
Figure 483762DEST_PATH_IMAGE017
And node
Figure 5879DEST_PATH_IMAGE025
The pore pressure of (a) is set,
Figure 380DEST_PATH_IMAGE028
is the viscosity coefficient of carbon dioxide;
step S602, obtaining a time-varying curve of the carbon dioxide flow rate of the carbon dioxide transport passage by using the time as a variable through the method of step S601.
FIG. 3 shows a calculation of the fluid conductance characteristics of the carbon dioxide transport channel according to an embodiment of the present invention
Figure 79194DEST_PATH_IMAGE053
Schematic diagram and node of upper left corner part of time plane subdivision method
Figure 718117DEST_PATH_IMAGE022
Is a node
Figure 20923DEST_PATH_IMAGE017
And 4 nearest neighbor nodes after the peripheral plane is split.
As shown in fig. 4, a graph of the carbon dioxide flow rate of the carbon dioxide transport channel with time is shown for the fluid conductance characteristic analysis method of the carbon dioxide transport channel according to an embodiment of the present invention, where H =0.5 is a graph of the carbon dioxide flow rate of the carbon dioxide transport channel with time when the self-affine coefficient is 0.5, H =1 is a graph of the carbon dioxide flow rate of the carbon dioxide transport channel with time when the self-affine coefficient is 1, and H =1.5 is a graph of the carbon dioxide flow rate of the carbon dioxide transport channel with time when the self-affine coefficient is 1.5.
As shown in fig. 5, a permeability variation curve of a carbon dioxide migration channel with different self-affine coefficients, height distribution mean values and standard deviations is obtained by using a local cubic law, which is an embodiment of the method for analyzing the fluid conductance characteristics of the carbon dioxide migration channel according to the present invention, and the permeability of the carbon dioxide migration channel is as shown in formulas (6) to (8):
Figure 892932DEST_PATH_IMAGE054
Figure 306596DEST_PATH_IMAGE055
Figure 65605DEST_PATH_IMAGE056
wherein the content of the first and second substances,
Figure 742574DEST_PATH_IMAGE032
is the permeability of the carbon dioxide migration channel,
Figure 711667DEST_PATH_IMAGE033
representing the pressure difference across the channel driving carbon dioxide transport,
Figure 53655DEST_PATH_IMAGE034
and
Figure DEST_PATH_IMAGE057
representing the length of the carbon dioxide migration passageAnd a width.
As shown in fig. 6, in an embodiment of the method for analyzing fluid transport characteristics of a carbon dioxide transport channel according to the present invention, based on a hooke's law, an elastic modulus variation curve of the carbon dioxide transport channel with different self-affine coefficients is obtained based on numerical simulation, a surface of the carbon dioxide transport channel is divided into fine cylindrical systems, and the elastic modulus of the carbon dioxide transport channel is calculated by using hooke's law, where the elastic modulus of the carbon dioxide transport channel is shown in formula (9) -formula (10):
Figure 463908DEST_PATH_IMAGE058
Figure DEST_PATH_IMAGE059
wherein the content of the first and second substances,
Figure 842937DEST_PATH_IMAGE038
is the elastic modulus of the carbon dioxide migration channel,
Figure 686609DEST_PATH_IMAGE039
in order to realize the maximum length in the cylindrical system when the plane of the carbon dioxide migration channel is split,
Figure 442076DEST_PATH_IMAGE040
is the first in a cylinder system
Figure 503573DEST_PATH_IMAGE017
The amount of contraction in the length of each cylinder,
Figure 928869DEST_PATH_IMAGE041
referring to fig. 7, a schematic diagram of a fracture surface being divided into cylinder systems according to an embodiment of the method for analyzing fluid transport characteristics of a carbon dioxide transport channel according to the present invention is shown, and the divided cylinder systems shown in fig. 6 can be seen.
And step S70, generating an analysis gauge plate for inverting the form of the carbon dioxide migration channel and the carbon dioxide migration capacity based on the elastic characteristics of the carbon dioxide migration channel based on the relation among the elastic modulus, the form parameters and the carbon dioxide transportation capacity of the carbon dioxide migration channel.
Step S71, generating a chart form based on the elastic modulus, the form parameters and the carbon dioxide transmission capacity of the carbon dioxide migration channel;
and step S72, reversely deducing the morphological parameters of the carbon dioxide migration channel and the carbon dioxide transportation capacity of the carbon dioxide migration channel based on the elastic modulus of the carbon dioxide migration channel according to the data trend in the chart form.
As shown in fig. 8, a finally obtained rock physical quantity plate including the elastic property of the migration channel and the carbon dioxide transportation capability is an embodiment of the method for analyzing the fluid transportation characteristics of the carbon dioxide migration channel according to the present invention, and an application process of the present invention is as follows: the method comprises the steps of obtaining form information of underground cracks through exploration means such as seismic exploration and the like, extracting self-affine coefficients, height distribution mean values and standard deviations of the surface forms of the cracks, representing the cracks, of the cracks according to the crack forms obtained through seismic exploration, obtaining crack form parameters, and evaluating the fluid conductance capacity of the cracks by using the crack form parameters, the crack permeability and the change curves of the crack fluid flow along with time, which are formed in the graph 4 and the graph 8.
Although the foregoing embodiments describe the steps in the above sequential order, those skilled in the art will understand that, in order to achieve the effect of the present embodiments, the steps may not be executed in such an order, and may be executed simultaneously (in parallel) or in an inverse order, and these simple variations are within the scope of the present invention.
The fluid conducting characteristic analysis system of the carbon dioxide transport channel of the second embodiment of the invention comprises the following modules:
the data acquisition module is configured to acquire morphological parameter characteristics forming a carbon dioxide migration channel based on crack outcrop, core observation data and logging data; the morphological parameter characteristics comprise fracture surface morphology, fracture length and fracture aperture;
the carbon dioxide migration channel feature extraction module is configured to perform statistical analysis on the crack surface morphology, the crack length and the crack aperture based on the morphological parameter features to obtain carbon dioxide migration channel features; the carbon dioxide migration channel characteristics comprise an affine function for representing the height level spatial distribution of the surface of the crack, a probability density function for representing that the crack has a rough surface and morphological parameters for representing the surface morphology of the crack;
the power spectrum function construction module is configured to construct a power spectrum function of a carbon dioxide migration channel form based on the carbon dioxide migration channel characteristics; wherein the morphological parameters comprise self-affine coefficients, mean and standard deviation of height distribution;
the crack surface generation module is configured to set different self-affine coefficients, mean values and standard deviations of height distribution for the power spectrum function of the carbon dioxide migration channel form, and generate different random crack surfaces through Fourier transformation;
the migration channel building module is configured to form a closed crack space by taking the highest plane of the crack surface as a top interface of the crack, and obtain the carbon dioxide migration channel with different self-affine coefficients, mean values of height distribution and standard deviation as a carbon dioxide migration channel;
the migration channel parameter calculation module is configured to obtain the elastic modulus and the carbon dioxide transmission capacity of the carbon dioxide migration channel through numerical simulation by utilizing a local cubic law, a Navigneaux equation and a Hooke law; the carbon dioxide transportation capacity comprises a curve of the permeability of a carbon dioxide migration channel and the change of carbon dioxide flow along with time;
and the migration channel analysis gauge plate generation module is configured to generate an analysis gauge plate for inverting the carbon dioxide migration channel form and the carbon dioxide migration capacity based on the elastic characteristics of the carbon dioxide migration channel based on the relation among the elastic modulus, the form parameters and the carbon dioxide transportation capacity of the carbon dioxide migration channel.
It can be clearly understood by those skilled in the art that, for convenience and brevity of description, the specific working process and related description of the system described above may refer to the corresponding process in the foregoing method embodiments, and will not be described herein again.
It should be noted that, the fluid transportation characteristic analysis system of the carbon dioxide transportation channel provided in the above embodiment is only illustrated by the division of the above functional modules, and in practical applications, the above functions may be distributed by different functional modules according to needs, that is, the modules or steps in the embodiment of the present invention are further decomposed or combined, for example, the modules in the above embodiment may be combined into one module, or may be further split into a plurality of sub-modules, so as to complete all or part of the above described functions. The names of the modules and steps involved in the embodiments of the present invention are only for distinguishing the modules or steps, and are not to be construed as unduly limiting the present invention.
An electronic apparatus according to a third embodiment of the present invention includes:
at least one processor; and
a memory communicatively coupled to at least one of the processors; wherein the content of the first and second substances,
the memory stores instructions executable by the processor for execution by the processor to implement the fluid conductance characterization method of the carbon dioxide migration channel described above.
A computer-readable storage medium of a fourth embodiment of the present invention stores computer instructions for execution by the computer to implement the fluid conductance characterization method of the carbon dioxide transport channel described above.
It can be clearly understood by those skilled in the art that, for convenience and brevity of description, the specific working processes and related descriptions of the storage device and the processing device described above may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
Those of skill in the art would appreciate that the various illustrative modules, method steps, and modules described in connection with the embodiments disclosed herein may be implemented as electronic hardware, computer software, or combinations of both, and that programs corresponding to the software modules, method steps may be located in Random Access Memory (RAM), memory, Read Only Memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art. To clearly illustrate this interchangeability of electronic hardware and software, various illustrative components and steps have been described above generally in terms of their functionality. Whether such functionality is implemented as electronic hardware or software depends upon the particular application and design constraints imposed on the solution. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
The terms "first," "second," and the like are used for distinguishing between similar elements and not necessarily for describing or implying a particular order or sequence.
The terms "comprises," "comprising," or any other similar term are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus.
So far, the technical solutions of the present invention have been described in connection with the preferred embodiments shown in the drawings, but it is easily understood by those skilled in the art that the scope of the present invention is obviously not limited to these specific embodiments. Equivalent changes or substitutions of related technical features can be made by those skilled in the art without departing from the principle of the invention, and the technical scheme after the changes or substitutions can fall into the protection scope of the invention.

Claims (9)

1. A method for analyzing fluid transport characteristics of a carbon dioxide transport channel, the method comprising:
step S10, acquiring morphological parameter characteristics forming a carbon dioxide migration channel based on crack outcrop, core observation data and logging data; the morphological parameter characteristics comprise fracture surface morphology, fracture length and fracture aperture;
step S20, based on the morphological parameter characteristics, performing statistical analysis on the crack surface morphology, the crack length and the crack aperture to obtain the carbon dioxide migration channel characteristics; the carbon dioxide migration channel characteristics comprise an affine function for representing the height level spatial distribution of the surface of the crack, a probability density function for representing that the crack has a rough surface and morphological parameters for representing the surface morphology of the crack; wherein the morphological parameters comprise self-affine coefficients, mean and standard deviation of height distribution;
step S30, constructing a power spectrum function of the carbon dioxide migration channel form based on the carbon dioxide migration channel characteristics:
Figure 285521DEST_PATH_IMAGE001
wherein the content of the first and second substances,
Figure 625367DEST_PATH_IMAGE002
as a function of the power spectrum of the carbon dioxide migration channel morphology,
Figure 559825DEST_PATH_IMAGE003
the wave number of the surface of the carbon dioxide migration channel in the wave number domain,
Figure 84347DEST_PATH_IMAGE004
is the self-affine coefficient of the carbon dioxide migration channel form,
Figure 979622DEST_PATH_IMAGE005
respectively representing the mean and standard deviation of the height distribution of the carbon dioxide migration channel morphology,
Figure 982213DEST_PATH_IMAGE006
is a dioxide of oxygenThe carbon migration channel morphology satisfies a probability density function,
Figure 36756DEST_PATH_IMAGE007
an auto-affine function representing the height level spatial distribution of the fracture surface;
step S40, setting different self-affine coefficients, mean values and standard deviations of height distribution for the power spectrum function of the carbon dioxide migration channel form, and generating different random crack surfaces through Fourier transformation;
step S50, taking the highest plane of the surface of the crack as the top interface of the crack to form a closed crack space as a carbon dioxide migration channel, and obtaining the carbon dioxide migration channel with different self-affine coefficients, average values of height distribution and standard deviations;
step S60, acquiring the elastic modulus and the carbon dioxide transportation capacity of the carbon dioxide migration channel through numerical simulation by using a local cubic law, a Navigneaux equation and a Hooke law; the carbon dioxide transportation capacity comprises a curve of the permeability of a carbon dioxide migration channel and the change of carbon dioxide flow along with time;
and step S70, generating an analysis gauge plate for inverting the form of the carbon dioxide migration channel and the carbon dioxide migration capacity based on the elastic characteristics of the carbon dioxide migration channel based on the relation among the elastic modulus, the form parameters and the carbon dioxide transportation capacity of the carbon dioxide migration channel.
2. The method for analyzing fluid conductance characteristics of a carbon dioxide transport channel according to claim 1, wherein the step S40 is performed by generating different random fracture surfaces through fourier transform by:
Figure 607546DEST_PATH_IMAGE008
wherein the content of the first and second substances,
Figure 849172DEST_PATH_IMAGE009
are data points
Figure 655454DEST_PATH_IMAGE010
The FFT represents the fourier transform.
3. The method for analyzing the fluid conductance characteristics of a carbon dioxide migration channel according to claim 1, wherein the curve of the change of the carbon dioxide flow rate with time is obtained by a random walk-around particle tracing method;
step S601, when the total number is
Figure 203984DEST_PATH_IMAGE011
When the carbon dioxide fluid unit passes through the carbon dioxide transport passage, the second
Figure 70309DEST_PATH_IMAGE012
The time required for each cell was:
Figure 64810DEST_PATH_IMAGE013
wherein the content of the first and second substances,
Figure 284570DEST_PATH_IMAGE014
is as follows
Figure 48127DEST_PATH_IMAGE012
A unit arrives at a node
Figure 819774DEST_PATH_IMAGE015
The time required for the operation of the apparatus,
Figure 301570DEST_PATH_IMAGE016
representing slave nodes
Figure 590600DEST_PATH_IMAGE015
Arriving node
Figure 208664DEST_PATH_IMAGE017
At any one point in time, the time of the day,
Figure 416791DEST_PATH_IMAGE018
is as follows
Figure 995671DEST_PATH_IMAGE012
A volume of carbon dioxide comprising a carbon dioxide fluid unit,
Figure 681867DEST_PATH_IMAGE019
as slave nodes in unit time
Figure 685595DEST_PATH_IMAGE017
Any point arrives at node
Figure 172946DEST_PATH_IMAGE015
The volume of carbon dioxide in the gas mixture,
Figure 363756DEST_PATH_IMAGE020
for the mesh width of the plane subdivision of the carbon dioxide migration channel,
Figure 853643DEST_PATH_IMAGE021
for nodes in carbon dioxide channels
Figure 446299DEST_PATH_IMAGE015
And node
Figure 871595DEST_PATH_IMAGE022
The geometric mean of the height cube,Lrepresentative node
Figure 549701DEST_PATH_IMAGE015
And node
Figure 843279DEST_PATH_IMAGE022
The distance between the two or more of the two or more,
Figure 165807DEST_PATH_IMAGE023
and
Figure 886639DEST_PATH_IMAGE024
is a node
Figure 52041DEST_PATH_IMAGE015
And node
Figure 883730DEST_PATH_IMAGE022
The pore pressure of (a) is set,
Figure 60765DEST_PATH_IMAGE025
is the viscosity coefficient of carbon dioxide;
step S602, obtaining a time-varying curve of the carbon dioxide flow rate of the carbon dioxide transport passage by using the time as a variable through the method of step S601.
4. The method of analyzing fluid transport properties of a carbon dioxide transport passageway according to claim 3, wherein the permeability of the carbon dioxide transport passageway is expressed as:
Figure 952498DEST_PATH_IMAGE026
wherein the content of the first and second substances,
Figure 870775DEST_PATH_IMAGE027
is the permeability of the carbon dioxide migration channel,
Figure 614478DEST_PATH_IMAGE028
representing the pressure difference across the channel driving carbon dioxide transport,
Figure 36232DEST_PATH_IMAGE029
and
Figure 98866DEST_PATH_IMAGE030
representing the length and width of the carbon dioxide migration channel.
5. The method for analyzing the fluid transport characteristics of a carbon dioxide transport passageway according to claim 4, wherein the surface of the carbon dioxide transport passageway is sectioned into a fine cylindrical system, and the elastic modulus of the carbon dioxide transport passageway is calculated by hooke's law, which is expressed as:
Figure 973281DEST_PATH_IMAGE031
wherein the content of the first and second substances,
Figure 287719DEST_PATH_IMAGE032
is the elastic modulus of the carbon dioxide migration channel,
Figure 829559DEST_PATH_IMAGE033
in order to realize the maximum length in the cylindrical system when the plane of the carbon dioxide migration channel is split,
Figure 63094DEST_PATH_IMAGE034
is the first in a cylinder system
Figure 300171DEST_PATH_IMAGE015
The amount of contraction in the length of each cylinder,
Figure 11775DEST_PATH_IMAGE035
the total number of cylinders in the cylinder system.
6. The method for analyzing the fluid conducting characteristics of a carbon dioxide transport passageway according to claim 1, wherein the step S70 includes:
step S71, generating a chart form based on the elastic modulus, the form parameters and the carbon dioxide transmission capacity of the carbon dioxide migration channel;
and step S72, reversely deducing the morphological parameters of the carbon dioxide migration channel and the carbon dioxide transportation capacity of the carbon dioxide migration channel based on the elastic modulus of the carbon dioxide migration channel according to the data trend in the chart form.
7. A fluid conductance characterization system for a carbon dioxide transport channel, the characterization system comprising the following modules:
the data acquisition module is configured to acquire morphological parameter characteristics forming a carbon dioxide migration channel based on crack outcrop, core observation data and logging data; the morphological parameter characteristics comprise fracture surface morphology, fracture length and fracture aperture;
the carbon dioxide migration channel feature extraction module is configured to perform statistical analysis on the crack surface morphology, the crack length and the crack aperture based on the morphological parameter features to obtain carbon dioxide migration channel features; the carbon dioxide migration channel characteristics comprise an affine function for representing the height level spatial distribution of the surface of the crack, a probability density function for representing that the crack has a rough surface and morphological parameters for representing the surface morphology of the crack; wherein the morphological parameters comprise self-affine coefficients, mean and standard deviation of height distribution;
a power spectrum function construction module configured to construct a power spectrum function of a carbon dioxide migration channel form based on the carbon dioxide migration channel characteristics:
Figure 408122DEST_PATH_IMAGE001
wherein the content of the first and second substances,
Figure 953503DEST_PATH_IMAGE002
as a function of the power spectrum of the carbon dioxide migration channel morphology,
Figure 802511DEST_PATH_IMAGE003
for transporting carbon dioxideThe wave number of the trace surface in the wavenumber domain,
Figure 317806DEST_PATH_IMAGE004
is the self-affine coefficient of the carbon dioxide migration channel form,
Figure 568659DEST_PATH_IMAGE005
respectively representing the mean and standard deviation of the height distribution of the carbon dioxide migration channel morphology,
Figure 517898DEST_PATH_IMAGE006
is a probability density function satisfied by the carbon dioxide migration channel shape,
Figure 854201DEST_PATH_IMAGE007
an auto-affine function representing the height level spatial distribution of the fracture surface;
the crack surface generation module is configured to set different self-affine coefficients, mean values and standard deviations of height distribution for the power spectrum function of the carbon dioxide migration channel form, and generate different random crack surfaces through Fourier transformation;
the migration channel building module is configured to form a closed crack space by taking the highest plane of the crack surface as a top interface of the crack, and obtain the carbon dioxide migration channel with different self-affine coefficients, mean values of height distribution and standard deviation as a carbon dioxide migration channel;
the migration channel parameter calculation module is configured to obtain the elastic modulus and the carbon dioxide transmission capacity of the carbon dioxide migration channel through numerical simulation by utilizing a local cubic law, a Navigneaux equation and a Hooke law; the carbon dioxide transportation capacity comprises a curve of the permeability of a carbon dioxide migration channel and the change of carbon dioxide flow along with time;
and the migration channel analysis gauge plate generation module is configured to generate an analysis gauge plate for inverting the carbon dioxide migration channel form and the carbon dioxide migration capacity based on the elastic characteristics of the carbon dioxide migration channel based on the relation among the elastic modulus, the form parameters and the carbon dioxide transportation capacity of the carbon dioxide migration channel.
8. An electronic device, comprising:
at least one processor; and
a memory communicatively coupled to at least one of the processors; wherein the content of the first and second substances,
the memory stores instructions executable by the processor for performing a method of fluid conductance profiling of a carbon dioxide transport channel according to any of claims 1-6.
9. A computer-readable storage medium storing computer instructions for execution by the computer to perform the method for fluid conductance characterization of a carbon dioxide transport channel of any one of claims 1-6.
CN202111620381.XA 2021-12-28 2021-12-28 Fluid conductance characteristic analysis method, system and equipment for carbon dioxide migration channel Active CN113987851B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111620381.XA CN113987851B (en) 2021-12-28 2021-12-28 Fluid conductance characteristic analysis method, system and equipment for carbon dioxide migration channel

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111620381.XA CN113987851B (en) 2021-12-28 2021-12-28 Fluid conductance characteristic analysis method, system and equipment for carbon dioxide migration channel

Publications (2)

Publication Number Publication Date
CN113987851A CN113987851A (en) 2022-01-28
CN113987851B true CN113987851B (en) 2022-04-08

Family

ID=79734704

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111620381.XA Active CN113987851B (en) 2021-12-28 2021-12-28 Fluid conductance characteristic analysis method, system and equipment for carbon dioxide migration channel

Country Status (1)

Country Link
CN (1) CN113987851B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107526909A (en) * 2016-06-21 2017-12-29 中国石油化工股份有限公司 Determine the method and system of fault transfer system permeability
CN109284571A (en) * 2018-10-19 2019-01-29 西南石油大学 A kind of multiple dimensioned multi- scenarios method seepage flow Mathematical Modeling Methods of carbon dioxide replacement shale gas
CN112800592A (en) * 2021-01-08 2021-05-14 中国石油大学(华东) Evaluation method for leakage risk in carbon dioxide geological storage body

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10571605B2 (en) * 2016-03-28 2020-02-25 Exxonmobil Upstream Research Company Method and system for modeling and simulating a fractured reservoir

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107526909A (en) * 2016-06-21 2017-12-29 中国石油化工股份有限公司 Determine the method and system of fault transfer system permeability
CN109284571A (en) * 2018-10-19 2019-01-29 西南石油大学 A kind of multiple dimensioned multi- scenarios method seepage flow Mathematical Modeling Methods of carbon dioxide replacement shale gas
CN112800592A (en) * 2021-01-08 2021-05-14 中国石油大学(华东) Evaluation method for leakage risk in carbon dioxide geological storage body

Also Published As

Publication number Publication date
CN113987851A (en) 2022-01-28

Similar Documents

Publication Publication Date Title
Buscombe Transferable wavelet method for grain‐size distribution from images of sediment surfaces and thin sections, and other natural granular patterns
CA2704547C (en) Fracture network characterization method
US7266041B1 (en) Multi-attribute background relative scanning of 3D geophysical datasets for locally anomaluous data points
Sebacher et al. A probabilistic parametrization for geological uncertainty estimation using the ensemble Kalman filter (EnKF)
Pérez et al. Verifying the high-order consistency of training images with data for multiple-point geostatistics
Park Modeling uncertainty in metric space
CN112462438B (en) Method, device and equipment for quantitatively evaluating formation fluid based on nuclear magnetic resonance logging
Yang et al. Noise‐assisted multivariate empirical mode decomposition of saturated hydraulic conductivity along a south‐north transect across the Loess Plateau of China
Zhang et al. Evaluation and error analysis: Kalman gain regularization versus covariance regularization
Liu et al. An integrated 3D fracture network reconstruction method based on microseismic events
Tzu-hao et al. Reservoir uncertainty quantification using probabilistic history matching workflow
Chang et al. History matching of stimulated reservoir volume of shale-gas reservoirs using an iterative ensemble smoother
Kang et al. Initial model selection for efficient history matching of channel reservoirs using Ensemble Smoother
CN114529110A (en) Lithofacies inversion method and system based on deep neural network model
CN106650042B (en) Method and device for determining fracture conductance capacity parameter
CN113987851B (en) Fluid conductance characteristic analysis method, system and equipment for carbon dioxide migration channel
Amanbek et al. Framework of fracture network modeling using conditioned data with sequential gaussian simulation
Udegbe et al. Big Data Analytics for Seismic Fracture Identification, Using Amplitude-Based Statistics
Basu et al. Multi-Start method for reservoir model uncertainty quantification with application to robust decision-making
Kracht et al. Spatial representativeness evaluation of air quality monitoring sites by point-centred variography
Barros et al. Automated clustering based scenario reduction to accelerate robust life-cycle optimization
US20210349238A1 (en) Systems and methods for generating vertical and lateral heterogeneity indices of reservoirs
Koppe et al. A chart for judging optimal sample spacing for ore grade estimation
Ahmadi et al. A sensitivity study of FILTERSIM algorithm when applied to DFN modeling
Nejadi et al. Ensemble Kalman filter predictor bias correction method for non-Gaussian geological facies detection

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant