CN108181645A - A kind of controlled source frequency dividing while exciting method, apparatus and system - Google Patents

A kind of controlled source frequency dividing while exciting method, apparatus and system Download PDF

Info

Publication number
CN108181645A
CN108181645A CN201711190612.1A CN201711190612A CN108181645A CN 108181645 A CN108181645 A CN 108181645A CN 201711190612 A CN201711190612 A CN 201711190612A CN 108181645 A CN108181645 A CN 108181645A
Authority
CN
China
Prior art keywords
sub
scanning signal
ramp function
controlled source
preset
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.)
Pending
Application number
CN201711190612.1A
Other languages
Chinese (zh)
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 National Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
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 National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN201711190612.1A priority Critical patent/CN108181645A/en
Publication of CN108181645A publication Critical patent/CN108181645A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/02Generating seismic energy
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/24Recording seismic data
    • G01V1/247Digital recording of seismic data, e.g. in acquisition units or nodes

Abstract

The embodiment of the present application discloses a kind of controlled source frequency dividing exciting method, apparatus and system simultaneously, the method includes targeted scans signal is divided based on preset division number, obtains sub- scanning signal;The ramp function of sub- scanning signal is built based on the first preset rules, first preset rules include:Adjacent sub- scanning signal overlapping position is complementary relationship in the ramp function of frequency domain;Processing is optimized to the slope length of sub- scanning signal based on the ramp function, the sub- scanning signal is made to meet the second preset rules, the sub- scanning signal after being optimized, wherein, second preset rules include:The energy fluctuation at adjacent sub- scanning signal overlapping position meets preset precision conditions;The controlled source of preset number is carried out at the same time seismic wave excitation based on the sub- scanning signal after the optimization.Using each embodiment of the application, can effectively ensure to obtain the hi-fi of seismic data using frequency dividing while the method for excitation.

Description

