CN106094020A - A kind of seismic inversion method and device - Google Patents
A kind of seismic inversion method and device Download PDFInfo
- Publication number
- CN106094020A CN106094020A CN201610375072.3A CN201610375072A CN106094020A CN 106094020 A CN106094020 A CN 106094020A CN 201610375072 A CN201610375072 A CN 201610375072A CN 106094020 A CN106094020 A CN 106094020A
- Authority
- CN
- China
- Prior art keywords
- subregion
- analysis window
- dispersion curve
- work area
- frequency
- 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 45
- 239000006185 dispersion Substances 0.000 claims abstract description 95
- 238000001228 spectrum Methods 0.000 claims abstract description 42
- 238000011161 development Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 20
- 239000000843 powder Substances 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000010363 phase shift Effects 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 229910052704 radon Inorganic materials 0.000 description 2
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 2
- 238000005096 rolling process Methods 0.000 description 2
- 238000011282 treatment Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 239000011295 pitch Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The embodiment of the present application discloses a kind of seismic inversion method and device.Described seismic inversion method includes: divides described work area according to the geological data in work area, obtains the subregion in described work area；Obtain the face ripple type of each subregion and single big gun geological data；Face based on described each subregion ripple type and single big gun geological data, determine each analysis window of this subregion；Obtain single big gun geological data of described each subregion frequency rate spectrum in each analysis window of this subregion；Compose based on the frequency rate in the described each analysis window of each subregion, generate the dispersion curve of this analysis window；Dispersion curve based on the described each analysis window of each subregion, carries out seismic inversion to described work area.The seismic inversion method of the embodiment of the present application and device, can improve the precision of seismic inversion.
Description
Technical field
The application belongs to technical field of geophysical exploration, particularly to a kind of seismic inversion method and device.
Background technology
Underground strata spatial structure, mainly by the geological data of observation, is become by seismic inversion with physical property
Picture.At present, the face ripple of utilization carries out the research of seismic inversion and has been carried out for a long time, and achieves preferable effect, its
In, described ripple can include rayleigh surface wave (Rayleigh Wave) and pressgang face ripple (Love Wave) etc..
In prior art, the methods such as phase shift method, f-k method and τ-p method that are based primarily upon extract surface wave dispersion from geological data
Curve, carries out inverting, wherein, described dispersion curve based on the surface wave frequency dispersion curve extracted to spatial structure formation and physical property
Refer generally to the phase velocity of seismic wave with frequency or the change curve of wavelength.
The process extracting surface wave frequency dispersion curve based on methods such as phase shift method, f-k method and τ-p methods from geological data is general such as
Under:
From the geological data in work area, obtain single big gun geological data, and use an analysis window to described single big gun earthquake number
According to carrying out Analysis of dispersion characteristics, obtain the frequency-normal-moveout spectrum of this list big gun geological data, from described frequency-normal-moveout spectrum, pick up face
Wave energy focuses on preferable point, is generated the dispersion curve of this list big gun geological data, the dispersion curve that will generate by the point of pickup
Surface wave frequency dispersion curve as described work area.
During realizing the application, inventor finds that in prior art, at least there are the following problems:
During the face ripple of utilization carries out seismic inversion, the precision of surface wave frequency dispersion curve typically precision to seismic inversion
There is large effect.Above-mentioned method of the prior art, is not the work area grown very much for face ripple, uses an analysis window
Often cannot obtain frequency-normal-moveout spectrum that precision is higher, the precision causing the dispersion curve of generation is poor, so that earthquake is anti-
The precision drilled is relatively low.
Summary of the invention
The purpose of the embodiment of the present application is to provide a kind of seismic inversion method and device, to improve the precision of seismic inversion.
For solving above-mentioned technical problem, the embodiment of the present application provides a kind of seismic inversion method and device to be realized in
:
A kind of seismic inversion method, including:
According to the geological data in work area, described work area is divided, obtain the subregion in described work area；
Obtain the face ripple type of each subregion and single big gun geological data；
Face based on described each subregion ripple type and single big gun geological data, determine each analysis window of this subregion；
Obtain single big gun geological data of described each subregion frequency-normal-moveout spectrum in each analysis window of this subregion；
Based on the frequency-normal-moveout spectrum in the described each analysis window of each subregion, generate the dispersion curve of this analysis window；
Dispersion curve based on the described each analysis window of each subregion, carries out seismic inversion to described work area.
A kind of seismic inversion device, including:
Division unit, divides described work area according to the geological data in work area, obtains the subregion in described work area；
First acquiring unit, for obtaining the face ripple type of each subregion and single big gun geological data；
Determine unit, for face based on described each subregion ripple type and single big gun geological data, determine each of this subregion
Analysis window；
Second acquisition unit, for obtaining single big gun geological data of described each subregion in each analysis window of this subregion
Frequency-normal-moveout spectrum；
Signal generating unit, for based on the frequency-normal-moveout spectrum in the described each analysis window of each subregion, generates this analysis window
Dispersion curve；
Inverting unit, for dispersion curve based on the described each analysis window of each subregion, carries out earthquake to described work area
Inverting.
The technical scheme provided from above the embodiment of the present application, the embodiment of the present application can be according to the earthquake number in work area
Divide according to described work area, obtain the subregion in described work area；Can face based on each subregion ripple type and single big gun earthquake
Data, determine each analysis window of this subregion；The single big gun geological data that can obtain described each subregion is respectively analyzed at this subregion
Time window in frequency-normal-moveout spectrum；This analysis can be generated based on the frequency-normal-moveout spectrum in the described each analysis window of each subregion
Time window dispersion curve；And based on the dispersion curve of the described each analysis window of each subregion, described work area can be carried out earthquake
Inverting.Compared with prior art, work area can be divided by the embodiment of the present application；Can use multiple for each subregion
Analysis window carries out Analysis of dispersion characteristics to single big gun geological data of this subregion, to obtain the frequency in each analysis window of this subregion
Rate-normal-moveout spectrum and dispersion curve.Therefore, the embodiment of the present application when work area is carried out seismic inversion, the analysis window of employing
Quantity is more, can obtain frequency-normal-moveout spectrum that precision is higher, such that it is able to improve the precision of work area dispersion curve, Jin Erke
To improve the precision of work area seismic inversion.
Accompanying drawing explanation
In order to be illustrated more clearly that the embodiment of the present application or technical scheme of the prior art, below will be to embodiment or existing
In having technology to describe, the required accompanying drawing used is briefly described, it should be apparent that, the accompanying drawing in describing below is only this
Some embodiments described in application, for those of ordinary skill in the art, in the premise not paying creative work
Under, it is also possible to other accompanying drawing is obtained according to these accompanying drawings.
Fig. 1 is the flow chart of a kind of seismic inversion method of the embodiment of the present application；
Fig. 2 is the face ripple type schematic diagram of a subregion in the embodiment of the present application work area；
Fig. 3 is the analysis window schematic diagram of another subregion in the embodiment of the present application work area；
Fig. 4 a is an analysis window schematic diagram of the embodiment of the present application face ripple type 201；
Fig. 4 b is the dispersion curve schematic diagram of analysis window shown in Fig. 4 a；
Fig. 5 a is another analysis window schematic diagram of the embodiment of the present application face ripple type 201；
Fig. 5 b is the dispersion curve schematic diagram of analysis window shown in Fig. 5 a；
Fig. 6 a is an analysis window schematic diagram of the embodiment of the present application face ripple type 202；
Fig. 6 b is the dispersion curve schematic diagram of analysis window shown in Fig. 6 a；
Fig. 7 a is another analysis window schematic diagram of the embodiment of the present application face ripple type 202；
Fig. 7 b is the dispersion curve schematic diagram of analysis window shown in Fig. 7 a；
Fig. 8 is the illustrative view of functional configuration of the embodiment of the present application a kind of seismic inversion device.
Detailed description of the invention
For the technical scheme making those skilled in the art be more fully understood that in the application, real below in conjunction with the application
Execute the accompanying drawing in example, the technical scheme in the embodiment of the present application is clearly and completely described, it is clear that described enforcement
Example is only some embodiments of the present application rather than whole embodiments.Based on the embodiment in the application, this area is common
The every other embodiment that technical staff is obtained under not making creative work premise, all should belong to the application protection
Scope.
An embodiment of the application seismic inversion method is described below, as it is shown in figure 1, this embodiment may include that
S101: according to the geological data in work area, described work area is divided, obtain the subregion in described work area.
Described work area refers generally to carry out the region of geological prospecting.Described geological data can include being measured by observation system
Data, it is also possible to include the data measured by logging method.The number of partitions in described work area can be 1, it is also possible to is
Multiple.
Specifically, according to the attribute data of face ripple in the geological data of work area, described work area can be divided.Described genus
Property data can include apparent velocity data, video frequency data and development degree data etc..Wherein, face ripple can generally be used
Energy represent the development degree of this face ripple.So, described development degree data can include energy datum.
In one embodiment, apparent velocity data and the video frequency number of face ripple can be obtained from the geological data in work area
According to；Can according to described apparent velocity data and video frequency data the regularity of distribution in described work area, described work area is carried out point
District processes.Such as, certain ripple apparent velocity in A district in work area is 500～700 meter per seconds, and video frequency is 10～12 hertz；
This face ripple apparent velocity in B district in work area is 850～1100 meter per seconds, and video frequency is 9～11 hertz.It is possible to according to being somebody's turn to do
The apparent velocity of face ripple and the video frequency difference rule in described work area, carries out multidomain treat-ment to described work area, thus by described
Work area is divided into A district and two, B district subregion.
In another embodiment, the apparent velocity data of face ripple, video frequency can be obtained from the geological data in work area
Data and development degree data；Can be according to described apparent velocity data, video frequency data and development degree data in institute
State the regularity of distribution in work area, described work area is carried out multidomain treat-ment.Specifically, can be according to described apparent velocity data and video
The rate data regularity of distribution in described work area, carries out Preliminary division to described work area, obtains the subregion after Preliminary division；Then
According to the described development degree data regularity of distribution in described work area, the subregion after Preliminary division can be carried out further
Divide, thus described work area is divided into one or more subregion.
S102: obtain the face ripple type of each subregion and single big gun geological data.
The geological data that described single big gun geological data produces when referring generally to excite a big gun.
Specifically, for each subregion, the attribute data of this subregion inner face ripple can be obtained from described geological data, and
The face ripple type of this subregion and single big gun geological data based on described attribute data, can be obtained.
The quantity of each subregion face ripple type can be 1, it is also possible to for multiple.So, for each subregion, Ke Yicong
Described geological data obtains apparent velocity data and the video frequency data of this subregion inner face ripple, it is possible to based on described apparent velocity number
According to video frequency data, the face ripple in this subregion is carried out Type division, obtains the face ripple type of this subregion.Such as, C in work area
The quantity of face, district ripple type can be 2, specifically can include face ripple type one dough-making powder ripple type two.Wherein, face ripple type one
Apparent velocity be 500～700 meter per seconds, video frequency is 10～12 hertz；The apparent velocity of face ripple type two is 900～1000 meters/
Second, video frequency is 7～9 hertz.
Fig. 2 is the face ripple type schematic diagram of a certain subregion in work area.In fig. 2, the quantity of this subregion face ripple type is 2
Individual, specifically include face ripple type 201 dough-making powder ripple type 202.Wherein, the apparent velocity of face ripple type is 512 meter per seconds, face ripple type
The apparent velocity of 202 is 1050 meter per seconds.
The quantity of each subregion list big gun geological data is generally 1.Specifically, for each subregion, can be from described
Shake data obtain the apparent velocity codomain of this subregion inner face ripple, it is possible to from described geological data, obtain apparent velocity regard described
The face ripple that the median of speed codomain is neighbouring and energy is stronger, the single big gun geological data corresponding to the ripple of face that can obtain
Single big gun geological data as this subregion.Such as, in work area, the apparent velocity codomain of A district face ripple is 500～700 meter per seconds, described regards
The median of speed codomain is 600 meter per seconds.It is possible to from work geological data obtain apparent velocity value near 600 meter per seconds,
And the face ripple that energy is stronger, it is possible to single big gun geological data corresponding described in the face ripple that will obtain is as single big gun earthquake in A district
Data.
It should be noted that the face ripple type of each subregion can be identical in work area, it is also possible to different, it is also possible to part phase
With.Such as, in work area, the face ripple type in A district can include face ripple type one dough-making powder ripple type two, the face ripple type in B district in work area
Face ripple type two can be included.
Also, it should be noted the same face ripple type can be identical at the attribute data of different subregions, it is also possible to different.Example
As, the face ripple type one apparent velocity in A district can be 500～700 meter per seconds, and video frequency can be 10～12 hertz；This face ripple class
The type one apparent velocity in B district can be 1400～1500 meter per seconds, and video frequency can be 8～10 hertz.
S103: face based on described each subregion ripple type and single big gun geological data, determine the analysis window of this subregion.
The time window that described analysis window is used when being generally face wave analysis.
Specifically, for described each subregion, can obtain corresponding to the quantity of this subregion face ripple type, face ripple type
Apparent velocity and the geophone offset codomain of single big gun geological data, can be corresponding to quantity based on described ripple type, face ripple type
Apparent velocity and the geophone offset codomain of single big gun geological data, determine the analysis window of this subregion.Such as, the face in A district in work area
The quantity of ripple type can be 2, specifically can include face ripple type three dough-making powder ripple type four.Regarding corresponding to face ripple type three
Speed is 500～700 meter per seconds, corresponding to face ripple type four apparent velocity be 900～1000 meter per seconds.A district list big gun earthquake number
According to geophone offset codomain be L01～L02.It is possible to determine that the quantity of A district analysis window is 6, specifically it is respectively With
In one embodiment, for described each subregion, when single big gun geological data of this subregion is split shooting ground
During shake data, can all set up analysis window in the both sides of this list big gun geological data correspondence shot point.Wherein, described split shooting ground
Shake data are generally positioned at the geological data that the shot point in the middle of arrangement produces.
In order to save workload, " roll into and roll out " technology can be used when multi-fold.Such as split shooting two
Duan Ge 32 road receives, and during beginning, when the first bombard is put at survey line top, shot point trumpet end 0 road receives, and shot point large size end 32 road connects
Receive；Arrangement subsequently is pushed ahead 4 track pitches and is put the 2nd big gun, and now shot point trumpet end increases totally 4 roads, 4 roads and receives, and shot point large size end is still
Right 32 roads receive；When putting 3 big gun, shot point trumpet end increases totally 8 roads, 4 roads again and receives, and shot point large size end remains as 32 roads and receives；With this
Analogizing, during to 9 big gun, shot point trumpet end and large size end are 32 roads and receive, and this is the process that arrangement " is rolled into "；Otherwise, arrive
Survey line end, shot point trumpet end remains that 32 roads receive, and shot point large size end receives number of channels and successively decreases by 4 roads, and this is arrangement " rolling
Go out " process.So, in another embodiment, for described each subregion, when single big gun geological data of this subregion is for using
When " rolling into and roll out " geological data that mode produces, can set up in the one or both sides of this list big gun geological data correspondence shot point and divide
Window during analysis.
Fig. 3 is the analysis window schematic diagram of a certain subregion in work area.In figure 3,301 is analysis window 1, and 302 for analyzing
Time window 2,303 be analysis window 3.
S104: obtain single big gun geological data of described each subregion frequency-normal-moveout spectrum in each analysis window of this subregion.
Described frequency-normal-moveout spectrum refers generally to the frequency variation relation relative to speed of geological data.
Specifically, for described each subregion, the geophone offset codomain of this subregion each analysis window inner face ripple can be obtained；Can
With geophone offset codomain based on described each analysis window inner face ripple, calculate single big gun geological data of this subregion when described each analysis
Frequency-normal-moveout spectrum in window.
Such as, for described each subregion, the primary earthquake road quantity of this subregion can be obtained；For this subregion each point
Window during analysis, can obtain the focus condition of wave energy in this analysis window based on described primary earthquake road quantity；The most permissible
Increasing seismic channel quantity, after increasing seismic channel quantity, in this analysis window, the focus condition of wave energy can make moderate progress；
Earthquake number of channels amount can be gradually increased, till the focus condition of wave energy no longer improves in this analysis window, now
Geophone offset codomain corresponding to seismic channel quantity can be as the geophone offset codomain of this analysis window.
Further, for described each subregion, the geophone offset codomain of this subregion each analysis window inner face ripple can be obtained
With apparent velocity codomain；For each analysis window of this subregion, can geophone offset codomain based on this analysis window inner face ripple and regarding
Speed codomain, calculates single big gun geological data of this subregion frequency-normal-moveout spectrum in this analysis window.For example, it is possible to pass through
Radon conversion (radon transform) calculates single big gun geological data frequency-normal-moveout spectrum in this analysis window.
S105: based on the frequency-normal-moveout spectrum in the described each analysis window of each subregion, the frequency dispersion generating this analysis window is bent
Line.
Described dispersion curve refers generally to the relation curve between the cycle of seismic wave and velocity of wave.The cycle tool of described seismic wave
Body can include wavelength and the frequency of seismic wave.
Specifically, for each analysis window of described each subregion, can be based on the frequency-normal-moveout spectrum in this analysis window
Obtain the dispersion curve control point of this analysis window, it is possible to generate the frequency of this analysis window based on described dispersion curve control point
Non-dramatic song line.
The described dispersion curve control point obtaining this analysis window based on the frequency-normal-moveout spectrum in this analysis window, permissible
Including: pick up the dispersion curve control point of this analysis window the frequency-normal-moveout spectrum in this analysis window.For example, it is possible to edge
Above frequency-normal-moveout spectrum, ripple focuses on preferable position and picks up the dispersion curve control point of this analysis window.
The described dispersion curve generating this analysis window based on described dispersion curve control point, may include that described frequency
Non-dramatic song line traffic control point carries out line, generates the dispersion curve of this analysis window.
In one embodiment, in the frequency dispersion obtaining this analysis window based on the frequency-normal-moveout spectrum in this analysis window
After curve control point, described method can also include: for each analysis window of described each subregion, obtains this analysis window
The band of position；Dispersion curve control point in the described band of position is processed.Correspondingly, described bent based on described frequency dispersion
Line traffic control point generates the dispersion curve of this analysis window, specifically includes: based on the dispersion curve control point after processing, generate this point
The dispersion curve of window during analysis.
The described band of position can include combination center region and the region, intermediate data road etc. of analysis window.
Described dispersion curve control point in the described band of position is processed, may include that the described band of position
In dispersion curve control point be ranked up, or, reject the abnormity point in dispersion curve control point in the described band of position.
In another embodiment, at the frequency obtaining this analysis window based on the frequency-normal-moveout spectrum in this analysis window
After dissipating curve control point, described method can also include: for each analysis window of described each subregion, when obtaining this analysis
The band of position of window；According to the development degree of this analysis window inner face ripple, by the dispersion curve control point in the described band of position
Put together, obtain the dispersion curve control point set in the described band of position.Correspondingly, described based on described dispersion curve control
System point generates the dispersion curve of this analysis window, specifically includes: gather based on the dispersion curve control point in the described band of position,
Generate the dispersion curve of this analysis window.
Fig. 4 a is an analysis window schematic diagram of face ripple type 201.In fig .4, white portion is analysis window.
Fig. 4 b is the dispersion curve schematic diagram of analysis window shown in Fig. 4 a.
Fig. 5 a is another analysis window schematic diagram of face ripple type 201.In fig 5 a, white portion is analysis window.
Fig. 5 b is the dispersion curve schematic diagram of analysis window shown in Fig. 5 a.
Fig. 6 a is an analysis window schematic diagram of face ripple type 202.In Fig. 6 a, white portion is analysis window.
Fig. 6 b is the dispersion curve schematic diagram of analysis window shown in Fig. 6 a.
Fig. 7 a is another analysis window schematic diagram of face ripple type 202.In figure 7 a, white portion is analysis window.
Fig. 7 b is the dispersion curve schematic diagram of analysis window shown in Fig. 7 a.
S105: dispersion curve based on the described each analysis window of each subregion, carries out seismic inversion to described work area.
Specifically, the dispersion curve of the described each analysis window of each subregion can be carried out exponent number demarcation, based on described respectively
The calibrated dispersion curve of analysis window exponent number, carries out seismic inversion to described work area.
Further, the described dispersion curve to the described each analysis window of each subregion carries out exponent number demarcation, may include that
For described each subregion, according to the position relationship of the dispersion curve of each analysis window of this subregion, time respectively to described each analysis
The dispersion curve of window carries out exponent number demarcation.For example, it is possible to the dispersion curve of window is demarcated as 1 when will be located in lowest frequency position analysis
Rank, will be demarcated as 2 rank with the dispersion curve of this analysis window that lowest frequency position is adjacent and frequency location is higher.And for example, exist
In Fig. 3, the dispersion curve of analysis window 301 can be demarcated as 1 rank, the dispersion curve of analysis window 303 is demarcated as 2 rank.
In one embodiment, before step S105, described method can also include: for described each subregion
Each analysis window, obtains the band of position of this analysis window；Dispersion curve in the described band of position is processed.Accordingly
Ground, described dispersion curve based on the described each analysis window of each subregion, described work area is carried out seismic inversion, specifically includes:
Dispersion curve after process based on the described each analysis window of each subregion, carries out seismic inversion to described work area.
Described dispersion curve in the described band of position is processed, may include that the frequency in the described band of position
Non-dramatic song line is ranked up processing, or, edit the dispersion curve in the described band of position.
Embodiment corresponding to Fig. 1, can divide described work area according to the geological data in work area, obtain described work
The subregion in district；Can face based on each subregion ripple type and single big gun geological data, determine each analysis window of this subregion；Permissible
Obtain single big gun geological data of described each subregion frequency-normal-moveout spectrum in each analysis window of this subregion；Can be based on described
Frequency-normal-moveout spectrum in each analysis window of each subregion, generates the dispersion curve of this analysis window；And can be based on described each
The dispersion curve of each analysis window of subregion, carries out seismic inversion to described work area.Compared with prior art, the embodiment of the present application can
So that work area is divided；Can use multiple analysis window that single big gun geological data of this subregion is carried out for each subregion
Analysis of dispersion characteristics, to obtain the frequency-normal-moveout spectrum in each analysis window of this subregion and dispersion curve.Therefore, corresponding to Fig. 1
Embodiment is when carrying out seismic inversion to work area, and the quantity of the analysis window of employing is more, can obtain the higher frequency of precision-
Normal-moveout spectrum, such that it is able to improve the precision of work area dispersion curve, and then can improve the precision of work area seismic inversion.
The embodiment of the present application also provides for a kind of seismic inversion device, and as shown in Figure 8, this device includes:
Division unit 801, divides described work area according to the geological data in work area, obtains the subregion in described work area；
First acquiring unit 802, for obtaining the face ripple type of each subregion and single big gun geological data；
Determine unit 803, for face based on described each subregion ripple type and single big gun geological data, determine this subregion
Each analysis window；
Second acquisition unit 804, for obtaining single big gun geological data of described each subregion at each analysis window of this subregion
Interior frequency-normal-moveout spectrum；
Signal generating unit 805, for based on the frequency-normal-moveout spectrum in the described each analysis window of each subregion, generates this analysis
Time window dispersion curve；
Inverting unit 806, for dispersion curve based on the described each analysis window of each subregion, carries out ground to described work area
Shake inverting.
It is also known in the art that in addition to realizing controller in pure computer readable program code mode, complete
Full controller can be made with gate, switch, special IC, able to programme by method step carries out programming in logic
The form of logic controller and embedding microcontroller etc. realizes identical function.The most this controller is considered one
Hardware component, and the structure in hardware component can also be considered as to include in it for the device realizing various function.Or
Even, can be considered as not only can being the software module of implementation method but also can be Hardware Subdivision by the device being used for realizing various function
Structure in part.
System, device, module or the unit that above-described embodiment illustrates, specifically can be realized by computer chip or entity,
Or realized by the product with certain function.
For convenience of description, it is divided into various unit to be respectively described with function when describing apparatus above.Certainly, this is being implemented
The function of each unit can be realized in same or multiple softwares and/or hardware during application.
Each embodiment in this specification all uses the mode gone forward one by one to describe, identical similar portion between each embodiment
Dividing and see mutually, what each embodiment stressed is the difference with other embodiments.Real especially for system
For executing example, owing to it is substantially similar to embodiment of the method, so describe is fairly simple, relevant part sees embodiment of the method
Part illustrate.
The application can be used in numerous general or special purpose computing system environments or configuration.Such as: personal computer, clothes
Business device computer, handheld device or portable set, laptop device, multicomputer system, system based on microprocessor, put
Top box, programmable consumer-elcetronics devices, network PC, minicomputer, mainframe computer, include any of the above system or equipment
Distributed computing environment etc..
The application can be described in the general context of computer executable instructions, such as program
Module.Usually, program module includes performing particular task or realizing the routine of particular abstract data type, program, object, group
Part, data structure etc..The application can also be put into practice in a distributed computing environment, in these distributed computing environment, by
The remote processing devices connected by communication network performs task.In a distributed computing environment, program module is permissible
It is positioned in the local and remote computer-readable storage medium of storage device.
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application have many deformation and
Change is without deviating from spirit herein, it is desirable to appended claim includes that these deformation and change are without deviating from the application's
Spirit.
Claims (10)
1. a seismic inversion method, it is characterised in that including:
According to the geological data in work area, described work area is divided, obtain the subregion in described work area；
Obtain the face ripple type of each subregion and single big gun geological data；
Face based on described each subregion ripple type and single big gun geological data, determine each analysis window of this subregion；
Obtain single big gun geological data of described each subregion frequency-normal-moveout spectrum in each analysis window of this subregion；
Based on the frequency-normal-moveout spectrum in the described each analysis window of each subregion, generate the dispersion curve of this analysis window；
Dispersion curve based on the described each analysis window of each subregion, carries out seismic inversion to described work area.
2. the method for claim 1, it is characterised in that described work area is carried out drawing by the described geological data according to work area
Point, specifically include:
According to the attribute data of face ripple in the geological data of work area, dividing described work area, described attribute data includes regarding speed
Degrees of data, video frequency data and development degree data.
3. the method for claim 1, it is characterised in that the face ripple type of each subregion of described acquisition and single big gun earthquake number
According to, specifically include:
For each subregion, from described geological data, obtain the attribute data of this subregion inner face ripple, wherein, described attribute data
Including apparent velocity data, video frequency data and development degree data；
Attribute data based on described ripple, obtains the face ripple type of this subregion and single big gun geological data.
4. the method for claim 1, it is characterised in that described face based on described each subregion ripple type and single big gun ground
Shake data, determine the analysis window of this subregion, specifically include:
For described each subregion, obtain the apparent velocity corresponding to the quantity of this subregion face ripple type, face ripple type and single big gun
The geophone offset codomain of geological data；
Apparent velocity corresponding to quantity based on described ripple type, face ripple type and the geophone offset value of single big gun geological data
Territory, determines the analysis window of this subregion.
5. the method for claim 1, it is characterised in that single big gun geological data of the described each subregion of described acquisition is at this
Frequency-normal-moveout spectrum in each analysis window of subregion, specifically includes:
For described each subregion, the geophone offset codomain of this subregion each analysis window inner face ripple can be obtained；
Geophone offset codomain based on described each analysis window inner face ripple, calculates single big gun geological data of this subregion in described each analysis
Time window in frequency-normal-moveout spectrum.
6. the method for claim 1, it is characterised in that described based on the frequency in the described each analysis window of each subregion
Rate-normal-moveout spectrum, generates the dispersion curve of this analysis window, specifically includes:
For each analysis window of described each subregion, obtain this analysis window based on the frequency-normal-moveout spectrum in this analysis window
Dispersion curve control point；
The dispersion curve of this analysis window is generated based on described dispersion curve control point.
7. method as claimed in claim 6, it is characterised in that described based on the frequency in this analysis window-normal-moveout spectrum acquisition
The dispersion curve control point of this analysis window, specifically includes:
The dispersion curve control point of this analysis window is picked up frequency-normal-moveout spectrum in this analysis window.
8. method as claimed in claim 6, it is characterised in that at each analysis window for described each subregion, based on this
After frequency-normal-moveout spectrum in analysis window obtains the dispersion curve control point of this analysis window, described method also includes:
For each analysis window of described each subregion, obtain the band of position of this analysis window；
Dispersion curve control point in the described band of position is processed；
Correspondingly, the described dispersion curve generating this analysis window based on described dispersion curve control point, specifically include:
Based on the dispersion curve control point after processing, generate the dispersion curve of this analysis window.
9. the method for claim 1, it is characterised in that based on the frequency in the described each analysis window of each subregion-
Normal-moveout spectrum, after generating the dispersion curve of this analysis window, described method also includes:
For each analysis window of described each subregion, obtain the band of position of this analysis window；
Dispersion curve in the described band of position is processed；
Correspondingly, described dispersion curve based on the described each analysis window of each subregion, described work area is carried out seismic inversion, tool
Body includes:
Dispersion curve after process based on the described each analysis window of each subregion, carries out seismic inversion to described work area.
10. a seismic inversion device, it is characterised in that including:
Division unit, divides described work area according to the geological data in work area, obtains the subregion in described work area；
First acquiring unit, for obtaining the face ripple type of each subregion and single big gun geological data；
Determine unit, for face based on described each subregion ripple type and single big gun geological data, determine each analysis of this subregion
Time window；
Second acquisition unit, is used for the single big gun geological data obtaining described each subregion frequency in each analysis window of this subregion
Rate-normal-moveout spectrum；
Signal generating unit, for based on the frequency-normal-moveout spectrum in the described each analysis window of each subregion, generates the frequency of this analysis window
Non-dramatic song line；
Inverting unit, for dispersion curve based on the described each analysis window of each subregion, carries out seismic inversion to described work area.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610375072.3A CN106094020B (en) | 2016-05-31 | 2016-05-31 | A kind of seismic inversion method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610375072.3A CN106094020B (en) | 2016-05-31 | 2016-05-31 | A kind of seismic inversion method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106094020A true CN106094020A (en) | 2016-11-09 |
| CN106094020B CN106094020B (en) | 2018-09-04 |
Family
ID=57229538
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610375072.3A Active CN106094020B (en) | 2016-05-31 | 2016-05-31 | A kind of seismic inversion method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106094020B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2019071515A1 (en) * | 2017-10-12 | 2019-04-18 | 南方科技大学 | Surface wave detection method and terminal device |
| CN111856555A (en) * | 2020-06-19 | 2020-10-30 | 同济大学 | Underground detection method based on surface wave multi-scale window analysis |
| CN112083487A (en) * | 2020-09-16 | 2020-12-15 | 中国科学技术大学 | Method and device for extracting broadband frequency dispersion curve |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110134722A1 (en) * | 2009-12-07 | 2011-06-09 | Mediwound Ltd. | Simultaneous Joint Inversion of Surface Wave and Refraction Data |
| CN103616714A (en) * | 2013-11-19 | 2014-03-05 | 长安大学 | Method for detecting underground cavity boundary through rayleigh surface waves |
| CN103792572A (en) * | 2012-10-29 | 2014-05-14 | 中国石油化工股份有限公司 | Method employing surface wave information in seismic records to investigate surface structure |
| CN104122589A (en) * | 2014-07-30 | 2014-10-29 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Seismic record broadband expanding method |
| CN104345341A (en) * | 2013-08-08 | 2015-02-11 | 中国石油化工股份有限公司 | Region constraint-based frequency band division energy seismic surface wave processing method |
-
2016
- 2016-05-31 CN CN201610375072.3A patent/CN106094020B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110134722A1 (en) * | 2009-12-07 | 2011-06-09 | Mediwound Ltd. | Simultaneous Joint Inversion of Surface Wave and Refraction Data |
| CN103792572A (en) * | 2012-10-29 | 2014-05-14 | 中国石油化工股份有限公司 | Method employing surface wave information in seismic records to investigate surface structure |
| CN104345341A (en) * | 2013-08-08 | 2015-02-11 | 中国石油化工股份有限公司 | Region constraint-based frequency band division energy seismic surface wave processing method |
| CN103616714A (en) * | 2013-11-19 | 2014-03-05 | 长安大学 | Method for detecting underground cavity boundary through rayleigh surface waves |
| CN104122589A (en) * | 2014-07-30 | 2014-10-29 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Seismic record broadband expanding method |
Non-Patent Citations (1)
| Title |
|---|
| JIANG FUHAO ETL.: "MASW Inverse Near Surface Structure Method Study in The Applications", 《SPG/SEG BEIJING 2016 INTERNATIONAL GEOPHYSICAL CONFERENCE》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2019071515A1 (en) * | 2017-10-12 | 2019-04-18 | 南方科技大学 | Surface wave detection method and terminal device |
| CN111856555A (en) * | 2020-06-19 | 2020-10-30 | 同济大学 | Underground detection method based on surface wave multi-scale window analysis |
| CN112083487A (en) * | 2020-09-16 | 2020-12-15 | 中国科学技术大学 | Method and device for extracting broadband frequency dispersion curve |
| CN112083487B (en) * | 2020-09-16 | 2021-12-14 | 中国科学技术大学 | Method and device for extracting broadband frequency dispersion curve |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106094020B (en) | 2018-09-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6945895B2 (en) | Surface wave exploration methods and terminal devices | |
| CN106646615B (en) | A kind of data processing method and device of surface wave frequency dispersion curve | |
| CN107894613B (en) | Elastic wave vector imaging method, device, storage medium and equipment | |
| CN107728210B (en) | Multiple instruments acquire the determination method and apparatus that road is lacked in data | |
| CN106855640B (en) | A kind of submarine cable seismic data first arrival calculation method and device | |
| CN104316958B (en) | A kind of Coherent processing method for identifying different scale stratum breaking | |
| CN104597499B (en) | Adjacent shot interference suppressing method and device for independent simultaneous sweeping seismic data of seismic vibrators | |
| CN106094020A (en) | A kind of seismic inversion method and device | |
| CN109100795A (en) | A kind of inspection point distribution method, the apparatus and system of face element | |
| CN109765615A (en) | Stratum quality factor inversion method and device | |
| CN109001813A (en) | Method, device and system for suppressing multiple waves | |
| CN106501841B (en) | A kind of optimum design method and device of two dimension wode line seismic observation system | |
| CN109188520A (en) | Thin reservoir thickness prediction method and device | |
| CN107561574A (en) | Shot position, the detection method of receiving point position and device | |
| CN104182651B (en) | For the automatic quality control method in micro-seismic event azimuth that three-component geophone is received | |
| CN105301638B (en) | A kind of method and apparatus in extraction base of weathering face | |
| CN109709602A (en) | Remote detection sound wave offset imaging method, device and system | |
| CN106814394B (en) | Station formula wave detector combines low frequency restoration method with conventional wave detector | |
| CN106066491B (en) | A kind of high density wide-azimuth 3-d seismic exploration observation system optimization method | |
| CN109856672B (en) | Transient wave packet extracting method, storage medium and terminal based on depth wave-number spectrum | |
| CN111123361B (en) | Method and device for regularly reconstructing seismic data of vertical seismic profile and storage medium | |
| Díaz et al. | On the ability of the Benford’s Law to detect earthquakes and discriminate seismic signals | |
| CN107422375B (en) | The determination method and apparatus of the CRP gather of subsea node | |
| CN105425290B (en) | A kind of method and device of pre-stack time migration | |
| CN107976714A (en) | A kind of channel set computational methods of complicated earth surface classification spatial distance weighting |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |