CN112444865A - Method, system and device for exploring slot hole reservoir and storage medium - Google Patents
Method, system and device for exploring slot hole reservoir and storage medium Download PDFInfo
- Publication number
- CN112444865A CN112444865A CN201910804363.3A CN201910804363A CN112444865A CN 112444865 A CN112444865 A CN 112444865A CN 201910804363 A CN201910804363 A CN 201910804363A CN 112444865 A CN112444865 A CN 112444865A
- Authority
- CN
- China
- Prior art keywords
- random
- function
- dimensional
- medium
- autocorrelation
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000003860 storage Methods 0.000 title claims abstract description 19
- 238000005311 autocorrelation function Methods 0.000 claims abstract description 42
- 238000001228 spectrum Methods 0.000 claims abstract description 40
- 230000006870 function Effects 0.000 claims description 77
- 230000008569 process Effects 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 10
- 238000010276 construction Methods 0.000 claims description 10
- 238000011160 research Methods 0.000 abstract description 9
- 238000004891 communication Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 230000004044 response Effects 0.000 description 7
- 230000005236 sound signal Effects 0.000 description 4
- 238000005553 drilling Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 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
- 238000012545 processing Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000004364 calculation method Methods 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
- 238000011161 development Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000010365 information processing Effects 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
- 230000003094 perturbing effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
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. analysis, for interpretation, for correction
- 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. analysis, for interpretation, for correction
- G01V1/282—Application of seismic models, synthetic seismograms
Abstract
The disclosure relates to the technical field of seismic exploration, in particular to a method and a device for exploring a slot-hole reservoir and a storage medium, which are used for solving the technical problem that three-dimensional random medium modeling with directional characteristics is lacked in seismic exploration research in the related art. The method for exploring a fracture-cavity reservoir comprises the following steps: constructing an exponential type elliptical autocorrelation function under a three-dimensional coordinate system; carrying out Fourier transform 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 random medium model according to the normalized random disturbance function and the background medium parameters of the random medium in the slot storage layer; inputting background medium parameters of random media in a slot and hole reservoir to be explored into the constructed three-dimensional random medium model to obtain a random medium model corresponding to the slot and hole reservoir to be explored; and (5) exploring the fracture-cavity reservoir stratum by using a random medium model.
Description
Technical Field
The present disclosure relates to the field of seismic exploration, and in particular, to a method, an apparatus, and a storage medium for exploring a fracture-cavity reservoir.
Background
The Reservoir (or oil Reservoir) Geophysics (Reservoir Geophysics) is brought forward from the last 80 th century, and is widely concerned by geophysicists at home and abroad, so that the Reservoir (or oil Reservoir) Geophysics becomes a hot research field. The geophysical foundation of the oil reservoir is a heterogeneous 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) and provide services for oil and gas exploration, particularly oil and gas development. Because the geophysical relationship of oil reservoirs and complex and fine heterogeneous media is large, complete description is difficult to realize by using a conventional method, a forward model which is flexible, convenient and capable of completely describing the heterogeneity of the oil reservoirs needs to be provided, and a random medium model theory is generated.
Disclosure of Invention
The disclosure provides a method and a device for exploring a fracture-cavity reservoir and a storage medium, which are used for solving the technical problem that seismic exploration research in the related art lacks a three-dimensional random medium modeling with directional characteristics.
To achieve the above object, in a first aspect of the embodiments of the present disclosure, there is provided a method of exploring a fracture-cavity reservoir, the method comprising:
constructing an exponential type elliptical autocorrelation function under a three-dimensional coordinate system;
carrying out Fourier transform 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 the background medium parameters of the random medium in the slot and cave storage layer;
inputting background medium parameters of random media in a slot and hole reservoir to be explored into the constructed three-dimensional random medium model to obtain a random medium model corresponding to the slot and hole reservoir to be explored;
and exploring the fracture-cavity reservoir stratum by using the random medium model.
Optionally, the constructing an exponential elliptical autocorrelation function in a 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 space offset between points in the three-dimensional local coordinate system and the three-dimensional global coordinate system;
and constructing an exponential elliptical autocorrelation function according to the three local autocorrelation lengths and the three rotation angles.
Alternatively, a (x, y, z), b (x, y, z), c (x, y, z) are respectively set as three local autocorrelation lengths of (x, y, z) points in the three-dimensional global coordinate system xyz; according to the three local autocorrelation lengths and the three rotation angles, the expression of the constructed exponential elliptical autocorrelation function is as follows:
in formula (1):
a (x, y, z), b (x, y, z) and c (x, y, z) are respectively a local autocorrelation function and a local autocorrelation length of an OXYZ point (x, y, z) of a three-dimensional global coordinate system; x'1,y′1,z′1For a three-dimensional local coordinate system O' X1'Y1'Z1'spatial offset from point (X, y, z), three-dimensional local coordinate system O' X1'Y1'Z1'is formed by a three-dimensional local coordinate system O' X1Y1 Z1The rotation angle α, β, γ in the counterclockwise direction is obtained, and the expression of the spatial offset is:
in formula (2): d is the spatial offset, α, β, γ are the rotation angles corresponding to a (x, y, z), b (x, y, z), c (x, y, z).
Optionally, the formula for performing fourier transform on the constructed exponential elliptic autocorrelation function is as follows:
in formula (3): Γ denotes the spatial Fourier transform operator, φ (k)x,ky,kz) And a power spectrum function representing a random disturbance function, wherein k is a space wave number, kx is the space wave number in the x coordinate axis direction, ky is the space wave number in the y coordinate axis direction, and kz is the space wave number in the z coordinate axis direction.
Optionally, the obtaining a random perturbation function according to the power spectrum function includes:
generating a corresponding random spectrum sequence by a power spectrum function according to a random process method;
and 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 comprises:
and normalizing the random disturbance function according to the requirements of the random disturbance mean value and the variance.
Optionally, according to the normalized random perturbation function and the background medium parameter of the random medium, the expression of the constructed three-dimensional random medium model is as follows:
in formula (3): rho is density, lambda and mu are Lame parameters, rho0、λ0、μ0Background medium parameters for random media; δ ρ, δ λ, δ μ are the non-uniform disturbance amounts.
In a second aspect of embodiments of the present disclosure, there is provided a system for exploring a fracture-cavity reservoir, the system comprising:
a first construction module configured to construct an exponential elliptical autocorrelation function in a three-dimensional coordinate system;
a Fourier transform module configured to obtain a power spectrum function of a random disturbance function by performing Fourier transform on the constructed exponential elliptic autocorrelation function;
a first obtaining module configured to obtain and normalize a random perturbation function according to the power spectrum 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 slot and hole storage layer;
and the second obtaining module is configured to input background medium parameters of random media in the to-be-explored slot and hole reservoir in the constructed three-dimensional random medium model so as to obtain a random medium model corresponding to the to-be-explored slot and hole reservoir, and therefore exploration of the slot and hole reservoir is carried out by using the random medium model.
In a third aspect of the embodiments of the present disclosure, a computer-readable storage medium is provided, on which a computer program is stored, which when executed by a processor implements the steps of the method of any one of the above first aspects.
In a fourth aspect of the disclosed embodiments, there is provided an apparatus for exploring a fracture-cavity 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 first aspects above.
By adopting the technical scheme, the following technical effects can be at least achieved:
the method is based on a random medium modeling theory, three-dimensional random medium model construction with directional characteristics is carried out, different angles of complex stratums in the stratums are represented by introducing autocorrelation angles, three-dimensional modeling is carried out on the fracture-cavity reservoir, the constructed model is enabled to be more consistent with the actual fracture-cavity reservoir geological conditions, model data are provided for researching fracture-cavity seismic response modes and mechanisms, calibration and identification are better carried out on the seismic response characteristics of the 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 production increase.
Additional features and advantages of the disclosure will be set forth in the detailed description which follows.
Drawings
The accompanying drawings, which are included to provide a further understanding of the disclosure and are incorporated in and constitute a part of this specification, illustrate embodiments of the disclosure and together with the description serve to explain the disclosure without limiting the disclosure. In the drawings:
FIG. 1 is a flow chart of a method of exploring a fracture-cavity reservoir, shown in an exemplary embodiment of the present disclosure.
FIG. 2 is a flow chart illustrating steps involved in a method of exploring a fracture-cavity reservoir in which an exponential elliptical autocorrelation function is constructed according to an exemplary embodiment of the present disclosure.
Fig. 3 is a schematic diagram illustrating a relationship between a three-dimensional local coordinate system and a three-dimensional global coordinate system according to an exemplary embodiment of the disclosure.
Fig. 4 is a schematic diagram illustrating a relationship between two three-dimensional local coordinate systems according to an exemplary embodiment of the present disclosure.
FIG. 5 is a flow chart illustrating steps involved in a method of exploring a fracture-cavity reservoir in which a stochastic perturbation function is obtained according to an exemplary embodiment of the disclosure.
FIG. 6 is a schematic diagram illustrating a three-dimensional argument according to an exemplary embodiment of the present disclosure.
FIG. 7 is a schematic diagram of a three-dimensional stochastic media model shown in an exemplary embodiment of the present disclosure.
FIG. 8 is a block diagram of a system for exploring a fracture-cavity reservoir, shown in an exemplary embodiment of the present disclosure.
Fig. 9 is a block diagram of an apparatus for exploring a fracture-cavity reservoir, shown in an exemplary embodiment of the present disclosure.
Detailed Description
Embodiments of the present disclosure will be described in detail with reference to the accompanying drawings and examples, so that how to apply technical means to solve technical problems and achieve the corresponding technical effects can be fully understood and implemented. The embodiments and various features in the embodiments of the present application can be combined with each other without conflict, and the formed technical solutions are all within the protection scope of the present disclosure.
The inventor of the present disclosure finds, through research, that in the related art, random medium modeling is mainly based on a two-dimensional situation, a fracture-cavity reservoir and seismic response characteristics should be reasonably represented under a three-dimensional condition, and research on a three-dimensional random medium modeling technology is less, and research on three-dimensional random medium modeling involving different dip angle factors of the fracture-cavity reservoir is not found, so that the present disclosure provides a method for exploring the fracture-cavity reservoir, and lays a foundation for research on seismic response of an actual reservoir.
Example one
Fig. 1 is a flow chart illustrating a method for exploring a fracture-cavity reservoir in order to solve the technical problem of lack of three-dimensional stochastic medium modeling with directional characteristics in seismic exploration research in the related art according to an exemplary embodiment of the present disclosure. As shown in fig. 1, the method for exploring a fracture-cavity reservoir shown in the present embodiment may include the steps of:
and S11, constructing an exponential elliptical autocorrelation function under the three-dimensional coordinate system.
S12, carrying out Fourier transform on the constructed exponential elliptical autocorrelation function to obtain a power spectrum function of the random disturbance function.
And S13, obtaining a random disturbance function according to the power spectrum function and normalizing the random disturbance function.
And S14, constructing a three-dimensional random medium model according to the normalized random disturbance function and the background medium parameters of the random medium in the slot storage layer.
S15, inputting background medium parameters of random media in the to-be-explored slot and hole reservoir layer in the constructed three-dimensional random medium model to obtain a random medium model corresponding to the to-be-explored slot and hole reservoir layer;
and S16, exploring the fracture-cavity reservoir by using the random medium model.
As shown in fig. 2, the step S11 of constructing the exponential elliptical autocorrelation function under the three-dimensional coordinate system may include the following steps:
and S111, selecting three local autocorrelation lengths of points in the three-dimensional global coordinate system.
And S112, selecting three rotation angles corresponding to the three local autocorrelation lengths according to the space offset between points in the three-dimensional local coordinate system and the three-dimensional global coordinate system.
S113, constructing an exponential type elliptical autocorrelation function according to the three local autocorrelation lengths and the three rotation angles.
In step S111, three local autocorrelation lengths of the midpoint (x, y, z) in the three-dimensional global coordinate system xyz are assumed to be a (x, y, z), b (x, y, z), and c (x, y, z). x'1,y′1,z′1For a three-dimensional local coordinate system O' X1'Y1'Z1'spatial offset from the (X, y, z) point, three-dimensional local coordinate system O' X1'Y1'Z1The relationship between' and the three-dimensional global coordinate system xyz 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. Further, the coordinate transformation matrix may be constructed by a coordinate rotation method, as shown in equation (2).
Wherein, as shown in FIG. 4, a three-dimensional local coordinate system O' X1'Y1'Z1'is formed by a three-dimensional local coordinate system O' X1Y1 Z1The rotation angles alpha, beta and gamma are obtained 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 the rotation angles α (x, y, z), β (x, y, z), γ (x, y, z), the expression of the constructed exponential elliptical autocorrelation function may be:
in formula (1):
a (x, y, z), b (x, y, z) and c (x, y, z) are the local autocorrelation function and the local autocorrelation length of the xyz point (x, y, z) of the three-dimensional global coordinate system.
After the exponential elliptical autocorrelation function is constructed, step S12 may be executed to obtain a power spectrum function of the random perturbation function by performing fourier transform on the constructed exponential elliptical autocorrelation function.
The constructed exponential elliptical autocorrelation function
The formula for performing the fourier transform may be:
in formula (3): Γ denotes the spatial Fourier transform operator, φ (k)x,ky,kz) And a power spectrum function representing a random disturbance function, wherein k is a space wave number, kx is the space wave number in the x coordinate axis direction, ky is the space wave number in the y coordinate axis direction, and kz is the space wave number in the z coordinate axis direction. Power spectrum function phi (k)x,ky,kz) I.e. the power spectrum of the random disturbance sigma (x, y, z).
After obtaining the power spectrum function, step S13 may be executed to obtain and normalize the random perturbation function according to the power spectrum function. As shown in fig. 5, the obtaining a random perturbation function according to the power spectrum function may include the following steps:
s131, from power spectrum function phi (k)x,ky,kz) And generating a corresponding random spectrum sequence according to a random process method.
S132, applying a spectrum expansion formula of a random process,obtaining the constructed exponential elliptical autocorrelation function
The random perturbation function σ (x, y, z) is described. Wherein, the random spectrum sequence and the random disturbance function sigma (x, y, z) are both 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 value and variance. The mean value of the normalized random perturbation function value is equal to 0, and the variance is a certain value (generally, a value is obtained by percentage of background value perturbation, for example, 0.01, then the set perturbation value satisfies that the mean value is zero, and the variance is 0.01), and the normalization process is a well-known technology in the aspect of mathematics, and is not specifically developed here for the simplicity of the specification.
After the obtained random disturbance function sigma (x, y, z) is normalized, step S14 is executed, and a random medium model is constructed according to the normalized random disturbance function and the background medium parameters of the random medium in the slot and hole storage layer.
In the related art, a random medium can be understood as a theoretical model distributed with small-scale non-uniform disturbance under 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 is0Represents a large-scale homogeneous medium; δ m (X) denotes a small-scale inhomogeneous perturbing medium; σ (X) is the random perturbation characteristic of the medium; m (X) is a random medium; x is (X, y, z) a spatial position vector.
The isotropic elastic medium may be determined by the density ρ and the Lame parameters λ, μ, and when m (X) contains ρ, λ, μ, equation (6) may be expressed as:
in formula (7): rho0、λ0、μ0Background medium parameters for random media, assumed to be constantNumber or slowly varying with spatial coordinates (x, y, z); δ ρ, δ λ, δ μ are the non-uniform perturbation quantities added to the above background, and assume a spatially stationary random process with a mean value of zero, a variance of a constant value, and some autocorrelation function.
Since the directional characteristics, such as different dip angle factors of the fracture-cavity reservoir, are not considered in the formula (6) and the formula (7), the stochastic medium model for the complex formation is not accurate. The method adds directional characteristics in the modeling process, and can obtain the random disturbance function sigma (x, y, z) to be constructed after normalizing the obtained random disturbance function sigma (x, y, z)
The expression of the constructed three-dimensional random medium model, which is taken as an autocorrelation function and has a specified mean value and variance, is as follows:
in formula (4): rho is density, lambda and mu are Lame parameters, rho0、λ0、μ0Background medium parameters of random media in the seam hole storage layer; δ ρ, δ λ, δ μ are the non-uniform disturbance amounts.
After the three-dimensional random medium model is built, step S15 is executed, background medium parameters of random media in the slot-hole reservoir to be explored are input into the built three-dimensional random medium model to obtain a random medium model corresponding to the slot-hole reservoir to be explored, and then the random medium model is used for exploring the slot-hole reservoir. Because the actual medium condition contains the inhomogeneous body, the speed of the underground medium is expressed as certain randomness, and the inhomogeneous body has various spatial distributions in the underground medium, the method increases the directionality of the random medium space, can reasonably construct random medium models under different conditions, performs wave field simulation through models of different random medium occurrence states, analyzes the profile wave field characteristics, guides the outdoor actual seismic exploration by utilizing indoor analysis results and experiences, calibrates and identifies the seismic response characteristics of the outdoor actual deep layer micro-amplitude fracture-cavity reservoir, improves the accuracy of outdoor drilling, and provides data support for oil and gas production increase.
The whole process of building the three-dimensional random model with the directional characteristic is shown below by an example. The following are the three-dimensional stochastic model parameters in the example:
size of the model: nx is 100, Ny is 100, Nz is 100;
grid spacing: dx is 1m, dy is 1m, dz is 1 m;
speed average value: v. of0＝5500m/s2；
Speed disturbance standard deviation: 10 percent;
autocorrelation length: a is 3m, b is 10m, and c is 20 m; the autocorrelation angle is: α is 0 °, β is 0 °, γ is 75 °;
the exponential elliptic autocorrelation function is carried out by adopting a formula (1), firstly, according to a given model range, the argument required by the Fourier transform of the exponential elliptic function is completed, the argument range is between 0 and 2 pi and is changed in a random characteristic manner, an intercepted two-dimensional section is shown in figure 6, then the realization of a disturbance power spectrum is completed by combining a transformation function, the calculation of a disturbance quantity is further completed according to a disturbance mean value and a variance value, and a final three-dimensional random medium model containing the autocorrelation angle 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 using fig. 7, the background medium parameters of the 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 is obtained.
The method is based on a random medium modeling theory, three-dimensional random medium model construction with directional characteristics is carried out, different angles of a fracture-cavity reservoir stratum in the stratum are represented by introducing autocorrelation angles, three-dimensional modeling is carried out on the fracture-cavity reservoir stratum, the constructed model is enabled to better accord with actual geological conditions, model data is provided for researching a fracture-cavity seismic response mode and mechanism, calibration and identification are better carried out on seismic response characteristics of an actual deep micro-amplitude fracture-cavity reservoir stratum in the field, the accuracy of field drilling is improved, and support is provided for oil and gas production increase.
It should be noted that the method embodiment shown in fig. 1 is described as a series of acts or combinations for simplicity of description, but it should be understood by those skilled in the art that the present disclosure is not limited by the order of acts or steps described. Further, those skilled in the art will appreciate that the embodiments described in the specification are presently preferred and that no particular act is required in order to implement the disclosure.
Example two
Fig. 8 is a system for exploring a fracture and cave reservoir, according to an exemplary embodiment of the present disclosure, and as shown in fig. 8, the system 300 for exploring a fracture and cave reservoir includes:
a first construction module 310 configured to construct an exponential elliptical 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 elliptic autocorrelation function;
a first obtaining module 330 configured to obtain and normalize a random perturbation function according to the power spectrum function;
a second constructing module 340 configured to construct a random medium model according to the normalized random perturbation function and the background medium parameters of the random medium in the slot and cave storage layer;
and a second obtaining module 350, configured to input the background medium parameters of the random medium in the fracture-cavity reservoir to be explored into the three-dimensional random medium model 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 space offset between points in the three-dimensional local coordinate system and the three-dimensional global coordinate system; and constructing an exponential elliptical autocorrelation function according to the three local autocorrelation lengths and the three rotation angles.
Alternatively, a (x, y, z), b (x, y, z), c (x, y, z) are respectively set as three local autocorrelation lengths of (x, y, z) points in the three-dimensional global coordinate system xyz; based on the three local autocorrelation lengths and the three rotation angles, the expression of the exponential elliptical autocorrelation function constructed by the first construction module 310 is:
in formula (1):
a (x, y, z), b (x, y, z) and c (x, y, z) are respectively a local autocorrelation function and a local autocorrelation length of an OXYZ point (x, y, z) of a three-dimensional global coordinate system; x'1,y′1,z′1As a three-dimensional local coordinate system O1X′1Y′1Z′1Spatial offset of the medium relative to the point (x, y, z), three-dimensional local coordinate system O1X′1Y′1Z′1Is composed of a three-dimensional local coordinate system O1X1Y1Z1The rotation angles α, β, γ in the counterclockwise direction along the respective coordinates are obtained, and the expression of the spatial offset is:
in formula (2): d is the spatial offset, α, β, γ are the rotation angles corresponding to a (x, y, z), b (x, y, z), c (x, y, z).
Optionally, the fourier transform module 320 performs fourier transform according to the following formula:
in formula (3): Γ denotes the spatial Fourier transform operator, φ (k)x,ky,kz) A power spectrum function representing a random perturbation function, k being the number of spatial waves,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.
Optionally, the obtaining module 340 is further configured to: generating a corresponding random spectrum sequence by a power spectrum function according to a random process method; and 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 value and the variance.
Optionally, the expression of the stochastic medium model constructed by the second construction module 350 is:
in formula (3): rho0、λ0、μ0Background medium parameters for random media; δ ρ, δ λ, δ μ are the non-uniform disturbance amounts.
With regard to the system in the above embodiment, the specific manner in which each module performs the operation has been described in detail in the embodiment related to the method, and will not be elaborated here.
EXAMPLE III
The present disclosure also provides a computer-readable storage medium having stored thereon a computer program which, when being executed by a processor, carries out the method steps of exploring a fracture-cavity reservoir as set forth in any of the above-mentioned alternative embodiments.
The method implemented when the computer program of the method for exploring a fracture-cavity reservoir run on the processor is executed may refer to the specific embodiment of the method for exploring a fracture-cavity reservoir of the present disclosure, and details are not repeated here.
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 (CPU), a Network Processor (NP), and the like.
Example four
The present disclosure also provides a device for exploring a fracture-cavity 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-cavity reservoir of any of the alternative embodiments described above.
FIG. 9 is a block diagram illustrating an apparatus 400 for exploring a fracture-cavity reservoir, according to an exemplary 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.
The processor 401 is configured to control the overall operation of the apparatus 400, so as to complete all or part of the steps in the above-mentioned method for constructing a three-dimensional random medium model. The memory 402 is used to store various types of data to support operation of the apparatus 400, and such data may include, for example, instructions for any application or method operating on the apparatus 400, as well as application-related data. The Memory 402 may be implemented by any type of volatile or non-volatile Memory device or combination thereof, such as Static Random Access Memory (SRAM), Electrically Erasable Programmable Read-Only Memory (EEPROM), Erasable Programmable Read-Only Memory (EPROM), Programmable Read-Only Memory (PROM), Read-Only Memory (ROM), magnetic Memory, flash Memory, magnetic disk or optical disk. The multimedia components 403 may include a screen and an audio component. Wherein the screen may be, for example, a touch screen and the audio component is used 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 further be stored in the memory 402 or transmitted through the communication component 405. The audio assembly also includes at least one speaker for outputting audio signals. The I/O interface 404 provides an interface between the processor 401 and other interface modules, such as 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), 2G, 3G, or 4G, or a combination of one or more of them, so that the corresponding Communication component 405 may include: 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 (DSPs), Digital Signal Processing Devices (DSPDs), Programmable Logic Devices (PLDs), Field Programmable Gate Arrays (FPGAs), controllers, microcontrollers, microprocessors, or other electronic components for performing the above-described method for exploring a fracture-cavity reservoir.
In another exemplary embodiment, a computer readable storage medium, such as a memory 402, comprising program instructions executable by a processor 401 of the apparatus 400 to perform the above-described method of exploring a fracture-cavity reservoir is also provided.
The preferred embodiments of the present disclosure are described in detail with reference to the accompanying drawings, however, the present disclosure is not limited to the specific details of the above embodiments, and various simple modifications may be made to the technical solution of the present disclosure within the technical idea of the present disclosure, and these simple modifications all belong to the protection scope of the present disclosure.
It should be noted that the various features described in the above embodiments may be combined in any suitable manner without departing from the scope of the invention. In order to avoid unnecessary repetition, various possible combinations will not be separately described in this disclosure.
In addition, any combination of various embodiments of the present disclosure may be made, and the same should be considered as the disclosure of the present disclosure, as long as it does not depart from the spirit of the present disclosure.
Claims (10)
1. A method of exploring a fracture-cavity reservoir, the method comprising:
constructing an exponential type elliptical autocorrelation function under a three-dimensional coordinate system;
carrying out Fourier transform 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 the background medium parameters of the random medium in the slot and cave storage layer;
inputting background medium parameters of random media in a slot and hole reservoir to be explored into the constructed three-dimensional random medium model to obtain a random medium model corresponding to the slot and hole reservoir to be explored;
and exploring the fracture-cavity reservoir stratum by using the random medium model.
2. The method of claim 1, wherein constructing the exponential elliptical autocorrelation function in a three-dimensional coordinate system comprises:
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 space offset between points in the three-dimensional local coordinate system and the three-dimensional global coordinate system;
and constructing an exponential elliptical autocorrelation function according to the three local autocorrelation lengths and the three rotation angles.
3. The method according to claim 2, wherein a (x, y, z), b (x, y, z), c (x, y, z) are set as three local autocorrelation lengths of (x, y, z) points in a three-dimensional global coordinate system, OXYZ, respectively; according to the three local autocorrelation lengths and the three rotation angles, the expression of the constructed exponential elliptical autocorrelation function is as follows:
in formula (1):
a (x, y, z), b (x, y, z) and c (x, y, z) are respectively a local autocorrelation function and a local autocorrelation length of an OXYZ point (x, y, z) of a three-dimensional global coordinate system; x'1,y′1,z′1Is a three-dimensional local coordinate system O 'X'1Y'1Z'1Spatial offset from point (X, y, z), three-dimensional local coordinate system O 'X'1Y'1Z'1Is composed of a three-dimensional local coordinate system O' X1 Y1 Z1The rotation angle α, β, γ in the counterclockwise direction is obtained, and the expression of the spatial offset is:
in formula (2): d is the spatial offset, α, β, γ are the rotation angles corresponding to a (x, y, z), b (x, y, z), c (x, y, z).
4. The method of claim 3, wherein the formula for the constructed exponential elliptical autocorrelation function to perform Fourier transform is:
in formula (3): Γ denotes the spatial Fourier transform operator, φ (k)x,ky,kz) A power spectrum function representing a random disturbance function, k is a space wave number, and kx is an x coordinate axis squareAnd the space wave number in the direction is ky, and kz, the space wave number in the direction of the y coordinate axis and the space wave number in the direction of the z coordinate axis.
5. The method of claim 4, wherein obtaining a stochastic perturbation function from the power spectrum function comprises:
generating a corresponding random spectrum sequence by a power spectrum function according to a random process method;
and applying a spectrum expansion formula of a random process to obtain a random disturbance function described by the constructed exponential elliptic autocorrelation function.
6. The method of claim 5, wherein normalizing the random perturbation function comprises:
and normalizing the random disturbance function according to the requirements of the random disturbance mean value and the variance.
7. The method according to claim 3, wherein the expression of the three-dimensional random medium model constructed according to the normalized random perturbation function and the background medium parameters of the random medium is as follows:
in formula (4): rho is density, lambda and mu are Lame parameters, rho0、λ0、μ0Background medium parameters of random media in the seam hole storage layer; δ ρ, δ λ, δ μ are the non-uniform disturbance amounts.
8. A system for exploring a fracture-cavity reservoir, comprising:
a first construction module configured to construct an exponential elliptical autocorrelation function in a three-dimensional coordinate system;
a Fourier transform module configured to obtain a power spectrum function of a random disturbance function by performing Fourier transform on the constructed exponential elliptic autocorrelation function;
a first obtaining module configured to obtain and normalize a random perturbation function according to the power spectrum 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 slot and hole storage layer;
and the second obtaining module is configured to input background medium parameters of random media in the to-be-explored slot and hole reservoir in the constructed three-dimensional random medium model so as to obtain a random medium model corresponding to the to-be-explored slot and hole reservoir, and therefore exploration of the slot and hole reservoir is carried out by using the random medium model.
9. An apparatus for exploring a fracture-cavity reservoir, comprising:
a memory having a computer program stored thereon; and
a processor for executing the computer program in the memory to carry out the steps of the method of any one of claims 1 to 7.
10. A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, carries out the steps of the method according to any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910804363.3A CN112444865A (en) | 2019-08-28 | 2019-08-28 | Method, system and device for exploring slot hole reservoir and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910804363.3A CN112444865A (en) | 2019-08-28 | 2019-08-28 | Method, system and device for exploring slot hole reservoir and storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112444865A true CN112444865A (en) | 2021-03-05 |
Family
ID=74741452
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910804363.3A Pending CN112444865A (en) | 2019-08-28 | 2019-08-28 | Method, system and device for exploring slot hole reservoir and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112444865A (en) |
Citations (14)
| 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 |
| US20060241920A1 (en) * | 2004-04-27 | 2006-10-26 | Mickaele Le Ravalec-Dupin | Method of reconstructing a stochastic model, representative of a porous heterogeneous medium, to improve its calibration by production data |
| US20070120055A1 (en) * | 2005-11-30 | 2007-05-31 | Jeol Ltd. | Method of measuring aberrations and correcting aberrations using ronchigram and electron microscope |
| 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 |
| 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 |
| US20190156206A1 (en) * | 2017-11-17 | 2019-05-23 | Facebook, Inc. | Analyzing Spatially-Sparse Data Based on Submanifold Sparse Convolutional Neural Networks |
-
2019
- 2019-08-28 CN CN201910804363.3A patent/CN112444865A/en active Pending
Patent Citations (14)
| 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 |
| US20060241920A1 (en) * | 2004-04-27 | 2006-10-26 | Mickaele Le Ravalec-Dupin | Method of reconstructing a stochastic model, representative of a porous heterogeneous medium, to improve its calibration by production data |
| US20070120055A1 (en) * | 2005-11-30 | 2007-05-31 | Jeol Ltd. | Method of measuring aberrations and correcting aberrations using ronchigram and electron microscope |
| 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 |
| 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 |
| US20190156206A1 (en) * | 2017-11-17 | 2019-05-23 | Facebook, Inc. | Analyzing Spatially-Sparse Data Based on Submanifold Sparse Convolutional Neural Networks |
| CN109490965A (en) * | 2018-10-15 | 2019-03-19 | 长江大学 | A kind of heteropical method and device in quantitative assessment stratum |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Higdon et al. | Computer model calibration using high-dimensional output | |
| Belson et al. | Algorithm 945: modred—a parallelized model reduction library | |
| RU2648391C2 (en) | Downhole gradiometric ranging utilising transmitters and receivers having magnetic dipoles | |
| Barker et al. | Using performance modeling to design large-scale systems | |
| CN104854480A (en) | Apparatus and methods to find a position in an underground formation | |
| KR20190017454A (en) | Device and method for generating location estimation model | |
| Liu et al. | Model updating of complex structures using the combination of component mode synthesis and Kriging predictor | |
| Zhu et al. | A MEMS gyroscope noise suppressing method using neural architecture search neural network | |
| Overstall et al. | acebayes: An R package for Bayesian optimal design of experiments via approximate coordinate exchange | |
| Farago et al. | On the convergence and local splitting error of different splitting schemes | |
| Qaddouri et al. | Optimized Schwarz methods with an overset grid for the shallow-water equations: preliminary results | |
| CN110850057B (en) | Reservoir fracture modeling method and system based on self-similarity theory | |
| Fillacier et al. | Calculating Prediction Uncertainty using Posterior Ensembles Generated from Proxy Models | |
| Chen et al. | Forecasting and remote sensing outer belt relativistic electrons from low Earth orbit | |
| CN112444865A (en) | Method, system and device for exploring slot hole reservoir and storage medium | |
| KR102243917B1 (en) | Device and method for generating geomagnetic sensor based location estimation model using artificial neural networks | |
| Sun et al. | Computational benefit of GPU optimization for the atmospheric chemistry modeling | |
| US11139569B2 (en) | Antenna system, method and computer program product, with real time axial ratio polarization correction | |
| Arroyo et al. | Algorithm 1013: An R implementation of a continuous spectral algorithm for simulating vector Gaussian random fields in Euclidean spaces | |
| 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 | |
| Tian et al. | Optimizing wifi ap placement for both localization and coverage | |
| CN114449439B (en) | Underground pipe gallery space positioning method and device | |
| CN112764105B (en) | HTI medium quasi-longitudinal wave forward simulation method and device, storage medium and processor | |
| US11681838B2 (en) | Distributed Sequential Gaussian Simulation |
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 |