A kind of controlled source frequency dividing while exciting method, apparatus and system
Technical field
The present invention relates to oil gas development technology fields, particularly, are related to a kind of controlled source frequency dividing exciting method, dress simultaneously It puts and system.
Background technology
In seismic prospecting, controlled source has been obtained for being widely applied as the main tool for generating seismic signal. How to improve the hi-fi of its production efficiency and data becomes the hot spot studied.Traditional focus is scanned, in order to disappear Except influence of the harmonic wave interference to seismic data, the time interval of twice sweep has to be larger than equal to sweep time and listens the time.It listens Time is critically important, it can ensure that the data of generation are from harmonic wave interference.But listen time restriction production efficiency.
The Chinese invention patent of application number ZL201310545312.6 discloses a kind of controlled source frequency dividing while the side of excitation Method, the invention have carried out frequency-division section to scanning signal and have been divided into required subsignal, then by coding, make the letter excited simultaneously It is number mutually incoherent in frequency domain, so as to solve the high separation demand to adjacent big gun data, while production efficiency is improved, subtract The hi-fi of data after the input of few equipment and separation.When being excited simultaneously using fraction frequency signal progress focus, frequency dividing is swept It retouches signal needs and acquisition earthquake record is reconstructed after SNR estimation and compensation, but usual adjacent sub- scanning signal overlapping part is folded Closing energy can be there are certain fluctuation, if overlapping energy fluctuation is larger, it will reduce the high-fidelity of collected seismic data Property.Therefore, how effectively to ensure the hi-fi of seismic data obtained using frequency dividing while the method for excitation, become urgently The technical issues of solution.
Invention content
The embodiment of the present application is designed to provide a kind of controlled source frequency dividing exciting method, apparatus and system simultaneously, can With the effective hi-fi for ensureing the method using frequency dividing excitation simultaneously and obtaining seismic data.
Exciting method, apparatus and system are by including in the following manner to a kind of controlled source frequency dividing that the application provides simultaneously It realizes:
A kind of controlled source frequency dividing while exciting method, including:
Targeted scans signal based on preset division number is divided, obtains sub- scanning signal;
The ramp function of sub- scanning signal is built based on the first preset rules, first preset rules include:Adjacent son It is complementary relationship that scanning signal, which overlaps position in the ramp function of frequency domain,;
Processing is optimized to the slope length of sub- scanning signal based on the ramp function, expires the sub- scanning signal The second preset rules of foot, the sub- scanning signal after being optimized, wherein, second preset rules include:Adjacent sub- scanning letter Number overlapping position energy fluctuation meet preset precision conditions;
The controlled source of preset number is carried out at the same time seismic wave excitation based on the sub- scanning signal after the optimization.
Exciting method, the first preset rules structure that is based on scan the controlled source frequency dividing of the embodiment of the present application simultaneously The ramp function of signal, including:
The ramp function of time-domain, the slope letter of the time-domain are obtained according to the ramp function of sub- scanning signal frequency domain Evolution of the number for frequency domain ramp function.
Exciting method, the first preset rules structure that is based on scan the controlled source frequency dividing of the embodiment of the present application simultaneously The ramp function of signal, including:
The energy fluctuation value at adjacent sub- scanning signal overlapping position is calculated based on a plurality of types of ramp functions, by minimum energy The corresponding ramp function type of amount undulating value is determined as the ramp function type of the sub- scanning signal;
The slope of sub- scanning signal is built based on the first preset rules according to determining sub- scanning signal ramp function type Function.
The controlled source frequency dividing while exciting method of the embodiment of the present application, the controlled source of the preset number is based on described Sub- scanning signal after optimization is carried out at the same time after seismic wave excitation, is further included:
The controlled source of preset number is carried out at the same time seismic wave excitation based on the sub- scanning signal after the optimization, obtains ground The female record data of shake;
Each sub- scanning signal of excitation simultaneously or corresponding force signal are recorded into data into line number with corresponding earthquake mother According to separating treatment, recorded in data from earthquake mother and isolate the corresponding earthquake subrecord data of each sub- scanning signal;
The corresponding earthquake subrecord data of complete scan sequence of single shot point are obtained, are swept according to sub- scanning signal in target The earthquake subrecord data of single shot point are sequentially overlapped by the sequence retouched in signal, obtain single-shot earthquake record.
The corresponding earthquake subrecord data of complete scan sequence of single shot point are obtained, are swept according to sub- scanning signal in target The earthquake subrecord data of single shot point are sequentially overlapped by the sequence retouched in signal, obtain single-shot earthquake record.
On the other hand, the embodiment of the present application also provides a kind of controlled source frequency dividing excitation apparatus simultaneously, including:
Frequency division module for targeted scans signal to be divided based on preset division number, obtains sub- scanning signal;
Ramp function builds module, for building the ramp function of sub- scanning signal based on the first preset rules, described the One preset rules include:Adjacent sub- scanning signal overlapping position is complementary relationship in the ramp function of frequency domain;
Optimization module for optimizing processing to the slope length of sub- scanning signal based on the ramp function, makes institute It states sub- scanning signal and meets the second preset rules, the sub- scanning signal after being optimized, wherein, the second preset rules packet It includes:The energy fluctuation at adjacent sub- scanning signal overlapping position meets preset precision conditions;
Excitation module is carried out at the same time earthquake for the controlled source of preset number based on the sub- scanning signal after the optimization Wave excites.
Excitation apparatus, the ramp function structure module include the controlled source frequency dividing of the embodiment of the present application simultaneously:
The ramp function of time-domain, the slope letter of the time-domain are obtained according to the ramp function of sub- scanning signal frequency domain Evolution of the number for frequency domain ramp function.
Excitation apparatus, the ramp function structure module include the controlled source frequency dividing of the embodiment of the present application simultaneously:
Ramp function type optimizes unit, and adjacent sub- scanning signal overlapping is calculated for being based on a plurality of types of ramp functions The corresponding ramp function type of least energy undulating value is determined as the slope of the sub- scanning signal by the energy fluctuation value at position Type function;
Ramp function construction unit, for being based on the first preset rules according to determining sub- scanning signal ramp function type Build the ramp function of sub- scanning signal.
Excitation apparatus, described device further include the controlled source frequency dividing of the embodiment of the present application simultaneously:
Data capture unit is carried out at the same time for the controlled source of preset number based on the sub- scanning signal after the optimization Seismic wave excites, and obtains earthquake mother and records data;
Data separating unit, for each sub- scanning signal that will be excited simultaneously or corresponding force signal and corresponding earthquake Mother's record data carry out data separating processing, are recorded from earthquake mother and the corresponding earthquake note of each sub- scanning signal is isolated in data Record data;
Data investigation unit, for obtaining the corresponding earthquake subrecord data of the complete scan sequence of single shot point, according to The earthquake subrecord data of single shot point are sequentially overlapped by sequence of the sub- scanning signal in targeted scans signal, are obtained single Big gun earthquake record.
The controlled source frequency dividing while excitation apparatus of the embodiment of the present application, can hold including processor and for storing processor The memory of row instruction is realized when described instruction is performed by the processor and is included the following steps:
Targeted scans signal based on preset division number is divided, obtains sub- scanning signal;
The ramp function of sub- scanning signal is built based on the first preset rules, first preset rules include:Adjacent son It is complementary relationship that scanning signal, which overlaps position in the ramp function of frequency domain,;
Processing is optimized to the slope length of sub- scanning signal based on the ramp function, expires the sub- scanning signal The second preset rules of foot, the sub- scanning signal after being optimized, wherein, second preset rules include:Adjacent sub- scanning letter Number overlapping position energy fluctuation meet preset precision conditions;
The controlled source of preset number is carried out at the same time seismic wave excitation based on the sub- scanning signal after the optimization.
On the other hand, the embodiment of the present application also provides a kind of controlled source frequency dividing activating system simultaneously, including it is multiple simultaneously Controlled source, wave detector, at least one processor and the memory for storing computer executable instructions of excitation, the processing The step of device realizes above-mentioned any one the method when performing described instruction.
A kind of controlled source frequency dividing while exciting method, device and be that this specification one or more embodiment provides System, the ramp function of sub- scanning signal that can be obtained by being primarily based on the first preset rules structure frequency dividing, described first is pre- If rule can include:Adjacent sub- scanning signal overlapping position is complementary relationship in the ramp function of frequency domain.Then, based on institute It states ramp function and processing is optimized to the slope length of sub- scanning signal, the sub- scanning signal is made to meet the second default rule Then, the sub- scanning signal after being optimized, wherein, second preset rules include:Adjacent sub- scanning signal overlapping position Energy fluctuation meets preset precision conditions.Multiple controlled sources complete frequency dividing simultaneously based on the sub- scanning signal after the optimization The implementation of excitation so as to reduce the energy fluctuation at adjacent sub- scanning signal overlapping position in restructuring procedure, ensures single-shot earthquake number According to hi-fi.
Description of the drawings
In order to illustrate the technical solutions in the embodiments of the present application or in the prior art more clearly, to embodiment or it will show below There is attached drawing needed in technology description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this Some embodiments described in application, for those of ordinary skill in the art, in the premise of not making the creative labor property Under, it can also be obtained according to these attached drawings other attached drawings.In the accompanying drawings:
Fig. 1 is a kind of controlled source frequency dividing that this specification provides while the flow diagram of exciting method embodiment;
Sub- scanning signal is divided in one embodiment that Fig. 2 is provided for this specification and its corresponding overlapping position slope shows It is intended to;
Fig. 3 is the analog result signal that ramp function is Blackman functions in one embodiment that this specification provides Figure;
Fig. 4 is the analog result schematic diagram that ramp function is Hanning functions in one embodiment that this specification provides;
Fig. 5 is the analog result schematic diagram that ramp function is Hamming functions in one embodiment that this specification provides;
Fig. 6 is the analog result schematic diagram that ramp function is Kaiser functions in one embodiment that this specification provides;
Fig. 7 is the analog result schematic diagram that ramp function is sin cos functions in one embodiment that this specification provides;
Fig. 8 is the analog result schematic diagram that ramp function is linear function in one embodiment that this specification provides;
Fig. 9 is the analog result schematic diagram that slope length is 300ms in another embodiment that this specification provides;
Figure 10 is the analog result schematic diagram that slope length is 500ms in another embodiment that this specification provides;
Figure 11 is another embodiment neutron scanning signal lamination portion slope design diagram that this specification provides;
Figure 12 is that another embodiment neutron scanning signal that this specification provides overlaps energy fluctuation schematic diagram;
Figure 13 is that single-shot excites complete scan sequence diagram in another embodiment that this specification provides;
Figure 14 is the firsthand information schematic diagram of frequency dividing excitation field acquisition in another embodiment that this specification provides;
Figure 15 is the modular structure signal of a kind of controlled source frequency dividing that this specification provides while excitation apparatus embodiment Figure.
Specific embodiment
In order to which those skilled in the art is made to more fully understand the technical solution in this specification, below in conjunction with this explanation Attached drawing in book one or more embodiment carries out the technical solution in this specification one or more embodiment clear, complete Site preparation describes, it is clear that described embodiment is only specification part of the embodiment, instead of all the embodiments.Based on saying Bright book one or more embodiment, those of ordinary skill in the art are obtained all without making creative work Other embodiment should all belong to the range of this specification embodiment scheme protection.
The conventional use of controlled source can include following steps:First, one or more controlled source is disposed On a shot point.Then start to be driven controlled source with a scanning signal, focus starts to vibrate, time of vibration one As for 10 to 16 seconds, scanning signal is typically the signal that the frequency of one changes over time.Then, wave detector is by a period of time Response data is received and is recorded, and this period listens the time equal to sweep time plus one.This time should include seismic wave from It generates, main purpose layer is reflected into the minimum required time of reception.Finally, it is mutual by being carried out to record data and reference scan Pass obtains earthquake record.Certainly it repeats the above steps and will be remembered several times on same shot point sometimes for increase signal strength Record, and be added together and to form final earthquake record.Focus can be carried out the vibration of next shot point after the completion.
The Chinese invention patent of application number ZL201310545312.6 discloses a kind of controlled source frequency dividing while the side of excitation Method, the invention have carried out frequency-division section to scanning signal and have been divided into required subsignal, then by coding, make the letter excited simultaneously It is number mutually incoherent in frequency domain, so as to solve the high separation demand to adjacent big gun data, while production efficiency is improved, subtract The hi-fi of data after the input of few equipment and separation.When being excited simultaneously using fraction frequency signal progress focus, frequency dividing is swept It retouches signal needs and acquisition earthquake record is reconstructed after SNR estimation and compensation, but in restructuring procedure, adjacent sub- scanning signal overlapping Usually there are certain fluctuations for the overlapping energy at position, if overlapping energy fluctuation is larger, can influence collected earthquake number According to hi-fi.
In view of the above problems, the application provides a kind of controlled source frequency dividing exciting method simultaneously.It can be first pre- based on first If the ramp function of the sub- scanning signal after rule structure frequency dividing, first preset rules can include:Adjacent sub- scanning letter It is complementary relationship that number ramp function, which meets in frequency domain,.Then, it is based on second in advance according to the sub- scanning signal ramp function of structure If rule optimizes processing to sub- scanning signal slope length, the sub- scanning signal is made to meet the second preset rules, is obtained Sub- scanning signal after optimization, wherein, two preset rules can include:The energy wave at adjacent sub- scanning signal overlapping position It is dynamic to meet preset precision conditions.The controlled source of preset number can be carried out at the same time according to the sub- scanning signal after the optimization Seismic wave excites, so as to which the frequency dividing for completing controlled source excites simultaneously.It, can be with using the above-mentioned frequency dividing of the application exciting method simultaneously The energy fluctuation at adjacent sub- scanning signal overlapping position in restructuring procedure is reduced, is further ensured that the high-fidelity of single-shot seismic data Property.
Fig. 1 is a kind of controlled source frequency dividing that this specification provides while exciting method embodiment flow diagram. Although present description provides such as following embodiments or method operating procedure shown in the drawings or apparatus structure, based on conventional or Person can include operating procedure less after more or part merges without performing creative labour in the method or device Or modular unit.In logicality there is no the step of necessary causality or in structure, the execution sequence or dress of these steps The modular structure put is not limited to this specification embodiment or execution shown in the drawings sequence or modular structure.The method or mould Device in practice, server or the end product of block structure are in application, can be according to embodiment or side shown in the drawings Method or the execution of modular structure carry out sequence be either parallel to be performed (such as the environment of parallel processor or multiple threads, even Implementation environment including distributed treatment, server cluster).
Specific one embodiment is as shown in Figure 1, a kind of controlled source frequency dividing while exciting method that this specification provides One embodiment in, the method may include:
S2, targeted scans signal based on preset division number is divided, obtains sub- scanning signal;
In the present embodiment, depth of stratum, harmonic wave interference etc. can be reflected by comprehensive analysis because usually determining targeted scans Signal.The targeted scans signal can be raising frequency setting or frequency reducing is set.It is then possible to according to preset division number Targeted scans signal in frequency domain is divided into several frequency bands, obtains sub- scanning signal, then the sub- scanning signal is in frequency Domain is orthogonal.Frequency dividing can be the frequency range even single-frequency of decile or random cutting.Assuming that the sub- scanning signal after frequency dividing Respectively s1,s2,...,sn, frequency band is respectively f1-f1+Δf1,f2-f2+Δf2,f3-f3+Δf3,...,fn-fn+Δfn, then s1 And s2、s2And s3For adjacent sub- scanning signal.
It, can be according to the noiseless length of record needs and simultaneously between explosive source in one embodiment of this specification Spatial distribution, determine preset division number.In one or more embodiment of this specification, it is further contemplated that focus is walked Time calculates preset division number, required divider ratio when being excited simultaneously so as to the focus for being further reduced similary number of units.Then such as Fruit preset division number is certain, then can realize that more focus excite simultaneously, so as to improve working efficiency.
S4, the ramp function that sub- scanning signal is built based on the first preset rules, first preset rules are included:It is adjacent Sub- scanning signal overlapping position is complementary relationship in the ramp function of frequency domain;
Usually by targeted scans signal divide when, adjacent sub- scanning signal can be overlapped position and be designed as slope Form, to reduce jeep this effect.In the present embodiment, the ramp function that position slope can be overlapped by structure is swept to describe son Retouch the feature on signal slope.As shown in Fig. 2, Fig. 2 represents to divide sub- scanning signal and its corresponding overlapping position slope schematic diagram. In one embodiment of this specification, the ramp function of sub- scanning signal can be built according to preset rules, to reduce earthquake note Record the energy fluctuation that position is overlapped in restructuring procedure.It, can be according to the needs of practical operation in one embodiment of this specification The ramp function type of sub- scanning signal is predefined, the ramp function can be any one letter in conventional ramp function Several classes of types, such as ramp function type can be Blacknam function, breathe out bright function, linear function etc..
It is then possible to build the ramp function of sub- scanning signal based on the first preset rules, first preset rules can To include:It is complementary relationship that adjacent sub- scanning signal, which overlaps the corresponding ramp function in position in frequency domain,.To ensure to sweep after dividing The amplitude that signal is retouched in frequency domain is consistent with former targeted scans signal, thereby may be ensured that in earthquake record restructuring procedure The energy fluctuation value for overlapping position is smaller.
As illustrated in fig. 2, it is assumed that former targeted scans amplitude is A, in order to ensure fraction frequency signal and former targeted scans signal Amplitude is consistent, then the amplitude at adjacent sub- scanning signal overlapping position should be made to be added also as A.If adjacent sub- scanning signal overlapping portion Position is a (f) in a ramp function of frequency domain, then another side can be A-a (f) in the ramp function b (f) of frequency domain.To A After being normalized, adjacent sub- scanning signal overlapping position is complementary relationship in the ramp function of frequency domain, i.e. b (f)=1-a (f)。
In one embodiment of this specification, the seismic data separation process of usually frequency dividing excitation simultaneously can be swept by son The auto-correlation of signal is retouched to realize, auto-correlation correspondence is the relationship of product in frequency domain.Correspondingly, the slope letter by frequency domain Number is transformed into time-domain, then two neighboring sub- scanning signal is respectively in the ramp function of time-domain And meet
In one embodiment of this specification, multiple types ramp function can be based on and calculate adjacent sub- scanning signal overlapping The corresponding ramp function type of least energy undulating value is determined as the slope of the sub- scanning signal by the energy fluctuation value at position Type function.It is then possible to the ramp function of sub- scanning signal is built according to determining sub- scanning signal ramp function type.No The ramp function of same type has different influences to the energy fluctuation at sub- scanning signal overlapping position, can be in the present embodiment By carrying out analysis comparison to the influence for overlapping position energy fluctuation to different ramp function types, optimal slope is determined Type function is used for the structure of sub- scanning signal ramp function, so as to further reduce ramp function to overlapping position energy The influence of fluctuation.In one or more embodiment of this specification, the analyses sides such as numerical simulation analysis, tabular analysis can be passed through Method determines the ramp function type of sub- scanning signal.
It, can be respectively with Blacknam (Blackman) function, peaceful (Hanning) letter of the Chinese in this specification one embodiment Number breathes out bright (Hamming) function, Kai Ze (Kaiser) function, sin cos functions and linear (triangle) function as phase The type function of the ramp function of adjacent sub- scanning signal carries out numerical simulation analysis, determines optimal ramp function type, uses In the structure of ramp function.
In the specific example that this specification provides, it is assumed that the sweep speed of frequency dividing excitation simultaneously is 5hz/s, subsignal Sweep length 2 seconds, adjacent scanning subsignal frequency band is respectively 6-16hz and 15-25hz, overlapping bands width 1hz, slope length 200ms.Then, numerical simulation analysis is carried out to above-mentioned ramp function type, determines optimal ramp function type, such as Fig. 3 extremely Shown in Fig. 8.Fig. 3 represents that the ramp function of adjacent sub- scanning signal is respectively radical sign (Sqrt) Blackman functions and radical sign 1- The analog result of Blackman functions;Fig. 4 represent the ramp function of adjacent sub- scanning signal be respectively radical sign Hanning functions with And the analog result of radical sign 1-Hanning functions;Fig. 5 represents that the ramp function of adjacent sub- scanning signal is respectively radical sign The analog result of Hamming functions and radical sign 1-Hamming functions;Fig. 6 represents the ramp function point of adjacent sub- scanning signal Not Wei radical sign Kaiser functions and radical sign 1-Kaiser functions analog result;Fig. 7 represents the slope of adjacent sub- scanning signal Function is respectively the analog result of sinusoidal (sin) function and cosine (cos) function;Fig. 8 shows the oblique of adjacent sub- scanning signal Slope function is respectively the analog result of radical sign linear function and radical sign 1- linear functions.At the band overlapping of Fig. 3 to Fig. 8 Energy fluctuation is analyzed it is found that the ramp function of adjacent sub- scanning signal is radical sign linear function and radical sign 1- linear functions Energy fluctuation value it is smaller.Therefore, in one or more embodiment of this specification, it is preferred that linear function can be made Excite the ramp function type of adjacent sub- scanning signal simultaneously for frequency dividing.
S6, processing is optimized to the slope length of sub- scanning signal based on the ramp function, makes the sub- scanning letter Number meet the second preset rules, the sub- scanning signal after being optimized, wherein, second preset rules include:Adjacent son is swept The energy fluctuation for retouching signal overlapping position meets preset precision conditions;
In the present embodiment, the slope length can be that adjacent sub- scanning signal overlaps the corresponding time span in position, right Should be overlapping frequency bandwidth in frequency domain.When it is implemented, it can be swept according to ramp function antithetical phrase determining through the above scheme The slope length for retouching signal optimizes processing so that sub- scanning signal meets the second preset rules, i.e., adjacent sub- scanning signal The energy fluctuation at overlapping position meets preset precision conditions, so as to the sub- scanning signal after being optimized.When it is implemented, tool The optimization means of body are unlimited, and the realization rate for meeting the second preset rules belongs to the mode that the present embodiment can be implemented.Such as The slope length for the sub- scanning signal for meeting condition can be determined by modes such as tabular analysis, constructor analyses.It is described Preset precision conditions can preset the required precision for splicing energy according in practical operation.This specification provides One or more embodiment in, can be according to preset precision conditions, by digital simulation to the slope of sub- scanning signal Length optimizes processing.
Assuming that the frequency starting point at sub- scanning signal overlapping position is f1, terminal f2, ramp function type is linear letter Number, then the ramp function of frequency domain can be expressed as:
Wherein, A is amplitude, and f is closed interval [f1, f2] on arbitrary value, f2-f1To be overlapped frequency bandwidth Δ f, when being transformed into Between domain can obtain:
Wherein, Δ t overlaps the corresponding slope length in position for adjacent sub- scanning signal, and k is sweep speed, then overlapping bands Width is one-to-one relationship with slope length.The ramp function of time-domain can be expressed asAssuming that preset energy wave Dynamic precision conditions are:Frequency overlapping amplitude fluctuation is less than 3dB.Sweep speed is 5Hz/s, and sub- sweep length is 2 seconds.Pass through number The slope length of sub- scanning signal is determined in value sunykatuib analysis, as shown in Figure 9 and Figure 10.Fig. 9 represents that slope length is 300ms, Corresponding overlaid frequency is the analog result of 1.5hz;Figure 10 represents that slope length is 500ms, and corresponding overlaid frequency is 2.5hz Analog result.By comparison diagram 9 and Figure 10 it is found that the overlapping energy fluctuation that adjacent sub- scanning signal overlaps position in Figure 10 is expired The preset precision conditions of foot.Therefore, the overlaid frequency of adjacent sub- scanning signal and slope length can be identified as 500ms、2.5hz。
In another specific example of this specification, as shown in figure 11, by sweep length be 20 seconds, frequency is 3-84Hz's Scanning signal is divided into two sub- scanning signals, and the sweep length of two sub- scanning signals is respectively 10s, and sweep speed is identical, overlapping Less than 3dB, ramp function is Blackman functions for amplitude fluctuation requirement at frequency band.Figure 12 represents simulation overlapping as a result, when weight When folded frequency bandwidth reaches 3Hz, amplitude fluctuation is less than 2.6dB, and the energy fluctuation at adjacent sub- scanning signal overlapping position meets in advance If precision conditions, then can determine adjacent sub- scanning signal overlapping frequency bandwidth be 3Hz.
S8, preset number controlled source based on the sub- scanning signal after the optimization be carried out at the same time seismic wave excitation.
The number that controlled source excites simultaneously can be preset according to practical operation, then can be according to described excellent Sub- scanning signal after change establishes frequency dividing data, divides sub- scanning signal by command centre's point dynamic, commander's controlled source is simultaneously Implement seismic wave excitation;Frequency range task pool can also be established by command centre, and establishes corresponding task mechanism, controlled source is certainly It is dynamic that the excitation of task-cycle seismic wave is obtained from task pool.
In another embodiment of this specification, the controlled source of preset number is obtained based on the son scanning after the optimization Signal is carried out at the same time seismic wave and corresponding earthquake mother is excited to record data.When it is implemented, consecutive shock record number can be acquired According to, according to the firing time of controlled source to the consecutive shock data carry out cutting, it is corresponding to obtain each firing time section Earthquake mother records data.Then, by the sub- scanning signal of each firing time section or corresponding force signal, with the time The earthquake mother of section records data and carries out related or convolution, realizes data separating processing, records in data and isolate from earthquake mother The corresponding earthquake subrecord data of each sub- scanning signal.
Sub- scanning signal between the controlled source that frequency dividing excites simultaneously is orthogonal, therefore, son scanning can be utilized to believe Number to the earthquake mother of acquisition is excited to record related or convolution the method for data progress simultaneously, realize the separation of seismic data, obtain Obtain the corresponding earthquake subrecord of each sub- scanning signal of focus.If it is considered that the response of Seismic Source System, can also utilize son scanning to believe It is related that number corresponding force signal and earthquake mother record data progress, realizes the separation of seismic data.
It a focus (or more source patterns) can be set must complete targeted scans signal frequency split in same shot point and obtain The excitation task for completing a shot point is just calculated in the seismic wave excitation of all sub- scanning signals obtained, i.e., one shot point of completion is complete The excitation of whole scanning sequence.The corresponding each earthquake subrecord data of complete scan sequence of single shot point are obtained, according to sub- scanning The earthquake subrecord data of single shot point are sequentially overlapped by signal putting in order in targeted scans signal, obtain single-shot Earthquake record.
Assuming that the corresponding sub- scanning signal of N platform focus excited simultaneously is respectively ss1,ss2...ssi...ssN, ss1, ss2...ssi...ssNBetween it is orthogonal;The seismic response function on stratum is h (t), and N platform focus is based respectively on respective son and sweeps Retouch signal ss1,ss2...ssi...ssNIt is S (t) to be carried out at the same time the earthquake mother that seismic wave excitation obtains and record data, and the earthquake is female Record data can be expressed as:
Then the corresponding earthquake subrecord of the sub- scanning signal of each focus can be expressed as:
If the preset division number is n, then the corresponding n earthquake subrecord of same shot point complete scan sequence is obtained, it will N earthquake subrecord of same shot point is overlapped successively in order, obtains the earthquake record r (t) of the shot point:
In another specific example of this specification, Figure 13 represents that the single-shot that the sub- scanning signal after optimization is formed completely is swept Signal sequence is retouched, which is 20 fraction frequency signals.The left hand view of Figure 14 represents that focus is based on 20 frequency dividing and sweeps The single-shot earthquake record that signal sequence is excited and obtained after reconstructing is retouched, the right part of flg of Figure 14 represents that the focus is based on 20 frequency dividing and sweeps Collected firsthand information spectrogram after signal sequence excites is retouched, can be seen that the frequency spectrum after overlapping does not have by Figure 14 left hand views Apparent amplitude fluctuation.
Each embodiment in this specification is described by the way of progressive, identical similar portion between each embodiment Point just to refer each other, and the highlights of each of the examples are difference from other examples.Specifically it is referred to The description of aforementioned relevant treatment related embodiment, does not do repeat one by one herein.
It is above-mentioned that this specification specific embodiment is described.Other embodiments are in the scope of the appended claims It is interior.In some cases, the action recorded in detail in the claims or step can be come according to different from the sequence in embodiment It performs and still can realize desired result.In addition, the process described in the accompanying drawings not necessarily require show it is specific suitable Sequence or consecutive order could realize desired result.In some embodiments, multitasking and parallel processing be also can With or it may be advantageous.
A kind of controlled source while exciting method that this specification one or more embodiment provides, can pass through base first The ramp function for the sub- scanning signal that frequency dividing obtains is built in the first preset rules, first preset rules can include:Phase Adjacent sub- scanning signal overlapping position is complementary relationship in the ramp function of frequency domain.Then, it is swept based on the ramp function antithetical phrase The slope length for retouching signal optimizes processing, and the sub- scanning signal is made to meet the second preset rules, the son after being optimized Scanning signal, wherein, second preset rules include:The energy fluctuation at adjacent sub- scanning signal overlapping position meets preset Precision conditions.Multiple controlled sources complete the implementation of frequency dividing excitation simultaneously based on the sub- scanning signal after the optimization, so as to subtract The energy fluctuation at adjacent sub- scanning signal overlapping position in few restructuring procedure ensures the hi-fi of single-shot seismic data.
Based on controlled source frequency dividing described above, exciting method, this specification one or more embodiment also provide simultaneously A kind of true controlled source frequency dividing while excitation apparatus.The device can include the use of this specification embodiment the method System, software (application), module, component, server etc. and combine the necessary device for implementing hardware.Based on same innovation structure Think, the device in one or more embodiments that this specification embodiment provides is as described in the following examples.Due to device solution Certainly the implementation of problem is similar to method, therefore the implementation of the specific device of this specification embodiment may refer to preceding method Implementation, overlaps will not be repeated.Used below, term " unit " or " module " can realize the soft of predetermined function The combination of part and/or hardware.Although following embodiment described device is preferably realized with software, hardware or The realization of the combination of software and hardware is also what may and be contemplated.Specifically, one kind that Figure 15, which is this specification, to be provided is controllable The modular structure schematic diagram of focus frequency dividing while excitation apparatus embodiment, as shown in figure 15, described device can include:
Frequency division module 102 can be used for dividing targeted scans signal based on preset division number, obtain son scanning letter Number;
Ramp function builds module 104, can be used for building the ramp function of sub- scanning signal based on the first preset rules, First preset rules include:Adjacent sub- scanning signal overlapping position is complementary relationship in the ramp function of frequency domain;
Optimization module 106 can be used for optimizing place to the slope length of sub- scanning signal based on the ramp function Reason makes the sub- scanning signal meet the second preset rules, the sub- scanning signal after being optimized, wherein, described second presets Rule includes:The energy fluctuation at adjacent sub- scanning signal overlapping position meets preset precision conditions;
Excitation module 108 can be used for the controlled source of preset number based on the sub- scanning signal after the optimization simultaneously Carry out seismic wave excitation.
Certainly, it is described with reference to preceding method embodiment, in the other embodiment of described device, the ramp function builds mould Block 104 can include:
The ramp function of time-domain, the slope letter of the time-domain are obtained according to the ramp function of sub- scanning signal frequency domain Evolution of the number for frequency domain ramp function.
In another embodiment for the described device that this specification provides, the ramp function structure module 104 can wrap Ramp function type optimization unit and ramp function construction unit are included, wherein,
The ramp function type optimizes unit, can be used for calculating adjacent sub- scanning based on a plurality of types of ramp functions Signal overlaps the energy fluctuation value at position, and the corresponding ramp function type of least energy undulating value is determined as the sub- scanning believes Number ramp function type;
The ramp function construction unit can be used for being based on first according to determining sub- scanning signal ramp function type Preset rules build the ramp function of sub- scanning signal.
Need what is illustrated, device described above can also include other embodiment party according to the description of embodiment of the method Formula.Concrete implementation mode is referred to the description of related method embodiment, does not repeat one by one herein.
A kind of controlled source frequency dividing while excitation apparatus that this specification one or more embodiment provides, can pass through head The ramp function for the sub- scanning signal that frequency dividing obtains first is built based on the first preset rules, first preset rules can wrap It includes:Adjacent sub- scanning signal overlapping position is complementary relationship in the ramp function of frequency domain.Then, based on the ramp function pair The slope length of sub- scanning signal optimizes processing, the sub- scanning signal is made to meet the second preset rules, after being optimized Sub- scanning signal, wherein, second preset rules include:The energy fluctuation at adjacent sub- scanning signal overlapping position meets pre- If precision conditions.Multiple controlled sources complete the implementation of frequency dividing excitation simultaneously based on the sub- scanning signal after the optimization, from And the energy fluctuation at adjacent sub- scanning signal overlapping position in restructuring procedure is reduced, ensure the hi-fi of single-shot seismic data.
The method or apparatus described in above-described embodiment that this specification provides can realize that business is patrolled by computer program It collects and records on a storage medium, the storage medium can be read and be performed with computer, realize this specification embodiment institute The effect of description scheme.Therefore, this specification also provides a kind of controlled source frequency dividing excitation apparatus simultaneously, including processor and deposits The memory of processor-executable instruction is stored up, realizes and includes the following steps when described instruction is performed by the processor:
Targeted scans signal based on preset division number is divided, obtains sub- scanning signal;
The ramp function of sub- scanning signal is built based on the first preset rules, first preset rules include:Adjacent son It is complementary relationship that scanning signal, which overlaps position in the ramp function of frequency domain,;
Processing is optimized to the slope length of sub- scanning signal based on the ramp function, expires the sub- scanning signal The second preset rules of foot, the sub- scanning signal after being optimized, wherein, second preset rules include:Adjacent sub- scanning letter Number overlapping position energy fluctuation meet preset precision conditions;
The controlled source of preset number is carried out at the same time seismic wave excitation based on the sub- scanning signal after the optimization.
The storage medium can be included for the physical unit that stores information, typically by after information digitalization again with profit The media of the modes such as electricity consumption, magnetic or optics are stored.The storage medium has and can include:Letter is stored using electric energy mode The device of breath such as, various memory, such as RAM, ROM;The device of information is stored using magnetic energy mode such as, hard disk, floppy disk, magnetic Band, core memory, magnetic bubble memory, USB flash disk;Using optical mode storage information device such as, CD or DVD.Certainly, also its The readable storage medium storing program for executing of his mode, such as quantum memory, graphene memory etc..
Need what is illustrated, device described above can also include other embodiment party according to the description of embodiment of the method Formula.Concrete implementation mode is referred to the description of related method embodiment, does not repeat one by one herein.
A kind of controlled source frequency dividing while excitation apparatus described in above-described embodiment, can be default by being primarily based on first The ramp function of sub- scanning signal that rule structure frequency dividing obtains, first preset rules can include:Adjacent sub- scanning letter Number overlapping position frequency domain ramp function be complementary relationship.Then, based on the ramp function to the oblique of sub- scanning signal Length of grade degree optimizes processing, and the sub- scanning signal is made to meet the second preset rules, the sub- scanning signal after being optimized, In, second preset rules include:The energy fluctuation at adjacent sub- scanning signal overlapping position meets preset precision conditions.It is more A controlled source completes the implementation of frequency dividing excitation simultaneously based on the sub- scanning signal after the optimization, so as to reduce in restructuring procedure The energy fluctuation at adjacent sub- scanning signal overlapping position ensures the hi-fi of single-shot seismic data.
This specification also provides a kind of controlled source frequency dividing activating system simultaneously, and the system can be individually controllable shake Source while activating system, can also apply in a plurality of types of system for acquiring seismic data.The system can be independent Server, one or more the methods of this specification or one or more embodiment device can also be include the use of Server cluster, system (including distributed system), software (application), practical operation device, logic gates device, quantum meter Calculation machine etc. simultaneously combines the necessary terminal installation for implementing hardware.The excitation simultaneously of controlled source frequency dividing can include it is multiple simultaneously Controlled source, wave detector, at least one processor and the memory for storing computer executable instructions of excitation, the processing The step of device realizes method described in any one above-mentioned embodiment when performing described instruction.
Need what is illustrated, system described above can also include others according to the description of method or device embodiment Embodiment, concrete implementation mode are referred to the description of related method embodiment, do not repeat one by one herein.
A kind of controlled source frequency dividing while activating system described in above-described embodiment, can be default by being primarily based on first The ramp function of sub- scanning signal that rule structure frequency dividing obtains, first preset rules can include:Adjacent sub- scanning letter Number overlapping position frequency domain ramp function be complementary relationship.Then, based on the ramp function to the oblique of sub- scanning signal Length of grade degree optimizes processing, and the sub- scanning signal is made to meet the second preset rules, the sub- scanning signal after being optimized, In, second preset rules include:The energy fluctuation at adjacent sub- scanning signal overlapping position meets preset precision conditions.It is more A controlled source completes the implementation of frequency dividing excitation simultaneously based on the sub- scanning signal after the optimization, so as to reduce in restructuring procedure The energy fluctuation at adjacent sub- scanning signal overlapping position ensures the hi-fi of single-shot seismic data.
It should be noted that this specification device described above or system according to the description of related method embodiment also It can include other embodiments, concrete implementation mode is referred to the description of embodiment of the method, does not go to live in the household of one's in-laws on getting married one by one herein It states.Each embodiment in this specification is described by the way of progressive, and identical similar part is mutual between each embodiment Mutually referring to the highlights of each of the examples are difference from other examples.Especially for hardware+program For class, storage medium+program embodiment, since it is substantially similar to embodiment of the method, so description is fairly simple, it is related Part illustrates referring to the part of embodiment of the method.
It is above-mentioned that this specification specific embodiment is described.Other embodiments are in the scope of the appended claims It is interior.In some cases, the action recorded in detail in the claims or step can be come according to different from the sequence in embodiment It performs and still can realize desired result.In addition, the process described in the accompanying drawings not necessarily require show it is specific suitable Sequence or consecutive order could realize desired result.In some embodiments, multitasking and parallel processing be also can With or it may be advantageous.
System, device, module or the unit that above-described embodiment illustrates can specifically realize by computer chip or entity, Or it is realized by having the function of certain product.A kind of typical realization equipment is computer.Specifically, computer for example may be used Think personal computer, laptop computer, vehicle-mounted human-computer interaction device, cellular phone, camera phone, smart phone, individual Digital assistants, media player, navigation equipment, electronic mail equipment, game console, tablet computer, wearable device or The combination of any equipment in these equipment of person.
For convenience of description, it is divided into various modules during description apparatus above with function to describe respectively.Certainly, implementing this The function of each module is realized can in the same or multiple software and or hardware during specification one or more, it can also The module for realizing same function is realized by the combination of multiple submodule or subelement etc..Device embodiment described above is only It is only illustrative, for example, in addition the division of the unit, only a kind of division of logic function can have in actual implementation Dividing mode, such as multiple units or component may be combined or can be integrated into another system or some features can be with Ignore or do not perform.Another point, shown or discussed mutual coupling, direct-coupling or communication connection can be logical Some interfaces are crossed, the INDIRECT COUPLING or communication connection of device or unit can be electrical, machinery or other forms.
It is also known in the art that other than realizing controller in a manner of pure computer readable program code, it is complete Entirely can by by method and step carry out programming in logic come controller with logic gate, switch, application-specific integrated circuit, may be programmed The form of logic controller and embedded microcontroller etc. realizes identical function.Therefore this controller is considered one kind Hardware component, and the structure that can also be considered as to the device for being used to implement various functions that its inside includes in hardware component.Or The device for being used to implement various functions even, can be considered as either the software module of implementation method can be hardware again by person Structure in component.
The present invention be with reference to according to the method for the embodiment of the present invention, the flow of equipment (system) and computer program product Figure and/or block diagram describe.It should be understood that it can be realized by computer program instructions every first-class in flowchart and/or the block diagram The combination of flow and/or box in journey and/or box and flowchart and/or the block diagram.These computer programs can be provided The processor of all-purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices is instructed to produce A raw machine so that the instruction performed by computer or the processor of other programmable data processing devices is generated for real The device of function specified in present one flow of flow chart or one box of multiple flows and/or block diagram or multiple boxes.
These computer program instructions, which may also be stored in, can guide computer or other programmable data processing devices with spy Determine in the computer-readable memory that mode works so that the instruction generation being stored in the computer-readable memory includes referring to Enable the manufacture of device, the command device realize in one flow of flow chart or multiple flows and/or one box of block diagram or The function of being specified in multiple boxes.
These computer program instructions can be also loaded into computer or other programmable data processing devices so that counted Series of operation steps are performed on calculation machine or other programmable devices to generate computer implemented processing, so as in computer or The instruction offer performed on other programmable devices is used to implement in one flow of flow chart or multiple flows and/or block diagram one The step of function of being specified in a box or multiple boxes.
In a typical configuration, computing device includes one or more processors (CPU), input/output interface, net Network interface and memory.
It should also be noted that, term " comprising ", "comprising" or its any other variant are intended to nonexcludability Comprising so that process, method, commodity or equipment including a series of elements are not only including those elements, but also wrap Include other elements that are not explicitly listed or further include for this process, method, commodity or equipment it is intrinsic will Element.In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that wanted including described Also there are other identical elements in the process of element, method or equipment.
It will be understood by those skilled in the art that this specification one or more embodiment can be provided as method, system or calculating Machine program product.Therefore, this specification one or more embodiment can be used complete hardware embodiment, complete software embodiment or With reference to the form of the embodiment in terms of software and hardware.Moreover, this specification one or more embodiment can be used at one or It is multiple wherein include computer usable program code computer-usable storage medium (include but not limited to magnetic disk storage, CD-ROM, optical memory etc.) on the form of computer program product implemented.
This specification one or more embodiment can computer executable instructions it is general on Described in hereafter, such as program module.Usually, program module includes performing particular task or realizes particular abstract data type Routine, program, object, component, data structure etc..This this specification one can also be put into practice in a distributed computing environment Or multiple embodiments, in these distributed computing environment, by being held by communication network and connected remote processing devices Row task.In a distributed computing environment, program module can be located at the local and remote computer including storage device In storage medium.
Each embodiment in this specification is described by the way of progressive, identical similar portion between each embodiment Point just to refer each other, and the highlights of each of the examples are difference from other examples.Especially for system reality For applying example, since it is substantially similar to embodiment of the method, so description is fairly simple, related part is referring to embodiment of the method Part explanation.In the description of this specification, reference term " one embodiment ", " some embodiments ", " example ", The description of " specific example " or " some examples " etc. means to combine specific features, structure, material that the embodiment or example describe Or feature is contained at least one embodiment or example of this specification.In the present specification, to the signal of above-mentioned term Property statement must not necessarily be directed to identical embodiment or example.Moreover, specific features, structure, material or the spy of description Point may be combined in any suitable manner in any one or more of the embodiments or examples.In addition, without conflicting with each other, Those skilled in the art can be by the different embodiments or examples described in this specification and different embodiments or examples Feature is combined.
The foregoing is merely embodiments herein, are not limited to the application.For those skilled in the art For, the application can have various modifications and variations.All any modifications made within spirit herein and principle are equal Replace, improve etc., it should be included within the scope of claims hereof.

Claims (10)

1. a kind of controlled source frequency dividing while exciting method, which is characterized in that including:
Targeted scans signal based on preset division number is divided, obtains sub- scanning signal;
The ramp function of sub- scanning signal is built based on the first preset rules, first preset rules include:Adjacent sub- scanning It is complementary relationship that signal, which overlaps position in the ramp function of frequency domain,;
Processing is optimized to the slope length of sub- scanning signal based on the ramp function, the sub- scanning signal is made to meet the Two preset rules, the sub- scanning signal after being optimized, wherein, second preset rules include:Adjacent sub- scanning signal is folded The energy fluctuation for closing position meets preset precision conditions;
The controlled source of preset number is carried out at the same time seismic wave excitation based on the sub- scanning signal after the optimization.
2. controlled source according to claim 1 frequency dividing exciting method simultaneously, which is characterized in that described to be based on first default Rule builds the ramp function of sub- scanning signal, including:
The ramp function of time-domain is obtained according to the ramp function of sub- scanning signal frequency domain, the ramp function of the time-domain is The evolution of frequency domain ramp function.
3. controlled source according to claim 1 frequency dividing exciting method simultaneously, which is characterized in that described to be based on first default Rule builds the ramp function of sub- scanning signal, including:
The energy fluctuation value at adjacent sub- scanning signal overlapping position is calculated based on a plurality of types of ramp functions, by least energy wave The dynamic corresponding ramp function type of value is determined as the ramp function type of the sub- scanning signal;
The ramp function of sub- scanning signal is built based on the first preset rules according to determining sub- scanning signal ramp function type.
4. controlled source according to claim 1 frequency dividing exciting method simultaneously, which is characterized in that the preset number can After control focus is carried out at the same time seismic wave excitation based on the sub- scanning signal after the optimization, further include:
The controlled source of preset number is carried out at the same time seismic wave excitation based on the sub- scanning signal after the optimization, and it is female to obtain earthquake Record data;
Each sub- scanning signal of excitation simultaneously or corresponding force signal are recorded into data with corresponding earthquake mother and carry out data point From processing, recorded in data from earthquake mother and isolate the corresponding earthquake subrecord data of each sub- scanning signal;
The corresponding earthquake subrecord data of complete scan sequence of single shot point are obtained, are believed according to sub- scanning signal in targeted scans The earthquake subrecord data of single shot point are sequentially overlapped by the sequence in number, obtain single-shot earthquake record.
5. a kind of controlled source frequency dividing while excitation apparatus, which is characterized in that including:
Frequency division module for targeted scans signal to be divided based on preset division number, obtains sub- scanning signal;
Ramp function builds module, and for building the ramp function of sub- scanning signal based on the first preset rules, described first is pre- If rule includes:Adjacent sub- scanning signal overlapping position is complementary relationship in the ramp function of frequency domain;
Optimization module for optimizing processing to the slope length of sub- scanning signal based on the ramp function, makes the son Scanning signal the second preset rules of satisfaction, the sub- scanning signal after being optimized, wherein, second preset rules include:Phase The energy fluctuation at adjacent sub- scanning signal overlapping position meets preset precision conditions;
Excitation module is carried out at the same time seismic wave based on the sub- scanning signal after the optimization and swashs for the controlled source of preset number Hair.
6. controlled source frequency dividing according to claim 5 while excitation apparatus, which is characterized in that the ramp function structure Module includes:
The ramp function of time-domain is obtained according to the ramp function of sub- scanning signal frequency domain, the ramp function of the time-domain is The evolution of frequency domain ramp function.
7. controlled source frequency dividing according to claim 5 while excitation apparatus, which is characterized in that the ramp function structure Module includes:
Ramp function type optimizes unit, and adjacent sub- scanning signal overlapping position is calculated for being based on a plurality of types of ramp functions Energy fluctuation value, the corresponding ramp function type of least energy undulating value is determined as to the ramp function of the sub- scanning signal Type;
Ramp function construction unit, for being built according to determining sub- scanning signal ramp function type based on the first preset rules The ramp function of sub- scanning signal.
8. controlled source frequency dividing according to claim 5 while excitation apparatus, which is characterized in that described device further includes:
Data capture unit is carried out at the same time earthquake for the controlled source of preset number based on the sub- scanning signal after the optimization Wave excites, and obtains earthquake mother and records data;
Data separating unit, for each sub- scanning signal of excitation or corresponding force signal simultaneously to be remembered with corresponding earthquake mother It records data and carries out data separating processing, recorded in data from earthquake mother and isolate the corresponding earthquake subrecord number of each sub- scanning signal According to;
Data investigation unit for obtaining the corresponding earthquake subrecord data of the complete scan sequence of single shot point, is swept according to son It retouches sequence of the signal in targeted scans signal to be sequentially overlapped the earthquake subrecord data of single shot point, with obtaining single-shot Shake record.
9. controlled source frequency dividing according to claim 5 while excitation apparatus, which is characterized in that including processor and be used for The memory of processor-executable instruction is stored, realizes and includes the following steps when described instruction is performed by the processor:
Targeted scans signal based on preset division number is divided, obtains sub- scanning signal;
The ramp function of sub- scanning signal is built based on the first preset rules, first preset rules include:Adjacent sub- scanning It is complementary relationship that signal, which overlaps position in the ramp function of frequency domain,;
Processing is optimized to the slope length of sub- scanning signal based on the ramp function, the sub- scanning signal is made to meet the Two preset rules, the sub- scanning signal after being optimized, wherein, second preset rules include:Adjacent sub- scanning signal is folded The energy fluctuation for closing position meets preset precision conditions;
The controlled source of preset number is carried out at the same time seismic wave excitation based on the sub- scanning signal after the optimization.
10. a kind of controlled source frequency dividing while activating system, which is characterized in that including multiple while excitation controlled source, inspection Wave device, at least one processor and the memory for storing computer executable instructions, when the processor performs described instruction The step of realizing any one the method in claim 1-4.
CN201711190612.1A 2017-11-24 2017-11-24 A kind of controlled source frequency dividing while exciting method, apparatus and system Pending CN108181645A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711190612.1A CN108181645A (en) 2017-11-24 2017-11-24 A kind of controlled source frequency dividing while exciting method, apparatus and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711190612.1A CN108181645A (en) 2017-11-24 2017-11-24 A kind of controlled source frequency dividing while exciting method, apparatus and system

Publications (1)

Publication Number Publication Date
CN108181645A true CN108181645A (en) 2018-06-19

Family

ID=62545220

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711190612.1A Pending CN108181645A (en) 2017-11-24 2017-11-24 A kind of controlled source frequency dividing while exciting method, apparatus and system

Country Status (1)

Country Link
CN (1) CN108181645A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112255681A (en) * 2020-10-26 2021-01-22 中国石油天然气集团有限公司 Vibroseis frequency reduction scanning data processing method and device
CN112415579A (en) * 2020-11-03 2021-02-26 中国石油天然气集团有限公司 Simultaneous source random excitation method, system and device
CN112859150A (en) * 2019-11-28 2021-05-28 中国石油天然气集团有限公司 Method and device for generating scanning signals of vibroseis

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120147699A1 (en) * 2010-12-10 2012-06-14 Bp Corporation North America Inc. Distance- and frequency-separated swept-frequency seismic sources
CN103605154A (en) * 2013-11-06 2014-02-26 中国石油集团东方地球物理勘探有限责任公司 Vibroseis frequency division simultaneous excitation method
CN104950326A (en) * 2015-03-04 2015-09-30 中石化石油工程地球物理有限公司胜利分公司 Target-stratum-frequency-spectrum-based design method of non-linear scanning signal of vibroseis
CN105204068A (en) * 2014-06-26 2015-12-30 中石化石油工程地球物理有限公司胜利分公司 Actual data spectrum characteristic-based nonlinear scanning signal design method
CN105572723A (en) * 2014-10-14 2016-05-11 中石化石油工程地球物理有限公司胜利分公司 Autocorrelation-wavelet-based design method for scanning signal of vibroseis

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120147699A1 (en) * 2010-12-10 2012-06-14 Bp Corporation North America Inc. Distance- and frequency-separated swept-frequency seismic sources
CN103605154A (en) * 2013-11-06 2014-02-26 中国石油集团东方地球物理勘探有限责任公司 Vibroseis frequency division simultaneous excitation method
CN105204068A (en) * 2014-06-26 2015-12-30 中石化石油工程地球物理有限公司胜利分公司 Actual data spectrum characteristic-based nonlinear scanning signal design method
CN105572723A (en) * 2014-10-14 2016-05-11 中石化石油工程地球物理有限公司胜利分公司 Autocorrelation-wavelet-based design method for scanning signal of vibroseis
CN104950326A (en) * 2015-03-04 2015-09-30 中石化石油工程地球物理有限公司胜利分公司 Target-stratum-frequency-spectrum-based design method of non-linear scanning signal of vibroseis

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
骆飞 等: "可控震源分频同时扫描技术", 《石油地球物理勘探》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112859150A (en) * 2019-11-28 2021-05-28 中国石油天然气集团有限公司 Method and device for generating scanning signals of vibroseis
CN112255681A (en) * 2020-10-26 2021-01-22 中国石油天然气集团有限公司 Vibroseis frequency reduction scanning data processing method and device
CN112415579A (en) * 2020-11-03 2021-02-26 中国石油天然气集团有限公司 Simultaneous source random excitation method, system and device

Similar Documents

Publication Publication Date Title
Savioja et al. Audio signal processing using graphics processing units
CN108205155B (en) A kind of controlled source alternating independent synchronization exciting method, apparatus and system
CN108181645A (en) A kind of controlled source frequency dividing while exciting method, apparatus and system
CN106855640B (en) A kind of submarine cable seismic data first arrival calculation method and device
CN108181646B (en) A kind of controlled source exciting method, apparatus and system simultaneously
Sadhu An integrated multivariate empirical mode decomposition method towards modal identification of structures
CN109471171A (en) A kind of method, apparatus and system of aliased seismic data separation
CN108387933A (en) A kind of method, apparatus and system of definitely interval quality factors
CN109100795A (en) A kind of inspection point distribution method, the apparatus and system of face element
CN108107471B (en) The acquisition methods and device of a kind of point of orientation first arrival data volume
Stott et al. Gravitational wave bursts from cosmic string cusps and pseudocusps
CN108828668A (en) A kind of pre-stack time migration data processing method and device
Khan et al. Recent open source wireless sensor network supporting simulators: A performance comparison
Dreger et al. Introduction to the focus section on validation of the SCEC broadband platform V14. 3 simulation methods
CN108108707A (en) Gesture processing method and processing device based on video data, computing device
Fichtner et al. Seismic noise correlation on heterogeneous supercomputers
CN106324666A (en) Extraction method and extraction device for extracting converted wave angle gather of transversely isotropic media
CN104181600A (en) Seismic data linear noise attenuation method and device
CN109709602A (en) A kind of remote detection sonic wave shift imaging method, apparatus and system
CN102651214B (en) Method, device and sound card for generating sounds in virtual world
CN108469633B (en) The calculation method and device of a kind of ground interval quality factors
Behar et al. Resonant Whistler‐Electron Interactions: MMS Observations Versus Test‐Particle Simulation
CN108828669A (en) A kind of two-dimensional intersection survey line static corrections processing method, apparatus and system
CN112649887B (en) Method and device for quantitatively dividing stratum sequence based on drilling data
CN106094020A (en) A kind of seismic inversion method and device

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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20180619