CN112051613A - Aliasing acquisition chatter time generation method and device - Google Patents
Aliasing acquisition chatter time generation method and device Download PDFInfo
- Publication number
- CN112051613A CN112051613A CN202010766404.7A CN202010766404A CN112051613A CN 112051613 A CN112051613 A CN 112051613A CN 202010766404 A CN202010766404 A CN 202010766404A CN 112051613 A CN112051613 A CN 112051613A
- Authority
- CN
- China
- Prior art keywords
- time
- value
- seismic source
- source excitation
- initial
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/37—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy specially adapted for seismic systems using continuous agitation of the ground, e.g. using pulse compression of frequency swept signals for enhancement of received signals
- G01V1/375—Correlating received seismic signals with the emitted source signal
Abstract
The invention provides a method and a device for generating aliasing acquisition chatter time, which relate to the technical field of geophysical exploration, and comprise the following steps: acquiring a value range parameter, a value precision parameter and an initial vibration time value; determining a plurality of to-be-selected flutter time values according to the value range parameter and the value precision parameter; determining the selection probability of each to-be-selected vibration time value corresponding to the target seismic source excitation time according to the value range parameter, the value precision parameter and the initial vibration time value; and generating aliasing acquisition chatter time corresponding to the target seismic source excitation time according to the selection probability. The method can improve the randomness of the adjacent flutter time difference corresponding to the target seismic source excitation time based on probability weighting, and provides support for obtaining better seismic data separation effect.
Description
Technical Field
The invention relates to the technical field of geophysical exploration, in particular to a method and a device for generating aliasing acquisition chatter time.
Background
The high-efficiency aliasing acquisition technology greatly improves the daily efficiency of seismic acquisition and is a technology widely adopted by the high-precision high-density seismic exploration at present. However, the high-efficiency mixed data has serious adjacent shot interference noise, and the mixed data needs to be separated subsequently. At present, the mainstream mixed data separation method is based on the characteristic that effective signals of mixed data are continuous in a certain domain and aliasing noise is randomly distributed. Therefore, the more the difference in firing time between the current shot and the adjacent shot tends to be uniformly randomly distributed, the more data separation is facilitated. In order to make the excitation time of each cannon generate irregular change, a method of adding random flutter time to the excitation time of each cannon is generally adopted at home and abroad at present, the method only aims at the flutter time, and although the single flutter time is uniform and random in a specified range, the difference is not uniform and random, and a better data separation effect cannot be ensured.
Disclosure of Invention
The invention provides a method and a device for generating aliasing acquisition chattering time, which can realize uniform random distribution of adjacent chattering time difference values and ensure a better data separation effect.
In a first aspect, an embodiment of the present invention provides an aliasing acquisition dither time generation method, where the method includes: acquiring a value range parameter, a value precision parameter and an initial vibration time value; determining a plurality of jitter time values to be selected according to the value range parameter and the value precision parameter; determining the selection probability of each to-be-selected flutter time value corresponding to the target seismic source excitation moment according to the value range parameter, the value precision parameter and the initial flutter time value; and generating aliasing acquisition flutter time corresponding to the target seismic source excitation time according to the selection probability.
In a second aspect, an embodiment of the present invention further provides an aliasing acquisition jitter time generation apparatus, including: the acquisition module is used for acquiring a value range parameter, a value precision parameter and an initial vibration time value; the determining module is used for determining a plurality of to-be-selected flutter time values according to the value range parameter and the value precision parameter; the probability module is used for determining the selection probability of each to-be-selected flutter time value corresponding to the target seismic source excitation moment according to the value range parameter, the value precision parameter and the initial flutter time value; and the generating module is used for generating aliasing acquisition chatter time corresponding to the target seismic source excitation time according to the selection probability.
In a third aspect, an embodiment of the present invention further provides a computer device, which includes a memory, a processor, and a computer program stored on the memory and executable on the processor, where the processor implements the above aliasing acquisition dithering time generation method when executing the computer program.
In a fourth aspect, an embodiment of the present invention further provides a computer-readable storage medium, where a computer program for executing the above aliasing acquisition dithering time generation method is stored.
The embodiment of the invention has the following beneficial effects: the embodiment of the invention provides an aliasing acquisition flutter time generation scheme, which comprises the steps of firstly obtaining a value range parameter, a value precision parameter and an initial flutter time value, determining a plurality of flutter time values to be selected according to the value range parameter and the value precision parameter, then determining the selection probability of each flutter time value to be selected corresponding to the target seismic source excitation time according to the value range parameter, the value precision parameter and the initial flutter time value, and then generating the aliasing acquisition flutter time corresponding to the target seismic source excitation time according to the selection probability. The method and the device can improve the randomness of the adjacent flutter time difference corresponding to the target seismic source excitation time based on probability weighting, and provide support for obtaining a better seismic data separation effect.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
In order to make the aforementioned and other objects, features and advantages of the present invention comprehensible, preferred embodiments accompanied with figures are described in detail below.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
Fig. 1 is a flowchart of an aliasing acquisition dither time generation method according to an embodiment of the present invention;
FIG. 2 is a graph of jitter time difference distribution for a conventional stochastic method according to an embodiment of the present invention;
fig. 3 is a schematic diagram illustrating a jitter time range and a combination of adjacent jitter time ranges according to an embodiment of the present invention;
fig. 4 is a schematic diagram of a value law (linear equation set) of adjacent flutter time difference provided in the embodiment of the present invention;
FIG. 5 is a general flowchart of a probability weighted random method according to an embodiment of the present invention;
fig. 6 is a histogram of jitter time when a probability weighted random method is adopted according to an embodiment of the present invention;
fig. 7 is a histogram of the jitter time difference when the probability weighted random method is adopted according to the embodiment of the present invention;
FIG. 8 is a graph illustrating the effect of the jitter free time provided by an embodiment of the present invention;
FIG. 9 is a diagram illustrating the effect of a conventional stochastic method according to an embodiment of the present invention;
FIG. 10 is a diagram illustrating the effect of the probability weighted random method according to the embodiment of the present invention;
fig. 11 is a block diagram of a structure of an aliasing acquisition jitter time generation apparatus according to an embodiment of the present invention;
fig. 12 is a block diagram of another aliasing acquisition dither time generation apparatus according to an embodiment of the present invention;
fig. 13 is a block diagram of a computer device according to an embodiment of the present invention.
Detailed Description
To make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is apparent that the described embodiments are some, but not all embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
When the marine seismic source ship operates, although the ship speed is changed due to the influence of the external environment, the mass of the hull of the seismic source ship is too large, and the external influence is relatively limited. Therefore, from the general trend, the ship speed is still relatively constant, and the time difference between shots also tends to a certain value (for example, the ship speed is stabilized at 2.5m/s, the distance between shots is 12.5m, and the time difference between shots is basically distributed near 5 s). The centralized distribution rule does not meet the requirement of the mixed data separation method on the random distribution of time difference between shots. Therefore, random flutter time needs to be added to the excitation time of each cannon, so that the excitation time of each cannon generates irregular change.
The fundamental purpose of the dither time is to excite the randomness of the time instants. Let TnIs the excitation time of the nth shot, dnFor the random jitter time of the nth shot, the inter-shot time difference formula is as follows:
ΔT＝(Tn+dn)-(Tn-1+dn-1)＝(Tn-Tn-1)+(dn-dn-1)
it can be seen that in the case of poor natural randomness, the difference between the original excitation times varies little, so the main factor affecting the randomness of the final excitation times is the difference between adjacent dither times (d)n-dn-1). At present, the flutter time random method commonly adopted at home and abroad only aims at the flutter time, although the single flutter time is uniform and random within a specified range, the difference is not uniform and random, but presents isosceles triangle distribution gradually decreasing from the middle to two sides, such as a flutter time difference distribution diagram in the conventional random method shown in fig. 2, wherein the flutter time value range is [ -500ms, 500ms]The sampling interval is 2ms, the number of sampling points is 10 tens of thousands, and the existing method cannot ensure the optimal data separation effect.
Based on this, the aliasing acquisition chatter time generation method and device provided by the embodiment of the invention can realize uniform random distribution of adjacent chatter time difference values.
For the convenience of understanding the present embodiment, a detailed description will be given first of all on an aliasing acquisition dithering time generation method disclosed in the present embodiment.
The embodiment of the invention provides an aliasing acquisition jittering time generation method, which is shown in a flow chart of the aliasing acquisition jittering time generation method shown in figure 1 and comprises the following steps:
and S102, acquiring a value range parameter, a value precision parameter and an initial flutter time value.
In the embodiment of the invention, the value range parameter is used for determining the time length range of the offset of the aliasing acquisition vibration time relative to the target seismic source excitation time. The numeric precision parameter is used to determine the time duration interval between optional aliased acquisition dither times. The initial dither time value may be generated in advance.
And step S104, determining a plurality of jitter time values to be selected according to the value range parameter and the value precision parameter.
In the embodiment of the invention, after the value range parameter is determined, the selectable range of the jitter time value to be selected can be determined, and a plurality of jitter time values to be selected can be determined in the selectable range according to the value precision parameter. For example, if the value range parameter includes plus or minus 500 milliseconds, the selectable range of the to-be-selected chattering time value is a value between minus 500 milliseconds and plus 500 milliseconds, and if the value precision parameter is 2 milliseconds, the to-be-selected chattering time value includes 500, 498, 496, … …, minus 498, and minus 500 milliseconds.
And S106, determining the selection probability of each to-be-selected vibration time value corresponding to the target seismic source excitation time according to the value range parameter, the value precision parameter and the initial vibration time value.
In the embodiment of the present invention, the target seismic source excitation time may include a plurality of seismic source excitation times, and for each seismic source excitation time, a group of data is determined according to a value range parameter, a value precision parameter, and an initial flutter time value, where the group of data includes: and each chattering time value to be selected and the selection probability corresponding to each chattering time value to be selected.
And S108, generating aliasing acquisition chatter time corresponding to the target seismic source excitation time according to the selection probability.
In the embodiment of the invention, the selection probability of the to-be-selected trembling time value is used as the selection weight, and for each target seismic source excitation time, one value is selected from a plurality of to-be-selected trembling time values corresponding to the target seismic source excitation time to obtain aliasing acquisition trembling time corresponding to the target seismic source excitation time, and further obtain a trembling time sequence corresponding to a plurality of target seismic sources.
The embodiment of the invention provides an aliasing acquisition flutter time generation scheme, which comprises the steps of firstly obtaining a value range parameter, a value precision parameter and an initial flutter time value, determining a plurality of flutter time values to be selected according to the value range parameter and the value precision parameter, then determining the selection probability of each flutter time value to be selected corresponding to the target seismic source excitation time according to the value range parameter, the value precision parameter and the initial flutter time value, and then generating the aliasing acquisition flutter time corresponding to the target seismic source excitation time according to the selection probability. The method and the device can improve the randomness of the adjacent flutter time difference corresponding to the target seismic source excitation time based on probability weighting, and provide support for obtaining a better seismic data separation effect.
The target seismic source excitation time comprises an initial seismic source excitation time, and the following steps can be executed before the value range parameter, the value precision parameter and the initial vibration time value are acquired in order to improve the calculation efficiency:
randomly generating a vibration time value corresponding to the initial seismic source excitation time according to the value range parameter; and taking the vibration time value corresponding to the initial seismic source excitation time as an initial vibration time value.
In the embodiment of the present invention, when acquiring the jitter time sequence, the first jitter time value is acquired by a conventional random method. The initial source activation time is an activation time corresponding to the source that is activated earliest among the plurality of target sources. And taking the vibration time value corresponding to the initial seismic source excitation time as an initial vibration time value.
Considering that in order to improve randomness, the selection probability of each to-be-selected dither time value corresponding to the target seismic source excitation time is determined according to the value range parameter, the value precision parameter and the initial dither time value, and the method can be executed according to the following steps.
If the target seismic source excitation time is a second seismic source excitation time adjacent to the initial seismic source excitation time, determining the selection probability of each to-be-selected vibration time value corresponding to the second seismic source excitation time according to the value range parameter, the value precision parameter and the initial vibration time value; and if the target seismic source excitation time is the target seismic source excitation time except the initial seismic source excitation time and the second seismic source excitation time, determining the selection probability of each to-be-selected vibration time value corresponding to the target seismic source excitation time according to the value range parameter, the value precision parameter and the aliasing acquisition vibration time corresponding to the last target seismic source excitation time.
The selection probability of each candidate tremor time value corresponding to the target seismic source excitation time can be determined according to the following formula:
wherein p isdIn order to select the probability, d is the chattering time value to be selected, d has S +1 possible values, S is a variable from 0 to S, dprevIs the initial dither time valueOr last dither time value, dstepTo obtain a value precision parameter, dmaxIs a value range parameter.
Referring to a general flow chart shown in fig. 5 when the probability weighted random method is adopted, in the embodiment of the present invention, when a tremor time sequence is obtained, a first tremor time value is obtained by a conventional random method, then all tremor times are obtained by a probability weighted random method, and each time, the last tremor time is used as a parameter to calculate a weight of each value of the current tremor time, so as to obtain a probability of each value. And finally, randomly sampling by a computer to obtain the current flutter time, and repeating the steps in a circulating way until the number of the flutter time sequences meets the requirement.
The method may further perform the steps of: and adjusting the target seismic source excitation time according to the aliasing acquisition vibration time so that the target seismic source is excited according to the adjusted target seismic source excitation time.
In the embodiment of the invention, the aliasing acquisition flutter time is a time offset, and the target seismic source excitation time is adjusted according to the aliasing acquisition flutter time, so that the adjusted target seismic source excitation time can be obtained. After obtaining the aliasing acquisition vibration time, for example, the aliasing acquisition vibration time corresponding to the target seismic source a, the target seismic source b and the target seismic source c is 5 seconds, 8 seconds and minus 2 seconds in sequence, the excitation time corresponding to the target seismic source a, the target seismic source b and the target seismic source c is 13 points, 14 points and 15 points in sequence, and then the excitation time of the post-seismic source is adjusted to 13 points for 5 seconds, 14 points for 8 seconds and 2 seconds for 15 points.
The present invention can be realized as follows.
1) Analyzing a flutter time value-taking rule: and analyzing the flutter time value-taking principle and a mathematical model.
2) Analyzing the value law of the flutter time difference: and analyzing the flutter time difference value-taking principle and a mathematical model.
3) And (3) analyzing the probability distribution of the flutter time difference: and determining the probability distribution rule of the vibration according to the distribution mathematical model of the vibration time difference.
4) And (3) calculating the weight distribution of the dithering time difference value: and adjusting the weight according to the probability distribution rule of the flutter time difference value, so that the equal probability of different values of the flutter time difference value is realized.
5) And (3) calculating the fibrillation time weighted probability: and converting the weight value of the single sampling point in the linear equation determined by the vibrating time difference value into the integral weight and probability when the vibrating time takes different values.
6) Cyclic acquisition of dithered time series: and acquiring all the tremor time sequences according to a set flow loop.
7) And (3) analyzing results of the flutter time and the flutter time difference value: the final tremor time series was analyzed and evaluated.
8) And (3) analyzing the time difference result between the guns: and analyzing and evaluating the final shot-to-shot time difference distribution.
The implementation process of the method is described below with reference to the drawings and formulas.
1) Analyzing a flutter time value-taking rule:
if the absolute maximum value of the flutter time is dmaxThe value precision is dstep。
The value range is-dmaxTo dmaxIs provided with
Taking the current chattering time as the vertical axis y and the adjacent previous chattering time as the horizontal axis x, see fig. 3 showing the combination diagram of chattering time value range and adjacent chattering time value, where d is an example shown in the diagrammax＝500，d step100. All points within the black square region represent all possible combinations of values for (x, y), which exhibit a uniform random distribution within the square.
2) Analyzing the value law of the flutter time difference:
and if the difference value between the current dithering time and the last dithering time is delta d, then: and d is y-x.
The two-dimensional coordinate axis is a straight line with a slope of 1, the intersection point of the straight line and the y axis is equal to deltad, and all grid points in the square area penetrated by the straight line are all possible (x, y) value combinations with the difference value equal to deltad.
Delta d has a value range of-2 dmaxTo 2dmaxThere are (2S +1) possible values. Referring to a schematic diagram of a value law (a linear equation set) of adjacent flutter time difference shown in fig. 4, as shown by oblique lines in the diagram, each value corresponds to a unique linear equation, and the set of all equations is:
3) and (3) analyzing the probability distribution of the flutter time difference:
from FIG. 4, the total number of all dots in the black square region is (S +1)2The grid points are uniformly and randomly distributed.
The number of grid points traversed by the straight line y + Δ d, i.e. all possible combinations of (x, y) values, is equal to:
probability of value of
Therefore, the method comprises the following steps:
when Δ d is 0, there are a total of S +1 possibilities with a probability of
Highest;
when Δ d is 2dmaxOr Δ d ═ 2dmaxWhen the probability is 1, the probability is
And the lowest.
4) And (3) calculating the weight distribution of the dithering time difference value:
due to the uniform random distribution of the grid points, the weight of the numeric probability of each grid point is equal to 1 by default.
The number of grid points that the straight line y + Δ d passes through varies with the value of Δ d, and all the straight lines are equal in probability of appearing and are connected toIndividual probability weights for passing grid points are set to
That is, the number is the inverse, the weight of the straight line is the weight of a single grid point multiplied by the number of grid points, the equal weight is 1, and the equal probability of all straight lines can be realized.
5) And (3) calculating the fibrillation time weighted probability:
when the last dithering time value x has been determined, there are (S +1) possibilities for the current dithering time value y, corresponding to (S +1) grid points on the two-dimensional coordinate axis. And combining the grid point single point weight formula calculated in the last step.
The value probability when y is any possible value can be obtained as follows:
replace x by the last dither time dprevAnd y is replaced by any possible value d of the current chattering time, and the value probability of any value d is as follows:
wherein-dmax≤d≤dmaxD has (S +1) possible values.
6) Cyclic acquisition of dithered time series:
when a flutter time sequence is obtained, the first flutter time value is still obtained by a conventional random method, then all the flutter times are obtained by a probability weighted random method, the last flutter time is taken as a parameter each time to calculate the weight of each value of the current flutter time, and then the probability of each value is obtained, finally, the current flutter time is obtained by computer random sampling, and the steps are repeated in this way until the number of the flutter time sequence meets the requirement. The flow is shown in fig. 5.
7) And analyzing results of the flutter time and the flutter time difference value.
The final flutter time distribution is a rule of more two sides and less middle, but the whole is symmetrical, and the average value approaches to 0.
The distribution of jitter time differences is consistent with a uniform random signature, which is the ultimate goal of the method. Fig. 6 and 7 are a histogram of the jitter time and a histogram of the jitter time difference when the probability weighted random method is used (in the figure, the jitter time ranges from-500 ms to 500ms, the sampling interval is 2ms, and the number of sampling points is 10 ten thousand).
8) And analyzing the time difference result between the shots.
Referring to the effect graph of the non-tremor time shown in fig. 8, the effect graph of the conventional random method shown in fig. 9, and the effect graph of the probability weighted random method shown in fig. 10, compared with the conventional random method, the probability weighted random method has a more obvious improvement effect on the time difference between shots, and the time difference between shots becomes more uniform and random, so that the distribution requirement of the mixed-mining seismic data separation on the time difference between shots is better met.
The embodiment of the invention provides a method and a device for generating aliasing acquisition flutter time, wherein the method can generate flutter time sequences required by aliasing acquisition of a marine seismic exploration multi-source ship, the flutter time sequences generated by the method can keep symmetrical and random distribution, and adjacent difference values of the flutter time sequences generated by the method can keep symmetrical, random and uniform distribution.
The embodiment of the present invention further provides an aliasing acquisition jitter time generation apparatus, as described in the following embodiments. Because the principle of the device for solving the problems is similar to the aliasing acquisition dithering time generation method, the implementation of the device can refer to the implementation of the aliasing acquisition dithering time generation method, and repeated parts are not described again. Referring to fig. 11, a block diagram of an aliasing acquisition dither time generation apparatus is shown, which includes:
an obtaining module 71, configured to obtain a value range parameter, a value precision parameter, and an initial jitter time value; the determining module 72 is configured to determine a plurality of to-be-selected jitter time values according to the value range parameter and the value precision parameter; the probability module 73 is configured to determine, according to the value range parameter, the value precision parameter, and the initial tremor time value, a selection probability of each to-be-selected tremor time value corresponding to the target seismic source excitation time; and the generating module 74 is configured to generate an aliasing acquisition tremor time corresponding to the target seismic source excitation time according to the selection probability.
In one embodiment, referring to another structural block diagram of an aliasing acquisition dither time generation apparatus shown in fig. 12, the target source firing time comprises an initial source firing time, and the apparatus further comprises a randomizer module 75 for: randomly generating a vibration time value corresponding to the initial seismic source excitation time according to the value range parameter; and taking the vibration time value corresponding to the initial seismic source excitation time as an initial vibration time value.
In one embodiment, the probability module is specifically configured to: if the target seismic source excitation time is a second seismic source excitation time adjacent to the initial seismic source excitation time, determining the selection probability of each to-be-selected vibration time value corresponding to the second seismic source excitation time according to the value range parameter, the value precision parameter and the initial vibration time value; and if the target seismic source excitation time is the target seismic source excitation time except the initial seismic source excitation time and the second seismic source excitation time, determining the selection probability of each to-be-selected vibration time value corresponding to the target seismic source excitation time according to the value range parameter, the value precision parameter and the aliasing acquisition vibration time corresponding to the last target seismic source excitation time.
In one embodiment, the probability module is specifically configured to: determining the selection probability of each candidate tremor time value corresponding to the target seismic source excitation time according to the following formula:
wherein p isdIn order to select the probability, d is the chattering time value to be selected, d has S +1 possible values, S is a variable from 0 to S, dprevIs the initial or last dither time value, dstepTo obtain a value precision parameter, dmaxIs a value range parameter.
In one embodiment, referring to another structural block diagram of the aliasing acquisition dither time generation apparatus shown in fig. 12, the apparatus further includes an adjusting module 76 for adjusting the target seismic source excitation time according to the aliasing acquisition dither time, so that the target seismic source is excited according to the adjusted target seismic source excitation time.
An embodiment of the present invention further provides a computer device, referring to the schematic block diagram of the structure of the computer device shown in fig. 13, the computer device includes a memory 81, a processor 82, and a computer program stored in the memory and executable on the processor, and when the processor executes the computer program, the processor implements any of the steps of the above-mentioned aliasing acquisition jitter time generation method.
It can be clearly understood by those skilled in the art that, for convenience and brevity of description, the specific working process of the computer device described above may refer to the corresponding process in the foregoing method embodiment, and is not described herein again.
An embodiment of the present invention further provides a computer-readable storage medium, in which a computer program for executing any one of the above aliasing acquisition jitter time generation methods is stored.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
Finally, it should be noted that: the above-mentioned embodiments are only specific embodiments of the present invention, which are used for illustrating the technical solutions of the present invention and not for limiting the same, and the protection scope of the present invention is not limited thereto, although the present invention is described in detail with reference to the foregoing embodiments, those skilled in the art should understand that: any person skilled in the art can modify or easily conceive the technical solutions described in the foregoing embodiments or equivalent substitutes for some technical features within the technical scope of the present disclosure; such modifications, changes or substitutions do not depart from the spirit and scope of the embodiments of the present invention, and they should be construed as being included therein. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (12)
1. An aliasing acquisition flutter time generation method is applied to a seismic exploration multi-source ship, and comprises the following steps:
acquiring a value range parameter, a value precision parameter and an initial vibration time value;
determining a plurality of jitter time values to be selected according to the value range parameter and the value precision parameter;
determining the selection probability of each to-be-selected flutter time value corresponding to the target seismic source excitation moment according to the value range parameter, the value precision parameter and the initial flutter time value;
and generating aliasing acquisition flutter time corresponding to the target seismic source excitation time according to the selection probability.
2. The method of claim 1, wherein the target source firing time comprises an initial source firing time;
before obtaining the value range parameter, the value precision parameter and the initial flutter time value, the method further comprises the following steps:
randomly generating a vibration time value corresponding to the initial seismic source excitation time according to the value range parameter;
and taking the vibration time value corresponding to the initial seismic source excitation time as an initial vibration time value.
3. The method according to claim 2, wherein determining a selection probability of each to-be-selected tremor time value corresponding to a target seismic source excitation time according to the value range parameter, the value precision parameter, and the initial tremor time value comprises:
if the target seismic source excitation time is a second seismic source excitation time adjacent to the initial seismic source excitation time, determining the selection probability of each to-be-selected vibration time value corresponding to the second seismic source excitation time according to the value range parameter, the value precision parameter and the initial vibration time value;
and if the target seismic source excitation time is the target seismic source excitation time except the initial seismic source excitation time and the second seismic source excitation time, determining the selection probability of each to-be-selected vibration time value corresponding to the target seismic source excitation time according to the value range parameter, the value precision parameter and the aliasing acquisition vibration time corresponding to the last target seismic source excitation time.
4. The method of claim 3, comprising determining a probability of selection for each of the candidate tremor time values corresponding to a target source firing time according to the following formula:
wherein p isdIn order to select the probability, d is the chattering time value to be selected, d has S +1 possible values, S is a variable from 0 to S, dprevIs the initial or last dither time value, dstepTo obtain a value precision parameter, dmaxIs a value range parameter.
5. The method according to any one of claims 1-4, further comprising:
and adjusting the target seismic source excitation time according to the aliasing acquisition vibration time so that the target seismic source is excited according to the adjusted target seismic source excitation time.
6. An aliasing acquisition flutter time generation device is applied to a seismic exploration multi-source ship and comprises the following components:
the acquisition module is used for acquiring a value range parameter, a value precision parameter and an initial vibration time value;
the determining module is used for determining a plurality of to-be-selected flutter time values according to the value range parameter and the value precision parameter;
the probability module is used for determining the selection probability of each to-be-selected flutter time value corresponding to the target seismic source excitation moment according to the value range parameter, the value precision parameter and the initial flutter time value;
and the generating module is used for generating aliasing acquisition chatter time corresponding to the target seismic source excitation time according to the selection probability.
7. The apparatus of claim 6, wherein the target source firing time comprises an initial source firing time; the apparatus further comprises a randomization module to:
randomly generating a vibration time value corresponding to the initial seismic source excitation time according to the value range parameter;
and taking the vibration time value corresponding to the initial seismic source excitation time as an initial vibration time value.
8. The apparatus of claim 7, wherein the probability module is specifically configured to:
if the target seismic source excitation time is a second seismic source excitation time adjacent to the initial seismic source excitation time, determining the selection probability of each to-be-selected vibration time value corresponding to the second seismic source excitation time according to the value range parameter, the value precision parameter and the initial vibration time value;
and if the target seismic source excitation time is the target seismic source excitation time except the initial seismic source excitation time and the second seismic source excitation time, determining the selection probability of each to-be-selected vibration time value corresponding to the target seismic source excitation time according to the value range parameter, the value precision parameter and the aliasing acquisition vibration time corresponding to the last target seismic source excitation time.
9. The apparatus of claim 8, wherein the probability module is specifically configured to: determining the selection probability of each candidate tremor time value corresponding to the target seismic source excitation time according to the following formula:
wherein p isdIn order to select the probability, d is the chattering time value to be selected, d has S +1 possible values, S is a variable from 0 to S, dprevIs the initial or last dither time value, dstepTo obtain a value precision parameter, dmaxIs a value range parameter.
10. The apparatus of any of claims 6-9, further comprising an adjustment module to adjust the target source firing time according to the aliased acquisition dither time such that the target source is fired according to the adjusted target source firing time.
11. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method of any one of claims 1 to 5 when executing the computer program.
12. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program for executing the method of any one of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010766404.7A CN112051613B (en) | 2020-08-03 | 2020-08-03 | Method and device for generating aliasing acquisition dithering time |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010766404.7A CN112051613B (en) | 2020-08-03 | 2020-08-03 | Method and device for generating aliasing acquisition dithering time |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112051613A true CN112051613A (en) | 2020-12-08 |
| CN112051613B CN112051613B (en) | 2023-09-26 |
Family
ID=73602582
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010766404.7A Active CN112051613B (en) | 2020-08-03 | 2020-08-03 | Method and device for generating aliasing acquisition dithering time |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112051613B (en) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120176861A1 (en) * | 2011-01-12 | 2012-07-12 | Bp Corporation North America Inc. | Shot scheduling limits for seismic acquisition with simultaneous source shooting |
| US8706536B1 (en) * | 2007-01-09 | 2014-04-22 | David Alan McPhetrige | Systems and methods for estimating safety stock levels |
| US20140372044A1 (en) * | 2013-06-17 | 2014-12-18 | Westerngeco L.L.C. | Seismic data processing |
| US20150247940A1 (en) * | 2014-03-03 | 2015-09-03 | King Abdullah University Of Science And Technology | History matching of time-lapse crosswell data using ensemble kalman filtering |
| CN107966728A (en) * | 2016-10-19 | 2018-04-27 | 中国石油化工股份有限公司 | The earthquake collection method and device of the more Seismic Source Systems of dynamite source |
| GB201816844D0 (en) * | 2017-10-19 | 2018-11-28 | Cgg Services Sas | System and method for generating dithering sequences for seismic exploration |
| CN109471171A (en) * | 2018-09-21 | 2019-03-15 | 中国石油天然气集团有限公司 | A kind of method, apparatus and system of aliased seismic data separation |
| CN109581484A (en) * | 2018-11-01 | 2019-04-05 | 中国石油天然气集团有限公司 | P wave component multiple wave general ambient light degree index analysis method and device |
| US20200072998A1 (en) * | 2018-09-04 | 2020-03-05 | Sercel | System and method for generating dithering sequences with minimum value for seismic exploration |
-
2020
- 2020-08-03 CN CN202010766404.7A patent/CN112051613B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8706536B1 (en) * | 2007-01-09 | 2014-04-22 | David Alan McPhetrige | Systems and methods for estimating safety stock levels |
| US20120176861A1 (en) * | 2011-01-12 | 2012-07-12 | Bp Corporation North America Inc. | Shot scheduling limits for seismic acquisition with simultaneous source shooting |
| US20140372044A1 (en) * | 2013-06-17 | 2014-12-18 | Westerngeco L.L.C. | Seismic data processing |
| US20150247940A1 (en) * | 2014-03-03 | 2015-09-03 | King Abdullah University Of Science And Technology | History matching of time-lapse crosswell data using ensemble kalman filtering |
| CN107966728A (en) * | 2016-10-19 | 2018-04-27 | 中国石油化工股份有限公司 | The earthquake collection method and device of the more Seismic Source Systems of dynamite source |
| GB201816844D0 (en) * | 2017-10-19 | 2018-11-28 | Cgg Services Sas | System and method for generating dithering sequences for seismic exploration |
| US20200072998A1 (en) * | 2018-09-04 | 2020-03-05 | Sercel | System and method for generating dithering sequences with minimum value for seismic exploration |
| CN109471171A (en) * | 2018-09-21 | 2019-03-15 | 中国石油天然气集团有限公司 | A kind of method, apparatus and system of aliased seismic data separation |
| CN109581484A (en) * | 2018-11-01 | 2019-04-05 | 中国石油天然气集团有限公司 | P wave component multiple wave general ambient light degree index analysis method and device |
Non-Patent Citations (3)
| Title |
|---|
| FÁTIMA GOUVEIA ，RUI CARRILHO GOMES ， ISABEL LOPES: "Shallow and in depth seismic testing in urban environment: A case study in Lisbon Miocene stiff soils using joint inversion of active and passive Rayleigh wave measurements", JOURNAL OF APPLIED GEOPHYSICS * |
| 韩志雄;倪宇东;余建鹏;马涛;卢秀丽;张华;: "强环噪地区可控震源高效混叠采集应用效果", 地球物理学进展, no. 01 * |
| 马竹，侯昆鹏，柳兴刚等: "多震源船混叠釆集导航模拟方法及应用", pages 180 - 183 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112051613B (en) | 2023-09-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR102206324B1 (en) | Methods for ascertaining a model of a starting variable of a technical system | |
| CN110632674B (en) | Weak information extraction method for aviation gravity measurement data | |
| CA2969101A1 (en) | Seismic elastic wave simulation for tilted transversely isotropic media using adaptive lebedev staggered grid | |
| CN106842306A (en) | The staggered-mesh finite difference analogy method and device of a kind of global optimization | |
| CN112163335A (en) | Training method, prediction method and device of NOx concentration prediction model | |
| CN114578417B (en) | Pulse type earthquake motion fitting method based on wavelet adjustment | |
| CN110058276A (en) | Abnormal point judgment method and device | |
| KR20200049422A (en) | Effective Network Compression using Simulation-guided Iterative Pruning | |
| CN112051613A (en) | Aliasing acquisition chatter time generation method and device | |
| CN111738400A (en) | Iteration method for optimization of variational modal decomposition parameters of non-stationary signals | |
| Kalandros et al. | Randomization and super-heuristics in choosing sensor sets for target tracking applications | |
| US10902428B1 (en) | Maintaining a risk model using feedback directed to other risk models | |
| CN111435174B (en) | Method and device for compensating amplitude of seismic data in strong reflection area | |
| CN113411140B (en) | Channel self-adaptive channel characteristic learning method, simulator thereof and satellite navigation system | |
| JP2020024119A5 (en) | ||
| CN110376642B (en) | Three-dimensional seismic velocity inversion method based on conical surface waves | |
| CN113239636A (en) | Performance evaluation method, system, equipment and storage medium of weapon equipment | |
| CN108805319B (en) | Method and system for determining optimal state number of wind power modeling | |
| RU2629641C1 (en) | Digital predictor | |
| CN113391354A (en) | Aliasing noise suppression method and device | |
| CN111856552A (en) | Method and device for generating scanning signals of vibroseis | |
| CN111856557A (en) | Method and device for making depth domain synthetic seismic record | |
| RU97113716A (en) | FIRE MANAGEMENT SYSTEM | |
| SU903873A1 (en) | Generator of random numbers for simulating general population by objects of a sample | |
| CN114814930A (en) | Method and device for determining distortion of vibroseis force signal |
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 |