CN106291678A - A kind of seismic data acquisition method and system - Google Patents
A kind of seismic data acquisition method and system Download PDFInfo
- Publication number
- CN106291678A CN106291678A CN201510280348.5A CN201510280348A CN106291678A CN 106291678 A CN106291678 A CN 106291678A CN 201510280348 A CN201510280348 A CN 201510280348A CN 106291678 A CN106291678 A CN 106291678A
- Authority
- CN
- China
- Prior art keywords
- observation system
- shot point
- data
- earthquake
- distribution
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a kind of seismic data acquisition method, said method comprising the steps of: creating compression and gather observation system, wherein, the geophone station that described compression gathers observation system is uniformly distributed, the sparse arrangement of shot point of corresponding described geophone station and skewness；Gather observation system acquisition based on described compression and meet the earthquake sampled data of specific distribution rule；It is reconstructed described earthquake sampled data to process and finally obtains complete geological data to obtain the geological data of the shot point lacked.The invention allows for a kind of system for acquiring seismic data.Compared with prior art, seismic data acquisition method of the present invention and system are greatly reduced shot point number on the basis of not reducing data precision, its process is the most convenient to carry out, the collecting efficiency of geological data can be effectively improved and reduce the acquisition cost of geological data, it is possible to providing a kind of means obtaining complete rule seismic data for can not obtain the area of rule big gun record due to objective factor.
Description
Technical field
The present invention relates to field of seismic exploration, in particular relate to a kind of seismic data acquisition method and system.
Background technology
During seismic prospecting, it is frequently used earthquake data acquisition observation system to carry out adopting of geological data
Collection.Earthquake data acquisition observation system shot position sends Artificial Seismic Wave, artificial by receiving and processing
Seismic wave passback seismic wave after the transmission reflection of stratum obtains corresponding formation data.In order to realize above-mentioned merit
Can, earthquake data acquisition observation system comprises the shot point generating Artificial Seismic Wave and the inspection receiving passback seismic wave
Wave point.
The most in theory, the number of shot point is the most, the highest to the superposition degree of covering of survey area, finally obtains
Signal to noise ratio and the lateral resolution of the seismic signal taken are the highest.It is accomplished by increasing to improve the precision of seismic signal
Add the number of shot point.But shot point increase the rising that certainly will bring earthquake data acquisition observation system cost of layout with
And the growth of the data volume of the later stage required seismic data processed.
Therefore, in actual geological prospecting, often according to concrete exploration demand design packet containing given number
Geophone station and the earthquake data acquisition observation system of shot point.To realize relatively optimum exploration effects and cost
Control.In existing earthquake data acquisition observation system, shot point-geophone station arrangement meets the rule of certain rule
Distribution.When gathering geological data, collected the geological data of each shot point respectively by all geophone stations.
But the number of the shot point that earthquake data acquisition observation system based on said method design is comprised is still
The most, the cost setting up whole system is the most higher.Meanwhile, by field nature acquisition condition (as
Landform, river, road, village, building or marine streamer turning etc.) limit, earthquake in some cases
Data acquisition observation system can not gather whole earthquake big gun record completely.
Therefore, for the shot point number needed for reducing further under the precision premise not affecting seismic signal, repair
The irregular disappearance big gun record that multiple field objective factor causes, needs a kind of new seismic data acquisition method and is
System.
Summary of the invention
For the shot point number needed for reducing further under the precision premise not affecting seismic signal, the present invention carries
Supply a kind of seismic data acquisition method, said method comprising the steps of:
Creating compression and gather observation system, wherein, the geophone station that described compression gathers observation system is uniformly distributed,
The sparse arrangement of shot point of corresponding described geophone station and skewness；
Gather observation system acquisition based on described compression and meet the earthquake sampled data of specific distribution rule；
It is reconstructed described earthquake sampled data to process and finally obtains to obtain the shot point geological data of disappearance
Take whole geological data.
In one embodiment, during obtaining described earthquake sampled data, gather observation according to described compression
System arranges corresponding described geophone station and described shot point and obtains described earthquake sampled data to gather.
In one embodiment, during gathering the described earthquake sampled data of acquisition, gather based on described compression
The shot point distribution of observation system obtains institute from the geological data that common seismic data acquisition observation system collects
State earthquake sampled data.
In one embodiment, during creating described compression collection observation system, according to common seismic data
Gather observation system and build described compression collection observation system, wherein:
Described compression gathers the geophone station distribution of observation system and described common seismic data acquisition observation system
Geophone station distribution keeps consistent；
The shot point of described common seismic data acquisition observation system is distributed and samples to obtain described compression and adopt
The shot point distribution of collection observation system.
In one embodiment, the optimization i.e. shake Bernoulli Jacob of Bernoulli process distribution is used to be distributed as described specific point
Cloth rule.
In one embodiment, it is distributed described common seismic data acquisition observation system based on described shake Bernoulli Jacob
Shot point distribution carry out sparse sampling with obtain described compression gather observation system shot point distribution.
In one embodiment, during described earthquake sampled data is reconstructed process:
Shot point distribution and described earthquake sampled data according to described compression earthquake data acquisition observation system are divided
Choosing generates Xu Bubao road collection；
Complete big gun road collection is reconstructed to obtain disappearance big gun road collection based on described Xu Bubao road collection；
Described complete geological data is generated based on described earthquake sampled data and described disappearance big gun road collection.
In one embodiment, use fast Fourier transform and spectrum Projected restructing algorithm to described earthquake hits
According to being reconstructed process.
The invention allows for a kind of system for acquiring seismic data, described system comprises uniformly divides according to ad hoc rule
The geophone station of cloth and the skewness corresponding with described geophone station and the shot point of sparse arrangement, wherein:
The distribution of described geophone station is consistent with the distribution of the geophone station of conventional earthquake data acquisition observation system；
The shot point of the earthquake data acquisition observation system being distributed as described routine of described shot point is distributed in specific point
Sampling under cloth rule.
Compared with prior art, seismic data acquisition method of the present invention and system are on the basis not reducing data precision
On be greatly reduced shot point number, its process is the most convenient to carry out, it is possible to be effectively improved the collecting efficiency of geological data
And reduce the acquisition cost of geological data, it is possible to carry for can not obtain the area of rule big gun record due to objective factor
For a kind of means obtaining complete rule seismic data.
The further feature of the present invention or advantage will illustrate in the following description.Further, the part of the present invention is special
Levy or advantage will be become apparent by description, or be appreciated that by implementing the present invention.The present invention
Purpose and certain advantages can be come by step specifically noted in description, claims and accompanying drawing
Realize or obtain.
Accompanying drawing explanation
Accompanying drawing is for providing a further understanding of the present invention, and constitutes a part for description, with the present invention
Embodiment be provided commonly for explain the present invention, be not intended that limitation of the present invention.In the accompanying drawings:
Fig. 1 is flow chart according to an embodiment of the invention；
Fig. 2 is seismic channel sparse sampling mode schematic diagram according to an embodiment of the invention；
Fig. 3 is the analogue simulation stratum physical model under a simulated environment and shot point, cymoscope laying schematic diagram；
Fig. 4 a is that the shot point laying of 2D physical modeling common seismic data acquisition observation system under a simulated environment is shown
It is intended to；
Fig. 4 b-4d is the 2D physical modeling of different acquisition rate according to an embodiment of the invention under a simulated environment
Compression gathers the shot point of observation system and lays schematic diagram；
Fig. 5 a is earthquake big gun and the trace record of the collection of common seismic data acquisition observation system under a simulated environment
Position view；
Fig. 5 b is the ground that the compression simulated according to an embodiment of the invention under a simulated environment gathers observation system
Shake big gun and trace record position view；
Fig. 6 a is that under a simulated environment, compression gathers the big gun note that observation system gathers according to an embodiment of the invention
Record data volume schematic diagram；
Fig. 6 b is the big gun record data volume schematic diagram after reconstructing according to one embodiment of the invention under a simulated environment；
Fig. 6 c is the big gun record data volume signal that under a simulated environment, common seismic data acquisition observation system gathers
Figure；
Fig. 7 a is the 100th big gun record signal that under a simulated environment, common seismic data acquisition observation system gathers
Figure；
Fig. 7 b is the 100th big gun record schematic diagram after reconstructing according to one embodiment of the invention under a simulated environment.
Detailed description of the invention
Embodiments of the present invention are described in detail, whereby the enforcement of the present invention below with reference to drawings and Examples
Personnel can fully understand how application technology means are to solve technical problem for the present invention, and reaches technique effect
Realize process and be embodied as the present invention according to the above-mentioned process that realizes.As long as it should be noted that do not constitute conflict,
Each embodiment in the present invention and each feature in each embodiment can be combined with each other, the technology formed
Scheme is all within protection scope of the present invention.
During seismic prospecting, it is frequently used earthquake data acquisition observation system to carry out adopting of geological data
Collection.Artificial Seismic Wave is sent, by receiving and processing artificially in earthquake data acquisition observation system shot position
Seismic wave passback seismic wave after the transmission reflection of stratum obtains corresponding formation data.In order to realize above-mentioned merit
Can, earthquake data acquisition observation system comprises the shot point generating Artificial Seismic Wave and the inspection receiving passback seismic wave
Wave point position.
The most in theory, the number of shot point is the most, the highest to the superposition degree of covering of survey area, finally obtains
Signal to noise ratio and the lateral resolution of the seismic signal taken are the highest.It is accomplished by increasing to improve the precision of seismic signal
Add the number of shot point.But shot point increase the rising that certainly will bring earthquake data acquisition observation system cost of layout with
And the growth of the data volume of the later stage required seismic data processed.
Therefore, in actual seismic prospecting, often according to concrete exploration demand and actual seismic geology feelings
Condition design packet contains certain number of geophone station and the earthquake data acquisition observation system of shot point.To realize
Excellent exploration effects and cost control.In existing common seismic data acquisition observation system, shot point-geophone station
Arrangement meets the regular distribution of certain rule.When gathering geological data, collected respectively often by all geophone stations
The geological data of individual shot point.
But the number of the shot point that common seismic data acquisition observation system based on said method design is comprised
The most, the cost setting up whole system is the most higher.Sometimes by field nature acquisition condition (as
Landform, river, road, village, building or marine streamer turning etc.) limit, it is impossible to gather all completely
Earthquake big gun record.For the shot point number needed for reducing further on the premise of the precision not affecting seismic signal
Mesh, the present invention proposes a kind of system for acquiring seismic data based on compressive sensing theory.For ease of describing, difference
The system for acquiring seismic data of prior art, in this manual, the earthquake data acquisition system that the present invention is proposed
Named compression of uniting gathers observation system, and system for acquiring seismic data of the prior art is referred to as common seismic number
According to gathering observation system.
Compressed sensing (Compressed sensing is called for short CS) theory is a kind of between mathematics and information science
Between brand-new signals collecting, encoding and decoding theoretical.The classical signal collection of this theoretical breakthrough must is fulfilled for Nai Kuisi
The restriction of spy-Shannon (Shannon/Nyquist) sampling thheorem, can be with needed for traditional sampling theorem
Sampling number or the points of measurement high-quality reconstruct original signal or the image meeting sampling thheorem.Utilize compressed sensing
Technology, it is possible to utilize on the specific a part of shot point in the common seismic all shot points of data acquisition observation system
Geological data reconstruct the geological data on all shot points.
So, in actual exploration process, only need to gather earthquake big gun/track data irregular, incomplete, then pass through
Later stage indoor process quasi-complement disappearance seismic channel/big gun, formation rule data volume.Just can not reduce seismic signal
The number of the required shot point gathering geological data is reduced on the premise of precision.
In the present invention, based on compressive sensing theory to common seismic data acquisition observation system of the prior art
The shot point of required framework is compressed creating the relative small number of compression of shot point and gathers observation system.This
Sample, compared to common seismic data collecting system, compression gathers the shot point number of observation system and greatly reduces.
Next the execution process of the method for the present invention and the system of the present invention is described in detail based on accompanying drawing
Construction process.The logical order of each step shown in the flow chart of accompanying drawing, but in some cases, can be with
It is different from the step shown or described by order execution herein.
As it is shown in figure 1, in the present embodiment, step S110 is first carried out and creates compression collection observation system.
For the shot point number needed for reducing further on the premise of the precision not affecting seismic signal, at the present embodiment
In, create compression based on compressive sensing theory and gather observation system.
When specifically performing, step S111 is first carried out, according to concrete exploration demand design packet containing certain number of
The common seismic data acquisition observation system of geophone station and shot point.In common seismic data acquisition observation system
In, shot point-geophone station arrangement meets the regular distribution of certain rule.In the present embodiment, by common seismic data
Shot point-geophone station the arrangement gathering observation system is called rule geophone station distribution and rule shot point distribution.
In step S111, perform step S1111 (determining rule geophone station distributing position) and step respectively
S1112 (determines rule shot point distributing position).Then perform step S112, determine that compression gathers observation system
Geophone station distribution.In step S112, compression gathers geophone station distribution and the common seismic number of observation system
Keep consistent according to the regular geophone station distribution gathering observation system.
Next it is accomplished by determining the shot point distribution that compression gathers observation system.Compressive sensing theory shows, if
Signal is sparse or compressible at certain transform domain, it is possible to utilize one and the conversion incoherent random survey of base
Conversion coefficient is projected on a lower dimensional space by moment matrix, further according to these a small amount of observations, excellent by solving
Change problem realizes the high-quality reconstruct of signal.Based on above-mentioned theory, utilize compression to gather observation system and finally obtain
The data of each shot point must be fulfilled for certain regularity of distribution and the most just can reconstruct complete big gun with high-quality
Point data.
Therefore, during determining the shot point distribution that compression gathers observation system, to common seismic data acquisition
The shot point distribution of collection observation system carries out sampling based on specific distribution rule and gathers observation system to obtain compression
Shot point distribution.
In the present embodiment, shake Bernoulli Jacob (Bernoulli) process distribution is used to advise as described specific distribution
Rule.Perform step S113, based on the shake Bernoulli process distribution big gun to common seismic data acquisition observation system
Point distributing position carries out sparse sampling and gathers the shot point distribution of observation system to obtain compression.Meet shake Bernoulli Jacob
The rare shot point distribution of the regularity of distribution can effectively control the distribution of actual shot point, it is to avoid disappearance emplacement is put excessively
Concentrate, improve missing data reconstruction quality.
Specific in an application examples, if P is a sparse sampling factor, it represents will need collection intactly
Shake road Nt every P road is divided into one group, and often group has 1 road to be sampled, and the random number being sampled seismic channel is 1, shake
Sparse sampling number of channels Ns is the cumulative of random number 1.Assume again that Nt is the multiple of P, to guarantee seismic channel of sampling
Number Ns=Nt/ P is integer, in addition sparse sampling rate η=Ns/Nt(representing with mark or percent).
Fig. 2 shows the seismic channel sparse sampling mode schematic diagram meeting the shake Bernoulli regularity of distribution.At Fig. 2
In, each triangle represents a seismic channel, has 20 seismic channels.The situation of 1/5 data is gathered in compression
Under, 20 seismic channels are divided into 4 groups (vertical line signal packet situations), and often group is made up of 5 seismic channels, often
Group gathers 1 road, amounts to and gather the 1/5 of total number of channels, the most only gather 4 seismic channels and (gather seismic channel real
Heart inverted triangle represents), then reconstruct complete 20 road earthquake records by these 4 seismic channel data gathered.
So, the geophone station that the final compression obtained gathers observation system is uniformly distributed (for common seismic data acquisition
The geophone station distribution of collection observation system), the sparse arrangement of shot point of corresponding described geophone station and skewness (phase
The relatively shot point with common seismic data acquisition observation system is distributed, and shot point number greatly reduces).
Compression gathers after observation system creates just can perform step S120, gathers earthquake data step, base
Gather observation system acquisition in compression and meet the earthquake sampled data of specific distribution rule.
In step S120 of the present embodiment, according to compression gather observation system arrange corresponding geophone station and
Shot point is to obtain earthquake sampled data.
In the present embodiment, in order to, compared with S133 reconstruct complete documentation, perform after step S111,
Perform S140 and directly gather the complete geological data of routine.Certain, compression geological data can also be based on existing
Have, the complete geological data of the common seismic data acquisition observation system being complete, be compressed sampling and obtain
Take.I.e. obtain during earthquake sampled data gathering, based on compression gather the shot point distribution of observation system from
The geological data that common seismic data acquisition observation system collects obtains earthquake sampled data.
In the wild in seismic data acquisition process, geological data in space acquisition often by nature pick-up slip
The restriction of part (such as landform, river, road, village, building or marine streamer turning etc.), can lack ground
Shake big gun/road, the data volume along cymoscope direction or direction of earthquake source collection is irregular or is unsatisfactory for Shannon-Nyquist
Sampling thheorem.In these cases, it is possible to gather shot point distribution and the common seismic of observation system based on compression
The corresponding relation of the shot point distribution of data acquisition observation system is to not using common seismic data acquisition observation system
The geological data obtained carries out sparse sampling to obtaining complete rule earthquake sampled data.
Step S130, reconstruct disappearance geological data step just can be performed after getting compression earthquake sampled data
Suddenly, earthquake sampled data is reconstructed process finally to have obtained to obtain the geological data of the shot point of disappearance
Whole geological data.In the present embodiment, employing fast Fourier transform (Fast Fourier Transformation,
FFT) and spectrum Projected restructing algorithm (SPGL1) earthquake sampled data is reconstructed process.
In step s 130, step S131 is first carried out, generates Xu Bubao road collection step, gather according to compression
In observation system, the sparse sampling earthquake emplacement of distribution is put the earthquake big gun data sorting obtained and is generated Xu Bubao road collection；
Then perform step S132, obtain complete big gun road collection step, to the collection application fast Fourier transform of Xu Bubao road
(Fast Fourier Transformation, FFT) and spectrum Projected restructing algorithm (SPGL1) lack
Big gun road collection data reconstruction processes, to obtain complete big gun road collection；
Finally perform step S133, obtain complete earthquake big gun recording step, generated based on the sorting of complete big gun road collection
Whole earthquake big gun record.
Next describe with a concrete analog simulation example method and system of the present invention application process and
Effect.In a simulated environment, Fig. 3 is physical model and shot point, cymoscope laying schematic diagram.In figure 3,
Coordinate in longitudinal direction is depth of stratum and the seismic wave spread speed signal on the stratum of different depth.The top
Lateral coordinates be that geophone station and shot point lay signal, wherein, geophone station fixes (in figure, upper triangle represents),
Shot point moves (in figure, lower triangle represents).
In order to compare conventional Christmas design and compression acquisition design effect, first against shown in Fig. 3
Simulation model design routine 2D physical modeling's common seismic data acquisition observation system, shot point is positioned at geophone station and divides
In the range of cloth, shot interval 10m, geophone interval 10m, gather 256 roads, 256 big guns, 1mm samples, and 5000
Individual sampling point.Shot position on Yi Dao such as Fig. 4 a, the point on abscissa represents the shot point position needing to gather on every one
Put.Shot point number nz=256, shot position is covered with abscissa so that being shown as a reality on the abscissa of Fig. 4 a
Heart black line.
Design 2D physical modeling compresses collection observation system, according to compressive sensing theory, chooses and meets shake
The compression of Bernoulli Random sparseness sample mode gathers the observation system of 1/2 big gun, and the shot position on per pass is such as
Fig. 4 b, shot point number nz=128；The observation system of compression collection 1/3 big gun, the such as Fig. 4 c of the shot position on per pass,
Shot point number nz=85；Compression gathers the observation system of 1/4 big gun, the such as Fig. 4 d of the shot position on per pass, shot point number
Nz=64.By sparse degree and the big gun of correspondence of the point of the shot position in lateral coordinates in comparison diagram 4a-4d
Counting and be not difficult to find out, in Fig. 4 a-4d, shot point number reduces successively.
Gather by Fig. 4 b compression, as a example by the 2D physical simulation experiment of half big gun record, quality reconstruction is described.This test
The earthquake big gun and the trace record position plane net trrellis diagram that gather show, as shown in Fig. 5 a and 5b, abscissa is
Seismic channel, vertical coordinate is earthquake big gun.Each mesh point correspond to corresponding earthquake big gun and trace record position.In figure
Black line instruction the 100th emplacement is put.Fig. 5 a display conventional Christmas observation system gathers the full big gun net of whole big gun records
Lattice, Fig. 5 b display gathers observation system compression with new compression and gathers the plane grid of conventional half big gun record,
Gather shot position and determined by meeting shake Bernoulli random distribution rule point.
Complex chart 4a-4d and Fig. 5 a, 5b are it is apparent that gather observation system according to the compression of the present invention
Required shot point number greatly reduces.
Data volume such as Fig. 6 a to the half big gun record (emplacement shown in Fig. 4 b is put) obtained with compression acquisition method
Shown in.Disappearance big gun record is reconstructed by application FFT and spectrum Projected restructing algorithm (SPGL1), reconstruct
After big gun record data volume as shown in Figure 6 b, Fig. 6 c is conventional vision systems collection (emplacement shown in Fig. 4 a is put)
Whole big gun data volumes.Fig. 6 b data volume relatively Fig. 6 a data volume imaging definition and the precision of reconstruct have had bigger
Improve and Fig. 6 b and Fig. 6 c graphic feature is about the same.Thus prove feasibility and the effectiveness of restructuring procedure.
Fig. 7 a is that the 100th big gun record of common seismic data acquisition observation system collection shows under above-mentioned simulated environment
Being intended to, Fig. 7 b is the 100th big gun record gathering 1/2 big gun data reconstruction under above-mentioned simulated environment with compression.Right
Ratio Fig. 7 a and Fig. 7 b, does not almost see the difference of the two.Prove the feasibility of restructuring procedure the most further
And effectiveness.
Complex chart 6a-6c and Fig. 7 a, 7b are this it appears that shot point data before comparing and reconstructing, according to this
The precision of shot point data, integrity degree after invention reconstruct are all greatly improved.Shot point number after reconstructing according to the present invention
According to precision and the integrity degree that can reach shot point data in prior art.Thus prove the restructuring procedure of the present invention
Feasibility and effectiveness, thus the method for the present invention and the feasibility of system and effectiveness are described.
To sum up, compared with prior art, seismic data acquisition method of the present invention and system are not reducing data precision
On the basis of be greatly reduced shot point number, its process is the most convenient to carry out, it is possible to be effectively improved the receipts of geological data
Collect efficiency and reduce the compiling costs of geological data, it is possible to for rule big gun record can not be obtained due to objective factor
Area provides a kind of means obtaining complete rule seismic data.
While it is disclosed that embodiment as above, but described content is only to facilitate understand the present invention
And the embodiment used, it is not limited to the present invention.Method of the present invention also can have other multiple realities
Execute example.Without departing from the spirit of the present invention, those of ordinary skill in the art are when making according to the present invention
Go out various corresponding change or deformation, but these change accordingly or deform the claim that all should belong to the present invention
Protection domain.
Claims (9)
1. a seismic data acquisition method, it is characterised in that said method comprising the steps of:
Creating compression and gather observation system, wherein, the geophone station that described compression gathers observation system is uniformly distributed,
The sparse arrangement of shot point of corresponding described geophone station and skewness；
Gather observation system acquisition based on described compression and meet the earthquake sampled data of specific distribution rule；
It is reconstructed described earthquake sampled data to process and finally obtains to obtain the shot point geological data of disappearance
Take whole geological data.
Method the most according to claim 1, it is characterised in that obtaining described earthquake sampled data
During, gather observation system according to described compression and corresponding described geophone station and described shot point are set to gather
Obtain described earthquake sampled data.
Method the most according to claim 1, it is characterised in that obtain described earthquake hits gathering
During according to, the shot point distribution gathering observation system based on described compression from common seismic data acquisition observation is
The geological data that system collects obtains described earthquake sampled data.
Method the most according to claim 1, it is characterised in that gather observation system creating described compression
During system, build described compression according to common seismic data acquisition observation system and gather observation system, wherein:
Described compression gathers the geophone station distribution of observation system and described common seismic data acquisition observation system
Geophone station distribution keeps consistent；
The shot point of described common seismic data acquisition observation system is distributed and samples to obtain described compression and adopt
The shot point distribution of collection observation system.
Method the most according to claim 4, it is characterised in that use and optimize Bernoulli process distribution i.e.
Shake Bernoulli Jacob be distributed as described specific distribution rule.
Method the most according to claim 5, it is characterised in that it is right to be distributed based on described shake Bernoulli Jacob
The shot point distribution of described common seismic data acquisition observation system carries out sparse sampling and gathers to obtain described compression
The shot point distribution of observation system.
7. according to the method according to any one of claim 1-6, it is characterised in that described earthquake is being sampled
During data are reconstructed process:
Shot point distribution and described earthquake sampled data according to described compression earthquake data acquisition observation system are divided
Choosing generates Xu Bubao road collection；
Complete big gun road collection is reconstructed to obtain disappearance big gun road collection based on described Xu Bubao road collection；
Described complete geological data is generated based on described earthquake sampled data and described disappearance big gun road collection.
Method the most according to claim 7, it is characterised in that use fast Fourier transform and spectrum projection
Gradient reconstruction algorithm is reconstructed process to described earthquake sampled data.
9. a system for acquiring seismic data, it is characterised in that described system comprises according to ad hoc rule uniform
The geophone station of distribution and the skewness corresponding with described geophone station and the shot point of sparse arrangement, wherein:
The distribution of described geophone station is consistent with the distribution of the geophone station of conventional earthquake data acquisition observation system；
The shot point of the earthquake data acquisition observation system being distributed as described routine of described shot point is distributed in specific point
Sampling under cloth rule.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510280348.5A CN106291678B (en) | 2015-05-27 | 2015-05-27 | A kind of seismic data acquisition method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510280348.5A CN106291678B (en) | 2015-05-27 | 2015-05-27 | A kind of seismic data acquisition method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106291678A true CN106291678A (en) | 2017-01-04 |
| CN106291678B CN106291678B (en) | 2018-11-20 |
Family
ID=57635586
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510280348.5A Active CN106291678B (en) | 2015-05-27 | 2015-05-27 | A kind of seismic data acquisition method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106291678B (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108181644A (en) * | 2017-12-06 | 2018-06-19 | 中国石油天然气集团公司 | A kind of method of sampling and device |
| CN108732619A (en) * | 2018-04-27 | 2018-11-02 | 中国科学院地质与地球物理研究所 | A kind of submarine geophysics collecting method |
| CN111399036A (en) * | 2020-04-03 | 2020-07-10 | 中国海洋石油集团有限公司 | Design method of seismic data efficient acquisition observation system |
| CN111474574A (en) * | 2019-01-23 | 2020-07-31 | 中国石油天然气集团有限公司 | Compressed sensing-based generation method and device for seismic acquisition observation system |
| CN112782758A (en) * | 2019-11-11 | 2021-05-11 | 中国石油天然气集团有限公司 | State determination method and device for sparse sampling observation system |
| CN112782753A (en) * | 2019-11-07 | 2021-05-11 | 中国石油天然气集团有限公司 | Method and device for generating compressed sensing observation system |
| CN114415235A (en) * | 2022-03-28 | 2022-04-29 | 中国科学院地质与地球物理研究所 | Irregular seismic data acquisition method based on forward modeling |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101625417A (en) * | 2008-07-08 | 2010-01-13 | 中国石油集团东方地球物理勘探有限责任公司 | Method for optimizing design of vertical seismic profile observation system |
| EP2383588A2 (en) * | 2010-04-30 | 2011-11-02 | ConocoPhillips Company | Caterpillar-style seismic data acquisition using autonomous, continuously recording seismic data recorders |
-
2015
- 2015-05-27 CN CN201510280348.5A patent/CN106291678B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101625417A (en) * | 2008-07-08 | 2010-01-13 | 中国石油集团东方地球物理勘探有限责任公司 | Method for optimizing design of vertical seismic profile observation system |
| EP2383588A2 (en) * | 2010-04-30 | 2011-11-02 | ConocoPhillips Company | Caterpillar-style seismic data acquisition using autonomous, continuously recording seismic data recorders |
Non-Patent Citations (3)
| Title |
|---|
| CAI RUI ET AL.: "Bernoulli-based random undersampling schemes for 2D seismic data regularization", 《APPLIED GEOPHYSICS》 * |
| 王汉闯等: "地震数据高效采集方法初步研究", 《中国地球科学联合学术年会 2014》 * |
| 陈生昌等: "稀疏性约束的地球物理数据高效采集方法初步研究", 《石油物探》 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108181644A (en) * | 2017-12-06 | 2018-06-19 | 中国石油天然气集团公司 | A kind of method of sampling and device |
| CN108732619A (en) * | 2018-04-27 | 2018-11-02 | 中国科学院地质与地球物理研究所 | A kind of submarine geophysics collecting method |
| CN108732619B (en) * | 2018-04-27 | 2021-02-26 | 中国科学院地质与地球物理研究所 | Seabed geophysical data acquisition method |
| CN111474574A (en) * | 2019-01-23 | 2020-07-31 | 中国石油天然气集团有限公司 | Compressed sensing-based generation method and device for seismic acquisition observation system |
| CN111474574B (en) * | 2019-01-23 | 2023-04-25 | 中国石油天然气集团有限公司 | Method and device for generating seismic acquisition observation system based on compressed sensing |
| CN112782753A (en) * | 2019-11-07 | 2021-05-11 | 中国石油天然气集团有限公司 | Method and device for generating compressed sensing observation system |
| CN112782758A (en) * | 2019-11-11 | 2021-05-11 | 中国石油天然气集团有限公司 | State determination method and device for sparse sampling observation system |
| CN111399036A (en) * | 2020-04-03 | 2020-07-10 | 中国海洋石油集团有限公司 | Design method of seismic data efficient acquisition observation system |
| CN111399036B (en) * | 2020-04-03 | 2021-07-06 | 中国海洋石油集团有限公司 | Design method of seismic data efficient acquisition observation system |
| CN114415235A (en) * | 2022-03-28 | 2022-04-29 | 中国科学院地质与地球物理研究所 | Irregular seismic data acquisition method based on forward modeling |
| CN114415235B (en) * | 2022-03-28 | 2022-06-17 | 中国科学院地质与地球物理研究所 | Irregular seismic data acquisition method based on forward modeling |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106291678B (en) | 2018-11-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106291678A (en) | A kind of seismic data acquisition method and system | |
| CN104142514B (en) | Three-dimensional earthquake observing system quantitative designing method | |
| CN109407143B (en) | Seismic exploration irregular observation system design method based on compressed sensing | |
| WO2020123084A1 (en) | Machine learning-augmented geophysical inversion | |
| CN104268412B (en) | A kind of angle gathers ray tomography migration velocity analysis method and device | |
| CN104280775B (en) | Microseism monitoring and positioning method based on full-waveform vector offset superposition | |
| CN108445538B (en) | The method and system of Depth Domain layer Q model is established based on reflected seismic information | |
| CN105659111A (en) | System and method for regularizing seismic data | |
| CN103630120B (en) | Martian surface linear array image core line method for resampling based on tight geometric model | |
| CN106249291A (en) | A kind of high precision seismic data re-establishing method based on two-dimentional non-homogeneous warp wavelet | |
| CN110031896A (en) | Earthquake stochastic inversion methods and device based on Multiple-Point Geostatistics prior information | |
| CN104345340A (en) | Jitter sparse sampling method in seismic data regularization | |
| CN109188506A (en) | A kind of pure earth's surface stereo observing system suitable for high-speed rail tunnel bottom earthquake CT | |
| CN105093319A (en) | Ground micro-seismic static correction method based on three-dimensional seismic data | |
| CN105301636A (en) | Method and device for establishing speed model | |
| WO2021127382A1 (en) | Full waveform inversion in the midpoint-offset domain | |
| CN104570116A (en) | Geological marker bed-based time difference analyzing and correcting method | |
| CN109655890B (en) | Depth domain shallow-medium-deep layer combined chromatography inversion speed modeling method and system | |
| CN111638551A (en) | Seismic first-motion wave travel time chromatography method and device | |
| CN111474574A (en) | Compressed sensing-based generation method and device for seismic acquisition observation system | |
| CN103592684A (en) | Massive seismic data compression method and device for preserving spatial attribute information | |
| CN104280774B (en) | Quantitive analysis method of single-frequency seismic scattering noise | |
| CN105676292A (en) | 3D earthquake data de-noising method based on 2D curvelet transform | |
| CN113534259A (en) | Vibroseis efficient acquisition real-time prestack time migration imaging method | |
| Qu et al. | 3-D Least-Squares Reverse Time Migration in Curvilinear-τ Domain |
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 |