CN110286410A - Crack inversion method and device based on diffraction wave energy - Google Patents
Crack inversion method and device based on diffraction wave energy Download PDFInfo
- Publication number
- CN110286410A CN110286410A CN201910720254.3A CN201910720254A CN110286410A CN 110286410 A CN110286410 A CN 110286410A CN 201910720254 A CN201910720254 A CN 201910720254A CN 110286410 A CN110286410 A CN 110286410A
- Authority
- CN
- China
- Prior art keywords
- wave energy
- crack
- diffraction wave
- initial model
- data
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 48
- 208000010392 Bone Fractures Diseases 0.000 claims description 21
- 206010017076 Fracture Diseases 0.000 claims description 21
- 230000015572 biosynthetic process Effects 0.000 claims description 17
- 238000000926 separation method Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000000605 extraction Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 3
- 230000006870 function Effects 0.000 description 29
- 238000012512 characterization method Methods 0.000 description 4
- 238000004590 computer program Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012897 Levenberg–Marquardt algorithm Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (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 present invention provides a kind of crack inversion methods and device based on diffraction wave energy, are related to the technical field of FRACTURE PREDICTION, including the isolated diffraction wave energy from seismic data, and diffraction wave energy includes prestack diffraction wave energy and poststack diffraction wave energy；The initial model of crack attribute is determined according to poststack diffraction wave energy；Using prestack diffraction wave energy and initial model as input data, and inverting is carried out as follows, until the value of inversion objective function meets expection, using initial model as the crack attribute after inverting: by input data input reverse objective function, obtaining the value of inversion objective function；If the value of inversion objective function does not meet expection, initial model data are adjusted, and using prestack diffraction wave energy and initial model adjusted as new input data, inverting is carried out by the diffraction wave energy in crack, crack property parameters are obtained, and then realize the purpose that fracture oily situation is predicted.
Description
Technical field
The present invention relates to FRACTURE PREDICTION technical fields, more particularly, to a kind of crack inversion method based on diffraction wave energy
And device.
Background technique
Fracture-type reservoir is the important goal of oil-gas exploration, the property in crack be exploration plan and control drilling well it is important because
Element.
Currently, diffracted wave is the effective means for portraying crack, but compared to back wave, the energy of diffracted wave is very weak,
So it is poor using prediction result of the common seismic method for crack attribute, and then crack can not be distinguished by crack attribute
The feature of contained fluid, it is more inconvenient for oil-gas exploration and production aspect.
Summary of the invention
The purpose of the present invention is to provide crack inversion methods and device based on diffraction wave energy, pass through the diffraction in crack
Wave energy carries out inverting, obtains crack property parameters, and then realizes the purpose that fracture oily situation is predicted.
In a first aspect, the embodiment of the present invention provides a kind of crack inversion method based on diffraction wave energy, comprising:
The isolated diffraction wave energy from seismic data, the diffraction wave energy include prestack diffraction wave energy and poststack
Diffraction wave energy；
The initial model of crack attribute is determined according to the poststack diffraction wave energy；
Using the initial model and the prestack diffraction wave energy as input data, and carry out as follows anti-
It drills, until the value of inversion objective function meets expection, using the initial model as the crack attribute after inverting:
By the input data input reverse objective function, the value of inversion objective function is obtained；
If the value of the inversion objective function does not meet expection, the initial model is adjusted, and will be described folded
Preceding diffraction wave energy and initial model adjusted are as new input data.
In alternative embodiments, the initial model of crack attribute is determined according to the poststack diffraction wave energy, comprising:
Equivalent fissure porosity model is constructed according to the poststack diffraction wave energy；
The initial model for characterizing crack attribute is calculated according to the equivalent fissure porosity model.
In alternative embodiments, equivalent fissure porosity model is constructed according to the poststack diffraction wave energy, comprising:
The obliquity information and azimuth information in crack are obtained according to Image Logging Data；
Intensity information is generated according to the obliquity information and azimuth information；
According to the poststack diffracted wave Energy extraction fracture development information, it is based on the fracture development information and the crack
Strength information constructs equivalent fissure porosity model in the way of Decided modelling and stochastic modeling respectively.
In alternative embodiments, the inversion objective function is determined according to angle gathers data and forward modeling data, wherein
It carries out extracting the angle gathers data from the prestack diffraction wave energy, be obtained according to formation of elastic parameters and the initial model
To the forward modeling data.
In alternative embodiments, the angle gathers data are extracted from the prestack diffraction wave energy, comprising:
The common midpoint gather of adjacent predetermined number in the prestack diffraction wave energy is overlapped, channel set is formed
Data；
Angle gathers data are extracted from the channel set data according to predetermined angle.
In alternative embodiments, the formation of elastic parameters further include seismic wave spread speed and stratum in the earth formation
Media density obtains forward modeling data according to formation of elastic parameters and the initial model, comprising:
Reflection coefficient is calculated by the initial model, the spread speed and the stratum media density；
Convolution calculating is carried out according to the reflection coefficient and seismic wavelet, obtains forward modeling data.
In alternative embodiments, the method also includes:
It is predicted according to the crack attribute fracture oily situation after the inverting.
It in alternative embodiments, further include similarity constraint item in the inversion objective function, the similitude is about
Beam item is used to calculate the average value of crack attribute after the inverting of adjacent seismic channel.
In alternative embodiments, the method also includes:
The inversion objective function is indicated according to the following formula:
Wherein, SobsFor angle gathers data, S is forward modeling data, and A and B indicate that likeness coefficient matrix, m indicate introductory die
Type, λ2And γ2Regularization parameter is represented, i is sampling number, and j is seismic channel number.
Second aspect, the embodiment of the present invention provide a kind of crack inverting device based on diffraction wave energy, comprising:
Separation module, for the isolated diffraction wave energy from seismic data, the diffraction wave energy include prestack around
Ejected wave energy and poststack diffraction wave energy；
Determining module, for determining the initial model of crack attribute according to the poststack diffraction wave energy；
Inverting module is used for using the initial model and the prestack diffraction wave energy as input data, and by such as
Lower step carries out inverting, until the value of inversion objective function meets expection, belongs to the initial model as the crack after inverting
Property:
By the input data input reverse objective function, the value of inversion objective function is obtained；
If the value of the inversion objective function does not meet expection, the initial model is adjusted, and will be described folded
Preceding diffraction wave energy and initial model adjusted are as new input data.
The embodiment of the invention provides a kind of crack inversion methods and device based on diffraction wave energy, are based on seismic data
Poststack diffraction wave energy obtained after separation determines the initial model of crack attribute, further according to inversion objective function, by introductory die
Type and prestack diffraction wave energy are updated, until inversion objective function is accorded with as input by being adjusted iteration to initial model
When closing expected value, the initial model of the crack attribute after obtaining characterization inverting, and then realize that fracture oily situation carries out
The purpose of prediction.
Detailed description of the invention
It, below will be to specific in order to illustrate more clearly of the specific embodiment of the invention or technical solution in the prior art
Embodiment or attached drawing needed to be used in the description of the prior art be briefly described, it should be apparent that, it is described below
Attached drawing is some embodiments of the present invention, for those of ordinary skill in the art, before not making the creative labor
It puts, is also possible to obtain other drawings based on these drawings.
Fig. 1 is a kind of crack inversion method flow chart based on diffraction wave energy provided in an embodiment of the present invention；
Fig. 2 is a kind of crack attribution inversion result schematic diagram provided in an embodiment of the present invention；
Fig. 3 is a kind of functional block diagram of the crack inverting device based on diffraction wave energy provided in an embodiment of the present invention.
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with attached drawing to the present invention
Technical solution be clearly and completely described, it is clear that described embodiments are some of the embodiments of the present invention, rather than
Whole embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art are not making creative work premise
Under every other embodiment obtained, shall fall within the protection scope of the present invention.
Early in the 1970s, people begin to the prediction for carrying out crack using seismic method, such as shear wave splitting
Theoretical (Crampin, 1978；Alford, 1986), the poststacks attribute such as curvature body and ant body, indication using prestack seismic amplitude with offset distance change
Change (Amplitude variation with offset, AVO), AVAZ, VVAZ isoinversion method, these methods are in certain journey
The qualitative description in crack is realized on degree.
Wherein, AVO theory fracture reservoir is described, and carries out AVO using the angle gathers data of earthquake reflected wave
Inverting, to obtain the elastic parameter and crack property parameters of reservoir.It is by formation impedance interface and crack according to seismic reflection
It is reflecting jointly as a result, having derived approximate AVO formula.But reflected energy caused by crack is far smaller than caused by impedance interface
Reflected energy carries out inverting by common seismic data, and Simultaneous Inversion goes out the parameters such as formation velocity, density and crack attribute, this
When crack response be usually submerged, and then can not Accurate Prediction crack attribute, therefore, when reservoir elastic parameter calculating miss
Difference can make the fracture parameters of prediction that huge variation occur error very little, while can also have the interference of noise, into
And the prediction result stability of crack attribute is poor.
To sum up, since response of the seismic reflection for crack is faint, so the side of crack attribute forecast can accurately be realized
Method is less.And the feature of fluid contained by crack is distinguished by crack attribute, there is weight again for oil-gas exploration and production
Big practical significance, therefore, the problem of crack attribute Accurate Prediction, is urgently to be resolved.
Based on this, a kind of crack inversion method and device based on diffraction wave energy provided in an embodiment of the present invention can be with
Inverting is carried out by the diffraction wave energy in crack, obtains crack property parameters, and then it is pre- to realize that fracture oily situation carries out
The purpose of survey.
To be based on diffraction wave energy to one kind disclosed in the embodiment of the present invention first convenient for understanding the present embodiment
Crack inversion method describe in detail.
Fig. 1 is a kind of crack inversion method flow diagram based on diffraction wave energy provided in an embodiment of the present invention.
Referring to Fig.1, the embodiment of the present invention is based on diffraction wave energy fracture attribute and carries out inverting, comprising the following steps:
Step S102, the isolated diffraction wave energy from seismic data, diffraction wave energy include prestack diffraction wave energy
With poststack diffraction wave energy.
It should be noted that when common seismic data carry out inverting, it also can inverting while attribute in inverting crack
The parameters such as interval velocity, density, since diffraction energy is very weak compared to reflected energy, so the response in crack is usually submerged in
In various inverted parameters, and then can not Accurate Prediction crack attribute.
Here, earthquake data before superposition includes the impedance interface (interfaces of the Parameters variations such as formation velocity or density) on stratum
The stratum reflected energy and crack diffraction wave energy of generation.By in earthquake data before superposition stratum reflected energy and crack around
Ejected wave energy is separated, and prestack diffraction wave energy is obtained.It is real using the methods of decomposition of plane wave, FK filtering during the separation process
Guarantor's width separation of existing diffracted wave field, retains effective diffraction information of the spokes such as fracture, similarly separate by post-stack seismic data
To poststack diffraction wave energy.
Step S104 determines the initial model of crack attribute according to poststack diffraction wave energy.
Wherein, the initial model proposed in the embodiment of the present invention is the set for characterizing the numerical matrix of crack attribute, is belonged to
The usual noun that those skilled in the art understand that, details are not described herein.Crack attribute includes normal direction flexibility SNWith it is tangential soft
Spend ST。
Step S106 using initial model and prestack diffraction wave energy as input data, and is carried out anti-as follows
It drills, until the value of inversion objective function meets expection, using initial model as the crack attribute after inverting:
By input data input reverse objective function, the value of inversion objective function is obtained；
If the value of inversion objective function does not meet expection, initial model is adjusted, and by prestack diffraction wave energy
And initial model adjusted is as new input data.
In the preferred embodiment of practical application, split based on seismic data poststack diffraction wave energy determination obtained after separation
The initial model for stitching attribute, further according to inversion objective function, using initial model and prestack diffraction wave energy as input, by right
Initial model is adjusted iteration update, until crack when inversion objective function meets expected value, after obtaining characterization inverting
The initial model of attribute, and then realize the purpose that fracture oily situation is predicted.
In alternative embodiments, above-mentioned steps S104 is further comprising the steps of:
1, equivalent fissure porosity model is constructed according to poststack diffraction wave energy.
Specifically, the obliquity information and azimuth information in crack are obtained according to Image Logging Data；According to obliquity information and
Azimuth information generates intensity information；According to poststack diffracted wave Energy extraction fracture development information, believed based on fracture development
Breath and intensity information construct equivalent fissure porosity model in the way of Decided modelling and stochastic modeling respectively.
Wherein, equivalent fissure porosity model can be by multiple Data Integrations such as earthquake, well logging, geology and engineering, this hair
The method that bright embodiment uses stochastic modeling, inclination angle and the azimuth information in crack are obtained by Image Logging Data, and are generated
Intensity curve by poststack diffracted wave Energy extraction fracture development information, and passes through earthquake using the method for stochastic modeling
More attribute fusions come restraint crack intensity and fracture development information, complete equivalent fissure porosity model.
2, the initial model for characterizing crack attribute is calculated according to equivalent fissure porosity model.
On the basis of equivalent fissure porosity model, the introductory die for being used for quantitatively characterizing crack attribute is obtained by calculating
Type.Characterization crack normal direction flexibility S is calculated by following formulaNWith tangential flexibility STInitial model.
SN=0.0013*ln (φf)+0.0223
ST=0.0107*ln (φf)+0.1783
Wherein, SNFor normal direction flexibility, STFor tangential flexibility, unit m/Gpa, φfIndicate porosity.
As a kind of optional embodiment, inversion objective function is determined according to angle gathers data and forward modeling data, wherein from
It carries out extracting angle gathers data in prestack diffraction wave energy, forward modeling data is obtained according to formation of elastic parameters and initial model.
Further, for the interference of Attenuating Random Noise, and then improve signal-to-noise ratio, the embodiment of the present invention first from prestack around
Angle gathers data are extracted in ejected wave energy, then construct inversion objective function by angle gathers data, comprising:
The common midpoint gather of adjacent predetermined number in prestack diffraction wave energy is overlapped, channel set number is formed
According to；Angle gathers data are extracted from channel set data according to predetermined angle.
Here, channel set is formed as a result, being overlapped using diffracted wave guarantor's width separation, then extracts and meets from channel set
The angle gathers data that predetermined angle requires.
Further, formation of elastic parameters include that spread speed and stratum media are close in the earth formation for crack attribute, seismic wave
Degree, obtains forward modeling data according to formation of elastic parameters and initial model, comprising:
1, reflection coefficient is calculated by initial model, spread speed and stratum media density.
It should be noted that reflection is the impedance interface (boundary of the Parameters variations such as formation velocity or density due to stratum
Face) generate, the embodiment of the present invention is based only upon diffracted wave data, is equivalently employed without stratum reflection, that is, does not have impedance interface, carries on the back
Scape medium is no resistance, i.e. P- and S-wave velocity variation is 0.Reflection coefficient can be asked using following formula in the calculating of forward modeling data
:
Rp(θ)=i ω STT(θ)+iωSNN(θ)
μ=ρ VS 2
Wherein, θ is incidence angle, and ω is angular frequency, and P is ray parameter, STFor tangential flexibility, SNFor normal direction flexibility, VPFor ground
Longitudinal wave propagation speed, V in seismic wave stratumSFor shear wave propagation speed in seismic wave stratum, ρ is density.
It is understood that μ, qp, T (θ), N (θ) can be described as middle coefficient, can be according to aforementioned with clear meaning
Parameters are indicated.
2, convolution calculating is carried out according to reflection coefficient and seismic wavelet, obtains forward modeling data.
Here, inversion objective function Q can be used for characterizing the error between angle gathers data and forward modeling data, as following formula carries out
It indicates:
Wherein, SobsFor angle gathers data, S is forward modeling data.
In order to improve the inversion accuracy of crack attribute, the embodiment of the present invention carries out inverting, inverting to multiple seismic channels simultaneously
Further include similarity constraint item in objective function, that is, constructs a coefficient matrix, similarity constraint item is for calculating adjacent earthquake
The average value of the inversion result (crack attribute after inverting) in road.
Inversion objective function Q and coefficient matrix are as follows:
Wherein, SobsFor angle gathers data, S is forward modeling data, and A and B indicate that likeness coefficient matrix, m indicate introductory die
Type, λ2And γ2Regularization parameter is represented, i is sampling number, and j is seismic channel number.
The embodiment of the present invention carries out inverting using the method for damped least squares, due to normal direction flexibility SNWith tangential flexibility ST
Value differ about two orders of magnitude, the influence degree both balanced in inverting, using the method for seeking logarithm, according to pre-
It is updated if step-length is adjusted iteration to initial model.
As a kind of possible application scenarios, the embodiment of the invention also includes: according to the crack attribute fracture after inverting
Situation predicted, the crack attribute normal direction flexibility S that method obtains through the embodiment of the present inventionNInversion result as shown in Fig. 2,
It is contemplated that the number of levels of the position in each crack and crack attribute, the crack attribute that will be obtained after inverting in three-dimensional grid
Numerical applications can be predicted current formation characteristic of crack into actual current formation model of place, i.e., the oil storage of stratum where,
Water storage, gas storage.
Crack attribute forecast method provided in an embodiment of the present invention based on diffraction wave energy, first by the ground in seismic data
Layer reflected energy and the separation of crack diffraction energy, obtain prestack diffraction wave energy and poststack diffraction wave energy, then utilize crack
Diffraction wave energy carry out AVO inverting, and then obtain characterization Cracks character property parameters.
Wherein, input data initial model includes the numerical matrix set of normal direction flexibility and tangential flexibility, is one rough
, data that indicate Long-term change trend, iteration by inverting updates, and the normal direction exported for characterizing the crack attribute after inverting is soft
The numerical matrix set of degree and tangential flexibility, data portrayed with details, accurate.
The embodiment of the present invention also provides a kind of crack inverting device based on diffraction wave energy, as shown in Figure 3, comprising:
Separation module, being used for isolated diffraction wave energy, diffraction wave energy from seismic data includes prestack diffracted wave
Energy and poststack diffraction wave energy；
Determining module, for determining the initial model of crack attribute according to poststack diffraction wave energy；
Inverting module, for using initial model and prestack diffraction wave energy as input data, and as follows into
Row inverting, until the value of inversion objective function meets expection, using initial model as the crack attribute after inverting:
By input data input reverse objective function, the value of inversion objective function is obtained；
If the value of inversion objective function does not meet expection, initial model is adjusted, and by prestack diffraction wave energy
And initial model adjusted is as new input data.
Crack inverting device provided in an embodiment of the present invention based on diffraction wave energy, is based on provided by the above embodiment
The crack inversion method technical characteristic having the same of diffraction wave energy reaches phase so also can solve identical technical problem
Same technical effect.
The computer program of crack inversion method and device based on diffraction wave energy provided by the embodiment of the present invention produces
Product, the computer readable storage medium including storing program code, before the instruction that said program code includes can be used for execution
Method described in the embodiment of the method for face, specific implementation can be found in embodiment of the method, and details are not described herein.
It is apparent to those skilled in the art that for convenience and simplicity of description, the system of foregoing description
It with the specific work process of device, can refer to corresponding processes in the foregoing method embodiment, details are not described herein.
In addition, in the description of the embodiment of the present invention unless specifically defined or limited otherwise, term " installation ", " phase
Even ", " connection " shall be understood in a broad sense, for example, it may be being fixedly connected, may be a detachable connection, or be integrally connected；It can
To be mechanical connection, it is also possible to be electrically connected；It can be directly connected, can also can be indirectly connected through an intermediary
Connection inside two elements.For the ordinary skill in the art, above-mentioned term can be understood at this with concrete condition
Concrete meaning in invention.
It, can be with if the function is realized in the form of SFU software functional unit and when sold or used as an independent product
It is stored in a computer readable storage medium.Based on this understanding, technical solution of the present invention is substantially in other words
The part of the part that contributes to existing technology or the technical solution can be embodied in the form of software products, the meter
Calculation machine software product is stored in a storage medium, including some instructions are used so that a computer equipment (can be a
People's computer, server or network equipment etc.) it performs all or part of the steps of the method described in the various embodiments of the present invention.
And storage medium above-mentioned includes: that USB flash disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), arbitrary access are deposited
The various media that can store program code such as reservoir (RAM, Random Access Memory), magnetic or disk.
In the description of the present invention, it should be noted that term " center ", "upper", "lower", "left", "right", "vertical",
The orientation or positional relationship of the instructions such as "horizontal", "inner", "outside" be based on the orientation or positional relationship shown in the drawings, merely to
Convenient for description the present invention and simplify description, rather than the device or element of indication or suggestion meaning must have a particular orientation,
It is constructed and operated in a specific orientation, therefore is not considered as limiting the invention.In addition, term " first ", " second ",
" third " is used for descriptive purposes only and cannot be understood as indicating or suggesting relative importance.
The embodiment of the present invention also provides a kind of electronic equipment, including memory, processor and storage are on a memory and can
The computer program run on a processor, processor are realized provided by the above embodiment based on diffraction when executing computer program
The step of crack inversion method of wave energy.
The embodiment of the present invention also provides a kind of computer readable storage medium, and meter is stored on computer readable storage medium
Calculation machine program executes the crack inversion method based on diffraction wave energy of above-described embodiment when computer program is run by processor
The step of.
Finally, it should be noted that embodiment described above, only a specific embodiment of the invention, to illustrate the present invention
Technical solution, rather than its limitations, scope of protection of the present invention is not limited thereto, although with reference to the foregoing embodiments to this hair
It is bright to be described in detail, those skilled in the art should understand that: anyone skilled in the art
In the technical scope disclosed by the present invention, it can still modify to technical solution documented by previous embodiment or can be light
It is readily conceivable that variation or equivalent replacement of some of the technical features；And these modifications, variation or replacement, do not make
The essence of corresponding technical solution is detached from the spirit and scope of technical solution of the embodiment of the present invention, should all cover in protection of the invention
Within the scope of.
Claims (10)
1. a kind of crack inversion method based on diffraction wave energy characterized by comprising
The isolated diffraction wave energy from seismic data, the diffraction wave energy include prestack diffraction wave energy and poststack diffraction
Wave energy；
The initial model of crack attribute is determined according to the poststack diffraction wave energy；
Using the initial model and the prestack diffraction wave energy as input data, and inverting is carried out as follows, directly
Value to inversion objective function meets expection, using the initial model as the crack attribute after inverting:
By the input data input reverse objective function, the value of inversion objective function is obtained；
If the value of the inversion objective function does not meet expection, the initial model is adjusted, and by the prestack around
Ejected wave energy and initial model adjusted are as new input data.
2. the crack inversion method according to claim 1 based on diffraction wave energy, which is characterized in that according to the poststack
Diffraction wave energy determines the initial model of crack attribute, comprising:
Equivalent fissure porosity model is constructed according to the poststack diffraction wave energy；
The initial model for characterizing crack attribute is calculated according to the equivalent fissure porosity model.
3. the crack inversion method according to claim 2 based on diffraction wave energy, which is characterized in that according to the poststack
Diffraction wave energy constructs equivalent fissure porosity model, comprising:
The obliquity information and azimuth information in crack are obtained according to Image Logging Data；
Intensity information is generated according to the obliquity information and azimuth information；
According to the poststack diffracted wave Energy extraction fracture development information, it is based on the fracture development information and the intensity
Information constructs equivalent fissure porosity model in the way of Decided modelling and stochastic modeling respectively.
4. the crack inversion method according to claim 1 based on diffraction wave energy, which is characterized in that the inverting target
Function is determined according to angle gathers data and forward modeling data, wherein carries out extracting the angle road from the prestack diffraction wave energy
Collect data, the forward modeling data are obtained according to formation of elastic parameters and the initial model.
5. the crack inversion method according to claim 4 based on diffraction wave energy, which is characterized in that from the prestack around
The angle gathers data are extracted in ejected wave energy, comprising:
The common midpoint gather of adjacent predetermined number in the prestack diffraction wave energy is overlapped, channel set number is formed
According to；
Angle gathers data are extracted from the channel set data according to predetermined angle.
6. the crack inversion method according to claim 4 based on diffraction wave energy, which is characterized in that the stratum elasticity
Parameter further includes seismic wave spread speed and stratum media density in the earth formation, according to formation of elastic parameters and the initial model
Obtain forward modeling data, comprising:
Reflection coefficient is calculated by the initial model, the spread speed and the stratum media density；
Convolution calculating is carried out according to the reflection coefficient and seismic wavelet, obtains forward modeling data.
7. the crack inversion method according to claim 1 based on diffraction wave energy, which is characterized in that the method is also wrapped
It includes:
It is predicted according to the crack attribute fracture oily situation after the inverting.
8. the crack inversion method according to claim 1 based on diffraction wave energy, which is characterized in that the inverting target
It further include similarity constraint item in function, crack attribute after inverting of the similarity constraint item for calculating adjacent seismic channel
Average value.
9. the crack inversion method according to claim 8 based on diffraction wave energy, which is characterized in that the method is also wrapped
It includes:
The inversion objective function is indicated according to the following formula:
Wherein, SobsFor angle gathers data, S is forward modeling data, and A and B indicate that likeness coefficient matrix, m indicate initial model, λ2With
γ2Regularization parameter is represented, i is sampling number, and j is seismic channel number.
10. a kind of crack inverting device based on diffraction wave energy characterized by comprising
Separation module, being used for isolated diffraction wave energy, the diffraction wave energy from seismic data includes prestack diffracted wave
Energy and poststack diffraction wave energy；
Determining module, for determining the initial model of crack attribute according to the poststack diffraction wave energy；
Inverting module is used for using the initial model and the prestack diffraction wave energy as input data, and by walking as follows
It is rapid to carry out inverting, until the value of inversion objective function meets expection, using the initial model as the crack attribute after inverting:
By the input data input reverse objective function, the value of inversion objective function is obtained；
If the value of the inversion objective function does not meet expection, the initial model is adjusted, and by the prestack around
Ejected wave energy and initial model adjusted are as new input data.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910720254.3A CN110286410B (en) | 2019-08-06 | 2019-08-06 | Fracture inversion method and device based on diffracted wave energy |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910720254.3A CN110286410B (en) | 2019-08-06 | 2019-08-06 | Fracture inversion method and device based on diffracted wave energy |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110286410A true CN110286410A (en) | 2019-09-27 |
| CN110286410B CN110286410B (en) | 2020-06-09 |
Family
ID=68024733
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910720254.3A Active CN110286410B (en) | 2019-08-06 | 2019-08-06 | Fracture inversion method and device based on diffracted wave energy |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110286410B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111025383A (en) * | 2019-11-21 | 2020-04-17 | 徐州工程学院 | Method for qualitatively judging water filling condition of tunnel front karst cave based on diffracted transverse waves |
| CN111273354A (en) * | 2020-03-24 | 2020-06-12 | 中国矿业大学(北京) | Method and system for analyzing underground water change caused by underground mining |
| CN111880219A (en) * | 2020-08-07 | 2020-11-03 | 中国科学院地质与地球物理研究所 | Diffracted wave imaging method and device based on azimuth-dip angle gather |
| CN113987797A (en) * | 2021-10-28 | 2022-01-28 | 中国矿业大学(北京) | Crack flexibility parameter obtaining method and device, electronic equipment and storage medium |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6687659B1 (en) * | 2000-03-24 | 2004-02-03 | Conocophillips Company | Method and apparatus for absorbing boundary conditions in numerical finite-difference acoustic applications |
| CN101907725A (en) * | 2010-06-12 | 2010-12-08 | 中国石油天然气股份有限公司 | Crack prediction method and device |
| CN103399345A (en) * | 2013-08-16 | 2013-11-20 | 中国科学院地质与地球物理研究所 | Survey method and device for buried-hill crevice distribution |
| CN106405644A (en) * | 2016-08-30 | 2017-02-15 | 中国石油天然气股份有限公司 | Crack determination method and apparatus |
| CN106597545A (en) * | 2016-12-07 | 2017-04-26 | 中国矿业大学(北京) | Seismic prestack inversion method and device for horizontal cracks |
| CN108957526A (en) * | 2017-05-17 | 2018-12-07 | 中国石油化工股份有限公司 | The method for obtaining crack fused data body |
| CN109387835A (en) * | 2018-08-28 | 2019-02-26 | 国家能源投资集团有限责任公司 | Ground Penetrating Radar diffracted wave Weak Signal Extraction and electronic equipment |
-
2019
- 2019-08-06 CN CN201910720254.3A patent/CN110286410B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6687659B1 (en) * | 2000-03-24 | 2004-02-03 | Conocophillips Company | Method and apparatus for absorbing boundary conditions in numerical finite-difference acoustic applications |
| CN101907725A (en) * | 2010-06-12 | 2010-12-08 | 中国石油天然气股份有限公司 | Crack prediction method and device |
| CN103399345A (en) * | 2013-08-16 | 2013-11-20 | 中国科学院地质与地球物理研究所 | Survey method and device for buried-hill crevice distribution |
| CN106405644A (en) * | 2016-08-30 | 2017-02-15 | 中国石油天然气股份有限公司 | Crack determination method and apparatus |
| CN106597545A (en) * | 2016-12-07 | 2017-04-26 | 中国矿业大学(北京) | Seismic prestack inversion method and device for horizontal cracks |
| CN108957526A (en) * | 2017-05-17 | 2018-12-07 | 中国石油化工股份有限公司 | The method for obtaining crack fused data body |
| CN109387835A (en) * | 2018-08-28 | 2019-02-26 | 国家能源投资集团有限责任公司 | Ground Penetrating Radar diffracted wave Weak Signal Extraction and electronic equipment |
Non-Patent Citations (1)
| Title |
|---|
| 舒梦: "潜山裂缝型油气藏绕射波属性预测方法", 《CT理论与应用研究》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111025383A (en) * | 2019-11-21 | 2020-04-17 | 徐州工程学院 | Method for qualitatively judging water filling condition of tunnel front karst cave based on diffracted transverse waves |
| CN111025383B (en) * | 2019-11-21 | 2021-09-24 | 徐州工程学院 | Method for qualitatively judging water filling condition of tunnel front karst cave based on diffracted transverse waves |
| CN111273354A (en) * | 2020-03-24 | 2020-06-12 | 中国矿业大学(北京) | Method and system for analyzing underground water change caused by underground mining |
| CN111880219A (en) * | 2020-08-07 | 2020-11-03 | 中国科学院地质与地球物理研究所 | Diffracted wave imaging method and device based on azimuth-dip angle gather |
| CN111880219B (en) * | 2020-08-07 | 2021-06-08 | 中国科学院地质与地球物理研究所 | Diffracted wave imaging method and device based on azimuth-dip angle gather |
| US11327193B2 (en) | 2020-08-07 | 2022-05-10 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Method and device for imaging diffracted waves based on azimuth-dip angle gathers, and storage medium |
| CN113987797A (en) * | 2021-10-28 | 2022-01-28 | 中国矿业大学(北京) | Crack flexibility parameter obtaining method and device, electronic equipment and storage medium |
| CN113987797B (en) * | 2021-10-28 | 2022-06-14 | 中国矿业大学(北京) | Crack flexibility parameter obtaining method and device, electronic equipment and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110286410B (en) | 2020-06-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2920499C (en) | Stratigraphic function | |
| EP3011366B1 (en) | Seismic data processing | |
| CN110286410A (en) | Crack inversion method and device based on diffraction wave energy | |
| EP3069171B1 (en) | Seismic data processing | |
| US9852373B2 (en) | Properties link for simultaneous joint inversion | |
| CN101932954B (en) | Subsurface prediction method and system | |
| US11042676B2 (en) | Representing structural uncertainty in a mesh representing a geological environment | |
| CN108139499A (en) | The full wave field inversion of Q- compensation | |
| CN109425896A (en) | Dolomite oil and gas reservoir distribution prediction method and device | |
| EP2548052B1 (en) | System and method of 3d salt flank vsp imaging with transmitted waves | |
| CN102053263B (en) | Method for inspecting surface structure | |
| CN101999086A (en) | Method for determining seismic data quality | |
| CN104975851B (en) | For amplitude with the reservoir model optimization method of geophone offset variation road set analysis | |
| CN117546051A (en) | Method and system for seismic imaging using an S-wave velocity model and machine learning | |
| US9921324B2 (en) | Systems and methods employing upward beam propagation for target-oriented seismic imaging | |
| CN106990436B (en) | The recognition methods of karst collapse col umn and device | |
| CN104199088B (en) | Incident angle gather extraction method and system | |
| EP3359983B1 (en) | Seismic polynomial filter | |
| Barison et al. | Processing and interpretation of seismic reflection data from the Los Humeros super-hot geothermal system | |
| WO2020205169A1 (en) | Low-frequency seismic survey design | |
| CN106990433B (en) | Identification method for micro erosion groove in bump area | |
| CN104316959A (en) | Fluid identification based on equivalent fluid acoustic wave impedance | |
| CN116699678A (en) | Speed modeling method and device for multi-information fusion | |
| Thomassen | Full-waveform inversion studies | |
| MX2013015014A (en) | Methods and systems for quality control of seismic illumination maps. |
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 |