CN109031417A - The bilateral synthesis beam offset imaging method and system picked up based on lineups - Google Patents

The bilateral synthesis beam offset imaging method and system picked up based on lineups Download PDF

Info

Publication number
CN109031417A
CN109031417A CN201810777046.2A CN201810777046A CN109031417A CN 109031417 A CN109031417 A CN 109031417A CN 201810777046 A CN201810777046 A CN 201810777046A CN 109031417 A CN109031417 A CN 109031417A
Authority
CN
China
Prior art keywords
plane wave
local plane
big gun
face element
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.)
Granted
Application number
CN201810777046.2A
Other languages
Chinese (zh)
Other versions
CN109031417B (en
Inventor
朱兆林
曹丹平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum East China
Original Assignee
China University of Petroleum East China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN201810777046.2A priority Critical patent/CN109031417B/en
Publication of CN109031417A publication Critical patent/CN109031417A/en
Application granted granted Critical
Publication of CN109031417B publication Critical patent/CN109031417B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30—Analysis
    • G01V1/301—Analysis for determining seismic cross-sections or geostructures
    • G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/50—Corrections or adjustments related to wave propagation
    • G01V2210/51—Migration

Abstract

The present invention discloses a kind of bilateral synthesis beam offset imaging method and system based on lineups pickup.The method and system divide big gun using second grade segmentation (rectangular mesh+auxiliary is oval) method and examine face element, solve the problems, such as the quick indexing of seismic channel in face element；It is sparse feature using Local plane wave data, effective data information is picked up using lineups pick-up method, solve the problems, such as effective expression, effective storage and fast excursion of higher-dimension (five dimensions) the plane wave conversion data in Shu Hecheng, it is stored again by compression, then fast excursion is carried out again using the lineups of effective small number, data volume to be treated in bilateral synthesis beam migration process is effectively reduced, the offset efficiency of seismic migration imaging is substantially increased.

Description

The bilateral synthesis beam offset imaging method and system picked up based on lineups
Technical field
The present invention relates to geophysical prospecting for oil Seismic Data Processing Technique fields, pick up more particularly to based on lineups The bilateral synthesis beam offset imaging method and system taken.
Background technique
Seismic data acquisition, processing and explanation are key step and content of the seismic technology in mine locating application, and earthquake is inclined The important component that imaging method is seism processing is moved, this method can totally be divided into integral class and wave equation side Method, beam offset imaging method are to belong to long-pending classification method, it is one integrated between wave equation approach between gram western Hough Kind is realized.Beam migration and imaging techniques have been extensively studied and have applied in oil seismic exploration field, Hill (1990) and Hale (1991) has delivered Gaussian beam post-stack migration method article respectively, elaborates basic principle, derivation and the realization of Gaussian beam； Sun et al (2000) has been delivered as a gram gram western hoeve's bundle offset method for western Hough offset fast algorithm, is illustrated The implementation and efficiency quantitative estimation of gram western hoeve's bundle offset method；Hill (2001) has further delivered common offset Gaussian beam Prestack migration method, on its basis, Gray (2007) and Gray (2009) have delivered total big gun and the altogether height of angle domain respectively This beam prestack migration method；Subsequent Geophysicist, which have delivered, numerous about Gaussian beam offset method and simplifies beam offset side The research paper of method.Existing beam offset method and its reduced form are all at unilateral domain (big gun domain or common offset altogether) It realizes, i.e., records the imaging domain finally intersected in shot point beam and geophone station beam after unilateral Local plane wave synthesis to seismic receiver Imaging, even if gram western hoeve's bundle offset method theory ray parameter of Sun et al (2000) can be the vector of higher-dimension, but still So emphasize the method that two vector parameters are reduced to a vector parameter, i.e., unilateral offset is realized.It is bilateral in practical application Local plane wave synthesis is just carried out to seismic data and then carries out beam decomposition migration imaging again, is the theoretical and ideal choosing of comparison It selects, direction be entered to available beam can select big gun inspection is bilateral simultaneously, but the data volume of bilateral synthesis is larger, particularly with For three dimensional seismic data, it is necessary first to which three dimensional seismic data is transformed into quintuple space, and (shot point side is three-dimensional, geophone station side three Dimension shares time dimension), binning and Local plane wave synthesis (having more bidimensional) then are carried out along plane to shot point side three-dimensional, Binning and Local plane wave synthesis (and having more bidimensional) then are carried out along plane to geophone station side three-dimensional again, finally carried out again Migration imaging.These processes, are difficult to realize in practical applications, and a large amount of data redundancy, operation efficiency are very low.
Summary of the invention
The object of the present invention is to provide it is a kind of based on lineups pick up bilateral synthesis beam offset imaging method and system, with Solve the problems, such as that existing bilateral synthesis beam offset imaging method data processing amount is big, operation efficiency is low.
To achieve the above object, the present invention provides following schemes:
The bilateral synthesis beam offset imaging method picked up based on lineups, which comprises
Obtain earthquake record data and external definition parameter；The earthquake record data include the time sampling of seismic data Number, time sampling interval, frequency range and big gun examine coordinate range；The external definition parameter includes face element defined parameters and speed Defined parameters: when the face element defined parameters include bin size, Local plane wave orientation angle interval and number and output Effective lineups pick up number；The speed defined parameters include earth's surface speed and average speed；
According to the earthquake record data and the face element defined parameters, big gun inspection face element is carried out using second grade segmentation method and is drawn Point, obtain multiple big gun inspection face elements；The big gun inspection face element includes shot point face element and geophone station face element；
The face element index of multiple big gun inspection face elements is established according to the earthquake record data；
The pickup of energy spectrum attribute is carried out to multiple big gun inspection face elements according to face element index and valid data store, it is raw At multiple Local plane wave compressed data blocks；It include plane wave data header file in each Local plane wave compressed data block And data file；
The offset distance of each Local plane wave is calculated according to the Local plane wave compressed data block；
The corresponding trace gather Local plane wave index of the Local plane wave is established according to the offset distance of the Local plane wave Number；
All imaging roads are completed according to the trace gather Local plane wave call number and the Local plane wave compressed data block Concentrate the migration imaging of all Local plane wave compressed data blocks.
Optionally, described according to the earthquake record data and the face element defined parameters, using second grade segmentation method into Row big gun examines binning, obtains multiple big gun inspection face elements, specifically includes:
Coordinate range, the bin size, plane wave orientation angle interval and number, which are examined, according to the big gun determines office The big gun of facial planes wave number examines plane sample space；
Big gun inspection plane sample space is divided into the rectangular mesh of rule according to the bin size, generates level-one big gun Examine face element；The level-one big gun inspection face element includes a shot point rectangle face element and a geophone station rectangle face element；
According to the frequency range of seismic data and the average speed, calculated using Gaussian beam offset parameter formula auxiliary Help elliptical radius；
Determine that auxiliary is oval according to the elliptical radius of auxiliary；
The big gun inspection face element is marked off according to the auxiliary ellipse；The big gun inspection face element includes shot point face element and geophone station face Member.
Optionally, described energy spectrum attributes are carried out to multiple big guns inspection face elements according to face element index to pick up and effectively Data storage, generates multiple Local plane wave compressed data blocks, specifically includes:
Search for the seismic channel for being included in i-th of big gun inspection face element；
The seismic channel data of the seismic channel is read according to face element index；
According to the seismic channel data, using the Local plane wave of big gun inspection face element described in i-th of transformation calculations of linear pulling Data；
Using the auto-correlation formula of Local plane wave data The Correlated Spectroscopy of Local plane wave data is calculated；
Number and the Correlated Spectroscopy are picked up according to effective lineups, using formulaPick up effective Local plane wave data；
The position of effective lineups in effective Local plane wave data is determined according to the amplitude power of the Correlated Spectroscopy And range；
According to the location and range of effective lineups, using compressor compresses and Local plane wave data head text is stored Part and data file generate i-th of Local plane wave compressed data block；The Local plane wave header file includes big gun Have described in the inspection coordinate information of face element, the location information of affiliated file, the Local plane wave data storage size of generation and one Imitate lineups；
Judge whether i is less than n, obtains the first judging result；N is the quantity that the big gun examines face element, 1≤i < n；
If first judging result is that i is less than n, i=i+1 is enabled, institute in big gun inspection face element described in i-th of described search of return The step of seismic channel for including；
If first judging result is that i is not less than n, output the 1st to the n Local plane wave compressed data blocks.
Optionally, the offset distance according to the Local plane wave establishes the corresponding trace gather part of the Local plane wave Plane wave call number, specifically includes:
According to the offset distance of the Local plane wave, the Local plane wave is established using formula Hi=(H1-H0)/DH+1 Corresponding trace gather Local plane wave call number；Wherein Hi indicates the corresponding trace gather Local plane wave index of i-th of Local plane wave Number；H1 indicates the offset distance of the Local plane wave；H0 indicates first offset distance in offset gather parameter information；DH table Show offset distance interval.
Optionally, described to be completed according to the trace gather Local plane wave call number and the Local plane wave compressed data block The migration imaging of all Local plane wave compressed data blocks in all imaging trace gathers, specifically includes:
Obtain i-th of Local plane wave compressed data block in k-th of trace gather；
Obtain the effective lineups stored in i-th of Local plane wave compressed data block spatial position and when meta position It sets；
The ray tracing of big gun inspection position is carried out according to the spatial position and the time location, and crossover location is carried out Migration imaging generates imaging temporary file；
Judge whether i is less than L, obtains the second judging result；Wherein L is Local plane wave data described in k-th of trace gather The quantity of compression blocks；
If second judging result is that i is less than L, i=i+1 is enabled, returns to acquisition i-th of part plan wave number According to the effective lineups stored in compression blocks spatial position and time location the step of；
If second judging result is that i is not less than L, judge whether k is less than K, obtains third judging result；Wherein K is The quantity of trace gather is imaged；
If the third judging result is that k is less than K, k=k+1 is enabled, returns to i-th innings obtained in k-th of trace gather The step of facial planes wave compressed data block；
If the third judging result is that k is not less than K, merges the 1st imaging temporary file for arriving K trace gather, generate imaging Trace gather file.
The invention also discloses a kind of bilateral synthesis beam migration imaging system picked up based on lineups, the system packets It includes:
Data and parameter acquisition module, for obtaining earthquake record data and external definition parameter；The earthquake record number Coordinate range is examined according to the time sampling number, time sampling interval, frequency range and the big gun that include seismic data；The external definition ginseng Number includes face element defined parameters and speed defined parameters: the face element defined parameters include bin size, Local plane wave direction Effective lineups pick up number when angle interval and number and output；The speed defined parameters include earth's surface speed and are averaged Speed；
Big gun examines binning module, is used for according to the earthquake record data and the face element defined parameters, using second level Dividing method carries out big gun and examines binning, obtains multiple big gun inspection face elements；The big gun inspection face element includes shot point face element and geophone station face Member；
Face element index establishes module, for establishing the face element rope of multiple big gun inspection face elements according to the earthquake record data Draw；
Local plane wave compressed data block generation module, for according to face element index to multiple big guns examine face elements into Row energy spectrum attribute picks up and valid data storage, generates multiple Local plane wave compressed data blocks；Each part plan It include plane wave data header file and data file in wave compressed data block；
Offset distance computing module, for calculating the inclined of each Local plane wave according to the Local plane wave compressed data block Move away from；
Trace gather Local plane wave call number establishes module, for establishing the office according to the offset distance of the Local plane wave The corresponding trace gather Local plane wave call number of facial planes wave；
Migration imaging module, for according to the trace gather Local plane wave call number and the Local plane wave data compression Block completes the migration imaging of all Local plane wave compressed data blocks in all imaging trace gathers.
Optionally, the big gun inspection binning module specifically includes:
Big gun examines plane sample space determination unit, for examining the coordinate range, bin size, described flat according to the big gun Surface wave orientation angle interval and number determine the big gun inspection plane sample space of part plan wave number；
Level-one big gun examines face element generation unit, for being divided into big gun inspection plane sample space according to the bin size The rectangular mesh of rule generates level-one big gun and examines face element；The level-one big gun inspection face element includes a shot point rectangle face element and an inspection Wave point rectangle face element；
Assist ellipse radii computing unit, for according to seismic data the frequency range and the average speed, adopt It is calculated with Gaussian beam offset parameter formula and assists elliptical radius；
Oval determination unit is assisted, for determining that auxiliary is oval according to the elliptical radius of auxiliary；
Big gun examines binning unit, examines face element for marking off the big gun according to the auxiliary ellipse；The big gun examines face element Including shot point face element and geophone station face element.
Optionally, the Local plane wave compressed data block generation module specifically includes:
Seismic channel search unit, for searching for the seismic channel for being included in i-th of big gun inspection face element；
Seismic channel data reading unit, for reading the seismic channel data of the seismic channel according to face element index；
Local plane wave Data Computation Unit, for pulling transformation calculations i-th using linear according to the seismic channel data The Local plane wave data of a big gun inspection face element；
Correlated Spectroscopy computing unit, for the auto-correlation formula using Local plane wave dataThe Correlated Spectroscopy of Local plane wave data is calculated；
Effective Local plane wave data pick-up unit, for according to effective lineups pick up number to it is described related Spectrum, using formulaPick up effective Local plane wave data；
Effective lineups position range determination unit, determines effective office for the amplitude power according to the Correlated Spectroscopy The location and range of facial planes wave number effective lineups in；
Local plane wave compressed data block generation unit is used for the location and range according to effective lineups Compressor compresses simultaneously store Local plane wave header file and data file, generate i-th of Local plane wave data pressure Contracting block；The Local plane wave header file includes the coordinate information of big gun inspection face element, the location information of affiliated file, generation Local plane wave data storage size and effective lineups；
First judging unit obtains the first judging result for judging whether i is less than n；N is the number that the big gun examines face element Amount, 1≤i < n；
First return unit, for enabling i=i+1, returning to the seismic channel when first judging result is that i is less than n Search unit；
Local plane wave compressed data block output unit, for exporting when first judging result is that i is not less than n 1st to the n Local plane wave compressed data block.
Optionally, the trace gather Local plane wave call number is established module and is specifically included:
Trace gather Local plane wave call number establishes unit, for the offset distance according to the Local plane wave, using formula Hi=(H1-H0)/DH+1 establishes the corresponding trace gather Local plane wave call number of the Local plane wave；Wherein Hi is indicated i-th The corresponding trace gather Local plane wave call number of Local plane wave；H1 indicates the offset distance of the Local plane wave；H0 indicates offset Away from first offset distance in trace gather parameter information；DH indicates offset distance interval.
Optionally, the migration imaging module specifically includes:
Local plane wave compressed data block acquiring unit, for obtaining i-th of Local plane wave data in k-th of trace gather Compression blocks；
Effective lineups position acquisition unit stores in i-th of Local plane wave compressed data block for obtaining The spatial position of effective lineups and time location；
Migration imaging unit, the ray for carrying out big gun inspection position according to the spatial position and the time location chase after Track, and migration imaging is carried out to crossover location, generate imaging temporary file；
Second judgment unit obtains the second judging result for judging whether i is less than L；Wherein L is institute in k-th of trace gather State the quantity of Local plane wave compressed data block；
Second return unit returns described effectively same for enabling i=i+1 when second judging result is that i is less than L Phase shaft position acquiring unit；
Third judging unit obtains for judging whether k is less than K when second judging result is that i is not less than L Three judging results；Wherein K is the quantity that trace gather is imaged；
Third return unit, for enabling k=k+1, it is flat returning to the part when the third judging result is that k is less than K Surface wave compressed data block acquiring unit；
Trace gather file generating unit is imaged, for merging the 1st to K when the third judging result is that k is not less than K The imaging temporary file of trace gather generates imaging trace gather file.
The specific embodiment provided according to the present invention, the invention discloses following technical effects:
The present invention provides a kind of bilateral synthesis beam offset imaging method picked up based on lineups and system, the methods And system divides big gun using second grade segmentation (rectangular mesh+auxiliary oval) method and examines face element, solves the fast of seismic channel in face element Fast index problem；It is sparse feature using Local plane wave data, effective data is picked up using lineups pick-up method and are believed Breath solves the problems, such as effective expression, effective storage and fast excursion of higher-dimension (five dimensions) the plane wave conversion data in Shu Hecheng, It is stored again by compression, then carries out fast excursion again using the lineups of effective small number, effectively reduce bilateral Data volume to be treated in beam migration process is synthesized, the offset efficiency of seismic migration imaging is substantially increased.
Detailed description of the invention
It in order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, below will be to institute in embodiment Attached drawing to be used is needed to be briefly described, it should be apparent that, the accompanying drawings in the following description is only some implementations of the invention Example, for those of ordinary skill in the art, without any creative labor, can also be according to these attached drawings Obtain other attached drawings.
Fig. 1 is the method flow diagram of the bilateral synthesis beam offset imaging method provided by the invention picked up based on lineups；
Fig. 2 is the schematic diagram that big gun inspection binning is carried out using second grade segmentation method provided by the invention；
Fig. 3 is that the big gun provided by the invention picked up based on lineups examines the bilateral synthetic method stream of domain Local plane wave data Cheng Tu；
Fig. 4 is the method that the big gun provided by the invention picked up based on lineups examines bilateral Local plane wave data-bias imaging Flow chart；
Fig. 5 is the system construction drawing of the bilateral synthesis beam migration imaging system provided by the invention picked up based on lineups；
Fig. 6 is Marmousi rate pattern schematic diagram provided in an embodiment of the present invention；
Fig. 7 is the schematic diagram of synthetic seismic record data provided in an embodiment of the present invention；
Fig. 8 is the Local plane wave schematic diagram data before progress lineups pickup provided in an embodiment of the present invention；
Fig. 9 is the Local plane wave schematic diagram data after progress lineups pickup provided in an embodiment of the present invention；
Figure 10 is Marmousi model composite traces provided in an embodiment of the present invention quickly bilateral synthesis beam offset method The schematic diagram of imaging results.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other Embodiment shall fall within the protection scope of the present invention.
The object of the present invention is to provide it is a kind of based on lineups pick up bilateral synthesis beam offset imaging method and system, with Solve the problems, such as that existing bilateral synthesis beam offset imaging method data processing amount is big, operation efficiency is low.
In order to make the foregoing objectives, features and advantages of the present invention clearer and more comprehensible, with reference to the accompanying drawing and specific real Applying mode, the present invention is described in further detail.
Unlike traditional Gaussian beam or gram western hoeve's bundle offset, bilateral synthesis beam offset method is simultaneously in big gun The synthesis that the bilateral face element in position carries out plane wave is examined, the data volume of bilateral synthesis is larger.For three dimensional seismic data, It is transformed into quintuple space (shot point side is three-dimensional, and geophone station side is three-dimensional, shares time dimension) firstly the need of three dimensional seismic data, then Binning and Local plane wave synthesis (having more bidimensional) are carried out along plane to shot point side three-dimensional, then again to geophone station side three-dimensional Binning and Local plane wave synthesis (and having more bidimensional) are carried out along plane, finally carries out migration imaging again.This process meeting There is a large amount of data redundancy, causes operation efficiency very low.
In order to solve this technical problem in the prior art, the invention proposes a kind of based on the quick of lineups pickup Bilateral synthesis beam deviates (Fast event-picking-based double-sided beam migration, FDBM) method And system stores effective plane wave data, realizes fast excursion in practical applications by lineups pick-up method.This hair The bilateral synthesis beam migration imaging formula of bright offer picked up based on lineups is expressed as follows:
Wherein,It is imaging point coordinate,WithIt is shot point and geophone station face element centre coordinate respectively,WithIt is seismic trace coordinates in shot point and geophone station face element respectively, when t is the longitudinal coordinate travelling of earthquake record；I is initial seismic image Variable, AiIt is the corresponding amplitude weight of i-th of face element, a shared n face element, WjIt is the corresponding weighted value in the road face element Zhongj, it should Face element one shares the road m；D is the seismic channel data for including in face element, and τ is the intercept of plane wave conversion,It is shot point respectively The plane wave direction pointed out with detection,Be respectively in face element seismic channel to face element big gun center and detection dot center Distance,Indicate the valid data that five dimensional plane wave spaces are picked up according to data attribute.The bilateral synthesis beam is inclined It moves imaging formula (1) and further convert from general seismic imaging formula and obtain, mainly joined the present invention " based on picking up together The core concept of phase axis ", byExpression, the inside D are the Local plane wave expression formulas of higher-dimension.
Fig. 1 is the method flow diagram of the bilateral synthesis beam offset imaging method provided by the invention picked up based on lineups. Referring to Fig. 1, the bilateral synthesis beam offset imaging method provided by the invention picked up based on lineups is specifically included:
Step 101: obtaining earthquake record data and external definition parameter.
The source of earthquake record data is the seismic data or laboratory model analog synthesis of field earthquake prospecting collection Seismic data.The earthquake record data include time sampling number, time sampling interval, frequency range and the big gun of seismic data Examine coordinate range.The external definition parameter includes face element defined parameters and speed defined parameters: the face element defined parameters packet Effective lineups pick up number when including bin size, Local plane wave orientation angle interval and number and output.The speed Defined parameters include earth's surface speed and average speed.
Step 102: according to the earthquake record data and the face element defined parameters, big gun being carried out using second grade segmentation method Binning is examined, multiple big gun inspection face elements are obtained.The big gun inspection face element includes shot point face element and geophone station face element.
Determined according to big gun inspection coordinate range, bin size, plane wave orientation angle interval and the number in earthquake record data The sample space of adopted part plan wave number, and big gun is carried out by second grade segmentation method and examines domain binning: the first order is according to face element Size divides rectangular mesh, and the second level is deviated according to earthquake data frequency range, average speed size using general Gaussian beam Parameter equation, which calculates, assists elliptical radius, divides second level face element by the way that auxiliary is oval.It divides first order face element (rectangle face element) After be arranged successively to obtain big gun inspection domain face element sequence number, facilitate earthquake record data to establish corresponding face element index.
Fig. 2 is the schematic diagram that big gun inspection binning is carried out using second grade segmentation method provided by the invention.S is indicated in Fig. 2 Shot point, G indicate that geophone station, I indicate imaging point；Rectangle frame indicates the segmentation of level-one face element, and oval frame is the face element section of auxiliary, claims For second level face element divide.
Three-dimensional seismic acquisition data are actually five dimension datas, and big gun inspection regional partition face element needs carry out on high dimensional plane, It may further be decomposed into from formula and carry out Gauss segmentation, Gaussian function numerical expression in the two-dimensional surface of shot point and the two-dimensional surface of geophone station It is as follows:
Wherein, (σx, σy) it is the horizontal width of Gauss clock in the x and y direction,It isOr PersonThe centre coordinate of face element, Gauss clock show as ellipse in the horizontal plane.Upper face element segmentation is such as Fig. 2 institute Show, second grade segmentation is illustrated in figure, i.e., carry out rectangle segmentation first, then carries out being limited with oval, be obtained from neighbouring face element Part seismic channel data.Big gun can be examined the rectangular mesh that plane space is divided into rule respectively according to the bin size of definition, As soon as all corresponding shot point rectangle face element of each track data of earthquake record data and geophone station rectangle face element, the two shape Face element is examined at a level-one big gun.Auxiliary is oval mainly using the seismic channel in the distance of major and minor axis and neighbouring face element into the face element The distance of heart position compares, and judges whether seismic channel has second level face element.
The step 102 according to the earthquake record data and the face element defined parameters, using second grade segmentation method into Row big gun examines binning, obtains multiple big gun inspection face elements, specifically includes:
Coordinate range, the bin size, plane wave orientation angle interval and number, which are examined, according to the big gun determines office The big gun of facial planes wave number examines plane sample space；
Big gun inspection plane sample space is divided into the rectangular mesh of rule according to the bin size, generates level-one big gun Examine face element；The level-one big gun inspection face element includes a shot point rectangle face element and a geophone station rectangle face element；
According to the frequency range of seismic data and the average speed, calculated using Gaussian beam offset parameter formula auxiliary Help elliptical radius；
Determine that auxiliary is oval according to the elliptical radius of auxiliary；
The big gun inspection face element is marked off according to the auxiliary ellipse；The big gun inspection face element includes shot point face element and geophone station face Member.
Step 103: being indexed according to the face element that the earthquake record data establish multiple big gun inspection face elements.
Earthquake record data are established according to the corresponding shot point face element of earthquake record data inspection point coordinate and geophone station face element Corresponding face element index.All seismic channels of earthquake record data include that face element call number, hereof position and big gun examine coordinate The structural body (present invention is referred to as index structure body) of information retains and uses when calculating in memory for Local plane wave.
Step 104: the pickup of energy spectrum attribute and significant figure are carried out to multiple big gun inspection face elements according to face element index According to storage, multiple Local plane wave compressed data blocks are generated.It include plane wave in each Local plane wave compressed data block Header file and data file.
Before the step 104, need first to examine face element Local plane wave calculating preparation memory and compressor for single big gun. Memory required for Local plane wave calculates is opened up according to the max number of channels that single big gun examines the seismic data that face element is included.It is described Compressor is defined by the inclusion of the index structure body of the spatial position of pickup lineups, time location and amplitude size, institute Compressor is stated to compress data using the threshold compression algorithm damaged.
The extraction of the effective lineups of Local plane wave, can method by linearly pulling transformation or Least squares inversion Obtain Local plane wave data.Either which kind of method, obtaining Local plane wave data has noise.If recording all parts The information of plane-wave field, data volume is huger, also unbearable in production even if being stored in disk.If directly according to vibration Width threshold value is extracted, and data volume reduction still not satisfies requirement, therefore the present invention passes through Correlated Spectroscopy formula (3) or other attributes extractions Effective part plan wave field information.Correlated Spectroscopy pickup can be obtained by the focusing of spectrum energy than sparse effective same phase Axis energy improves decrement, reduces the memory space of Local plane wave, mainly describes the effective lineups of Local plane wave here Pick process method, energy spectrum is not uniquely to select.
The Correlated Spectroscopy calculating formula that the present invention uses is as follows:
Wherein,It is Local plane wave dataCorrelated Spectroscopy, here n be Refer to the sampling point number of Local plane wave data in sliding window function.The Local plane wave defined by parameter list pick up block number andEffective part plan data are picked up, shown in following formula:
Wherein,It is effective Local plane wave data,It is above by Correlated Spectroscopy Pick up the process of effective Local plane wave data.
Fig. 3 is that the big gun provided by the invention picked up based on lineups examines the bilateral synthetic method stream of domain Local plane wave data Cheng Tu.Referring to Fig. 3, step 104 provided by the invention indexes according to the face element and carries out energy spectrum to multiple big gun inspection face elements Attribute picks up and valid data storage, generates multiple Local plane wave compressed data blocks, specifically includes:
The seismic channel for being included in step S31: i-th of the search big gun inspection face element.
Step S31 is mainly that subsequent step provides the input data that Local plane wave calculates, and is stored in from step 103 interior All seismic channel relevant informations of i-th of big gun inspection face element are extracted in the index structure body deposited, comprising: position hereof It sets, geophone offset coordinate information and the range information from face element center, and further positions and read from earthquake record data file Relative earthquake track data.
Step S32: the seismic channel data of the seismic channel is read according to face element index.
Step S33: according to the seismic channel data, using the office of big gun inspection face element described in i-th of transformation calculations of linear pulling Facial planes wave number evidence.
The Local plane wave data of face element are examined using the single big gun of linear pulling transformation calculations according to the seismic channel data, And energy spectrum is obtained using the auto-correlation formula of Local plane wave data and picks up attribute.
Step S 34: using the auto-correlation formula of Local plane wave dataThe Correlated Spectroscopy of Local plane wave data is calculated.
Step S35: number and the Correlated Spectroscopy are picked up according to effective lineups, using formulaPick up effective Local plane wave data.
Step S36: effective same phase in effective Local plane wave data is determined according to the amplitude power of the Correlated Spectroscopy The location and range of axis.
The vibration of the Local plane wave energy spectrum of number and step S4 calculating is picked up according to effective lineups that step 101 obtains Width power determines the location and range of effective lineups in Local plane wave data.Based on this location and range, step 104 is used Ready compressor compresses and store Local plane wave header file and data file before.The header file data include Big gun examines coordinate information, the location information of affiliated file and the Local plane wave data storage size of generation of face element, and has There are an effective lineups, the present invention is referred to as a compression blocks of Local plane wave data, i.e. Local plane wave data compression Block.
Step S37: according to the location and range of effective lineups, using compressor compresses and Local plane wave is stored Header file and data file generate i-th of Local plane wave compressed data block.The Local plane wave data head text Part include big gun inspection the coordinate information of face element, the location information of affiliated file, generation Local plane wave data storage size and one A effective lineups.
Step S38: judging whether i is less than n, obtains the first judging result.N is the quantity that the big gun examines face element, 1≤i < n.
Judge whether big gun inspection face element number calculated reaches the maximum value of the inspection domain face element sequence number of big gun obtained in step 103 N, if not having (i < n), return step S31 enters next big gun inspection face element and calculates；Otherwise (i >=n) moves on to step S39.
Step S39: output the 1st to the n Local plane wave compressed data blocks.
Step S310: header file and data file to the Local plane wave compression blocks of all outputs carry out correctness inspection It looks into, and output error big gun inspection face element number provides quality control information for user.Test mode are as follows: re-read Local plane wave Header file extracts relevant information from Local plane wave compressed data block according to header file information, and whether see can correctly read, It is otherwise exactly the output of mistake if the Local plane wave data for being judged as correctly exporting can correctly be read.
Step 105: the offset distance of each Local plane wave is calculated according to the Local plane wave compressed data block.
Before the step 105, first have to obtain Local plane wave compressed data block and external definition parameter.The office Facial planes wave compressed data block is obtained from step 104 Local plane wave Data Synthesis result.The external definition parameter is main It is conventional integration method migration imaging parameter, rate pattern parameter and ray tracing parameter.Step 106 mainly uses outside First offset distance size H0, offset distance gap size DH and offset number NH in the offset gather parameter information of definition.
Between the time sampling for reading original seismic data data according to the header file of the Local plane wave compressed data block Every, time sampling number, while counting according to the data file of Local plane wave compressed data block the number of part plan wave number evidence With the offset distance H1 of each Local plane wave, and each trace gather is calculated according to the offset gather parameter information of external definition Number (1,2,3 ...) corresponding offset distance size H2.
Step 106: it is locally flat that the corresponding trace gather of the Local plane wave being established according to the offset distance of the Local plane wave Surface wave call number.
The corresponding offset distance size H2 of each offset gather number and Local plane wave obtained according to step 105 is corresponding Offset distance size H1 compares the two and can establish the corresponding trace gather Local plane wave call number (Hi=(H1- of Local plane wave H0)/DH+1)。
Specifically, the offset distance according to the Local plane wave establishes the corresponding trace gather part of the Local plane wave Plane wave call number, specifically includes:
According to the offset distance of the Local plane wave, the Local plane wave is established using formula Hi=(H1-H0)/DH+1 Corresponding trace gather Local plane wave call number；Wherein Hi indicates the corresponding trace gather Local plane wave index of i-th of Local plane wave Number；H1 indicates the offset distance of the Local plane wave；H0 indicates first offset distance in offset gather parameter information；DH table Show offset distance interval.
Step 107: institute is completed according to the trace gather Local plane wave call number and the Local plane wave compressed data block There is the migration imaging of all Local plane wave compressed data blocks in imaging trace gather.
Fig. 4 is the method that the big gun provided by the invention picked up based on lineups examines bilateral Local plane wave data-bias imaging Flow chart.The Local plane wave data of storage are compressed after big gun examines bilateral part plan wave conversion and effective lineups pick up, Referred to as effective Local plane wave data, the lineups of pickup are referred to as Local plane wave valid data compression blocks.The present invention Beam shifted versions that beam offset based on Local plane wave valid data compression blocks is deviated with previous Gaussian beam and it simplifies are not It is with place, input data of the invention is Local plane wave valid data compression blocks, each compressed data block contains Directional information, temporal information and amplitude information independently can be deviated once.It is flat that process shown in Fig. 4 is namely based on part The beams of surface wave valid data compression blocks deviates process, including outer circulation of the outer layer based on offset gather number and internal layer with every The interior circulation for the valid data compression block number that one offset distance includes, the as nested double-deck circulation.The fineness ratio of outer loop Larger, the granularity of task parallelism easy to accomplish, interior loop is moderate, thread loops easy to accomplish.In this way, the base that the present invention uses It can be realized quickly and efficiently operation in the beam offset of Local plane wave valid data compression blocks.
Referring to fig. 4, the step 107 is according to the trace gather Local plane wave call number and the Local plane wave data pressure Contracting block completes the migration imaging of all Local plane wave compressed data blocks in all imaging trace gathers, specifically includes:
Step S41: i-th of Local plane wave compressed data block in k-th of trace gather is obtained.
According to the trace gather Local plane wave call number established in step 103, office all in single offset gather is searched for Facial planes wave number evidence, then count each Local plane wave data includes how many a compression blocks (effective lineups data).To Obtain Local plane wave compressed data block all in K trace gather.
Step S42: the spatial position of the effective lineups stored in i-th of Local plane wave compressed data block is obtained And time location.
According to the storage information in the compression blocks counted in step S41 obtain effective lineups spatial position (direction) and Time location.
Step S43: the ray tracing of big gun inspection position is carried out according to the spatial position and the time location, and to intersection Position carries out migration imaging, generates imaging temporary file.
According to the direction of effective lineups and temporal information, carry out the ray tracing of big gun inspection position, and to crossover location into Row imaging.
Step S44: judging whether i is less than L, obtains the second judging result；Wherein L is locally to put down described in k-th of trace gather The quantity of surface wave compressed data block.
Judge whether that all compression blocks for completing the related all Local plane wave data of k-th of offset gather are (effective Lineups data) migration imaging, if not, enable i=i+1, return step S32 obtains i-th of Local plane wave data The step of spatial position of the effective lineups stored in compression blocks and time location, into the circulation of next compression blocks；Such as Fruit is, into next step S35.
Step S45: judging whether k is less than K, obtains third judging result；Wherein K is the quantity that trace gather is imaged.
Judge whether to complete all migration imagings that trace gathers are imaged, if not, k=k+1 is enabled, i=1, return step S41 The step of obtaining i-th of Local plane wave compressed data block in k-th of trace gather, into next offset gather migration imaging Circulation；If so, the imaging temporary file that this migration imaging obtains is exported, into next step S46.
Step S46: merging the 1st imaging temporary file for arriving K trace gather, generates imaging trace gather file.
All single offset distance imaging temporary files that step S45 is obtained are merged into final imaging results: if pressed Offset distance rearrangement then exports imaging trace gather file；If all temporary offsets are superimposed away from imaging data, obtain Imaging section file.
As it can be seen that the bilateral synthesis beam offset imaging method provided by the invention picked up based on lineups, using second grade segmentation Method (rectangular mesh+auxiliary is oval) divides face element, solves the problems, such as seismic channel quick indexing in face element.Existing beam migration technology It is all deviated from Gaussian beam, the segmentation of space face element is all overlapping, and usual Gaussian beam offset is all to realize process in algorithm It is middle to realize that the face element of overlapping is divided in moving process with a Gaussian window, it is moved to a new face element position, is then searched for Included in all seismic channels here.And the present invention carries out binning (rectangle to fast implement beam migration technology in advance Division+auxiliary is oval), rectangle binning is mutually indepedent, can establish the unique face element index of all seismic channels, weight originally Folded part is obtained from neighbouring face element with auxiliary ellipse.Therefore the second level face element dividing method that the present invention uses being capable of quick rope Draw the seismic channel in face element.
The present invention also uses lineups pick-up method to solve higher-dimension (five dimensions) plane wave conversion data having in Shu Hecheng Effect expression, effectively storage and fast excursion problem.The earthquake beam data of five dimensions divides face element, each face element transformation in the planes To plane wave space, become nine dimensions, data volume is especially huge, bilateral synthesis beam offset is not available in actual production, also Say it is not to be available bilateral synthesis beam to carry out beam offset, and the present invention is picked up using lineups, is exactly in the prior art It is sparse feature using Local plane wave data, picks up effective data information, stored again by compression, then utilizes these The lineups of effective small number carry out fast excursion again, realize effective table of the plane wave conversion data in Shu Hecheng Reach, effectively storage and fast excursion.
The invention also discloses a kind of bilateral synthesis beam migration imaging systems picked up based on lineups.Fig. 5 is the present invention The system construction drawing of the bilateral synthesis beam migration imaging system picked up in lineups provided.Referring to Fig. 5, the system comprises:
Data and parameter acquisition module 501, for obtaining earthquake record data and external definition parameter；The earthquake record Data include time sampling number, time sampling interval, frequency range and the big gun inspection coordinate range of seismic data；The external definition Parameter includes face element defined parameters and speed defined parameters: the face element defined parameters include bin size, Local plane wave side Effective lineups pick up number when to angle interval and number and exporting；The speed defined parameters include earth's surface speed peace Equal speed；
Big gun examines binning module 502, is used for according to the earthquake record data and the face element defined parameters, using two Grade dividing method carries out big gun and examines binning, obtains multiple big gun inspection face elements；The big gun inspection face element includes shot point face element and geophone station Face element；
Face element index establishes module 503, for establishing the face of multiple big gun inspection face elements according to the earthquake record data Member index；
Local plane wave compressed data block generation module 504, for examining face to multiple big guns according to face element index Member carries out the pickup of energy spectrum attribute and valid data storage, generates multiple Local plane wave compressed data blocks；Each part It include plane wave data header file and data file in plane wave data compression blocks；
Offset distance computing module 505, for calculating each Local plane wave according to the Local plane wave compressed data block Offset distance；
Trace gather Local plane wave call number establishes module 506, for establishing institute according to the offset distance of the Local plane wave State the corresponding trace gather Local plane wave call number of Local plane wave；
Migration imaging module 507, for according to the trace gather Local plane wave call number and the Local plane wave data Compression blocks complete the migration imaging of all Local plane wave compressed data blocks in all imaging trace gathers.
The big gun inspection binning module 502 specifically includes:
Big gun examines plane sample space determination unit, for examining the coordinate range, bin size, described flat according to the big gun Surface wave orientation angle interval and number determine the big gun inspection plane sample space of part plan wave number；
Level-one big gun examines face element generation unit, for being divided into big gun inspection plane sample space according to the bin size The rectangular mesh of rule generates level-one big gun and examines face element；The level-one big gun inspection face element includes a shot point rectangle face element and an inspection Wave point rectangle face element；
Assist ellipse radii computing unit, for according to seismic data the frequency range and the average speed, adopt It is calculated with Gaussian beam offset parameter formula and assists elliptical radius；
Oval determination unit is assisted, for determining that auxiliary is oval according to the elliptical radius of auxiliary；
Big gun examines binning unit, examines face element for marking off the big gun according to the auxiliary ellipse；The big gun examines face element Including shot point face element and geophone station face element.
The Local plane wave compressed data block generation module 504 specifically includes:
Seismic channel search unit, for searching for the seismic channel for being included in i-th of big gun inspection face element；
Seismic channel data reading unit, for reading the seismic channel data of the seismic channel according to face element index；
Local plane wave Data Computation Unit, for pulling transformation calculations i-th using linear according to the seismic channel data The Local plane wave data of a big gun inspection face element；
Correlated Spectroscopy computing unit, for the auto-correlation formula using Local plane wave dataThe Correlated Spectroscopy of Local plane wave data is calculated；
Effective Local plane wave data pick-up unit, for according to effective lineups pick up number to it is described related Spectrum, using formulaPick up effective Local plane wave data；
Effective lineups position range determination unit, determines effective office for the amplitude power according to the Correlated Spectroscopy The location and range of facial planes wave number effective lineups in；
Local plane wave compressed data block generation unit is used for the location and range according to effective lineups Compressor compresses simultaneously store Local plane wave header file and data file, generate i-th of Local plane wave data pressure Contracting block；The Local plane wave header file includes the coordinate information of big gun inspection face element, the location information of affiliated file, generation Local plane wave data storage size and effective lineups；
First judging unit obtains the first judging result for judging whether i is less than n；N is the number that the big gun examines face element Amount, 1≤i < n；
First return unit, for enabling i=i+1, returning to the seismic channel when first judging result is that i is less than n Search unit；
Local plane wave compressed data block output unit, for exporting when first judging result is that i is not less than n 1st to the n Local plane wave compressed data block.
The trace gather Local plane wave call number is established module 506 and is specifically included:
Trace gather Local plane wave call number establishes unit, for the offset distance according to the Local plane wave, using formula Hi=(H1-H0)/DH+1 establishes the corresponding trace gather Local plane wave call number of the Local plane wave；Wherein Hi is indicated i-th The corresponding trace gather Local plane wave call number of Local plane wave；H1 indicates the offset distance of the Local plane wave；H0 indicates offset Away from first offset distance in trace gather parameter information；DH indicates offset distance interval.
The migration imaging module 507 specifically includes:
Local plane wave compressed data block acquiring unit, for obtaining i-th of Local plane wave data in k-th of trace gather Compression blocks；
Effective lineups position acquisition unit stores in i-th of Local plane wave compressed data block for obtaining The spatial position of effective lineups and time location；
Migration imaging unit, the ray for carrying out big gun inspection position according to the spatial position and the time location chase after Track, and migration imaging is carried out to crossover location, generate imaging temporary file；
Second judgment unit obtains the second judging result for judging whether i is less than L；Wherein L is institute in k-th of trace gather State the quantity of Local plane wave compressed data block；
Second return unit returns described effectively same for enabling i=i+1 when second judging result is that i is less than L Phase shaft position acquiring unit；
Third judging unit obtains for judging whether k is less than K when second judging result is that i is not less than L Three judging results；Wherein K is the quantity that trace gather is imaged；
Third return unit, for enabling k=k+1, it is flat returning to the part when the third judging result is that k is less than K Surface wave compressed data block acquiring unit；
Trace gather file generating unit is imaged, for merging the 1st to K when the third judging result is that k is not less than K The imaging temporary file of trace gather generates imaging trace gather file.
Below by Marmousi model test, the present invention is based on the bilateral synthesis beam offset imaging methods that lineups pick up And the technical effect of system.2.5D Marmousi model adapts to 3D algoritic module.
Fig. 6 is Marmousi rate pattern schematic diagram provided in an embodiment of the present invention.Lateral coordinates indicate horizontal blank in Fig. 6 Between the direction x coordinate, longitudinal to indicate the direction depth z coordinate, vertically and horizontally unit is all rice (m).Referring to Fig. 6, Marmousi speed mould The sampling grid of type is 3x369x750, respectively indicates the direction y, the sampling number in the direction x and the direction z, and the sampling interval is 25mx25mx4m。
Fig. 7 is the schematic diagram of synthetic seismic record data provided in an embodiment of the present invention.Lateral coordinates indicate focus in Fig. 7 Number (dimensionless), longitudinal coordinate indicate to receive the time coordinate of earthquake record, and unit is all millisecond (ms).Referring to Fig. 7, earthquake is closed At record data (big gun number 119,120 and 121), longitudinal time sampling interval is 4ms, and 240 big gun, inspection point are all located at second altogether Bar line (y=25m), first shot position are (25m, 3000m), and last shot position is (25m, 8975m), are spaced 25m；It is right While blowing out, 96 wave detectors of every big gun are received, sampling interval 25m.
Local plane wave synthetic parameters, bin size 200mx200m, inclination maximum are 88 degree, and 12 degree of intervals, three dimensional angular is adopted It is all 225 that it is identical, which to examine side, for sample big gun, i.e., the Local plane wave original angle sampling sum of one face element is 225*255.According to journey Sequence programming count, effective face element number are 402.If the whole Local plane wave data of output need 61.4G completely, and this hair The method and system of bright offer are picked up by lineups, and compression storage only needs 1G~150M (lineups number parameter 5000- 10000) data volume that, can probably reduce 60~400 times, greatly reduces redundant data.
Fig. 8 is the Local plane wave schematic diagram data before progress lineups pickup provided in an embodiment of the present invention.Fig. 9 is this The Local plane wave schematic diagram data carried out after lineups pickup that inventive embodiments provide.Fig. 8 and Fig. 9 is shown the 19th The part plan wave conversion data of face element, take shot point side azimuth and inclination angle is all 0, the plane wave that geophone station side azimuth is 0 Data, abscissa are inclination angle (unit Deg), and ordinate is time (unit ms).Big gun side inclination angle and azimuth be all in Fig. 8 and Fig. 9 It is 0, that is, takes 113 position of intermediate point of 225 points in big gun side, then geophone station side azimuth is 0.Part referring to Fig. 8, before pickup Plane wave data out-focus, including linearly pull energy leakage caused by algorithm and some noises.And the part after Fig. 9 is picked up is flat Surface wave data eliminate the data and noise of energy leakage part, are capable of the signal-to-noise ratio of final migration imaging result.Base of the present invention The compression ratio of data can be improved by the pickup of effective lineups in the bilateral beam synthetic method that lineups pick up, solved double Side synthesizes the problem of huge data output, and improving one kind for practical application can apply and effective method.
Synthetic seismic record shown in Fig. 7 is linearly cut off by direct wave, then using the present invention is based on lineups pickups Quick bilateral synthesis beam offset method obtain imaging results.Figure 10 is Marmousi model provided in an embodiment of the present invention synthesis Record the schematic diagram of the quickly imaging results of bilateral synthesis beam offset method.Lateral coordinates indicate the direction horizontal space x in Figure 10 Coordinate, longitudinal to indicate the direction depth z coordinate, vertically and horizontally unit is all rice (m).Figure 10 middle-shallow layer top missing is because direct wave is thick Cut through the data that journey has lost part superficial reflex.The result of Figure 10 illustrates the bilateral beam offset picked up the present invention is based on lineups Method can use in practice.As said before, by the definition of different lineups numbers, 5000-10000 is tested here The range of lineups data, so that it may can probably reduce by 60~400 times of Local plane wave data volume, so so big pressure Demagnification number is enough quickly to handle the bigger data of this method.
Each embodiment in this specification is described in a progressive manner, the highlights of each of the examples are with other The difference of embodiment, the same or similar parts in each embodiment may refer to each other.For system disclosed in embodiment For, since it is corresponded to the methods disclosed in the examples, so being described relatively simple, related place is said referring to method part It is bright.
Used herein a specific example illustrates the principle and implementation of the invention, and above embodiments are said It is bright to be merely used to help understand method and its core concept of the invention；At the same time, for those skilled in the art, foundation Thought of the invention, there will be changes in the specific implementation manner and application range.In conclusion the content of the present specification is not It is interpreted as limitation of the present invention.

Claims (10)

1. the bilateral synthesis beam offset imaging method picked up based on lineups, which is characterized in that the described method includes:
Obtain earthquake record data and external definition parameter；The earthquake record data include seismic data time sampling number, Time sampling interval, frequency range and big gun examine coordinate range；The external definition parameter includes that face element defined parameters and speed are fixed Adopted parameter: the face element defined parameters have when including bin size, Local plane wave orientation angle interval and number and output It imitates lineups and picks up number；The speed defined parameters include earth's surface speed and average speed；
According to the earthquake record data and the face element defined parameters, big gun is carried out using second grade segmentation method and examines binning, Obtain multiple big gun inspection face elements；The big gun inspection face element includes shot point face element and geophone station face element；
The face element index of multiple big gun inspection face elements is established according to the earthquake record data；
The pickup of energy spectrum attribute is carried out to multiple big gun inspection face elements according to face element index and valid data store, is generated more A Local plane wave compressed data block；It include plane wave data header file sum number in each Local plane wave compressed data block According to file；
The offset distance of each Local plane wave is calculated according to the Local plane wave compressed data block；
The corresponding trace gather Local plane wave call number of the Local plane wave is established according to the offset distance of the Local plane wave；
It is completed in all imaging trace gathers according to the trace gather Local plane wave call number and the Local plane wave compressed data block The migration imaging of all Local plane wave compressed data blocks.
2. bilateral synthesis beam offset imaging method according to claim 1, which is characterized in that described to be remembered according to the earthquake Data and the face element defined parameters are recorded, big gun is carried out using second grade segmentation method and examines binning, obtain multiple big gun inspection face elements, tool Body includes:
Coordinate range, the bin size, plane wave orientation angle interval and number, which are examined, according to the big gun determines that part is flat The big gun of surface wave number examines plane sample space；
Big gun inspection plane sample space is divided into the rectangular mesh of rule according to the bin size, level-one big gun is generated and examines face Member；The level-one big gun inspection face element includes a shot point rectangle face element and a geophone station rectangle face element；
According to the frequency range of seismic data and the average speed, it is ellipse that auxiliary is calculated using Gaussian beam offset parameter formula Round radius；
Determine that auxiliary is oval according to the elliptical radius of auxiliary；
The big gun inspection face element is marked off according to the auxiliary ellipse；The big gun inspection face element includes shot point face element and geophone station face element.
3. bilateral synthesis beam offset imaging method according to claim 2, which is characterized in that described according to the face element rope Draw and the pickup of energy spectrum attribute and valid data storage are carried out to multiple big gun inspection face elements, generates multiple Local plane wave data pressures Contracting block, specifically includes:
Search for the seismic channel for being included in i-th of big gun inspection face element；
The seismic channel data of the seismic channel is read according to face element index；
According to the seismic channel data, using the Local plane wave data of big gun inspection face element described in i-th of transformation calculations of linear pulling；
Using the auto-correlation formula of Local plane wave data The Correlated Spectroscopy of Local plane wave data is calculated；
Number and the Correlated Spectroscopy are picked up according to effective lineups, using formulaPick up effective Local plane wave data；
The position of effective lineups and model in effective Local plane wave data are determined according to the amplitude power of the Correlated Spectroscopy It encloses；
According to the location and range of effective lineups, using compressor compresses and store Local plane wave header file and Data file generates i-th of Local plane wave compressed data block；The Local plane wave header file includes big gun inspection face The coordinate information of member, the location information of affiliated file, the Local plane wave data storage size of generation and one are described effectively same Phase axis；
Judge whether i is less than n, obtains the first judging result；N is the quantity that the big gun examines face element, 1≤i < n；
If first judging result is that i is less than n, i=i+1 is enabled, is included in big gun inspection face element described in i-th of described search of return Seismic channel the step of；
If first judging result is that i is not less than n, output the 1st to the n Local plane wave compressed data blocks.
4. bilateral synthesis beam offset imaging method according to claim 3, which is characterized in that described flat according to the part The offset distance of surface wave establishes the corresponding trace gather Local plane wave call number of the Local plane wave, specifically includes:
According to the offset distance of the Local plane wave, it is corresponding that the Local plane wave is established using formula Hi=(H1-H0)/DH+1 Trace gather Local plane wave call number；Wherein Hi indicates the corresponding trace gather Local plane wave call number of i-th of Local plane wave；H1 Indicate the offset distance of the Local plane wave；H0 indicates first offset distance in offset gather parameter information；DH indicates inclined It moves away from interval.
5. bilateral synthesis beam offset imaging method according to claim 4, which is characterized in that described according to the trace gather office Facial planes wave call number and the Local plane wave compressed data block complete all Local plane wave data in all imaging trace gathers The migration imaging of compression blocks, specifically includes:
Obtain i-th of Local plane wave compressed data block in k-th of trace gather；
Obtain spatial position and the time location of the effective lineups stored in i-th of Local plane wave compressed data block；
The ray tracing of big gun inspection position is carried out according to the spatial position and the time location, and crossover location is deviated Imaging generates imaging temporary file；
Judge whether i is less than L, obtains the second judging result；Wherein L is Local plane wave data compression described in k-th of trace gather The quantity of block；
If second judging result is that i is less than L, i=i+1 is enabled, returns to acquisition i-th of Local plane wave data pressure The step of spatial position of the effective lineups stored in contracting block and time location；
If second judging result is that i is not less than L, judge whether k is less than K, obtains third judging result；Wherein K is imaging The quantity of trace gather；
If the third judging result is that k is less than K, k=k+1 is enabled, it is flat to return to i-th of the part obtained in k-th of trace gather The step of surface wave compressed data block；
If the third judging result is that k is not less than K, merge the 1st imaging temporary file for arriving K trace gather, generates imaging trace gather File.
6. the bilateral synthesis beam migration imaging system picked up based on lineups, which is characterized in that the system comprises:
Data and parameter acquisition module, for obtaining earthquake record data and external definition parameter；The earthquake record data packet Include time sampling number, time sampling interval, frequency range and the big gun inspection coordinate range of seismic data；The external definition parameter packet Include face element defined parameters and speed defined parameters: the face element defined parameters include bin size, Local plane wave orientation angle Effective lineups pick up number when interval and number and output；The speed defined parameters include earth's surface speed peace average rate Degree；
Big gun examines binning module, is used for according to the earthquake record data and the face element defined parameters, using second grade segmentation Method carries out big gun and examines binning, obtains multiple big gun inspection face elements；The big gun inspection face element includes shot point face element and geophone station face element；
Face element index establishes module, and the face element for establishing multiple big gun inspection face elements according to the earthquake record data indexes；
Local plane wave compressed data block generation module, for carrying out energy to multiple big gun inspection face elements according to face element index Amount spectrum attribute picks up and valid data storage, generates multiple Local plane wave compressed data blocks；Each part plan wave number According in compression blocks include plane wave data header file and data file；
Offset distance computing module, for calculating the offset of each Local plane wave according to the Local plane wave compressed data block Away from；
Trace gather Local plane wave call number establishes module, flat for establishing the part according to the offset distance of the Local plane wave The corresponding trace gather Local plane wave call number of surface wave；
Migration imaging module, for complete according to the trace gather Local plane wave call number and the Local plane wave compressed data block At the migration imaging of all Local plane wave compressed data blocks in all imaging trace gathers.
7. bilateral synthesis beam migration imaging system according to claim 6, which is characterized in that the big gun examines binning mould Block specifically includes:
Big gun examines plane sample space determination unit, for examining coordinate range, the bin size, the plane wave according to the big gun Orientation angle interval and number determine the big gun inspection plane sample space of part plan wave number；
Level-one big gun examines face element generation unit, for big gun inspection plane sample space to be divided into rule according to the bin size Rectangular mesh, generate level-one big gun examine face element；The level-one big gun inspection face element includes a shot point rectangle face element and a geophone station Rectangle face element；
Assist ellipse radii computing unit, for according to seismic data the frequency range and the average speed, using height This beam offset parameter formula, which calculates, assists elliptical radius；
Oval determination unit is assisted, for determining that auxiliary is oval according to the elliptical radius of auxiliary；
Big gun examines binning unit, examines face element for marking off the big gun according to the auxiliary ellipse；The big gun examines face element Shot point face element and geophone station face element.
8. bilateral synthesis beam migration imaging system according to claim 7, which is characterized in that the Local plane wave data Compression blocks generation module specifically includes:
Seismic channel search unit, for searching for the seismic channel for being included in i-th of big gun inspection face element；
Seismic channel data reading unit, for reading the seismic channel data of the seismic channel according to face element index；
Local plane wave Data Computation Unit, for pulling i-th of institute of transformation calculations using linear according to the seismic channel data State the Local plane wave data of big gun inspection face element；
Correlated Spectroscopy computing unit, for the auto-correlation formula using Local plane wave dataThe Correlated Spectroscopy of Local plane wave data is calculated；
Effective Local plane wave data pick-up unit is adopted for picking up number and the Correlated Spectroscopy according to effective lineups Use formulaPick up effective Local plane wave data；
Effective lineups position range determination unit determines that effective part is flat for the amplitude power according to the Correlated Spectroscopy The location and range of effective lineups in surface wave data；
Local plane wave compressed data block generation unit, for the location and range according to effective lineups, using compression Device compresses and stores Local plane wave header file and data file, generates i-th of Local plane wave compressed data block； The Local plane wave header file include big gun inspection the coordinate information of face element, the location information of affiliated file, generation part Plane wave data storage size and effective lineups；
First judging unit obtains the first judging result for judging whether i is less than n；N be the big gun examine face element quantity, 1 ≤i<n；
First return unit returns to the seismic channel search for enabling i=i+1 when first judging result is that i is less than n Unit；
Local plane wave compressed data block output unit, for when first judging result is that i is not less than n, output the 1st to be arrived The n Local plane wave compressed data blocks.
9. bilateral synthesis beam migration imaging system according to claim 8, which is characterized in that the trace gather Local plane wave Call number is established module and is specifically included:
Trace gather Local plane wave call number establishes unit, for the offset distance according to the Local plane wave, using formula Hi= (H1-H0)/DH+1 establishes the corresponding trace gather Local plane wave call number of the Local plane wave；Wherein Hi indicates i-th of part The corresponding trace gather Local plane wave call number of plane wave；H1 indicates the offset distance of the Local plane wave；H0 indicates offset distance road Collect first offset distance in parameter information；DH indicates offset distance interval.
10. bilateral synthesis beam migration imaging system according to claim 9, which is characterized in that the migration imaging module It specifically includes:
Local plane wave compressed data block acquiring unit, for obtaining i-th of Local plane wave data compression in k-th of trace gather Block；
Effective lineups position acquisition unit, for obtain stored in i-th of Local plane wave compressed data block it is effective The spatial position of lineups and time location；
Migration imaging unit, for carrying out the ray tracing of big gun inspection position according to the spatial position and the time location, and Migration imaging is carried out to crossover location, generates imaging temporary file；
Second judgment unit obtains the second judging result for judging whether i is less than L；Wherein L is office described in k-th of trace gather The quantity of facial planes wave compressed data block；
Second return unit, for enabling i=i+1, returning to effective lineups when second judging result is that i is less than L Position acquisition unit；
Third judging unit obtains third and sentences for judging whether k is less than K when second judging result is that i is not less than L Disconnected result；Wherein K is the quantity that trace gather is imaged；
Third return unit, for enabling k=k+1, returning to the Local plane wave when the third judging result is that k is less than K Compressed data block acquiring unit；
Trace gather file generating unit is imaged, for merging the 1st and arriving K trace gather when the third judging result is that k is not less than K Imaging temporary file, generate imaging trace gather file.
CN201810777046.2A 2018-07-16 2018-07-16 The bilateral synthesis beam offset imaging method and system picked up based on lineups Active CN109031417B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810777046.2A CN109031417B (en) 2018-07-16 2018-07-16 The bilateral synthesis beam offset imaging method and system picked up based on lineups

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810777046.2A CN109031417B (en) 2018-07-16 2018-07-16 The bilateral synthesis beam offset imaging method and system picked up based on lineups

Publications (2)

Publication Number Publication Date
CN109031417A true CN109031417A (en) 2018-12-18
CN109031417B CN109031417B (en) 2019-09-10

Family

ID=64642574

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810777046.2A Active CN109031417B (en) 2018-07-16 2018-07-16 The bilateral synthesis beam offset imaging method and system picked up based on lineups

Country Status (1)

Country Link
CN (1) CN109031417B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102667529A (en) * 2009-10-27 2012-09-12 雪佛龙美国公司 Method and system for seismic imaging and earth modeling using beam tomography
US20140293744A1 (en) * 2013-04-02 2014-10-02 Bp Corporation North America Inc. Specular filter (sf) and dip oriented partial imaging (dopi) seismic migration
US20150331123A1 (en) * 2014-05-15 2015-11-19 Acoustic Zoom, Inc. System and Method for Adaptive Classification and Filtering for Imaging Diffused Energy from the Earths Subsurface
CN106483559A (en) * 2015-08-24 2017-03-08 中国石油化工股份有限公司 A kind of construction method of subsurface velocity model
CN106855635A (en) * 2016-12-09 2017-06-16 中国石油天然气集团公司 A kind of method and apparatus of the bin attributes information for determining orthogonal observation system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102667529A (en) * 2009-10-27 2012-09-12 雪佛龙美国公司 Method and system for seismic imaging and earth modeling using beam tomography
US20140293744A1 (en) * 2013-04-02 2014-10-02 Bp Corporation North America Inc. Specular filter (sf) and dip oriented partial imaging (dopi) seismic migration
US20150331123A1 (en) * 2014-05-15 2015-11-19 Acoustic Zoom, Inc. System and Method for Adaptive Classification and Filtering for Imaging Diffused Energy from the Earths Subsurface
CN106483559A (en) * 2015-08-24 2017-03-08 中国石油化工股份有限公司 A kind of construction method of subsurface velocity model
CN106855635A (en) * 2016-12-09 2017-06-16 中国石油天然气集团公司 A kind of method and apparatus of the bin attributes information for determining orthogonal observation system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
甘利灯等: "测井-地震-油藏模拟一体化技术及其在老油田挖潜中的应用", 《石油物探》 *

Also Published As

Publication number Publication date
CN109031417B (en) 2019-09-10

Similar Documents

Publication Publication Date Title
CN108037526B (en) Reverse-time migration method based on all-wave wave field VSP/RVSP seismic data
CN104268412B (en) A kind of angle gathers ray tomography migration velocity analysis method and device
CN107607992B (en) Multi-wave matching method based on convolutional neural network
CN103645503B (en) A kind of three-dimensional time territory illumination analysis and vibration amplitude compensation method
CN102066980A (en) Seismic horizon skeletonization
CN110609320A (en) Pre-stack seismic reflection pattern recognition method based on multi-scale feature fusion
CN103675897A (en) Seismic diffracted wave separating and imaging method
US20040068375A1 (en) Petroleum exploration and prediction apparatus and method
CN103969683B (en) Method for picking position faces in batched mode based on constraint in three-dimensional seismic interpretation
CN109001813A (en) A kind of method, apparatus and system of multiple suppression
CN106443775A (en) High resolution converted wave crack prediction method
CN104199103B (en) Static correction method and static correction device for converted wave
CN109188506A (en) A kind of pure earth&#39;s surface stereo observing system suitable for high-speed rail tunnel bottom earthquake CT
CN103149592A (en) Method for separating variable offset vertical seismic profile (VSP) wave fields
CN106970417B (en) Ellipse expansion converted-wave velocity analysis method and system
CN110031898A (en) Data optimization methods and Kichhoff integral pre-stack depth migration method
CN101131435A (en) NMO correction stacking method preventing NMO stretching
CN109031417B (en) The bilateral synthesis beam offset imaging method and system picked up based on lineups
CN108226997A (en) A kind of seismic facies analysis method based on earthquake data before superposition
CN106990434B (en) Ellipse expansion converted wave imaging method and system
CN109901222A (en) Frequency method, storage medium and equipment are opened up in earthquake based on dynamic dictionary match tracing
CN104133244B (en) A kind of near surface static corrections processing method of focus shear wave exploration
CN112180433B (en) Method and device for picking up first arrival wave of earthquake
US20040114464A1 (en) Petroleum exploration and prediction apparatus and method
CN106226815A (en) The acquisition methods of converted wave velocity equivalent comparison and converted-wave velocity analysis method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant