CN106291678B - A kind of seismic data acquisition method and system - Google Patents
A kind of seismic data acquisition method and system Download PDFInfo
- Publication number
- CN106291678B CN106291678B CN201510280348.5A CN201510280348A CN106291678B CN 106291678 B CN106291678 B CN 106291678B CN 201510280348 A CN201510280348 A CN 201510280348A CN 106291678 B CN106291678 B CN 106291678B
- Authority
- CN
- China
- Prior art keywords
- observation system
- seismic data
- shot point
- acquisition
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a kind of seismic data acquisition methods, the described method comprises the following steps:Acquisition observation system is compressed in creation, wherein the geophone station of the compression acquisition observation system is uniformly distributed, and the shot point of the corresponding geophone station is sparse to be arranged and be unevenly distributed;The earthquake sampled data for meeting specific distribution rule is obtained based on the compression acquisition observation system;Processing is reconstructed to obtain the seismic data of the shot point of missing to finally obtain complete seismic data to the earthquake sampled data.The invention also provides a kind of system for acquiring seismic data.Compared with prior art, seismic data acquisition method of the present invention and system greatly reduce shot point number on the basis of not reducing data precision, its process is simply convenient to carry out, the collecting efficiency of seismic data can be effectively improved and reduce the acquisition cost of seismic data, a kind of means obtaining complete regular seismic data can also be provided for the area that cannot obtain regular big gun record due to objective factor.
Description
Technical field
The present invention relates to field of seismic exploration, a kind of seismic data acquisition method and system are in particulard relate to.
Background technique
During seismic prospecting, the acquisition of seismic data is usually carried out using earthquake data acquisition observation system.?
Earthquake data acquisition observation system sends Artificial Seismic Wave on shot position, is passed by receiving and processing Artificial Seismic Wave through stratum
Passback seismic wave after defeated reflection obtains corresponding formation data.In order to realize the above functions, earthquake data acquisition observation system
System is comprising generating the shot point of Artificial Seismic Wave and receiving the geophone station of passback seismic wave.
Usually theoretically, the number of shot point is more, the ground that finally obtains higher to the superposition degree of covering of survey area
Signal-to-noise ratio and the lateral resolution for shaking signal are also higher.Just need to increase the number of shot point to improve the precision of seismic signal.
But shot point increases the earthquake handled needed for the rising that certainly will bring earthquake data acquisition observation system cost of layout and later period
The growth of the data volume of wave number evidence.
It therefore, include often certain number of inspection according to the design of specific exploration demand in actual geological prospecting
The earthquake data acquisition observation system of wave point and shot point.To realize relatively optimal exploration effects and cost control.Existing
Earthquake data acquisition observation system in, the arrangement of shot point-geophone station meets the regular distribution of certain rule.In acquisition seismic data
When, collect the seismic data of each shot point respectively by all geophone stations.
But the number of the earthquake data acquisition observation system shot point that is included based on above method design is still opposite
More, the cost for setting up whole system is still relatively higher.Meanwhile by field nature acquisition condition (such as landform, river, road,
Village, building or marine streamer turning etc.) limitation, earthquake data acquisition observation system cannot acquire completely in some cases
Whole earthquake big gun records.
Therefore, it in order to further decrease required shot point number under the premise of not influencing the precision of seismic signal, repairs wild
Big gun record is irregularly lacked caused by outer objective factor, needs a kind of new seismic data acquisition method and system.
Summary of the invention
In order to further decrease required shot point number under the premise of not influencing the precision of seismic signal, the present invention provides
A kind of seismic data acquisition method, the described method comprises the following steps:
Creation compression acquisition observation system, wherein the geophone station of the compression acquisition observation system is uniformly distributed, corresponding institute
It states the sparse arrangement of shot point of geophone station and is unevenly distributed;
The earthquake sampled data for meeting specific distribution rule is obtained based on the compression acquisition observation system;
Processing is reconstructed to the earthquake sampled data to obtain the shot point seismic data of missing to finally obtain
Whole seismic data.
In one embodiment, during obtaining the earthquake sampled data, observation system is acquired according to the compression
The corresponding geophone station and the shot point are set and obtain the earthquake sampled data to acquire.
In one embodiment, during acquisition obtains the earthquake sampled data, based on compression acquisition observation
The shot point distribution of system obtains the earthquake hits from the common seismic data acquisition collected seismic data of observation system
According to.
In one embodiment, it during creating the compression acquisition observation system, is acquired according to common seismic data
Observation system constructs the compression and acquires observation system, wherein:
The detection of the geophone station distribution of the compression acquisition observation system and common seismic data acquisition observation system
Point distribution is consistent;
The shot point distribution of common seismic data acquisition observation system is sampled to obtain the compression acquisition and see
The shot point of examining system is distributed.
In one embodiment, it is distributed using optimization Bernoulli process distribution i.e. shake Bernoulli Jacob and is advised as the specific distribution
Rule.
In one embodiment, based on the shake Bernoulli Jacob distribution to the big gun of common seismic data acquisition observation system
Point distribution carries out sparse sampling to obtain the shot point distribution of the compression acquisition observation system.
In one embodiment, during processing is reconstructed to the earthquake sampled data:
According to the shot point distribution of the compression earthquake data acquisition observation system and earthquake sampled data sorting life
At big gun trace gather need to be mended;
Complete big gun trace gather is reconstructed based on the big gun trace gather that need to mend to obtain missing big gun trace gather;
The complete seismic data is generated based on the earthquake sampled data and the missing big gun trace gather.
In one embodiment, using fast Fourier transform and spectrum Projected restructing algorithm to the earthquake sampled data into
Row reconstruction processing.
The invention also provides a kind of system for acquiring seismic data, the system includes equally distributed according to ad hoc rules
Geophone station and it is corresponding with the geophone station be unevenly distributed and the shot point of sparse arrangement, wherein:
The distribution of the geophone station is consistent with the geophone station distribution of conventional earthquake data acquisition observation system;
The shot point for being distributed as the conventional earthquake data acquisition observation system of the shot point is distributed in specific distribution rule
Sampling under rule.
Compared with prior art, seismic data acquisition method of the present invention and system are big on the basis of not reducing data precision
Big to reduce shot point number, process is simply convenient to carry out, can effectively improve the collecting efficiency of seismic data and reduce earthquake number
According to acquisition cost, can also be provided for the area that cannot obtain regular big gun record due to objective factor it is a kind of obtain it is complete regularly
The means of shake data.
Other feature or advantage of the invention will illustrate in the following description.Also, Partial Feature of the invention or
Advantage will be become apparent by specification, or be appreciated that by implementing the present invention.The purpose of the present invention and part
Advantage can be realized or be obtained by step specifically noted in the specification, claims and drawings.
Detailed description of the invention
Attached drawing is used to provide further understanding of the present invention, and constitutes part of specification, with reality of the invention
It applies example and is used together to explain the present invention, be not construed as limiting the invention.In the accompanying drawings:
Fig. 1 is flow chart according to an embodiment of the invention;
Fig. 2 is seismic channel sparse sampling schematic diagram according to an embodiment of the invention;
Fig. 3 is analogue simulation stratum physical model and shot point, wave detector layout diagram under a simulated environment;
Fig. 4 a is the shot point layout diagram of 2D physical analogy common seismic data acquisition observation system under a simulated environment;
Fig. 4 b-4d is that the 2D physical analogy compression of different acquisition rate according to an embodiment of the invention under a simulated environment is adopted
Collect the shot point layout diagram of observation system;
Fig. 5 a is the earthquake big gun of the acquisition of common seismic data acquisition observation system and trace record position under a simulated environment
Schematic diagram;
Fig. 5 b is the earthquake big gun of the compression acquisition observation system of simulation according to an embodiment of the invention under a simulated environment
With trace record position view;
Fig. 6 a is the big gun record number according to an embodiment of the invention for compressing acquisition observation system acquisition under a simulated environment
According to body schematic diagram;
Fig. 6 b be under a simulated environment according to the present invention an embodiment reconstruct after big gun record data volume schematic diagram;
Fig. 6 c is the big gun record data volume schematic diagram of common seismic data acquisition observation system acquisition under a simulated environment;
Fig. 7 a is the 100th big gun record schematic diagram of common seismic data acquisition observation system acquisition under a simulated environment;
Fig. 7 b be under a simulated environment according to the present invention an embodiment reconstruct after the 100th big gun record schematic diagram.
Specific embodiment
Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings and examples, implementation personnel of the invention whereby
Can fully understand that how the invention applies technical means to solve technical problems, and reach technical effect realization process and according to
The present invention is embodied according to above-mentioned realization process.As long as each embodiment it should be noted that do not constitute conflict, in the present invention
And each feature in each embodiment can be combined with each other, be formed by technical solution protection scope of the present invention it
It is interior.
During seismic prospecting, the acquisition of seismic data is usually carried out using earthquake data acquisition observation system.?
Earthquake data acquisition observation system shot position sends Artificial Seismic Wave, is transmitted by receiving and processing Artificial Seismic Wave through stratum
Passback seismic wave after reflection obtains corresponding formation data.In order to realize the above functions, earthquake data acquisition observation system
Comprising generating the shot point of Artificial Seismic Wave and receiving the geophone station position of passback seismic wave.
Usually theoretically, the number of shot point is more, the ground that finally obtains higher to the superposition degree of covering of survey area
Signal-to-noise ratio and the lateral resolution for shaking signal are also higher.Just need to increase the number of shot point to improve the precision of seismic signal.
But shot point increases the earthquake handled needed for the rising that certainly will bring earthquake data acquisition observation system cost of layout and later period
The growth of the data volume of wave number evidence.
Therefore, it in actual seismic prospecting, is often set according to specific exploration demand and actual seismic geological condition
Earthquake data acquisition observation system of the meter comprising certain number of geophone station and shot point.To realize relatively optimal exploration effects
And cost control.In existing common seismic data acquisition observation system, shot point-geophone station arrangement meets the rule of certain rule
Then it is distributed.When acquiring seismic data, the seismic data of each shot point is collected respectively by all geophone stations.
But the common seismic data based on above method design acquire the number for the shot point that observation system is included still
Relatively more, the cost for setting up whole system is still relatively higher.Sometimes by field nature acquisition condition (such as landform, river, road
Road, village, building or marine streamer turning etc.) limitation, whole earthquake big gun records cannot be acquired completely.In order to not influence ground
Required shot point number is further decreased under the premise of the precision for shaking signal, the present invention is based on compressive sensing theories to propose one kind
System for acquiring seismic data.For ease of description, the system for acquiring seismic data of the difference prior art incite somebody to action this in the present specification
The system for acquiring seismic data that invention proposes is named as compression acquisition observation system, by earthquake data acquisition system in the prior art
It is referred to as common seismic data acquisition observation system.
Compressed sensing (Compressed sensing, abbreviation CS) theory is a kind of between mathematics and information science
Completely new signal acquisition, encoding and decoding are theoretical.Theoretical breakthrough classical signal acquisition must satisfy Nyquist-Shaimon
(Shannon/Nyquist) limitation of sampling thheorem, can be with far fewer than sampling number needed for traditional sampling theorem or observation
Points high quality reconstructs the original signal or image for meeting sampling thheorem.Utilize compressed sensing technology, so that it may using routinely
The seismic data on specific a part of shot point in shake data acquisition all shot points of observation system reconstructs on all shot points
Seismic data.
In this way, need to only acquire irregular, incomplete earthquake big gun/track data, then pass through later period room in practical exploration process
Interior processing quasi-complement lacks seismic channel/big gun, formation rule data volume.It can subtract under the premise of not reducing the precision of seismic signal
The number of the shot point of acquisition seismic data needed for few.
In the present invention, common seismic data in the prior art are acquired needed for observation system based on compressive sensing theory
The shot point of framework is wanted to be compressed to create the relative small number of compression acquisition observation system of shot point.In this way, compared to routine
The shot point number of system for acquiring seismic data, compression acquisition observation system greatly reduces.
Next based on attached drawing come the construction of the implementation procedure of detailed description of the present invention method and system of the invention
Process.The logical order of each step is shown, but in some cases in the flow chart of attached drawing, it can be to be different from herein suitable
Sequence executes shown or described step.
As shown in Figure 1, in the present embodiment, step S110 creation compression acquisition observation system is first carried out.In order to not
Required shot point number is further decreased under the premise of the precision for influencing seismic signal, in the present embodiment, is based on compressed sensing
Theory creation compression acquisition observation system.
When specific execution, step S111 is first carried out, includes certain number of detection according to the design of specific exploration demand
The common seismic data of point and shot point acquire observation system.In common seismic data acquisition observation system, shot point-geophone station
Arrangement meets the regular distribution of certain rule.In the present embodiment, by shot point-detection of common seismic data acquisition observation system
Point arrangement is referred to as regular geophone station distribution and regular shot point distribution.
In step S111, it is (true with step S1112 that step S1111 (determining rule geophone station distributing position) is executed respectively
Set pattern then shot point distributing position).Then step S112 is executed, determines the geophone station distribution of compression acquisition observation system.In step
In S112, the geophone station distribution of compression acquisition observation system and the regular geophone station of common seismic data acquisition observation system are distributed
It is consistent.
Next just it needs to be determined that the shot point distribution of compression acquisition observation system.Compressive sensing theory shows if signal
It is sparse or compressible in some transform domain, so that it may will be converted using one with the transformation incoherent random measurement matrix of base
On factor projection a to lower dimensional space, further according to these a small amount of observations, the high-quality of signal is realized by solving optimization problem
Amount reconstruct.Based on above-mentioned theory, the data of each shot point finally obtained using compression acquisition observation system be must satisfy centainly
The regularity of distribution complete shot point data finally can be just reconstructed with high quality.
Therefore, it during determining the shot point distribution for compressing acquisition observation system, is seen being acquired to common seismic data
The shot point distribution of examining system carries out the sampling based on specific distribution rule to obtain the shot point distribution of compression acquisition observation system.
In the present embodiment, using the distribution of shake Bernoulli Jacob's (Bernoulli) process as the specific distribution rule.It holds
Row step S113 is carried out based on shot point distributing position of the shake Bernoulli process distribution to common seismic data acquisition observation system
Sparse sampling is distributed with the shot point for obtaining compression acquisition observation system.Meet the rare shot point distribution of shake Bernoulli Jacob's regularity of distribution
The distribution that practical shot point can be effectively controlled avoids missing emplacement from setting and excessively concentrates, improves missing data reconstruction quality.
Specific in an application examples, if P is a sparse sampling factor, it indicates the complete seismic channel Nt for acquiring needs
Every road P is divided into one group, and every group has 1 to be sampled, by the random number of sampling seismic channel be 1, shake sparse sampling road number Ns be with
Machine number 1 adds up.Assume again that Nt is the multiple of P, to ensure to sample seismic channel number Ns=Nt/ P is integer, furthermore sparse sampling rate η
=Ns/Nt(being indicated with score or percentage).
Fig. 2 display meets the seismic channel sparse sampling schematic diagram of the shake Bernoulli regularity of distribution.In Fig. 2, often
A triangle represents a seismic channel, shares 20 seismic channels.In the case where compression acquires 1/5 data, 20 seismic channels point
At 4 groups (vertical line signal grouping situations), every group is made of 5 seismic channels, acquires 1 in each group, amounts to and acquires total road number
1/5, i.e., it only acquires 4 seismic channels (acquisition seismic channel is indicated with filled inverted triangles), then with this 4 seismic channel data weights acquired
Complete 20 earthquake record of structure.
In this way, the geophone station of the compression acquisition observation system finally obtained is uniformly distributed and (acquires and see for common seismic data
The geophone station of examining system is distributed), it corresponds to the sparse arrangement of shot point of the geophone station and is unevenly distributed and (compare and common seismic number
According to the shot point distribution of acquisition observation system, shot point number greatly reduces).
Step S120 can be executed after compression acquisition observation system creation, earthquake data step is acquired, based on pressure
Contracting acquisition observation system obtains the earthquake sampled data for meeting specific distribution rule.
In the step S120 of the present embodiment, according to compression acquisition observation system be arranged corresponding geophone station and shot point with
Obtain earthquake sampled data.
In the present embodiment, in order to which S140 can be executed after executing step S111 compared with S133 reconstruct complete documentation
The directly conventional complete seismic data of acquisition.Of course, compression seismic data can also be based on existing, the routine completed
The complete seismic data of earthquake data acquisition observation system carries out compression sampling acquisition.Earthquake sampled data is obtained in acquisition
During, the shot point distribution based on compression acquisition observation system acquires the collected earthquake of observation system from common seismic data
Earthquake sampled data is obtained in data.
In outer seismic data acquisition process out of office, seismic data is in space acquisition often by natural acquisition condition
The limitation at (such as landform, river, road, village, building or marine streamer turning), can lack earthquake big gun/road, along wave detector
The data volume of direction or direction of earthquake source acquisition is irregular or is unsatisfactory for Shannon-Nyquist sampling thheorem.In these cases,
Pair with the shot point distribution of common seismic data acquisition observation system can be distributed based on the shot point of compression acquisition observation system
It should be related to and sparse sampling is carried out to the seismic data that cannot acquire observation system acquisition with common seismic data to obtain completely
Regular earthquake sampled data.
Step S130 can be executed after getting compression earthquake sampled data, reconstruct missing seismic data step is right
Processing is reconstructed to obtain the seismic data of the shot point of missing to finally obtain complete seismic data in earthquake sampled data.
In the present embodiment, using fast Fourier transform (Fast Fourier Transformation, FFT) and spectrum Projected weight
Processing is reconstructed to earthquake sampled data in structure algorithm (SPGL1).
In step s 130, step S131 is first carried out, generation need to mend big gun trace gather step, acquire observation system according to compression
The earthquake big gun data sorting that the sparse sampling earthquake emplacement of upper distribution is set, which generates, need to mend big gun trace gather;
Then step S132 is executed, complete big gun trace gather step is obtained, to big gun trace gather application fast Fourier transform (Fast need to be mended
Fourier Transformation, FFT) and spectrum Projected restructing algorithm (SPGL1) carry out missing big gun trace gather data reconstruction
Processing, to obtain complete big gun trace gather;
Step S133 is finally executed, complete earthquake big gun recording step is obtained, is generated based on the sorting of complete big gun trace gather complete
Earthquake big gun record.
Next the application process and effect of method and system of the invention are described with a specific analog simulation example.
In a simulated environment, Fig. 3 is physical model and shot point, wave detector layout diagram.In Fig. 3, longitudinal coordinate is ground
The spread speed of layer depth and seismic wave on the stratum of different depth is illustrated.The lateral coordinates of the top are geophone station and big gun
Point lays signal, wherein geophone station is fixed (upper triangle indicates in figure), and shot point is mobile (lower triangle indicates in figure).
In order to compare conventional Christmas design and compression acquisition design effect, first against emulation shown in Fig. 3
Modelling routine 2D physical analogy common seismic data acquire observation system, and shot point is located in the range of geophone station distribution, big gun
Spacing 10m, geophone interval 10m acquire 256,256 big guns, 1mm sampling, 5000 sampling points.Shot position on together is as schemed
4a, the point on abscissa are indicated per the shot position for needing acquisition on together.Shot point number nz=256, shot position is covered with abscissa
So that being shown as a solid black line on the abscissa of Fig. 4 a.
2D physical analogy compression acquisition observation system is designed to choose according to compressive sensing theory and meet shake Bernoulli
The compression of Random sparseness sample mode acquires the observation system of 1/2 big gun, shot position such as Fig. 4 b, shot point number nz=on per pass
128;Compression acquires the observation system of 1/3 big gun, shot position such as Fig. 4 c, shot point number nz=85 on per pass;Compression 1/4 big gun of acquisition
Observation system, shot position such as Fig. 4 d, shot point number nz=64 on per pass.By comparing in lateral coordinates in Fig. 4 a-4d
The sparse degree of the point of shot position and corresponding shot point number are not difficult to find out that shot point number is successively reduced in Fig. 4 a-4d.
Illustrate quality reconstruction so that Fig. 4 b compresses the 2D physical simulation experiment of acquisition half big gun record as an example.This test acquisition
Earthquake big gun and trace record position are shown with plane net trrellis diagram, and as shown in Fig. 5 a and 5b, abscissa is seismic channel, and ordinate is ground
Shake big gun.Each mesh point corresponds to corresponding earthquake big gun and trace record position.Black line indicates that the 100th emplacement is set in figure.Fig. 5 a is aobvious
Show the full big gun grid that whole big gun records are acquired with conventional Christmas observation system, Fig. 5 b display acquires observation system pressure with new compression
Contracting acquires the plane grid of conventional half big gun record, and Bernoulli random distribution rule point is shaken by meeting in acquisition shot position
To determine.
Complex chart 4a-4d and Fig. 5 a, 5b are, it is apparent that needed for compression according to the present invention acquisition observation system
Shot point number greatly reduces.
It is as shown in Figure 6 a to the data volume of the half big gun record (setting in emplacement shown in Fig. 4 b) obtained with compression acquisition method.It answers
Missing big gun record is reconstructed with FFT and spectrum Projected restructing algorithm (SPGL1), the big gun record data volume after reconstruct is as schemed
Shown in 6b, Fig. 6 c is whole big gun data volumes that conventional vision systems acquire (setting in emplacement shown in Fig. 4 a).Fig. 6 b data volume of reconstruct
There is larger improvement compared with Fig. 6 a data volume imaging definition and precision and Fig. 6 b and Fig. 6 c graphic feature is about the same.Thus it proves
The feasibility and validity of restructuring procedure.
Fig. 7 a is that the 100th big gun that common seismic data acquire observation system acquisition under above-mentioned simulated environment records signal
Figure, Fig. 7 b are under above-mentioned simulated environment with the 100th big gun record of compression 1/2 big gun data reconstruction of acquisition.Comparison diagram 7a and figure
7b does not almost see the difference of the two.Thus the feasibility and validity of restructuring procedure are further proved.
Complex chart 6a-6c and Fig. 7 a, 7b this it appears that compare with reconstruct before shot point data, weigh according to the present invention
The precision of shot point data after structure, integrity degree all greatly improve.Shot point data after reconstructing according to the present invention can achieve existing
The precision and integrity degree of shot point data in technology.Thus the feasibility and validity of restructuring procedure of the invention are proved, thus
Illustrate the feasibility and validity of method and system of the invention.
To sum up, compared with prior art, seismic data acquisition method of the present invention and system are in the base for not reducing data precision
Shot point number is greatly reduced on plinth, process is simply convenient to carry out, can effectively improve the collection efficiency of seismic data and reduction
The compiling costs of seismic data can also provide a kind of obtain completely for the area that cannot obtain regular big gun record due to objective factor
The means of regular seismic data.
While it is disclosed that embodiment content as above but described only to facilitate understanding the present invention and adopting
Embodiment is not intended to limit the invention.Method of the present invention can also have other various embodiments.Without departing substantially from
In the case where essence of the present invention, those skilled in the art make various corresponding changes or change in accordance with the present invention
Shape, but these corresponding changes or deformation all should belong to scope of protection of the claims of the invention.
Claims (7)
1. a kind of seismic data acquisition method, which is characterized in that the described method comprises the following steps:
Creation compression acquisition observation system, wherein the geophone station of the compression acquisition observation system is uniformly distributed, the corresponding inspection
The shot point of wave point is sparse to be arranged and is unevenly distributed;
The earthquake sampled data for meeting specific distribution rule is obtained based on the compression acquisition observation system;
It is complete to finally obtain to obtain the shot point seismic data of missing that processing is reconstructed to the earthquake sampled data
Seismic data,
Wherein, during processing is reconstructed to the earthquake sampled data:
Big gun trace gather need to be mended by being generated according to the shot point distribution of the compression acquisition observation system and earthquake sampled data sorting;
Complete big gun trace gather is reconstructed based on the big gun trace gather that need to mend to obtain missing big gun trace gather;
The complete seismic data is generated based on the earthquake sampled data and the missing big gun trace gather.
2. the method according to claim 1, wherein during obtaining the earthquake sampled data, according to
The compression acquisition observation system is arranged the corresponding geophone station and the shot point and obtains the earthquake hits to acquire
According to.
3. the method according to claim 1, wherein acquisition obtain the earthquake sampled data during,
Shot point distribution based on the compression acquisition observation system acquires the collected seismic data of observation system from common seismic data
It is middle to obtain the earthquake sampled data.
4. the method according to claim 1, wherein during creating the compression and acquiring observation system,
Observation system, which is acquired, according to common seismic data constructs the compression acquisition observation system, wherein:
The geophone station point of the geophone station distribution of the compression acquisition observation system and common seismic data acquisition observation system
Cloth is consistent;
The shot point distribution of common seismic data acquisition observation system is sampled to obtain the compression acquisition observation system
The shot point of system is distributed.
5. according to the method described in claim 4, it is characterized in that, using optimization Bernoulli process distribution i.e. shake Bernoulli Jacob point
Cloth is as the specific distribution rule.
6. according to the method described in claim 5, it is characterized in that, being distributed based on the shake Bernoulli Jacob to the common seismic
The shot point distribution of data acquisition observation system carries out sparse sampling to obtain the shot point distribution of the compression acquisition observation system.
7. method according to claim 1 to 6, which is characterized in that using fast Fourier transform and spectrum projection ladder
Processing is reconstructed to the earthquake sampled data in degree restructing algorithm.
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 CN106291678A (en) | 2017-01-04 |
CN106291678B true 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) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108181644B (en) * | 2017-12-06 | 2020-01-07 | 中国石油天然气集团公司 | Sampling method and device |
CN108732619B (en) * | 2018-04-27 | 2021-02-26 | 中国科学院地质与地球物理研究所 | Seabed geophysical data acquisition method |
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 |
CN112782758B (en) * | 2019-11-11 | 2024-06-25 | 中国石油天然气集团有限公司 | State determination method and device for sparse sampling observation system |
CN111399036B (en) * | 2020-04-03 | 2021-07-06 | 中国海洋石油集团有限公司 | Design method of seismic data efficient acquisition observation system |
CN114415235B (en) * | 2022-03-28 | 2022-06-17 | 中国科学院地质与地球物理研究所 | Irregular seismic data acquisition method based on forward modeling |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101625417B (en) * | 2008-07-08 | 2011-09-07 | 中国石油集团东方地球物理勘探有限责任公司 | Method for optimizing design of vertical seismic profile observation system |
EP2383588B1 (en) * | 2010-04-30 | 2018-09-05 | Conoco Phillips 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
Also Published As
Publication number | Publication date |
---|---|
CN106291678A (en) | 2017-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106291678B (en) | A kind of seismic data acquisition method and system | |
Li et al. | Refraction corrected transmission ultrasound computed tomography for application in breast imaging | |
US8000513B2 (en) | System and method for 3D time of flight PET forward projection based on an exact axial inverse rebinning relation in fourier space | |
CN107837090A (en) | Scattered ray correction based on sinogram in computer tomography | |
CN108577876A (en) | A kind of static CT of polygon and its working method | |
TWI517093B (en) | Computer tomography reconstruction method | |
Zhou et al. | Three-dimensional Vp and Vp/Vs structure in the Longtan reservoir area by local earthquake tomography | |
Colombi et al. | Seismic waveform inversion for core–mantle boundary topography | |
CN111830557B (en) | Artificial fracture complexity index obtaining method and system based on fracturing microseism | |
CN111474574A (en) | Compressed sensing-based generation method and device for seismic acquisition observation system | |
CN108364326B (en) | CT imaging method | |
Qu et al. | 3-D Least-Squares Reverse Time Migration in Curvilinear-τ Domain | |
Haggag et al. | Seismicity and 3D velocity structure in the Aswan Reservoir Lake area, Egypt | |
CN110687598A (en) | Method and device for accelerating microseismic numerical simulation | |
CN115754007A (en) | Damage detection method based on acoustic emission technology and tomography technology | |
Diez Zaldivar et al. | Evaluation of the event detection level of the cuban seismic network | |
Vales et al. | Evidence of high lateral variations of coda Q from local earthquakes in Western Iberia and its SW offshore area | |
Pulido et al. | Rupture process of the 2007 Notohanto earthquake by using an isochrones back-projection method and K-NET/KiK-net data | |
CN101375182B (en) | An improved method of interpolation between a plurality of observed tensors | |
Lehmann et al. | Synthetic ground motions in heterogeneous geologies: the HEMEW-3D dataset for scientific machine learning | |
Moldoveanu et al. | Quad coil towed streamer marine acquisition | |
CN114442146A (en) | Seismic exploration data acquisition method and device | |
CN112782758A (en) | State determination method and device for sparse sampling observation system | |
CN105317433B (en) | Granulation mass microcosmic duct based on level-set function extracting method | |
CN111965694B (en) | Method and device for determining position of physical point of earthquake and earthquake observation system |
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 |