CN106569262B - Background velocity model reconstruction method under low frequency seismic data missing - Google Patents
Background velocity model reconstruction method under low frequency seismic data missing Download PDFInfo
- Publication number
- CN106569262B CN106569262B CN201510655283.8A CN201510655283A CN106569262B CN 106569262 B CN106569262 B CN 106569262B CN 201510655283 A CN201510655283 A CN 201510655283A CN 106569262 B CN106569262 B CN 106569262B
- Authority
- CN
- China
- Prior art keywords
- seismic data
- envelope
- velocity model
- low frequency
- inverting
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The background velocity model reconstruction method being related under a kind of low frequency seismic data missing, this approach includes the following steps:1) envelope of the simulation seismic data according to the envelope of observation seismic data and based on initial velocity model, the phase-fitting mutually related objects functional of structure adding window normalization nonzero-lag;2) it is based on the mutually related objects functional and builds corresponding inverse time propagation with focus;3) it is solved with wave field using the adjoint focus backpropagation, and carries out envelope inverting, to which the initial velocity model be reconstructed, obtain the rate pattern after being reconstructed under low-frequency data deletion condition;4) precision for verifying envelope inverting carries out traditional wave equation chromatography in the case where precision is met the requirements using the rate pattern after reconstruct as initial velocity model.
Description
Technical field
The present invention relates to oil seismic exploration field, the background velocity under being lacked more particularly, to low frequency seismic data
Model reconstruction method.
Background technology
Traditional Full wave shape chromatographic technique based on wave equation utilizes earthquake data before superposition kinematics and power due to it
Information is learned, waveform, amplitude and the phase information of seismic wave are matched in refutation process, thus there is higher resolution ratio, is recognized
To be the highest method of velocity modeling precision.Due to all information of matching seismic wave, the earthquake number recorded using complex model
When according to carrying out inverting, which is the nonlinear problem of a height, cost functional often there is a large amount of local minimum,
The precision of inverting depends critically upon the precision (Bunk etc. 1995) of initial model.In order to weaken the non-linear property of inverting, Bunks etc.
(1995), Sigure and Pratt (2004) give from the angle in time and space domain and frequency space domain from low frequency to height respectively
The Multi-scale inversion strategy of frequency.Baeten etc. (2013) points out the data of 1.5Hz to the 2.0Hz in actual seismic data for extensive
Multiple underground medium macroscopic model plays a crucial role.However often to lack 5Hz below low for actual seismic data
Frequency information or low-frequency range unreliable information, thus multiple dimensioned Full wave shape chromatographic technique is carried out using the seismic data of low frequency missing
Also inversion result can not be made to converge to global minimum.
Wu and Luo (2013) goes to solution Conventional Time domain in field of seismic exploration using the envelope inverting of data residual error form
Inversion result instability problem caused by wave equation chromatography is lacked due to seismic data low frequency.In seismic data acquisition process
In, due to by source level, focus response, surface layer attenuation by absorption, various noise jamming, mechanism of transmission, geologic(al) factor and detection
The influences such as device coupling effect, land seismic data focus and geophone station energy unbalanced phenomena are extremely serious, even if using earth's surface
Consistency amplitude compensation, the energy of different focus is also up to an order of magnitude.Envelope using various amplitude amount of energy grade is residual
Poor inverting can make the solution of envelope inverting unstable, and when the time delay of envelope is more than the half period of envelope signal, envelope
Inverting is similarly faced with the risk for being absorbed in local extremum.
Invention content
For this purpose, the present invention provides a kind of envelope inverting based on normalization nonzero-lag phase-fitting, to reconstruct initially
Macroscopic rate pattern weakens traditional wave equation and chromatographs the defect relied on low frequency height, and keeps the technology more suitable
For the inverting of land seismic data, its robustness is improved.
On the one hand propose the background velocity model reconstruction method under a kind of low frequency seismic data missing, this method include with
Lower step:1) envelope of the simulation seismic data according to the envelope of observation seismic data and based on initial velocity model, structure add
Window normalizes the phase-fitting mutually related objects functional of nonzero-lag;2) it is corresponding inverse to be based on mutually related objects functional structure
When propagate with focus;3) it is solved with wave field using the adjoint focus backpropagation, and carries out envelope inverting, to institute
It states initial velocity model to be reconstructed, obtains the rate pattern after being reconstructed under low-frequency data deletion condition;4) envelope inverting is verified
Precision carry out traditional fluctuation using the rate pattern after reconstruct as initial velocity model in the case where precision is met the requirements
Equation chromatographs.
Steady macroscopic speed that the invention belongs to oil seismic exploration fields under low frequency seismic data deletion condition
Modeling method.The heavy dependence in practical land data high precision velocity modeling process is chromatographed for Conventional temporal domain wave equation
Seismic data low-frequency information, and due to earth's surface is inconsistent, focus and wave detector energy it is unbalanced caused by inverting instability problem,
The present invention provides the envelope data inversion method that nonzero-lag phase-fitting is normalized based on adding window, weakens traditional wave equation layer
The defect relied on low frequency height is analysed, and the technology is made to be more applicable for the inverting of land seismic data, improves its robustness.
Description of the drawings
Disclosure illustrative embodiments are described in more detail in conjunction with the accompanying drawings, the disclosure above-mentioned and its
Its purpose, feature and advantage will be apparent, wherein in disclosure illustrative embodiments, identical reference label
Typically represent same parts.
Fig. 1 shows the background velocity model reconstruction side under a kind of low frequency seismic data missing according to the ... of the embodiment of the present invention
The flow chart of method;
Fig. 2 shows low frequency missing Ricker wavelet and its schematic diagrames of envelope;
Fig. 3 shows that low frequency lacks wavelet forward modeling big gun collection;
Fig. 4 shows big gun collection envelope natural logrithm;
Fig. 5 shows true Marmousi models;
Fig. 6 shows normal gradient initial model;
Fig. 7 shows traditional L2 norms wave equation tomographic inversion result;
Fig. 8 shows that adding window nonzero-lag normalizes phase-fitting envelope inversion result;
Fig. 9 shows envelope inverting+tradition L2 norm wave equation tomographic inversion results;
Figure 10 shows observation seismic data envelope and observes the normalized autocorrelation of seismic data envelope;
Figure 11 shows the normalized crosscorrelation of observation seismic data envelope and initial model analogue data envelope;
Figure 12 shows the normalized crosscorrelation of observation seismic data envelope and envelope inverse model analogue data envelope.
Specific implementation mode
The preferred embodiment of the disclosure is more fully described below with reference to accompanying drawings.Although showing the disclosure in attached drawing
Preferred embodiment, however, it is to be appreciated that may be realized in various forms the disclosure without the embodiment party that should be illustrated here
Formula is limited.On the contrary, these embodiments are provided so that the disclosure is more thorough and complete, and can be by the disclosure
Range is completely communicated to those skilled in the art.
Fig. 1 shows the background velocity model reconstruction side under a kind of low frequency seismic data missing according to the ... of the embodiment of the present invention
The flow chart of method, this approach includes the following steps.
Step 101, the envelope of the simulation seismic data according to the envelope of observation seismic data and based on initial velocity model,
Build the phase-fitting mutually related objects functional of adding window normalization nonzero-lag:
In formula, χXcorNenlogCan be adding window normalizing for the normalized crosscorrelation cost functional of earthquake envelope natural logrithm
Change the phase-fitting mutually related objects functional of nonzero-lag, XcorNen1og (τ, h) is observation seismic data envelope natural logrithm
With the Normalized Cross Correlation Function of simulation seismic data envelope natural logrithm.P (τ) is weighting time histories sample, and τ is time shift amount, u
(t, h) is observation seismic data, and y (t, h) is simulation seismic data, uH(t, h) is the Hilbert transform for observing seismic data,
yH(t, h) is the Hilbert transform for simulating seismic data, eobs(t, h) is the envelope for observing seismic data, esyn(t, h) is mould
The envelope of quasi- seismic data, c are in order to avoid the constant that natural logrithm antilog is zero introducing, and h indicates the serial number of wave detector.
Envelope delay time using the permissible observation data of the cost functional and analogue data is more than half of envelope period, anti-to make
It is more steady to drill, and avoids the requirement to phase delay in half period based on envelope residual error.
Step 102, corresponding inverse time propagation is built with focus according to the mutually related objects functional that step 101 provides.
According to the research of Tromp etc. (2005) et al., for different Data Matchings, no matter being to solve for the residual of two signals
The solution expression formula of the similarity degree of difference or two signals, gradient is identical, is all to use forward-propagating wave field and inverse time
The adjoint wave field cross-correlation superposition propagated acquires.Unlike unique, for different Signal Matching modes, adjoint equation right end
Adjoint focus item differ.The present embodiment can realize that the sound wave of use is closed based on the ACOUSTIC WAVE EQUATION in isotropic medium
It is in the pressure gradient expression formula of rate pattern:
Wherein Ns indicates the number of focus, vpFor rate pattern, ρ is underground medium model,For the particle of forward-propagating
Vibration velocity wavefield component, Q are the adjoint pressure wavefield component propagated the inverse time, and is indicates focus index, xiIndicate x, y, the directions z
Coordinate, T indicate the maximum moment propagated.
In order to accurately solve gradient of the normalization nonzero-lag envelope mutually related objects functional about rate pattern, need
It accurately solves with wave field.And solved by adjoint equation with wave field, it is consistent with the ACOUSTIC WAVE EQUATION of forward-propagating, only
The focus item of right end unlike one.Therefore it accurately derives corresponding very necessary with focus.It can be by target with focus
Functional is solved about the gradient of simulation wave field:
It can be obtained by a series of mathematical derivations:
WhereinTo observe the packet of seismic data
The cross-correlation of the envelope of network and simulation seismic data.Hilbert indicates the Hilbert transform of signal.
It is corresponding to be with focus:
Step 103, it is solved with wave field using the adjoint focus backpropagation of structure, and carries out envelope inverting, to right
Initial velocity model is reconstructed, and obtains the rate pattern after being reconstructed under low-frequency data deletion condition;
Ghost is indicated using the Acoustic Wave-equation in isotropic medium, wherein u (x, t) when solving with wave field
:
In formula:
Wherein b is the inverse of density, and κ is bulk modulus, u (x, t)=[vx(x,t),vy(x,t),vz(x,t),p(x,t)
]T, s (x, t)=[0,0,0, Adjsource]T。Indicate the derivative to space,It is the derivative to the directions x,It is to the side y
To derivative,Indicate the derivative to the directions z.vx(x,t),vy(x,t),vz(x, t) is x, y, the particle vibrations in the directions z respectively
For speed with wave field, p (x, t) is with stress wave field.For with the time-derivative (abbreviated form) of wave field.
The inverting flow process that envelope inverting of the present invention uses is consistent with traditional wave equation waveform tomography flow, and envelope inverting can
Using fore condition LBFGS optimization algorithms, the fore condition PLBFGS algorithms of quasi- newton can also be used to improve the rate of convergence of inverting,
To weaken spherical diffusion effect, balance model deep layer and shallow-layer speed renewal amount, the precondition operator that the present embodiment uses can be for:
Z indicates depth in formula.
Step 104, the precision for verifying envelope inverting, in the case where precision is met the requirements, the rate pattern after reconstruct
Traditional wave equation chromatography is carried out as initial velocity model.
In one example, simulation earthquake can be obtained to the rate pattern forward simulation after reconstruct on the basis of step 103
Data, and its envelope is solved, carry out adding window normalization zero with the envelope of the envelope of observation seismic data and the simulation seismic data
Postpone cross-correlation, and is compared with the adding window normalization zero-lag auto-correlation of observation seismic data envelope, it is anti-to verify envelope
The precision drilled, with the quality for the model quality that quality monitoring is reconstructed through envelope inverting.Result of the comparison is closer, and it is anti-to represent envelope
The precision drilled is higher.
In another example, the rate pattern after reconstruct can be used as initial velocity model and carry out routine fluctuations equation
After chromatography, if tomographic results do not have noise, then it represents that the precision of envelope inverting is met the requirements, and is unsatisfactory for wanting if there is noise
It asks.Any means known can be used to carry out the judgement of noise in those skilled in the art.
In one example, when the precision of envelope inverting is unsatisfactory for requiring, modification initial velocity model re-executes step
Rapid 101, or modification inverted parameters, re-execute step 103;When the precision of envelope inverting is met the requirements, envelope inverting knot
Fruit carries out traditional wave equation chromatography as initial velocity model, completes final velocity modeling.It is chromatographed in traditional wave equation
When fore condition PLBFGS optimization algorithms equally can be used.
Due to the limitation of acquisition technique, seismic prospecting is difficult that 5Hz low-frequency information components below are recorded, and is lacked using low frequency
When the seismic data of mistake carries out traditional time-domain wave equation waveform tomography, inversion result can be absorbed in local extremum, can not obtain
To accurate inversion speed model.In addition actual seismic data shot point geophone station energy is unbalanced serious, and the earthquake envelope time
When delay is more than envelope half period, it can also become unstable using the envelope inverting of envelope residual error, land seismic can not be adapted to
The macro-scale velocity modeling of data.The present invention provides a kind of envelope inverting based on normalization nonzero-lag phase-fitting thus
Method, macromodel structure when this method is lacked for seismic data low frequency weaken traditional wave equation chromatography to low frequency height
The defect relied on is spent, and the time delay that the introducing for normalizing nonzero-lag phase-fitting makes the technology allow envelope is more than packet
The half period of network signal, making the present invention, it is more steady effective.The present invention can be the land data wave equation chromatography of tradition
Technology provides one and accurately includes the initial velocity model of lower wave number information, and then it is complete so that traditional wave equation chromatography is converged to
Office's minimum solution.
Compared with existing disclosed technology, the beneficial effects of the invention are as follows:Use the phase of adding window normalization nonzero-lag
Position fitting mutually related objects functional, simulation earthquake envelope can be handled be more than with the delay for observing data envelopment half envelope period this
Kind situation, improves the robustness of envelope inverting, and the introducing of normalized crosscorrelation can effectively suppress real data focus and inspection
Inverting noise, can be through the invention conventional wave under real data low frequency deletion condition caused by wave device envelope energy is unbalanced
Dynamic equation tomographic inversion structure includes the initial background rate pattern of accurate lower wave number information, is the high precision velocity of real data
Modeling lays the foundation.
Using example
The scheme and its effect of the embodiment of the present invention for ease of understanding provides a concrete application below in conjunction with Fig. 2-12 and shows
Example.It will be understood by those skilled in the art that the example is of the invention only for the purposes of understanding, any detail is not intended to appoint
Where formula limitation is of the invention.
This is using example by taking low frequency lacks Ricker wavelet as an example, and Fig. 2 shows low frequency missing Ricker wavelet and its envelopes
Schematic diagram;Fig. 3 shows that low frequency lacks wavelet forward modeling big gun collection;Fig. 4 shows big gun collection envelope natural logrithm;Fig. 5 is shown really
Marmousi models;Fig. 6 shows normal gradient initial model;Fig. 7 shows traditional L2 norms wave equation tomographic inversion result;
Fig. 8 shows that adding window nonzero-lag normalizes phase-fitting envelope inversion result;Fig. 9 shows envelope inverting+tradition L2 norms
Wave equation tomographic inversion result;Figure 10 shows observation seismic data envelope and observes the normalization of seismic data envelope from phase
It closes;Figure 11 shows the normalized crosscorrelation of observation seismic data envelope and initial model analogue data envelope;Figure 12 is shown
Observe the normalized crosscorrelation of seismic data envelope and envelope inverse model analogue data envelope.
Specifically, Fig. 7 be using Fig. 6 normal gradient initial model carry out routine fluctuations equation tomographic inversion as a result,
It can be seen from figure 7 that when low-frequency data lacks, there is noise in inversion result, influence the precision of inverting.And this application is shown
Example uses the normal gradient initial model of Fig. 6 to carry out adding window nonzero-lag and normalizes phase-fitting envelope inverting, obtains underground first
Medium macromodel (Fig. 8) uses model shown in Fig. 8 as initial model and then carries out traditional wave equation chromatography instead later
Drill, obtain Fig. 9, from Fig. 9 and Fig. 7 Comparative results can be seen that the noise of Fig. 9 substantially without, and with Fig. 5 (true model)
That matches is more preferable.
The disclosure can be system, method and/or computer program product.Computer program product may include computer
Readable storage medium storing program for executing, containing for making processor realize the computer-readable program instructions of various aspects of the disclosure.
Computer readable storage medium can be can keep and store the instruction used by instruction execution equipment tangible
Equipment.Computer readable storage medium for example can be-- but be not limited to-- storage device electric, magnetic storage apparatus, optical storage
Equipment, electromagnetism storage device, semiconductor memory apparatus or above-mentioned any appropriate combination.Computer readable storage medium
More specific example (non exhaustive list) includes:Portable computer diskette, random access memory (RAM), read-only is deposited hard disk
It is reservoir (ROM), erasable programmable read only memory (EPROM or flash memory), static RAM (SRAM), portable
Compact disk read-only memory (CD-ROM), digital versatile disc (DVD), memory stick, floppy disk, mechanical coding equipment, for example thereon
It is stored with punch card or groove internal projection structure and the above-mentioned any appropriate combination of instruction.Calculating used herein above
Machine readable storage medium storing program for executing is not interpreted that instantaneous signal itself, the electromagnetic wave of such as radio wave or other Free propagations lead to
It crosses the electromagnetic wave (for example, the light pulse for passing through fiber optic cables) of waveguide or the propagation of other transmission mediums or is transmitted by electric wire
Electric signal.
Computer-readable program instructions as described herein can be downloaded to from computer readable storage medium it is each calculate/
Processing equipment, or outer computer or outer is downloaded to by network, such as internet, LAN, wide area network and/or wireless network
Portion's storage device.Network may include copper transmission cable, optical fiber transmission, wireless transmission, router, fire wall, interchanger, gateway
Computer and/or Edge Server.Adapter or network interface in each calculating/processing equipment are received from network to be counted
Calculation machine readable program instructions, and the computer-readable program instructions are forwarded, for the meter being stored in each calculating/processing equipment
In calculation machine readable storage medium storing program for executing.
For execute the disclosure operation computer program instructions can be assembly instruction, instruction set architecture (ISA) instruction,
Machine instruction, machine-dependent instructions, microcode, firmware instructions, condition setup data or with one or more programming languages
Arbitrarily combine the source code or object code write, the programming language include the programming language-of object-oriented such as
Smalltalk, C++ etc., and conventional procedural programming languages-such as " C " language or similar programming language.Computer
Readable program instructions can be executed fully, partly execute on the user computer, is only as one on the user computer
Vertical software package executes, part executes or on the remote computer completely in remote computer on the user computer for part
Or it is executed on server.In situations involving remote computers, remote computer can pass through network-packet of any kind
It includes LAN (LAN) or wide area network (WAN)-is connected to subscriber computer, or, it may be connected to outer computer (such as profit
It is connected by internet with ISP).In some embodiments, by using computer-readable program instructions
Status information carry out personalized customization electronic circuit, such as programmable logic circuit, field programmable gate array (FPGA) or can
Programmed logic array (PLA) (PLA), the electronic circuit can execute computer-readable program instructions, to realize each side of the disclosure
Face.
Referring herein to according to the flow chart of the method, apparatus (system) of the embodiment of the present disclosure and computer program product and/
Or block diagram describes various aspects of the disclosure.It should be appreciated that flowchart and or block diagram each box and flow chart and/
Or in block diagram each box combination, can be realized by computer-readable program instructions.
These computer-readable program instructions can be supplied to all-purpose computer, special purpose computer or other programmable datas
The processor of processing unit, to produce a kind of machine so that these instructions are passing through computer or other programmable datas
When the processor of processing unit executes, work(specified in one or more of implementation flow chart and/or block diagram box is produced
The device of energy/action.These computer-readable program instructions can also be stored in a computer-readable storage medium, these refer to
It enables so that computer, programmable data processing unit and/or other equipment work in a specific way, to be stored with instruction
Computer-readable medium includes then a manufacture comprising in one or more of implementation flow chart and/or block diagram box
The instruction of the various aspects of defined function action.
Computer-readable program instructions can also be loaded into computer, other programmable data processing units or other
In equipment so that series of operation steps are executed on computer, other programmable data processing units or miscellaneous equipment, with production
Raw computer implemented process, so that executed on computer, other programmable data processing units or miscellaneous equipment
Instruct function action specified in one or more of implementation flow chart and/or block diagram box.
Flow chart and block diagram in attached drawing show the system, method and computer journey of multiple embodiments according to the disclosure
The architecture, function and operation in the cards of sequence product.In this regard, each box in flowchart or block diagram can generation
One module of table, program segment or a part for instruction, the module, program segment or a part for instruction include one or more use
The executable instruction of the logic function as defined in realization.In some implementations as replacements, the function of being marked in box
It can occur in a different order than that indicated in the drawings.For example, two continuous boxes can essentially be held substantially in parallel
Row, they can also be executed in the opposite order sometimes, this is depended on the functions involved.It is also noted that block diagram and/or
The combination of each box in flow chart and the box in block diagram and or flow chart can use function or dynamic as defined in executing
The dedicated hardware based system made is realized, or can be realized using a combination of dedicated hardware and computer instructions.
The presently disclosed embodiments is described above, above description is exemplary, and non-exclusive, and
It is not limited to disclosed each embodiment.Without departing from the scope and spirit of illustrated each embodiment, for this skill
Many modifications and changes will be apparent from for the those of ordinary skill in art field.The selection of term used herein, purport
In the principle, practical application or technological improvement to the technology in market for best explaining each embodiment, or this technology is made to lead
Other those of ordinary skill in domain can understand each embodiment disclosed herein.
Claims (9)
1. a kind of background velocity model reconstruction method under low frequency seismic data missing, this approach includes the following steps:
1) envelope of the simulation seismic data according to the envelope of observation seismic data and based on initial velocity model, structure adding window are returned
One changes the phase-fitting mutually related objects functional of nonzero-lag;
2) it is based on the mutually related objects functional and builds corresponding inverse time propagation with focus;
3) it is solved with wave field using the adjoint focus backpropagation, and carries out envelope inverting, to the initial velocity
Model is reconstructed, and obtains the rate pattern after being reconstructed under low-frequency data deletion condition;
4) precision for verifying envelope inverting, in the case where precision is met the requirements, using the rate pattern after reconstruct as initial speed
It spends model and carries out traditional wave equation chromatography;
Wherein, using the Acoustic Wave-equation in following isotropic medium when solving the adjoint wave field:
In formula:
Wherein b is the inverse of density, and κ is bulk modulus, u (x, t)=[vx(x,t),vy(x,t),vz(x,t),p(x,t)]T, s
(x, t)=[0,0,0, Adjsource]T, the adjoint wave field of u (x, t) expressions,Indicate the derivative to space,It is to the directions x
Derivative,It is the derivative to the directions y,Indicate the derivative to the directions z, vx(x,t),vy(x,t),vz(x, t) is x, y, z respectively
For the particle vibration velocity in direction with wave field, p (x, t) is adjoint stress wave field,To be led with the time of wave field
Number.
2. the background velocity model reconstruction method under low frequency seismic data missing according to claim 1, wherein adding window is returned
One changes the phase-fitting mutually related objects functional χ of nonzero-lagXcorNenlogIt is expressed as:
Wherein, XcorNen1og (τ, h) is observation seismic data envelope natural logrithm and simulation seismic data envelope natural logrithm
Normalized Cross Correlation Function, P (τ) is weighting time histories sample, and τ is time shift amount, and u (t, h) is observation seismic data, and y (t, h) is
Simulate seismic data, uH(t, h) is the Hilbert transform for observing seismic data, yH(t, h) is the Martin Hilb for simulating seismic data
Spy's transformation, eobs(t, h) is the envelope for observing seismic data, esyn(t, h) be simulate seismic data envelope, c be in order to avoid
Natural logrithm antilog is a constant of zero introducing, and h indicates the serial number of wave detector.
3. the background velocity model reconstruction method under low frequency seismic data missing according to claim 2, wherein with shake
Source Adjsource is expressed as:
WhereinFor observe seismic data envelope and
The cross-correlation of the envelope of seismic data is simulated, Hilbert indicates the Hilbert transform of signal.
4. the background velocity model reconstruction method under low frequency seismic data missing according to claim 1, wherein envelope is anti-
Drill the fore condition PLBFGS algorithms using fore condition LBFGS optimization algorithms or quasi- newton.
5. the background velocity model reconstruction method under low frequency seismic data missing according to claim 4, wherein use
Precondition operator is:
Z indicates depth in formula.
6. the background velocity model reconstruction method under low frequency seismic data missing according to claim 1, wherein verification packet
The precision of network inverting includes:
Forward simulation is carried out to the rate pattern after reconstruct, the simulation seismic data of the rate pattern after obtaining based on reconstruct is simultaneously asked
Its envelope is solved, it is mutual to carry out adding window normalization zero-lag with the envelope of the envelope of observation seismic data and the simulation seismic data
It closes, and is compared with the adding window normalization zero-lag auto-correlation of observation seismic data envelope, to verify the precision of envelope inverting.
7. the background velocity model reconstruction method under low frequency seismic data missing according to claim 1, wherein verification packet
The precision of network inverting includes:
The rate pattern after reconstruct is used as initial velocity model and carries out routine fluctuations equation chromatography, if tomographic results are not made an uproar
Sound, then it represents that the precision of envelope inverting is met the requirements, and is unsatisfactory for requiring if there is noise.
8. the background velocity model reconstruction method under low frequency seismic data missing according to claim 1, further includes:
In the case where the precision of envelope inverting is unsatisfactory for requirement, initial velocity model is changed, step 1) is re-executed.
9. the background velocity model reconstruction method under low frequency seismic data missing according to claim 1, further includes:It is wrapping
In the case that the precision of network inverting is unsatisfactory for requirement, inverted parameters are changed, step 3) is re-executed.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510655283.8A CN106569262B (en) | 2015-10-12 | 2015-10-12 | Background velocity model reconstruction method under low frequency seismic data missing |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510655283.8A CN106569262B (en) | 2015-10-12 | 2015-10-12 | Background velocity model reconstruction method under low frequency seismic data missing |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106569262A CN106569262A (en) | 2017-04-19 |
CN106569262B true CN106569262B (en) | 2018-10-02 |
Family
ID=58506307
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510655283.8A Active CN106569262B (en) | 2015-10-12 | 2015-10-12 | Background velocity model reconstruction method under low frequency seismic data missing |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106569262B (en) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109212602B (en) * | 2018-09-05 | 2019-11-08 | 湖南科技大学 | A kind of reflection coefficient inversion method improving seismic data resolution |
CN110007340B (en) * | 2019-02-01 | 2020-09-25 | 西安理工大学 | Salt dome velocity density estimation method based on angle domain direct envelope inversion |
CN110244351A (en) * | 2019-04-22 | 2019-09-17 | 西安石油大学 | A kind of Uniform Construction inversion method of different constraint Geophysical Inverse Problems |
CN112241024B (en) * | 2019-07-18 | 2024-04-09 | 中国石油化工股份有限公司 | Method for improving signal-to-noise ratio of seismic data, computer storage medium and system |
CN111007565B (en) * | 2019-12-24 | 2020-12-11 | 清华大学 | Three-dimensional frequency domain full-acoustic wave imaging method and device |
CN114240788B (en) * | 2021-12-21 | 2023-09-08 | 西南石油大学 | Complex scene-oriented robustness and adaptive background restoration method |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101630018B (en) * | 2008-07-16 | 2011-12-07 | 中国石油天然气股份有限公司 | Seismic exploration data processing method for controlling full acoustic wave equation inversion process |
CN103163554A (en) * | 2013-02-04 | 2013-06-19 | 西安交通大学 | Self-adapting wave form retrieval method through utilization of zero offset vertical seismic profile (VSP) data to estimate speed and Q value |
WO2013093468A2 (en) * | 2011-12-20 | 2013-06-27 | Shah Nikhil Koolesh | Full waveform inversion quality control method |
CN103513277A (en) * | 2013-09-27 | 2014-01-15 | 中国石油天然气股份有限公司 | Seismic stratum fracture crack density inversion method and system |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9857488B2 (en) * | 2012-11-20 | 2018-01-02 | International Business Machines Corporation | Efficient wavefield compression in seismic imaging |
-
2015
- 2015-10-12 CN CN201510655283.8A patent/CN106569262B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101630018B (en) * | 2008-07-16 | 2011-12-07 | 中国石油天然气股份有限公司 | Seismic exploration data processing method for controlling full acoustic wave equation inversion process |
WO2013093468A2 (en) * | 2011-12-20 | 2013-06-27 | Shah Nikhil Koolesh | Full waveform inversion quality control method |
CN103163554A (en) * | 2013-02-04 | 2013-06-19 | 西安交通大学 | Self-adapting wave form retrieval method through utilization of zero offset vertical seismic profile (VSP) data to estimate speed and Q value |
CN103513277A (en) * | 2013-09-27 | 2014-01-15 | 中国石油天然气股份有限公司 | Seismic stratum fracture crack density inversion method and system |
Non-Patent Citations (3)
Title |
---|
Velocity model building from seismic reflection data by full-waveform inversion;Romain Brossier 等;《Geophysical Prospecting》;20141111;第63卷(第2期);第354-367页 * |
低频缺失情况下的弹性波包络反演;黄超 等;《应用地球物理》;20150930;第12卷(第3期);第362-377页 * |
声波全波形反演目标函数性态;董良国 等;《地球物理学报》;20131231;第56卷(第10期);第3445-3460页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106569262A (en) | 2017-04-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106569262B (en) | Background velocity model reconstruction method under low frequency seismic data missing | |
CN106526674B (en) | Three-dimensional full waveform inversion energy weighting gradient preprocessing method | |
US11609352B2 (en) | Machine learning-augmented geophysical inversion | |
CN111239802B (en) | Deep learning speed modeling method based on seismic reflection waveform and velocity spectrum | |
Chen et al. | Elastic least-squares reverse time migration via linearized elastic full-waveform inversion with pseudo-Hessian preconditioning | |
CA2690373C (en) | Method for velocity analysis using waveform inversion in laplace domain for geophysical imaging | |
CN102892972B (en) | Artefact in the iterative inversion of geophysical data is reduced | |
CN105467444B (en) | A kind of elastic wave full waveform inversion method and device | |
CA2947410A1 (en) | Fast viscoacoustic and viscoelastic full-wavefield inversion | |
CN110058302A (en) | A kind of full waveform inversion method based on pre-conditional conjugate gradient accelerating algorithm | |
CN104965223B (en) | Method and device for inverting full waveform of viscoacoustic wave | |
CN105093278B (en) | Full waveform inversion gradient operator extracting method based on the main energy-optimised algorithm of excitation | |
CN104237937B (en) | Pre-stack seismic inversion method and system thereof | |
CN107894618B (en) | A kind of full waveform inversion gradient preprocess method based on model smoothing algorithm | |
CN104965222B (en) | Three-dimensional longitudinal wave impedance full-waveform inversion method and device | |
CN113341455B (en) | Viscous anisotropic medium seismic wave numerical simulation method, device and equipment | |
KR101820850B1 (en) | Seismic imaging apparatus and method using iterative direct waveform inversion | |
US9921324B2 (en) | Systems and methods employing upward beam propagation for target-oriented seismic imaging | |
Ivanov et al. | Traveltime parameters in tilted orthorhombic medium | |
Toxopeus et al. | Simulating migrated and inverted seismic data by filtering a geologic model | |
CN108680968B (en) | Evaluation method and device for seismic exploration data acquisition observation system in complex structural area | |
Ha et al. | 3D Laplace-domain waveform inversion using a low-frequency time-domain modeling algorithm | |
CN107179547A (en) | A kind of question of seismic wave impedance inversion low frequency model method for building up | |
CN109738944B (en) | Wide-angle reflection-based seismic acquisition parameter determination method and device | |
Jaimes-Osorio et al. | Amplitude variation with offset inversion using acoustic-elastic local solver |
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 |