CN107942374B - Diffracted wave field extracting method and device - Google Patents
Diffracted wave field extracting method and device Download PDFInfo
- Publication number
- CN107942374B CN107942374B CN201711133853.2A CN201711133853A CN107942374B CN 107942374 B CN107942374 B CN 107942374B CN 201711133853 A CN201711133853 A CN 201711133853A CN 107942374 B CN107942374 B CN 107942374B
- Authority
- CN
- China
- Prior art keywords
- wave
- indicate
- object function
- inclination angle
- reflection slope
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
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
Abstract
The present invention provides a kind of diffracted wave field extracting method and devices, it is related to the technical field of diffracted wave field extraction, this method comprises: obtaining the prestack common offset trace gather data in pending area, wherein, prestack common offset trace gather data are to carry the data of bed boundary information in pending area；On the basis of carrying out decomposition of plane wave to prestack common offset trace gather data, L is used after converting by warp wavelet to back wave part inclination angle0Norm carries out regularization constraint, obtains first object function of the diffracted wave field to be extracted about back wave part inclination angle；Target is solved using Trust Region Algorithm and reflects slope of wave surface, wherein target reflection slope of wave surface is back wave part inclination angle when first object function reaches minimum value；Combining target reflection slope of wave surface, prestack common offset trace gather data and first object function determine diffracted wave field to be extracted.The present invention alleviates the poor technical problem of diffracted wave precision that traditional diffracted wave extracting method is extracted.
Description
Technical field
The present invention relates to the technical fields that diffracted wave field extracts, more particularly, to a kind of diffracted wave field extracting method and dress
It sets.
Background technique
Conventional seismic exploration technique when seismic wave is when the underground space is propagated, is encountered based on theoretical reflection
Lithology interface generates back wave, is received by earth's surface wave detector, therefore back wave is with carrying large scale according to Snell law
Plastid information can portray stratigraphic structure and trend in detail.Simultaneously as the interaction of ground interlayer crustal stress, the underground space
There is also the small scale geologic body such as tomography, crack, theoretical reflection cannot accurately portray small scale geologic body, but such geology
Body is often closely bound up with the energy such as petroleum, natural gas, coal bed gas, thus the exploration of small scale geologic body is most important.
Diffracted wave is the seismic response of small scale geologic body, and diffracted wave carries the tectonic information of small scale geologic body, because
And diffracted wave can be used to be accurately positioned discontinuously plastid heterogeneous, provide stronger underground space illuminance.But
It is that diffraction wave energy is lower, is almost submerged in back wave, it is difficult to identification and utilization, therefore diffracted wave separation is diffracted wave application
An essential step.In existing method, decomposition of plane wave method is a kind of common diffracted wave separation method, anti-by estimating
Ejected wave part inclination angle removes back wave.However, decomposition of plane wave method is when estimating back wave part inclination angle, often with solution
Unstability influence back wave pressing result so that the estimation of back wave part inclination angle is inaccurate so that the diffracted wave extracted
Precision is poor.
Summary of the invention
In view of this, the purpose of the present invention is to provide a kind of diffracted wave field extracting method and device, with alleviate tradition around
The poor technical problem of the diffracted wave precision that ejected wave extracting method is extracted.
In a first aspect, the embodiment of the invention provides a kind of diffracted wave field extracting methods, comprising:
Obtain the prestack common offset trace gather data in pending area, wherein the prestack common offset trace gather data
For the data for carrying bed boundary information in the pending area；
On the basis of carrying out decomposition of plane wave to the prestack common offset trace gather data, by warp wavelet to reflection
Wave part inclination angle uses L after being converted0Norm carries out regularization constraint, obtains diffracted wave field to be extracted about back wave part
The first object function at inclination angle；
It solves target and reflects slope of wave surface, wherein the target reflection slope of wave surface is that the first object function reaches minimum
Back wave part inclination angle when value；
It is determined in conjunction with target reflection slope of wave surface, the prestack common offset trace gather data and the first object function
The diffracted wave field to be extracted.
With reference to first aspect, the embodiment of the invention provides the first possible embodiments of first aspect, wherein
On the basis of carrying out decomposition of plane wave to the prestack common offset trace gather data, by warp wavelet to back wave part inclination angle
L is used after carrying out sparse transformation0Norm carries out regularization constraint, obtains diffracted wave field to be extracted about back wave part inclination angle
First object function, comprising:
Data fit term is established according to the prestack common offset trace gather data and back wave part inclination angle；
By warp wavelet and then L is passed through to the reflection slope of wave surface0The result that norm obtains the warp wavelet
It carries out regularization constraint and establishes regularization term；
The result that the data fit term is added with the regularization term, and will add up is as the first object letter
Number.
The possible embodiment of with reference to first aspect the first, the embodiment of the invention provides second of first aspect
Possible embodiment, wherein the formula of the first object function are as follows:
Wherein, r indicates the diffracted wave field to be extracted,Indicate the data fit term, α | | Ω (ρ) | |0Table
Show the regularization term, C indicates that decomposition of plane wave filter operator, d indicate the prestack common offset trace gather data, and Ω is indicated
Warp wavelet operator, ρ indicate back wave part inclination angle, | | | |2Indicate L2Norm, | | | |0Indicate L0Norm, α indicate tuning ginseng
Number.
With reference to first aspect, the embodiment of the invention provides the third possible embodiments of first aspect, wherein asks
It solves target and reflects slope of wave surface, comprising:
The Trust-region subproblem for constructing the first object function obtains the of the update variable quantity about reflection slope of wave surface
Two objective functions and constraint condition, wherein second objective function are as follows:It is described
Constraint condition is | | ξk||2≤Δk, whereinIndicate second objective function, ξkIndicate the solution of the Trust-region subproblem,
gkIndicate the gradient of the first object function,Indicate gkTransposition, BkIndicate the Hessian matrix of the first object function,
ΔkIndicate the Trust Region Radius of kth time iteration；
The Trust-region subproblem is solved, obtains the update variable quantity of the reflection slope of wave surface, and become based on the update
Change amount determines the target reflection slope of wave surface.
The third possible embodiment with reference to first aspect, the embodiment of the invention provides the 4th kind of first aspect
Possible embodiment, wherein solve the Trust-region subproblem, obtain the update variable quantity of the reflection slope of wave surface, comprising:
The Trust-region subproblem is solved using truncation conjugate gradient method, obtains the more new change of the reflection slope of wave surface
Amount.
The third possible embodiment with reference to first aspect, the embodiment of the invention provides the 5th kind of first aspect
Possible embodiment, wherein the method also includes:
In spatio-temporal domain, mean value smoothing processing is carried out to target reflection slope of wave surface, and with after smoothing processing
The target reflection slope of wave surface replacement smoothing processing foregoing description target reflect slope of wave surface.
Second aspect, the embodiment of the present invention also provide a kind of diffracted wave field extraction element, comprising: module are obtained, for obtaining
Take the prestack common offset trace gather data in pending area, wherein the prestack common offset trace gather data are to carry
State the data of bed boundary information in pending area；
Module is constructed, for passing through on the basis of carrying out decomposition of plane wave to the prestack common offset trace gather data
Warp wavelet uses L after converting to back wave part inclination angle0Norm carries out regularization constraint, obtains diffracted wave field to be extracted
First object function about back wave part inclination angle；
Module is solved, for solving target reflection slope of wave surface, wherein the target reflection slope of wave surface is the first object
Function reaches back wave part inclination angle when minimum value；
Determining module, in conjunction with target reflection slope of wave surface, the prestack common offset trace gather data and described the
One objective function determines the diffracted wave field to be extracted.
The embodiment of the present invention bring it is following the utility model has the advantages that
The diffracted wave field extracting method includes: the prestack common offset trace gather data obtained in pending area, wherein folded
Preceding common offset trace gather data are to carry the data of bed boundary information in pending area；To the prestack common offset
On the basis of trace gather data carry out decomposition of plane wave, used after carrying out sparse transformation to back wave part inclination angle by warp wavelet
L0Norm carries out regularization constraint, obtains first object function of the diffracted wave field to be extracted about back wave part inclination angle；It solves
Target reflects slope of wave surface, wherein target reflection slope of wave surface is back wave part inclination angle when first object function reaches minimum value；
Combining target reflection slope of wave surface, prestack common offset trace gather data and first object function determine diffracted wave field to be extracted.
The diffracted wave field extracting method uses L on the basis of decomposition of plane wave0Norm and warp wavelet are to back wave part
Inclination angle carries out sparsity constraints, wherein warp wavelet is as a kind of sparse transformation method, by the sparse earth's surface in back wave part inclination angle
It shows and；L0Norm directly describes the sparsity of data after sparse variation as the method for calculating nonzero element number in data,
Pass through L0Norm carries out regularization constraint to data after sparse variation, so that the solution that first object function obtains is more stable, from
And reach depressor reflex wave, the purpose of diffracted wave field is extracted, the diffracted wave precision that traditional diffracted wave extracting method is extracted is alleviated
Poor technical problem.
Other features and advantages of the present invention will illustrate in the following description, also, partly become from specification
It obtains it is clear that understand through the implementation of the invention.The objectives and other advantages of the invention are in specification, claims
And specifically noted structure is achieved and obtained in attached drawing.
To enable the above objects, features and advantages of the present invention to be clearer and more comprehensible, preferred embodiment is cited below particularly, and cooperate
Appended attached drawing, is described in detail below.
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 flow chart for diffracted wave field extracting method that the embodiment of the present invention one provides；
Fig. 2 is the flow chart that a kind of solution target that the embodiment of the present invention one provides reflects slope of wave surface method；
Fig. 3 is a kind of structural block diagram of diffracted wave field extraction element provided by Embodiment 2 of the present invention；
Fig. 4 is the structural block diagram of another diffracted wave field extraction element provided by Embodiment 2 of the present invention.
Icon: 100- obtains module；200- constructs module；300- solves module；400- determining module；The smooth mould of 500-
Block.
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.
In existing method, decomposition of plane wave method is a kind of common diffracted wave separation method, by estimating back wave office
Portion inclination angle, to remove back wave.However, decomposition of plane wave method is when estimating back wave part inclination angle, often not with solution
Stability influences back wave pressing result so that the estimation of back wave part inclination angle is inaccurate, so that the diffracted wave precision extracted
It is poor.Based on this, a kind of diffracted wave field extracting method provided in an embodiment of the present invention and device are based on L0Norm and bent wave become
Back wave part inclination angle can be effectively estimated in the regularization constraint scheme changed, and reach depressor reflex and extract the beneficial of diffraction
Effect alleviates the poor technical problem of the diffracted wave precision of traditional diffracted wave extracting method extraction.For convenient for the present embodiment
Understood, is described in detail first to a kind of diffracted wave field extracting method disclosed in the embodiment of the present invention.
Embodiment one
A kind of diffracted wave field extracting method provided in an embodiment of the present invention, as shown in Figure 1, comprising:
Step S102 obtains the prestack common offset trace gather data in pending area, wherein prestack common offset trace gather
Data are to carry the data of bed boundary information in pending area.
In embodiments of the present invention, common offset seismic channel set data are the seismic wave detected by acquiring unit, tool
Body process are as follows: seismic wave is excited in shot point using manual method, seismic wave is to all the winds propagated, and is encountering underground not
With lithology decomposition face when produce back wave and diffracted wave, these back waves and diffracted wave return to ground again and cause ground
Then vibration detects the seismic wave for causing ground vibration along the place setting acquiring unit equidistant with shot point.
It should be noted that since the seismic wave detected is to receive the transformation of subsurface formations medium, ground stratum boundary
Face information includes at least geological structure and formation lithology, and the data for carrying bed boundary information include time, speed, frequency.
Step S104 passes through warp wavelet on the basis of carrying out decomposition of plane wave to prestack common offset trace gather data
L is used after converting to back wave part inclination angle0Norm carries out regularization constraint, obtains diffracted wave field to be extracted about reflection
The first object function at wave part inclination angle.
Specifically, current decomposition of plane wave method removes back wave simply by estimation back wave part inclination angle；This
The diffracted wave field extracting method provided in inventive embodiments is carrying out the same of decomposition of plane wave to prestack common offset trace gather data
When, sparse transformation is carried out to back wave part inclination angle by warp wavelet, and warp wavelet result is used into L0Norm carries out just
Then change constraint, to construct first object function of the diffracted wave field to be extracted about back wave part inclination angle.
Step S106 solves target and reflects slope of wave surface, wherein target reflection slope of wave surface is that first object function reaches minimum
Back wave part inclination angle when value.
Step S108, combining target reflects slope of wave surface, prestack common offset trace gather data and first object function determine to
Extract diffracted wave field.
Specifically, target reflection slope of wave surface is substituted into first object function, with determination diffracted wave field to be extracted.
In embodiments of the present invention, warp wavelet is as a kind of sparse transformation method, sparsely by back wave part inclination angle
It shows；L0Norm directly describes the sparse of data after sparse variation as the method for calculating nonzero element number in data
Property, pass through L0Norm carries out regularization constraint to data after warp wavelet, so that the solution that first object function obtains can be more steady
It is fixed, alleviate the poor technical problem of the diffracted wave precision of traditional diffracted wave extracting method extraction.
Wherein, on the basis of carrying out decomposition of plane wave to prestack common offset trace gather data, by warp wavelet to anti-
Ejected wave part inclination angle uses L after carrying out sparse transformation0Norm carries out regularization constraint, obtains diffracted wave field to be extracted about reflection
The first object function at wave part inclination angle, comprising:
Data fit term is established according to prestack common offset trace gather data and back wave part inclination angle；
By warp wavelet and then L is passed through to reflection slope of wave surface0Norm carries out canonical to the result that warp wavelet obtains
Change constraint and establishes regularization term；
The result that data fit term is added with regularization term, and will add up is as first object function.
Specifically, the formula for the first object function established are as follows:
Wherein, r indicates diffracted wave field to be extracted,Indicate data fit term, α | | Ω (ρ) | |0Indicate regularization
, C indicates that decomposition of plane wave filter operator, d indicate prestack common offset trace gather data, and Ω indicates warp wavelet operator, ρ table
Show back wave part inclination angle, | | | |2Indicate L2Norm, | | | |0Indicate L0Norm, α indicate tuner parameters.
In another optional embodiment of the embodiment of the present invention, as shown in Fig. 2, solving target reflects slope of wave surface, comprising:
Step S201 constructs the Trust-region subproblem of first object function, obtains the more new change about reflection slope of wave surface
The second objective function and constraint condition of amount, wherein the second objective function are as follows:About
Beam condition is | | ξk||2≤Δk, whereinIndicate the second objective function, ξkIndicate the solution of Trust-region subproblem, gkIndicate first
The gradient of objective function,Indicate gkTransposition, BkIndicate the Hessian matrix of first object function, ΔkIndicate kth time iteration
Trust Region Radius.
It should be noted that gkThe gradient for indicating first object function, that is, refer to: gkIndicate first object function to back wave
The first derivative at inclination angle；BkThe Hessian matrix for indicating first object function, that is, refer to: BkIndicate that first object function inclines to back wave
The second dervative at angle.
Step S202 solves Trust-region subproblem, obtains the update variable quantity of reflection slope of wave surface, and based on update variable quantity
Determine that target reflects slope of wave surface.
Specifically, the process of above-mentioned solution target reflection slope of wave surface uses Trust Region Algorithm, Trust Region Algorithm it is specific
Steps are as follows for execution:
A, reflection slope of wave surface initial value and Trust Region Radius initial value are obtained；
B, from reflection slope of wave surface initial value ρ1With Trust Region Radius initial value Δ1Start iteration, calculates first object function pair
The first derivative and second dervative at inclination angle, i.e. gradient gkWith Hessian matrix Bk, and it is based on gk、BkWithConstruct the second objective function
C, in conjunction with constraint condition | | ξk||2≤ΔkSolve the solution ξ of the second objective functionk；
D, at the kth iteration, it is based on ρkCalculate rk, and by ρkAnd ξkSum as first object argument of function come
Calculate rk+1；
E, it is based on ξkIt calculates
F, it is based on rk、rk+1WithMetric parameter p is calculated using following formulak:
G, according to pkThe size of value, according to reflection slope of wave surface field more new formula and Trust Region Radius more new formula to trusted zones
Radius and reflection slope of wave surface field are updated, and obtain the reflection slope of wave surface updated value ρ exported after kth time iterative calculationk+1, letter
Rely domain radius updated value Δk+1, wherein reflection slope of wave surface field more new formula are as follows:
Trust Region Radius more new formula are as follows:
η1Indicate the first preset threshold, η2Indicate the second preset threshold, σ1Indicate the first adjustment parameter, σ2Indicate that second adjusts
Save parameter, σ3Indicate third adjustment parameter, ΔmaxIndicate the default maximum of Trust Region Radius,
Wherein, successively value is 1 to K to k, and K indicates that the process of Trust Region Algorithm solution target reflection slope of wave surface reaches convergence
When the number of iterations, the ρ acquired after kth iterationK+1As target reflects slope of wave surface.
In the embodiment of the present invention, reflection slope of wave surface is calculated using Trust Region Algorithm, thus preferably depressor reflex wave, into one
Step is conducive to the high-fidelity of the diffracted wave field extracted.
In another optional embodiment of the embodiment of the present invention, trusted zones is solved by using truncation conjugate gradient method
Problem obtains the update variable quantity of reflection slope of wave surface.
Specifically, that is, above-mentioned steps c, in conjunction with constraint condition | | ξk||2≤ΔkSolve the solution ξ of the second objective functionkTo adopt
It is solved with truncation conjugate gradient method, following formula can be used and carry out detailed solution:
Wherein, ξk(0)For the initial value of the update variable quantity of preset reflection slope of wave surface, q(j)For more new direction；It indicates
The gradient of second objective function；δjAnd βjIndicate step-length regulatory factor；T indicates that iteration front and back is twice in truncation conjugate gradient method
Error,
Wherein, successively value is 1 to J to j, and when J is that t reaches t < σ (σ is default error threshold) for the first time, conjugate gradient is truncated
The number of iterations of method；The update variable quantity ξ for the reflection slope of wave surface that the J times iteration of truncation conjugate gradient method obtainsk(J+1)As believe
Rely the solution ξ of the second objective function in the kth time iteration of domain algorithmk。
In another optional embodiment of the embodiment of the present invention, diffracted wave field extracting method further include:
In spatio-temporal domain, mean value smoothing processing is carried out to target reflection slope of wave surface, and with the mesh after smoothing processing
Target reflects slope of wave surface before mark reflection slope of wave surface replacement smoothing processing.
In the embodiment of the present invention, target reflection slope of wave surface is subjected to spatio-temporal domain mean value smoothing, then by smoothing processing
Target reflection slope of wave surface afterwards substitutes into objective function, and then eliminates back wave, extracts diffracted wave field, above-mentioned time-space domain smoothing processing
Target reflection slope of wave surface after smoothing processing is substituted into objective function and carries out diffraction by the consecutive variations property for being conducive to keep wave field
Wave field extracts, and help to obtain the diffracted wave field of high-fidelity.
Embodiment two
As shown in Fig. 3 to Fig. 4, a kind of diffracted wave field extraction element provided in an embodiment of the present invention.
Referring to Fig. 3, which includes:
Module 100 is obtained, for obtaining the prestack common offset trace gather data in pending area, wherein prestack is inclined altogether
Moving away from trace gather data is to carry the data of bed boundary information in pending area；
Module 200 is constructed, for passing through song on the basis of carrying out decomposition of plane wave to prestack common offset trace gather data
Wave conversion uses L after converting to back wave part inclination angle0Norm carries out regularization constraint, obtains diffracted wave field to be extracted and closes
First object function in back wave part inclination angle；
Module 300 is solved, for solving target reflection slope of wave surface, wherein it is first object function that target, which reflects slope of wave surface,
Reach back wave part inclination angle when minimum value；
Determining module 400, for combining target reflection slope of wave surface, prestack common offset trace gather data and first object function
Determine diffracted wave field to be extracted.
In embodiments of the present invention, it obtains module 100 and obtains prestack common offset trace gather data, building module 200 uses
L0Norm and warp wavelet carry out sparsity constraints to back wave part inclination angle, wherein warp wavelet is as a kind of sparse transformation
Method sparsely shows at back wave part inclination angle；L0Norm is as the method for calculating nonzero element number in data, directly
The sparsity for connecing data after describing sparse variation, passes through L0Norm carries out regularization constraint to data after sparse variation, so that the
One objective function is more stable by solving the obtained solution of module 300, then 400 combining target of determining module reflection slope of wave surface,
The diffracted wave field precision to be extracted that prestack common offset trace gather data and first object function are determined is preferable, alleviate tradition around
The poor technical problem of the diffracted wave precision that ejected wave extracting method is extracted.
In one optional embodiment of the embodiment of the present invention, building module 200 is used for:
Data fit term is established according to prestack common offset trace gather data and back wave part inclination angle；
By warp wavelet and then L is passed through to reflection slope of wave surface0Norm carries out canonical to the result that warp wavelet obtains
Change constraint and establishes regularization term；
The result that data fit term is added with regularization term, and will add up is as first object function.
In another optional embodiment of the embodiment of the present invention, the formula of first object function are as follows:
Wherein, r indicates diffracted wave field to be extracted,Indicate data fit term, α | | Ω (ρ) | |0Indicate regularization
, C indicates that decomposition of plane wave filter operator, d indicate prestack common offset trace gather data, and Ω indicates warp wavelet operator, ρ table
Show back wave part inclination angle, | | | |2Indicate L2Norm, | | | |0Indicate L0Norm, α indicate tuner parameters.
In another optional embodiment of the embodiment of the present invention, module 300 is solved, is used for:
The Trust-region subproblem for constructing the first object function obtains the of the update variable quantity about reflection slope of wave surface
Two objective functions and constraint condition, wherein the second objective function are as follows:Constraint condition
For | | ξk||2≤Δk, whereinIndicate the second objective function, ξkIndicate the solution of Trust-region subproblem, gkIndicate first object letter
Several gradients,Indicate gkTransposition, BkIndicate the Hessian matrix of first object function, ΔkIndicate the trusted zones of kth time iteration
Radius；
Trust-region subproblem is solved, obtains the update variable quantity of reflection slope of wave surface, and determine target based on variable quantity is updated
Reflect slope of wave surface.
In another optional embodiment of the embodiment of the present invention, solves module 300 and is used for:
Trust-region subproblem is solved using truncation conjugate gradient method, obtains the update variable quantity of reflection slope of wave surface.
In another optional embodiment of the embodiment of the present invention, diffracted wave field extraction element further include:
Leveling Block 500, for carrying out mean value smoothing processing to target reflection slope of wave surface, being used in combination in spatio-temporal domain
Target reflects slope of wave surface before target reflection slope of wave surface replacement smoothing processing after smoothing processing.
The computer program product of diffracted wave field extracting method and device provided by the embodiment of the present invention, including store
The computer readable storage medium of program code, the instruction that program code includes can be used for executing the side in previous methods embodiment
Method, 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 " are used for description purposes only, it is not understood to indicate or imply phase
To importance.
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.Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (6)
1. a kind of diffracted wave field extracting method characterized by comprising
Obtain the prestack common offset trace gather data in pending area, wherein the prestack common offset trace gather data are to take
Data with bed boundary information in the pending area；
On the basis of carrying out decomposition of plane wave to the prestack common offset trace gather data, by warp wavelet to back wave office
Portion inclination angle uses L after being converted0Norm carries out regularization constraint, obtains diffracted wave field to be extracted about back wave part inclination angle
First object function；
It solves target and reflects slope of wave surface, wherein when the target reflection slope of wave surface is that the first object function reaches minimum value
Back wave part inclination angle；
In conjunction with described in target reflection slope of wave surface, the prestack common offset trace gather data and first object function determination
Diffracted wave field to be extracted；
Wherein, on the basis of carrying out decomposition of plane wave to the prestack common offset trace gather data, by warp wavelet to anti-
Ejected wave part inclination angle uses L after carrying out sparse transformation0Norm carries out regularization constraint, obtains diffracted wave field to be extracted about reflection
The first object function at wave part inclination angle, comprising:
Data fit term is established according to the prestack common offset trace gather data and back wave part inclination angle；
By warp wavelet and then L is passed through to the reflection slope of wave surface0Norm carries out the result that the warp wavelet obtains
Regularization constraint and establish regularization term；
The result that the data fit term is added with the regularization term, and will add up is as the first object function；
Wherein, the formula of the first object function are as follows:
Wherein, r indicates the diffracted wave field to be extracted,Indicate the data fit term, α | | Ω (ρ) | |0Indicate institute
Regularization term is stated, C indicates that decomposition of plane wave filter operator, d indicate the prestack common offset trace gather data, and Ω indicates Qu Bo
Transformation operator, ρ indicate back wave part inclination angle, | | | |2Indicate L2Norm, | | | |0Indicate L0Norm, α indicate tuner parameters；
Wherein, target reflection slope of wave surface is substituted into first object function, with determination diffracted wave field to be extracted.
2. the method according to claim 1, wherein solving target reflects slope of wave surface, comprising:
The Trust-region subproblem for constructing the first object function obtains the second mesh of the update variable quantity about reflection slope of wave surface
Scalar functions and constraint condition, wherein second objective function are as follows:The constraint
Condition is | | ξk||2≤Δk, whereinIndicate second objective function, ξkIndicate the solution of the Trust-region subproblem, gkTable
Show the gradient of the first object function,Indicate gkTransposition, BkIndicate the Hessian matrix of the first object function, Δk
Indicate the Trust Region Radius of kth time iteration；
The Trust-region subproblem is solved, obtains the update variable quantity of the reflection slope of wave surface, and be based on the update variable quantity
Determine the target reflection slope of wave surface.
3. according to the method described in claim 2, obtaining the back wave it is characterized in that, solve the Trust-region subproblem
The update variable quantity at inclination angle, comprising:
The Trust-region subproblem is solved using truncation conjugate gradient method, obtains the update variable quantity of the reflection slope of wave surface.
4. according to the method described in claim 2, it is characterized in that, the method also includes:
In spatio-temporal domain, mean value smoothing processing is carried out to target reflection slope of wave surface, and with the institute after smoothing processing
State target reflection slope of wave surface replacement smoothing processing foregoing description target reflection slope of wave surface.
5. a kind of diffracted wave field extraction element characterized by comprising
Module is obtained, for obtaining the prestack common offset trace gather data in pending area, wherein the prestack common offset
Trace gather data are to carry the data of bed boundary information in the pending area；
Module is constructed, for passing through Qu Bo on the basis of carrying out decomposition of plane wave to the prestack common offset trace gather data
Transformation uses L after converting to back wave part inclination angle0Norm carry out regularization constraint, obtain diffracted wave field to be extracted about
The first object function at back wave part inclination angle；
Module is solved, for solving target reflection slope of wave surface, wherein the target reflection slope of wave surface is the first object function
Reach back wave part inclination angle when minimum value；
Determining module, in conjunction with target reflection slope of wave surface, the prestack common offset trace gather data and first mesh
Scalar functions determine the diffracted wave field to be extracted；
Wherein, the building module is used for:
Data fit term is established according to the prestack common offset trace gather data and back wave part inclination angle；
By warp wavelet and then L is passed through to the reflection slope of wave surface0Norm carries out the result that the warp wavelet obtains
Regularization constraint and establish regularization term；
The result that the data fit term is added with the regularization term, and will add up is as the first object function；
Wherein, the formula of the first object function are as follows:
Wherein, r indicates the diffracted wave field to be extracted,Indicate the data fit term, α | | Ω (ρ) | |0Indicate institute
Regularization term is stated, C indicates that decomposition of plane wave filter operator, d indicate the prestack common offset trace gather data, and Ω indicates Qu Bo
Transformation operator, ρ indicate back wave part inclination angle, | | | |2Indicate L2Norm, | | | |0Indicate L0Norm, α indicate tuner parameters；
Wherein, target reflection slope of wave surface is substituted into first object function, with determination diffracted wave field to be extracted.
6. device according to claim 5, which is characterized in that the solution module is used for:
The Trust-region subproblem for constructing the first object function obtains the second mesh of the update variable quantity about reflection slope of wave surface
Scalar functions and constraint condition, wherein second objective function are as follows:The constraint
Condition is | | ξk||2≤Δk, whereinIndicate second objective function, ξkIndicate the solution of the Trust-region subproblem, gkTable
Show the gradient of the first object function,Indicate gkTransposition, BkIndicate the Hessian matrix of the first object function, ΔkTable
Show the Trust Region Radius of kth time iteration；
The Trust-region subproblem is solved, obtains the update variable quantity of the reflection slope of wave surface, and be based on the update variable quantity
Determine the target reflection slope of wave surface.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711133853.2A CN107942374B (en) | 2017-11-15 | 2017-11-15 | Diffracted wave field extracting method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711133853.2A CN107942374B (en) | 2017-11-15 | 2017-11-15 | Diffracted wave field extracting method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107942374A CN107942374A (en) | 2018-04-20 |
| CN107942374B true CN107942374B (en) | 2019-03-15 |
Family
ID=61932466
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711133853.2A Active CN107942374B (en) | 2017-11-15 | 2017-11-15 | Diffracted wave field extracting method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107942374B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109490951B (en) * | 2018-12-03 | 2019-11-15 | 中国矿业大学(北京) | Diffraction wave imaging method, device and electronic equipment |
| CN111929729B (en) | 2020-08-20 | 2021-04-06 | 中国矿业大学(北京) | Diffracted wave imaging method and device and electronic equipment |
| CN114114420B (en) * | 2020-09-01 | 2024-02-23 | 中国石油化工股份有限公司 | Diffraction identification imaging method, diffraction identification imaging device, electronic equipment and medium |
| CN113640872B (en) * | 2021-08-12 | 2022-03-08 | 中国矿业大学(北京) | Diffracted wave separation method and device and electronic equipment |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102455439B (en) * | 2010-11-02 | 2013-10-23 | 中国石油大学(北京) | Diffracted wave field separation method based on Kirchhoff integral method |
| CN102778693B (en) * | 2011-05-13 | 2014-09-10 | 中国石油化工股份有限公司 | Diffracted wave separation processing method based on reflection wave layer leveling extraction and elimination |
| CN102520444B (en) * | 2011-12-13 | 2014-10-01 | 中国科学院地质与地球物理研究所 | Diffraction wave information extraction method in post-stack seismic wave |
| CN103675897B (en) * | 2012-08-30 | 2017-10-27 | 中国石油化工股份有限公司 | A kind of earthquake diffracted wave method for separate imaging |
| CN106443785B (en) * | 2016-11-03 | 2018-05-29 | 中国矿业大学(北京) | Diffracted wave field acquisition methods and device |
-
2017
- 2017-11-15 CN CN201711133853.2A patent/CN107942374B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN107942374A (en) | 2018-04-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107942374B (en) | Diffracted wave field extracting method and device | |
| CN105182444B (en) | A kind of Inversion of geostatistics of High-Resolution Sequence Stratigraphic Framework constraint | |
| CN106443785B (en) | Diffracted wave field acquisition methods and device | |
| CN107272062B (en) | A kind of Q estimation methods of underground medium of data-driven | |
| SA520412321B1 (en) | Determining Sand-Dune Velocity Variations | |
| CN109709607B (en) | Method and device for predicting thickness of thin sandstone reservoir | |
| CN104459782B (en) | Horizon chromatography inversion earthquake speed modeling method with thin layers used and modeling unit | |
| CN110579802B (en) | High-precision inversion method for physical property parameters of natural gas hydrate reservoir | |
| AU2019340410B2 (en) | Reservoir characterization utilizing inversion of resampled seismic data | |
| EP2810101A1 (en) | Improving efficiency of pixel-based inversion algorithms | |
| CN109425900B (en) | Seismic reservoir prediction method | |
| CN101201409B (en) | Method for revising earthquake data phase | |
| CN113740901B (en) | Land seismic data full-waveform inversion method and device based on complex undulating surface | |
| CN105549084B (en) | A kind of three-dimensional high-precision velocity modeling method and system | |
| CN106291677A (en) | A kind of poststack sound impedance inversion method based on match tracing method | |
| CN103048690B (en) | Based on the stratum reflection pickup technology of the Rapid matching Projective decomposition of optimum seismic wavelet extraction | |
| CN102879814A (en) | Accurate depth domain layer speed updating method | |
| CN103439740A (en) | Method and device for predicting relative impedance based on dipole seismic wavelet multiple integrals | |
| CN109001840B (en) | A kind of mud shale genetic analysis method and device | |
| CN106291682A (en) | A kind of poststack sound impedance inversion method based on base method for tracing | |
| Cabello et al. | An outcrop-based comparison of facies modelling strategies in fan-delta reservoir analogues from the Eocene Sant Llorenç del Munt fan-delta (NE Spain) | |
| CN106526679A (en) | Method for recovering non-uniform erosion thickness through utilizing cycle analysis method | |
| CN106650192A (en) | Volcanic type uranium ore deposit magnetic interface retrieval method | |
| Cheng et al. | Compressible streamlines and three-phase history matching | |
| CN105301638B (en) | A kind of method and apparatus in extraction base of weathering face |
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 |