CN112444865B - Method, system, device and storage medium for exploring fracture-cave reservoir - Google Patents
Method, system, device and storage medium for exploring fracture-cave reservoir Download PDFInfo
- Publication number
- CN112444865B CN112444865B CN201910804363.3A CN201910804363A CN112444865B CN 112444865 B CN112444865 B CN 112444865B CN 201910804363 A CN201910804363 A CN 201910804363A CN 112444865 B CN112444865 B CN 112444865B
- Authority
- CN
- China
- Prior art keywords
- random
- function
- dimensional
- fracture
- coordinate system
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 60
- 238000005311 autocorrelation function Methods 0.000 claims abstract description 44
- 238000001228 spectrum Methods 0.000 claims abstract description 39
- 230000009466 transformation Effects 0.000 claims abstract description 7
- 230000006870 function Effects 0.000 claims description 77
- 230000008569 process Effects 0.000 claims description 13
- 238000010276 construction Methods 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 10
- 230000003595 spectral effect Effects 0.000 claims 1
- 238000011160 research Methods 0.000 abstract description 9
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000004891 communication Methods 0.000 description 7
- 230000004044 response Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 4
- 230000005236 sound signal Effects 0.000 description 4
- 238000005553 drilling Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000009918 complex formation Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 201000001098 delayed sleep phase syndrome Diseases 0.000 description 1
- 208000033921 delayed sleep phase type circadian rhythm sleep disease Diseases 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 235000019800 disodium phosphate Nutrition 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000004549 pulsed laser deposition Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The disclosure relates to the technical field of seismic exploration, in particular to a method, a device and a storage medium for exploration of a fracture-cave reservoir, which are used for solving the technical problem that the seismic exploration research lacks three-dimensional random medium modeling with directivity characteristics in the related technology. The method for exploring the fracture-cavity reservoir comprises the following steps: constructing an exponential elliptic autocorrelation function under a three-dimensional coordinate system; carrying out Fourier transformation on the constructed exponential elliptic autocorrelation function to obtain a power spectrum function of a random disturbance function; according to the power spectrum function, obtaining a random disturbance function and normalizing; constructing a random medium model according to the normalized random disturbance function and the background medium parameters of the random medium in the fracture-cavity reservoir; inputting background medium parameters of random mediums in a fracture-cavity reservoir to be explored into the constructed three-dimensional random medium model to obtain a random medium model corresponding to the fracture-cavity reservoir to be explored; and (5) exploration of the fracture-cavity reservoir stratum is carried out by using the random medium model.
Description
Technical Field
The present disclosure relates to the field of seismic exploration, and in particular, to a method, apparatus, and storage medium for exploration of a fracture-cave reservoir.
Background
The oil reservoir (or oil reservoir) geophysics (Reservoir Geophysics) is widely paid attention to by geophysicists at home and abroad from the 80 s of the last century, and becomes a hot research field. The geophysical basis of the oil reservoir is a non-uniform medium model, and the main content is to finely research the characteristics, parameters and oil and gas reservoirs of the oil reservoir (or oil reservoir), thereby serving the oil and gas exploration, in particular the oil and gas development. Because the reservoir geophysics are greatly connected with complex and fine heterogeneous media, the reservoir geophysics are difficult to completely describe by a conventional method, and a flexible, convenient forward model capable of completely describing the reservoir heterogeneity is proposed, so that a random medium model theory is generated.
Disclosure of Invention
The disclosure provides a method, a device and a storage medium for exploring a fracture-cave reservoir, which are used for solving the technical problem that the seismic exploration research lacks three-dimensional random medium modeling with directivity characteristics in the related technology.
To achieve the above object, a first aspect of embodiments of the present disclosure provides a method of exploring a fracture-cave reservoir, the method comprising:
Constructing an exponential elliptic autocorrelation function under a three-dimensional coordinate system;
Carrying out Fourier transformation on the constructed exponential elliptic autocorrelation function to obtain a power spectrum function of a random disturbance function;
According to the power spectrum function, a random disturbance function is obtained and normalized;
constructing a three-dimensional random medium model according to the normalized random disturbance function and background medium parameters of random medium in the fracture-cavity reservoir;
Inputting background medium parameters of random media in a fracture-cavity reservoir to be explored into the constructed three-dimensional random medium model to obtain a random medium model corresponding to the fracture-cavity reservoir to be explored;
and exploration of the fracture-cavity reservoir stratum is conducted by utilizing the random medium model.
Optionally, the constructing an exponential elliptic autocorrelation function in the three-dimensional coordinate system includes:
Selecting three local autocorrelation lengths of points in a three-dimensional global coordinate system;
Selecting three rotation angles corresponding to the three local autocorrelation lengths according to the three-dimensional local coordinate system and the space offset between points in the three-dimensional global coordinate system;
And constructing an exponential elliptic autocorrelation function according to the three local autocorrelation lengths and the three rotation angles.
Optionally, setting a (x, y, z), b (x, y, z), c (x, y, z) as three local autocorrelation lengths of (x, y, z) points in the three-dimensional global coordinate system ozz, respectively; according to the three local autocorrelation lengths and the three rotation angles, the constructed exponential elliptic autocorrelation function has the following expression:
In the formula (1): a (x, y, z), b (x, y, z), c (x, y, z) are the local autocorrelation function and the local autocorrelation length of the three-dimensional global coordinate system ozz points (x, y, z), respectively; x '1,y′1,z′1 is the spatial offset relative to the point (X, y, z) in the three-dimensional local coordinate system O' X 1'Y1'Z1 ', the three-dimensional local coordinate system O' X 1'Y1'Z1 'is obtained by rotating the three-dimensional local coordinate system O' X 1Y1 Z1 by an angle α, β, γ in the counterclockwise direction, and the expression of the spatial offset is:
in the formula (2): d is a space offset, and alpha, beta and gamma are rotation angles corresponding to a (x, y, z), b (x, y, z) and c (x, y, z).
Optionally, the formula for performing fourier transform on the constructed exponential elliptic autocorrelation function is:
In the formula (3): Γ denotes a spatial fourier transform operator, phi (k x,ky,kz) denotes a power spectrum function of a random disturbance function, k is a spatial wavenumber, kx is a spatial wavenumber in the x-axis direction, ky is a spatial wavenumber in the y-axis direction, and kz is a spatial wavenumber in the z-axis direction.
Optionally, the obtaining a random disturbance function according to the power spectrum function includes:
generating a corresponding random spectrum sequence according to a random process method by a power spectrum function;
and (3) applying a spectrum expansion formula of a random process to obtain a random disturbance function described by the constructed exponential elliptic autocorrelation function.
Optionally, normalizing the random perturbation function includes:
And normalizing the random disturbance function according to the requirements of the random disturbance mean and variance.
Optionally, according to the normalized random disturbance function and the background medium parameters of the random medium, the expression of the constructed three-dimensional random medium model is as follows:
in the formula (3): ρ is density, λ, μ are pull Mei Canshu, ρ 0、λ0、μ0 is the background medium parameter of the random medium; δρ, δλ, δμ are the non-uniform disturbance amounts.
In a second aspect of embodiments of the present disclosure, there is provided a system for exploration of a fracture-cave reservoir, the system comprising:
A first construction module configured to construct an exponential elliptic autocorrelation function in a three-dimensional coordinate system;
The Fourier transform module is configured to obtain a power spectrum function of a random disturbance function by carrying out Fourier transform on the constructed exponential elliptic autocorrelation function;
The first obtaining module is configured to obtain a random disturbance function according to the power spectrum function and normalize the random disturbance function;
The second construction module is configured to construct a three-dimensional random medium model according to the normalized random disturbance function and background medium parameters of random media in the fracture-cavity reservoir;
the second obtaining module is configured to input background medium parameters of random mediums in the fracture-cavity reservoir to be explored in the constructed three-dimensional random medium model so as to obtain a random medium model corresponding to the fracture-cavity reservoir to be explored, and exploration of the fracture-cavity reservoir is conducted by using the random medium model.
A third aspect of the disclosed embodiments provides a computer readable storage medium having stored thereon a computer program which when executed by a processor performs the steps of the method of any of the first aspects above.
In a fourth aspect of embodiments of the present disclosure, there is provided an apparatus for exploration of a fracture-cave reservoir, comprising:
A memory having a computer program stored thereon; and
A processor for executing the computer program in the memory to implement the steps of the method of any of the above first aspects.
By adopting the technical scheme, at least the following technical effects can be achieved:
The method and the device develop three-dimensional random medium model construction with direction characteristics based on a random medium modeling theory, represent different angles of complex stratum in stratum by introducing autocorrelation angles, and perform three-dimensional modeling on the fracture-cavity reservoir, so that the constructed model is more in line with actual fracture-cavity reservoir geology, model data are provided for researching a fracture-cavity seismic response mode and mechanism, calibration and identification are better performed on the seismic response characteristics of the actual deep micro-amplitude fracture-cavity reservoir in the wild, the accuracy of field drilling is improved, and support is provided for oil and gas yield increase.
Additional features and advantages of the present disclosure will be set forth in the detailed description which follows.
Drawings
The accompanying drawings are included to provide a further understanding of the disclosure, and are incorporated in and constitute a part of this specification, illustrate the disclosure and together with the description serve to explain, but do not limit the disclosure. In the drawings:
FIG. 1 is a flow chart of a method of exploring a fracture-cave reservoir, as shown in an exemplary embodiment of the present disclosure.
FIG. 2 is a flow chart illustrating a method of exploring a fracture-cave reservoir including constructing an exponential elliptic autocorrelation function in steps according to an exemplary embodiment of the present disclosure.
FIG. 3 is a schematic diagram of a relationship between a three-dimensional local coordinate system and a three-dimensional global coordinate system, as shown in an exemplary embodiment of the present disclosure.
Fig. 4 is a schematic diagram of a relationship between two three-dimensional local coordinate systems shown in an exemplary embodiment of the present disclosure.
FIG. 5 is a flowchart showing a method of exploring a fracture-cave reservoir including obtaining a random disturbance function in steps according to an exemplary embodiment of the present disclosure.
Fig. 6 is a schematic view of a three-dimensional argument shown in an exemplary embodiment of the present disclosure.
FIG. 7 is a schematic diagram of a three-dimensional random media model according to an exemplary embodiment of the present disclosure.
FIG. 8 is a system block diagram of an exploration fracture-cave reservoir as shown in an exemplary embodiment of the present disclosure.
FIG. 9 is a block diagram of an apparatus for exploration of a fracture-cave reservoir, as shown in an exemplary embodiment of the present disclosure.
Detailed Description
The embodiments of the present disclosure will be described in detail below with reference to the drawings and examples, so as to solve the technical problem by applying technical means to the present disclosure, and the implementation process for achieving the corresponding technical effects can be fully understood and implemented accordingly. The embodiment of the application and the features in the embodiment can be mutually combined on the premise of no conflict, and the formed technical scheme is within the protection scope of the disclosure.
The inventor of the present disclosure has found through researches that, in the related art, random medium modeling is mainly based on two-dimensional conditions, to reasonably characterize a fracture-cavity reservoir and seismic response characteristics should be performed under three-dimensional conditions, researches on three-dimensional random medium modeling techniques are relatively few, and researches on three-dimensional random medium modeling involving different inclination factors of the fracture-cavity reservoir are not found, so the present disclosure provides a method for exploring the fracture-cavity reservoir, and lays a foundation for researches on seismic response of an actual reservoir.
Example 1
FIG. 1 is a flow chart of a method for exploration of a fracture-cave reservoir shown in an exemplary embodiment of the present disclosure to solve the technical problem of lack of three-dimensional stochastic media modeling with directional characteristics in seismic exploration research in the related art. As shown in fig. 1, the method for exploring a fracture-cave reservoir shown in this embodiment may include the following steps:
S11, constructing an exponential elliptic autocorrelation function under a three-dimensional coordinate system.
S12, carrying out Fourier transformation on the constructed exponential elliptic autocorrelation function to obtain a power spectrum function of the random disturbance function.
S13, according to the power spectrum function, obtaining a random disturbance function and normalizing.
S14, constructing a three-dimensional random medium model according to the normalized random disturbance function and background medium parameters of random media in the fracture-cavity reservoir.
S15, inputting background medium parameters of random media in a fracture-cavity reservoir to be explored into the constructed three-dimensional random medium model to obtain a random medium model corresponding to the fracture-cavity reservoir to be explored;
s16, exploration of the fracture-cavity reservoir stratum is conducted by using the random medium model.
As shown in fig. 2, step S11, i.e. constructing an exponential elliptic autocorrelation function in a three-dimensional coordinate system, may include the steps of:
s111, selecting three local autocorrelation lengths of points in a three-dimensional global coordinate system.
S112, selecting three rotation angles corresponding to the three local autocorrelation lengths according to the three-dimensional local coordinate system and the space offset between points in the three-dimensional global coordinate system.
S113, constructing an exponential elliptic autocorrelation function according to the three local autocorrelation lengths and the three rotation angles.
In step S111, it is assumed that three local autocorrelation lengths at the midpoint (x, y, z) of the three-dimensional global coordinate system ozz are a (x, y, z), b (x, y, z), c (x, y, z). X ' 1,y′1,z′1 is the spatial offset relative to the (X, y, z) point in the three-dimensional local coordinate system O ' X 1'Y1'Z1 ', and the relationship between the three-dimensional local coordinate system O ' X 1'Y1'Z1 ' and the three-dimensional global coordinate system ozz is shown in fig. 3, and the relationship between the two can be obtained according to the coordinate conversion formula (5).
Equation (5) D is the spatial offset. The coordinate transformation matrix can be further constructed by a coordinate rotation method as shown in formula (2).
As shown in fig. 4, the three-dimensional local coordinate system O ' X 1'Y1'Z1 ' is obtained by rotating the three-dimensional local coordinate system O ' X 1Y1 Z1 by an angle α, β, γ along the counterclockwise direction of the respective coordinates.
Further, according to the selected autocorrelation lengths a (x, y, z), b (x, y, z), c (x, y, z), and rotation angles α (x, y, z), β (x, y, z), γ (x, y, z), the expression of the constructed exponential elliptic autocorrelation function may be:
In the formula (1): a (x, y, z), b (x, y, z), c (x, y, z) are the local autocorrelation function and the local autocorrelation length of the three-dimensional global coordinate system ozz points (x, y, z), respectively.
After the exponential elliptic autocorrelation function is constructed, step S12 may be executed to obtain a power spectrum function of the random disturbance function by performing fourier transform on the constructed exponential elliptic autocorrelation function.
The constructed exponential elliptic autocorrelation functionThe formula for performing the fourier transform may be:
In the formula (3): Γ denotes a spatial fourier transform operator, phi (k x,ky,kz) denotes a power spectrum function of a random disturbance function, k is a spatial wavenumber, kx is a spatial wavenumber in the x-axis direction, ky is a spatial wavenumber in the y-axis direction, and kz is a spatial wavenumber in the z-axis direction. The power spectrum function phi (k x,ky,kz) is the power spectrum of the random disturbance sigma (x, y, z).
After the power spectrum function is obtained, step S13 may be executed, and according to the power spectrum function, a random disturbance function is obtained and normalized. As shown in fig. 5, the obtaining a random disturbance function according to the power spectrum function may include the following steps:
S131, generating a corresponding random spectrum sequence according to a random process method by a power spectrum function phi (k x,ky,kz).
S132, applying a spectrum expansion formula of a random process to obtain the constructed exponential-elliptical autocorrelation functionThe random perturbation function σ (x, y, z) is described. Wherein, the random spectrum sequence and the random disturbance function sigma (x, y, z) are discrete data.
After the random disturbance function sigma (x, y, z) is obtained, the random disturbance function sigma (x, y, z) is normalized according to the requirements of the random disturbance mean and variance. The mean value of the normalized random disturbance function value is equal to 0, and the variance thereof is a certain value (generally, the disturbance value is obtained by a percentage of the background value disturbance, for example, 0.01, then the set disturbance value satisfies the mean value being zero, and the variance is 0.01), and the normalized process is a known technology in mathematical aspect, and is not specifically developed herein for the sake of brevity of the description.
After normalizing the obtained random disturbance function sigma (x, y, z), executing step S14, and constructing a random medium model according to the normalized random disturbance function and the background medium parameters of the random medium in the fracture-cavity reservoir.
In the related art, a random medium can be understood as a theoretical model with small-scale non-uniform disturbance distributed in a large-scale uniform background, and the construction theory can be described as the following expression:
m(X)＝m0+δm(X)＝m0(1+σ(X)) (6)
In formula (6): m 0 represents a large-scale homogeneous medium; δm (X) represents a small-scale non-uniform perturbation medium; sigma (X) is the random disturbance characteristic of the medium; m (X) is a random medium; x= (X, y, z) is a spatial position vector.
The isotropic elastic medium may be determined by the density ρ and the pull Mei Canshu λ, μ, when m (X) contains ρ, λ, μ, formula (6) may be expressed as:
in the formula (7): ρ 0、λ0、μ0 is the background medium parameter of the random medium, assumed to be constant or slowly varying with the spatial coordinates (x, y, z); δρ, δλ, δμ are the amount of non-uniform disturbance added to the background, and assume a spatially smooth random process with zero mean, constant variance, and some autocorrelation function.
Since equations (6) and (7) do not take into account directional characteristics such as different dip factors for fracture-cave reservoirs, the random medium model for complex formation formations is not accurate. The present disclosure adds directional characteristics in the modeling process, and after normalization of the obtained random disturbance function sigma (x, y, z), the desired structure can be obtainedAs an autocorrelation function, a random medium model with a specified mean and variance, the expression of the constructed three-dimensional random medium model is:
In the formula (4): ρ is density, λ, μ are pull Mei Canshu, ρ 0、λ0、μ0 is the background medium parameter of the random medium in the fracture-cavity reservoir; δρ, δλ, δμ are the non-uniform disturbance amounts.
After the three-dimensional random medium model is built, step S15 is executed, and background medium parameters of random medium in the fracture-cavity reservoir to be explored are input into the built three-dimensional random medium model so as to obtain a random medium model corresponding to the fracture-cavity reservoir to be explored, and further, the fracture-cavity reservoir is explored by utilizing the random medium model. The method has the advantages that the actual medium condition contains the non-uniform mass, so that the speed of the underground medium is represented as certain randomness, the non-uniform mass is distributed in various spaces of the underground medium, the random medium model under different conditions can be reasonably constructed, the wave field simulation is carried out through the models of different random medium shapes, the section wave field characteristics are analyzed, the field actual seismic exploration is guided by utilizing indoor analysis results and experience, the calibration and identification are carried out on the field actual deep micro-amplitude fracture-cavity reservoir seismic response characteristics, the accuracy of field drilling is improved, and the data support is provided for oil and gas production.
The entire set-up process of the three-dimensional stochastic model with directional features is shown below by way of an example. The following are three-dimensional stochastic model parameters in the example:
model size: nx=100, ny=100, nz=100;
grid spacing: dx=1 m, dy=1 m, dz=1 m;
Speed average: v 0＝5500m/s2;
standard deviation of velocity disturbance: 10%;
Autocorrelation length: a=3m, b=10m, c=20m; autocorrelation angle: α=0°, β=0°, γ=75°;
The exponential elliptic autocorrelation function is carried out by adopting a formula (1), firstly, according to the given model range, the amplitude angle required by the Fourier transformation of the exponential elliptic function is completed, the amplitude angle range is between 0 pi and 2 pi and is in random characteristic change, a two-dimensional section display of the cut-out is shown as shown in figure 6, the realization of a disturbance power spectrum is completed by combining a transformation function, the calculation of disturbance quantity is further completed according to a disturbance mean value and a variance value, and a final three-dimensional random medium model containing autocorrelation angles is completed according to a formula (4) by combining a background velocity field, as shown in figure 7. By using the three-dimensional random medium model constructed by the method shown in the figure 7, the background medium parameters of random medium in the fracture-cavity reservoir to be explored can be input into the model, and the random medium model corresponding to the fracture-cavity reservoir to be explored can be obtained.
According to the method, based on a random medium modeling theory, three-dimensional random medium model construction with direction characteristics is developed, different angles of a fracture-cavity reservoir in a stratum are represented by introducing an autocorrelation angle, and the fracture-cavity reservoir is subjected to three-dimensional modeling, so that the constructed model is more in line with actual geological conditions, model data are provided for researching a fracture-cavity seismic response mode and mechanism, calibration and identification are better carried out on the seismic response characteristics of an actual deep micro-amplitude fracture-cavity reservoir in the field, the accuracy of field drilling is improved, and support is provided for oil and gas yield increase.
It should be noted that, for simplicity of description, the method embodiment shown in fig. 1 is depicted as a series of acts, but it should be understood by those skilled in the art that the present disclosure is not limited by the order of acts depicted. Further, those skilled in the art will also appreciate that the embodiments described in the specification are all preferred embodiments and that the acts referred to are not necessarily required by the present disclosure.
Example two
FIG. 8 is a system for exploration of a fracture-cave reservoir, as shown in FIG. 8, according to an exemplary embodiment of the present disclosure, the system 300 for exploration of a fracture-cave reservoir comprising:
a first construction module 310 configured to construct an exponential elliptic autocorrelation function in a three-dimensional coordinate system;
A fourier transform module 320 configured to obtain a power spectrum function of a random disturbance function by performing fourier transform on the constructed exponential-elliptical autocorrelation function;
A first obtaining module 330 configured to obtain a random disturbance function and normalize according to the power spectrum function;
a second construction module 340, configured to construct a random medium model according to the normalized random disturbance function and the background medium parameters of the random medium in the fracture-cavity reservoir;
A second obtaining module 350, configured to input background medium parameters of random media in the fracture-cavity reservoir to be explored in the constructed three-dimensional random medium model, so as to obtain a random medium model corresponding to the fracture-cavity reservoir to be explored, so as to explore the fracture-cavity reservoir by using the random medium model.
Optionally, the first building module 310 is further configured to: selecting three local autocorrelation lengths of points in a three-dimensional global coordinate system; selecting three rotation angles corresponding to the three local autocorrelation lengths according to the three-dimensional local coordinate system and the space offset between points in the three-dimensional global coordinate system; and constructing an exponential elliptic autocorrelation function according to the three local autocorrelation lengths and the three rotation angles.
Optionally, setting a (x, y, z), b (x, y, z), c (x, y, z) as three local autocorrelation lengths of (x, y, z) points in the three-dimensional global coordinate system ozz, respectively; the expression of the exponential elliptic autocorrelation function constructed by the first construction module 310 is as follows:
In the formula (1): a (x, y, z), b (x, y, z), c (x, y, z) are the local autocorrelation function and the local autocorrelation length of the three-dimensional global coordinate system ozz points (x, y, z), respectively; x' 1,y′1,z′1 is the spatial offset relative to the point (x, y, z) in the three-dimensional local coordinate system O 1X′1Y′1Z′1, the three-dimensional local coordinate system O 1X′1Y′1Z′1 is obtained by rotating the three-dimensional local coordinate system O 1X1Y1Z1 by an angle α, β, γ along the counterclockwise direction of the respective coordinates, and the expression of the spatial offset is:
in the formula (2): d is a space offset, and alpha, beta and gamma are rotation angles corresponding to a (x, y, z), b (x, y, z) and c (x, y, z).
Alternatively, the fourier transform module 320 performs fourier transform as follows:
In the formula (3): Γ denotes a spatial fourier transform operator, phi (k x,ky,kz) denotes a power spectrum function of a random disturbance function, k is a spatial wavenumber, kx is a spatial wavenumber in the x-axis direction, ky is a spatial wavenumber in the y-axis direction, and kz is a spatial wavenumber in the z-axis direction.
Optionally, the obtaining module 340 is further configured to: generating a corresponding random spectrum sequence according to a random process method by a power spectrum function; and (3) applying a spectrum expansion formula of a random process to obtain a random disturbance function described by the constructed exponential elliptic autocorrelation function.
Optionally, the obtaining module 340 is further configured to include: and normalizing the random disturbance function according to the requirements of the random disturbance mean and variance.
Optionally, the expression of the random media model constructed by the second construction module 350 is:
In the formula (3): ρ 0、λ0、μ0 is the background medium parameter of the random medium; δρ, δλ, δμ are the non-uniform disturbance amounts.
The specific manner in which the various modules perform the operations in relation to the systems of the above embodiments have been described in detail in relation to the embodiments of the method and will not be described in detail herein.
Example III
The present disclosure also provides a computer readable storage medium having stored thereon a computer program which when executed by a processor performs the method steps of exploring a fracture-cave reservoir as set forth in any one of the alternative embodiments described above.
The method implemented when the computer program of the method for exploring a fracture-cave reservoir running on the processor is executed may refer to a specific embodiment of the method for exploring a fracture-cave reservoir of the present disclosure, and will not be described herein.
The processor may be an integrated circuit chip having information processing capabilities. The processor may be a general purpose processor including a central processing unit (Central Processing Unit, CPU), a network processor (Network Processor, NP), etc.
Example IV
The present disclosure also provides an apparatus for exploration of a fracture-cave reservoir, comprising:
A memory having a computer program stored thereon; and
A processor for executing the computer program in the memory to perform the method steps of exploring a fracture-cave reservoir as set forth in any one of the alternative embodiments above.
FIG. 9 is a block diagram illustrating an apparatus 400 for exploration of a fracture-cave reservoir, according to an example embodiment. As shown in fig. 9, the apparatus 400 may include: a processor 401, a memory 402, a multimedia component 403, an input/output (I/O) interface 404, and a communication component 405.
Wherein the processor 401 is configured to control the overall operation of the apparatus 400 to perform all or part of the steps in the method for constructing a three-dimensional random media model described above. The memory 402 is used to store various types of data to support operations at the device 400, which may include, for example, instructions for any application or method operating on the device 400, as well as application-related data. The Memory 402 may be implemented by any type or combination of volatile or non-volatile Memory devices, such as static random access Memory (Static Random Access Memory, SRAM for short), electrically erasable programmable Read-Only Memory (ELECTRICALLY ERASABLE PROGRAMMABLE READ-Only Memory, EEPROM for short), erasable programmable Read-Only Memory (Erasable Programmable Read-Only Memory, EPROM for short), programmable Read-Only Memory (Programmable Read-Only Memory, PROM for short), read-Only Memory (ROM for short), magnetic Memory, flash Memory, magnetic disk, or optical disk. The multimedia component 403 may include a screen and an audio component. Wherein the screen may be, for example, a touch screen, the audio component being for outputting and/or inputting audio signals. For example, the audio component may include a microphone for receiving external audio signals. The received audio signal may be further stored in the memory 402 or transmitted through the communication component 405. The audio assembly further comprises at least one speaker for outputting audio signals. The I/O interface 404 provides an interface between the processor 401 and other interface modules, which may be a keyboard, mouse, buttons, etc. These buttons may be virtual buttons or physical buttons. The communication component 405 is used for wired or wireless communication between the apparatus 400 and other devices. Wireless Communication, such as Wi-Fi, bluetooth, near field Communication (NFC for short), 2G, 3G, or 4G, or a combination of one or more thereof, and accordingly the Communication component 405 may comprise: wi-Fi module, bluetooth module, NFC module.
In an exemplary embodiment, the apparatus 400 may be implemented by one or more Application-specific integrated circuits (ASICs), digital signal processors (DIGITAL SIGNAL processors, DSPs), digital signal processing devices (DIGITAL SIGNAL Processing Device, DSPDs), programmable logic devices (Programmable Logic Device, PLDs), field programmable gate arrays (Field Programmable GATE ARRAY, FPGAs), controllers, microcontrollers, microprocessors, or other electronic components for performing the methods of exploring hole reservoirs described above.
In another exemplary embodiment, a computer readable storage medium is also provided, such as a memory 402, comprising program instructions executable by the processor 401 of the apparatus 400 to perform the method of exploring a hole reservoir described above.
The preferred embodiments of the present disclosure have been described in detail above with reference to the accompanying drawings, but the present disclosure is not limited to the specific details of the embodiments described above, and various simple modifications may be made to the technical solutions of the present disclosure within the scope of the technical concept of the present disclosure, and all the simple modifications belong to the protection scope of the present disclosure.
In addition, the specific features described in the above embodiments may be combined in any suitable manner without contradiction. The various possible combinations are not described further in this disclosure in order to avoid unnecessary repetition.
Moreover, any combination between the various embodiments of the present disclosure is possible as long as it does not depart from the spirit of the present disclosure, which should also be construed as the disclosure of the present disclosure.
Claims (8)
1. A method of exploring a fracture-cave reservoir, the method comprising:
Constructing an exponential elliptic autocorrelation function under a three-dimensional coordinate system; the method specifically comprises the following steps: selecting three local autocorrelation lengths of points in a three-dimensional global coordinate system; selecting three rotation angles corresponding to the three local autocorrelation lengths according to the three-dimensional local coordinate system and the space offset between points in the three-dimensional global coordinate system; constructing an exponential elliptic autocorrelation function according to the three local autocorrelation lengths and the three rotation angles; setting up Three-dimensional global coordinate system/>, respectivelyMiddle/>Three local autocorrelation lengths of a point; according to the three local autocorrelation lengths and the three rotation angles, the constructed exponential elliptic autocorrelation function has the following expression:
（1）
In the formula (1): ，/> Respectively three-dimensional global coordinate system Point/>A local autocorrelation function and a local autocorrelation length; /(I)Is a three-dimensional local coordinate systemMiddle relative to the point/>Spatial offset of (3) three-dimensional local coordinate system/>Is composed of three-dimensional local coordinate systemRotation angle in counter clockwise direction/>The expression of the spatial offset is:
（2）
in the formula (2): d is the amount of the spatial offset, For/>A corresponding rotation angle;
Carrying out Fourier transformation on the constructed exponential elliptic autocorrelation function to obtain a power spectrum function of a random disturbance function;
According to the power spectrum function, a random disturbance function is obtained and normalized;
constructing a three-dimensional random medium model according to the normalized random disturbance function and background medium parameters of random medium in the fracture-cavity reservoir;
Inputting background medium parameters of random media in a fracture-cavity reservoir to be explored into the constructed three-dimensional random medium model to obtain a random medium model corresponding to the fracture-cavity reservoir to be explored;
and exploration of the fracture-cavity reservoir stratum is conducted by utilizing the random medium model.
2. The method of claim 1, wherein the constructed exponential elliptic autocorrelation function is fourier transformed according to the formula:
（3）
In the formula (3): representing a spatial Fourier transform operator,/> The power spectrum function representing the random disturbance function, k is the spatial wave number, kx is the spatial wave number in the x-coordinate axis direction, ky is the spatial wave number in the y-coordinate axis direction, and kz is the spatial wave number in the z-coordinate axis direction.
3. The method of claim 2, wherein said obtaining a random perturbation function from said power spectral function comprises:
generating a corresponding random spectrum sequence according to a random process method by a power spectrum function;
and (3) applying a spectrum expansion formula of a random process to obtain a random disturbance function described by the constructed exponential elliptic autocorrelation function.
4. A method according to claim 3, wherein normalizing the random perturbation function comprises:
And normalizing the random disturbance function according to the requirements of the random disturbance mean and variance.
5. The method of claim 1, wherein the expression of the three-dimensional random medium model constructed according to the normalized random disturbance function and the background medium parameters of the random medium is:
（4）
In the formula (4): for density/> For pull Mei Canshu,/>Background medium parameters which are random media in a fracture-cavity reservoir; /(I)Is the non-uniform disturbance quantity.
6. A system for exploration of a fracture-cave reservoir, comprising:
A first construction module configured to construct an exponential elliptic autocorrelation function in a three-dimensional coordinate system; the method specifically comprises the following steps: selecting three local autocorrelation lengths of points in a three-dimensional global coordinate system; selecting three rotation angles corresponding to the three local autocorrelation lengths according to the three-dimensional local coordinate system and the space offset between points in the three-dimensional global coordinate system; constructing an exponential elliptic autocorrelation function according to the three local autocorrelation lengths and the three rotation angles; setting up Three-dimensional global coordinate system/>, respectivelyMiddle/>Three local autocorrelation lengths of a point; according to the three local autocorrelation lengths and the three rotation angles, the constructed exponential elliptic autocorrelation function has the following expression:
（1）
In the formula (1): ，/> Respectively three-dimensional global coordinate system Point/>A local autocorrelation function and a local autocorrelation length; /(I)Is a three-dimensional local coordinate systemMiddle relative to the point/>Spatial offset of (3) three-dimensional local coordinate system/>Is composed of three-dimensional local coordinate systemRotation angle in counter clockwise direction/>The expression of the spatial offset is:
（2）
in the formula (2): d is the amount of the spatial offset, For/>A corresponding rotation angle;
The Fourier transform module is configured to obtain a power spectrum function of a random disturbance function by carrying out Fourier transform on the constructed exponential elliptic autocorrelation function;
The first obtaining module is configured to obtain a random disturbance function according to the power spectrum function and normalize the random disturbance function;
The second construction module is configured to construct a three-dimensional random medium model according to the normalized random disturbance function and background medium parameters of random media in the fracture-cavity reservoir;
the second obtaining module is configured to input background medium parameters of random mediums in the fracture-cavity reservoir to be explored in the constructed three-dimensional random medium model so as to obtain a random medium model corresponding to the fracture-cavity reservoir to be explored, and exploration of the fracture-cavity reservoir is conducted by using the random medium model.
7. An apparatus for exploration of a fracture-cave reservoir, comprising:
A memory having a computer program stored thereon; and
A processor for executing the computer program in the memory to implement the steps of the method of any one of claims 1 to 5.
8. A computer readable storage medium, on which a computer program is stored, characterized in that the program, when being executed by a processor, implements the steps of the method according to any one of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910804363.3A CN112444865B (en) | 2019-08-28 | 2019-08-28 | Method, system, device and storage medium for exploring fracture-cave reservoir |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910804363.3A CN112444865B (en) | 2019-08-28 | 2019-08-28 | Method, system, device and storage medium for exploring fracture-cave reservoir |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112444865A CN112444865A (en) | 2021-03-05 |
| CN112444865B true CN112444865B (en) | 2024-05-17 |
Family
ID=74741452
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910804363.3A Active CN112444865B (en) | 2019-08-28 | 2019-08-28 | Method, system, device and storage medium for exploring fracture-cave reservoir |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112444865B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115267892A (en) * | 2022-07-21 | 2022-11-01 | 北京化工大学 | Organic-matter-rich shale seismic response numerical simulation method |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004090575A1 (en) * | 2003-04-09 | 2004-10-21 | Norsar | Method for simulating local prestack depth migrated seismic images |
| EA200702257A2 (en) * | 2007-11-15 | 2008-02-28 | Блаас Холдингс Лимитед | METHOD OF RADAR SENSING OF GROUND NEDR AND A DEVICE FOR ITS IMPLEMENTATION - A COMPLEX OF GEORADICAL LOCATION EXPLORATION |
| CN102253415A (en) * | 2011-04-19 | 2011-11-23 | 中国石油大学(华东) | Method for establishing earthquake response mode based on fracture equivalent medium model |
| WO2014003597A1 (en) * | 2012-06-26 | 2014-01-03 | Schlumberger, Holding Limited | A method for determining pore volume characteristics and porous materials' matrix thermal conductivity |
| CN103792573A (en) * | 2012-10-26 | 2014-05-14 | 中国石油化工股份有限公司 | Seismic wave impedance inversion method based on frequency spectrum fusion |
| CN104077738A (en) * | 2013-12-30 | 2014-10-01 | 辽宁师范大学 | Color image watermarking method based on local histogram characteristics |
| CN104977606A (en) * | 2014-04-02 | 2015-10-14 | 中国石油化工股份有限公司 | Method for establishing fracture-vuggy reservoir seismic numerical model |
| CN104992029A (en) * | 2015-07-20 | 2015-10-21 | 中国科学院国家天文台 | Modeling method for multi-scale non-uniform discrete random medium in lunar soil layer |
| CN105354421A (en) * | 2015-11-06 | 2016-02-24 | 吉林大学 | Magnetotelluric meshless numerical simulation method for random conductive medium model |
| CN107957594A (en) * | 2017-11-15 | 2018-04-24 | 中国石油集团东方地球物理勘探有限责任公司 | The oval bearing calibration of seismic data, dynamic bearing calibration and normal-moveout spectrum computational methods |
| CN109143362A (en) * | 2017-06-28 | 2019-01-04 | 中国石油化工股份有限公司 | Scattered wave separation method based on total scattering angle gathers |
| CN109490965A (en) * | 2018-10-15 | 2019-03-19 | 长江大学 | A kind of heteropical method and device in quantitative assessment stratum |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2869421B1 (en) * | 2004-04-27 | 2006-06-02 | Inst Francais Du Petrole | METHOD FOR RECONSTRUCTING A STOCHASTIC MODEL, REPRESENTATIVE OF A POROUS HETEROGENEOUS MEDIUM, FOR IMPROVING ITS SETTING BY PRODUCTION DATA |
| US7619220B2 (en) * | 2005-11-30 | 2009-11-17 | Jeol Ltd. | Method of measuring aberrations and correcting aberrations using Ronchigram and electron microscope |
| CN111615706A (en) * | 2017-11-17 | 2020-09-01 | 脸谱公司 | Analysis of spatial sparse data based on sub-manifold sparse convolutional neural network |
-
2019
- 2019-08-28 CN CN201910804363.3A patent/CN112444865B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004090575A1 (en) * | 2003-04-09 | 2004-10-21 | Norsar | Method for simulating local prestack depth migrated seismic images |
| EA200702257A2 (en) * | 2007-11-15 | 2008-02-28 | Блаас Холдингс Лимитед | METHOD OF RADAR SENSING OF GROUND NEDR AND A DEVICE FOR ITS IMPLEMENTATION - A COMPLEX OF GEORADICAL LOCATION EXPLORATION |
| CN102253415A (en) * | 2011-04-19 | 2011-11-23 | 中国石油大学(华东) | Method for establishing earthquake response mode based on fracture equivalent medium model |
| WO2014003597A1 (en) * | 2012-06-26 | 2014-01-03 | Schlumberger, Holding Limited | A method for determining pore volume characteristics and porous materials' matrix thermal conductivity |
| CN103792573A (en) * | 2012-10-26 | 2014-05-14 | 中国石油化工股份有限公司 | Seismic wave impedance inversion method based on frequency spectrum fusion |
| CN104077738A (en) * | 2013-12-30 | 2014-10-01 | 辽宁师范大学 | Color image watermarking method based on local histogram characteristics |
| CN104977606A (en) * | 2014-04-02 | 2015-10-14 | 中国石油化工股份有限公司 | Method for establishing fracture-vuggy reservoir seismic numerical model |
| CN104992029A (en) * | 2015-07-20 | 2015-10-21 | 中国科学院国家天文台 | Modeling method for multi-scale non-uniform discrete random medium in lunar soil layer |
| CN105354421A (en) * | 2015-11-06 | 2016-02-24 | 吉林大学 | Magnetotelluric meshless numerical simulation method for random conductive medium model |
| CN109143362A (en) * | 2017-06-28 | 2019-01-04 | 中国石油化工股份有限公司 | Scattered wave separation method based on total scattering angle gathers |
| CN107957594A (en) * | 2017-11-15 | 2018-04-24 | 中国石油集团东方地球物理勘探有限责任公司 | The oval bearing calibration of seismic data, dynamic bearing calibration and normal-moveout spectrum computational methods |
| CN109490965A (en) * | 2018-10-15 | 2019-03-19 | 长江大学 | A kind of heteropical method and device in quantitative assessment stratum |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112444865A (en) | 2021-03-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2879773C (en) | Multi-level reservoir history matching | |
| EP3555759A1 (en) | Systems and methods for generating, deploying, discovering, and managing machine learning model packages | |
| KR20190017454A (en) | Device and method for generating location estimation model | |
| MX2015004001A (en) | Propagating fracture plane updates. | |
| JP2014535044A5 (en) | ||
| Amornwongpaibun et al. | Scattering of anti-plane (SH) waves by a shallow semi-elliptical hill with a concentric elliptical tunnel | |
| CN111967169B (en) | Two-degree body weight abnormal product decomposition numerical simulation method and device | |
| CN112444865B (en) | Method, system, device and storage medium for exploring fracture-cave reservoir | |
| EP2975438B1 (en) | Multiscale method for reservoir models | |
| US20140156246A1 (en) | System for automated identification of surfaces for building of geologic hydrodynamic model of oil and gas deposit by seismic data | |
| CN113608262B (en) | Seismic data processing method and device for calculating rotation component by using translation component | |
| Liang et al. | Uncertainty quantification of geologic model parameters in 3D gravity inversion by Hessian-informed Markov chain Monte Carlo | |
| KR102243917B1 (en) | Device and method for generating geomagnetic sensor based location estimation model using artificial neural networks | |
| Gao et al. | A new code for calculating post-seismic displacements as well as Geoid and gravity changes on a layered visco-elastic spherical earth | |
| CN110083851B (en) | Method and device for determining bottom hole pressure of gas well and storage medium | |
| CN113221409B (en) | Two-dimensional numerical simulation method and device for acoustic waves with coupled finite elements and boundary elements | |
| US20150205003A1 (en) | Method of determining channelway trajectories | |
| US11681838B2 (en) | Distributed Sequential Gaussian Simulation | |
| CN114063161A (en) | Seismic inversion method, device and system | |
| US8494778B2 (en) | Variable grid for finite difference computation | |
| CN114706126B (en) | Quantitative carving method and system for seismic geologic body under geological awareness constraint | |
| WO2024125606A1 (en) | Seismic construction and interpretation method and apparatus based on multi-source multi-task deep learning | |
| CN114359491B (en) | Karst three-dimensional modeling method and device, electronic equipment and storage medium | |
| Sheng et al. | Analysis of 3D formation pressure based on logging data | |
| CN118826933A (en) | Wireless channel response prediction method, device, equipment and storage medium |
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 |