CN112051613B - Method and device for generating aliasing acquisition dithering time - Google Patents
Method and device for generating aliasing acquisition dithering time Download PDFInfo
- Publication number
- CN112051613B CN112051613B CN202010766404.7A CN202010766404A CN112051613B CN 112051613 B CN112051613 B CN 112051613B CN 202010766404 A CN202010766404 A CN 202010766404A CN 112051613 B CN112051613 B CN 112051613B
- Authority
- CN
- China
- Prior art keywords
- time
- value
- source excitation
- initial
- parameter
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
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 tremble time, which relate to the technical field of geophysical exploration, and the method comprises the following steps: acquiring a value range parameter, a value precision parameter and an initial tremor time value; determining a plurality of to-be-selected tremor time values according to the value range parameters and the value precision parameters; determining the selection probability of each vibration time value to be selected corresponding to the excitation moment of the target seismic source according to the value range parameter, the value precision parameter and the initial vibration time value; and generating aliasing acquisition dithering time corresponding to the target source excitation moment according to the selection probability. The method can improve the randomness of adjacent tremor time difference values 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 tremble time.
Background
The efficient aliasing acquisition technology greatly improves the earthquake acquisition efficiency, and is a technology widely adopted in the current high-precision high-density earthquake exploration. However, serious adjacent gun interference noise exists in the high-efficiency mixed data, and subsequent mixed data separation is needed. At present, the mainstream mixed data separation method is based on the characteristic that effective signals of mixed data in a certain domain are continuous, and mixed noise is randomly distributed. Therefore, the more the difference between the excitation moments of the current cannon and the adjacent cannons is close to uniform random distribution, the more data separation is facilitated. In order to make the excitation time of each gun produce irregular change, at present, a method of adding random trembling time at the excitation time of each gun is commonly adopted at home and abroad, but the method only aims at trembling time, although the independent trembling time is uniform and random in a designated range, the difference value is not uniform and random, and better data separation effect cannot be ensured.
Disclosure of Invention
The invention provides a method and a device for generating aliasing acquisition dithering time, which can realize uniform and random distribution of adjacent dithering time difference values and ensure a better data separation effect.
In a first aspect, an embodiment of the present invention provides a method for generating aliasing acquisition dithering time, where the method includes: acquiring a value range parameter, a value precision parameter and an initial tremor time value; determining a plurality of to-be-selected tremor time values according to the value range parameter and the value precision parameter; determining the selection probability of each vibration time value to be selected corresponding to the target source excitation moment according to the value range parameter, the value precision parameter and the initial vibration time value; generating aliasing acquisition trembling time corresponding to the target seismic source excitation moment according to the selection probability;
the determining the selection probability of each to-be-selected tremble time value corresponding to the target source excitation moment according to the value range parameter, the value precision parameter and the initial tremble time value comprises the following steps:
if the target source excitation time is a second source excitation time adjacent to the initial source excitation time, determining the selection probability of each to-be-selected tremble time value corresponding to the second source excitation time according to the value range parameter, the value precision parameter and the initial tremble time value; if the target source excitation time is the target source excitation time except the initial source excitation time and the second source excitation time, determining the selection probability of each to-be-selected vibration time value corresponding to the target source excitation time according to the value range parameter, the value precision parameter and the aliasing acquisition vibration time corresponding to the last target source excitation time.
In a second aspect, an embodiment of the present invention further provides an apparatus for generating an aliasing acquisition dithering time, where the apparatus includes: the acquisition module is used for acquiring the value range parameter, the value precision parameter and the initial tremor time value; the determining module is used for determining a plurality of to-be-selected tremor 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 tremble time value corresponding to the target source excitation moment according to the value range parameter, the value precision parameter and the initial tremble time value; the generation module is used for generating aliasing acquisition trembling time corresponding to the target seismic source excitation moment according to the selection probability;
wherein, probability module is specifically used for:
if the target source excitation time is a second source excitation time adjacent to the initial source excitation time, determining the selection probability of each to-be-selected tremble time value corresponding to the second source excitation time according to the value range parameter, the value precision parameter and the initial tremble time value; if the target source excitation time is the target source excitation time except the initial source excitation time and the second source excitation time, determining the selection probability of each to-be-selected vibration time value corresponding to the target source excitation time according to the value range parameter, the value precision parameter and the aliasing acquisition vibration time corresponding to the last target source excitation time.
In a third aspect, an embodiment of the present invention further provides a computer device, including a memory, a processor, and a computer program stored in the memory and capable of running on the processor, where the processor implements the method for generating aliasing acquisition dithering time when executing the computer program.
In a fourth aspect, an embodiment of the present invention further provides a computer readable storage medium storing a computer program for executing the above-described aliasing acquisition dither time generation method.
The embodiment of the invention has the following beneficial effects: the embodiment of the invention provides an aliasing acquisition trembling time generation scheme, which comprises the steps of firstly acquiring a value range parameter, a value precision parameter and an initial trembling time value, determining a plurality of to-be-selected trembling time values according to the value range parameter and the value precision parameter, then determining the selection probability of each to-be-selected trembling time value corresponding to a target source excitation moment according to the value range parameter, the value precision parameter and the initial trembling time value, and generating the aliasing acquisition trembling time corresponding to the target source excitation moment according to the selection probability. The embodiment of the invention can improve the randomness of the adjacent tremor time difference value corresponding to the target seismic source excitation moment based on probability weighting, and provides support for obtaining 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 above objects, features and advantages of the present invention more 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 that are needed in the description of the embodiments or the prior art will be briefly described, and it is obvious that the drawings in the description below are some embodiments of the present invention, and other drawings can be obtained according to the drawings without inventive effort for a person skilled in the art.
Fig. 1 is a flowchart of a method for generating aliasing acquisition dithering time according to an embodiment of the present invention;
FIG. 2 is a graph showing a distribution of jitter time difference values in a conventional random method according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of a combination of values of a dithering time range and values of adjacent dithering times according to an embodiment of the present invention;
fig. 4 is a schematic diagram of a rule (a set of linear equations) of values of adjacent tremor time difference values provided in an embodiment of the present invention;
FIG. 5 is a general flow chart of the probability weighted random method according to the embodiment of the present invention;
FIG. 6 is a graph of a square distribution of dithering time using a probability weighted random method in accordance with an embodiment of the present invention;
FIG. 7 is a graph of a square distribution of the difference in jitter time using a probability weighted random method according to an embodiment of the present invention;
FIG. 8 is a graph showing the effect of no jitter time provided by an embodiment of the present invention;
FIG. 9 is a graph showing the effect of a conventional random method according to an embodiment of the present invention;
FIG. 10 is a graph showing 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 dithering time generating device according to an embodiment of the present invention;
fig. 12 is a block diagram of another aliasing acquisition dithering 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
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, 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 embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
When the offshore seismic source ship works, although the ship speed is changed due to the influence of external environment, the external influence is relatively limited due to the huge mass of the hull of the seismic source ship. Thus, in general, the ship speed is still relatively constant, and the inter-gun time difference tends to be a certain value (e.g., the ship speed is stabilized at 2.5m/s, the inter-gun time difference is 12.5m, and the inter-gun time difference is basically distributed around 5 s). The rule of centralized distribution does not meet the requirement of the mixed data separation method on random distribution of time differences among cannons. Therefore, random dithering time is added to the firing time of each gun to cause irregular variation of the firing time.
The fundamental purpose of the dithering time is to excite the randomness of the moments. Let T be n D is the excitation time of the nth shot point n For the random trembling time of the nth shot, the formula of the inter-shot time difference is as follows:
ΔT＝(T n +d n )-(T n-1 +d n-1 )＝(T n -T n-1 )+(d n -d n-1 )
it can be seen that in the case of poor natural randomness, the difference between the original excitation timings varies little, so that the main factor affecting randomness of the final excitation timings is the difference between adjacent dither times (d n -d n-1 ). The vibration time random method commonly adopted at home and abroad at present only aims at the vibration time, although the independent vibration time is uniformly random within a specified range, the difference is not uniformly random, but isosceles triangle distribution gradually decreasing from the middle to two sides is presented, as shown in a vibration time difference distribution diagram of a conventional random method shown in fig. 2, the vibration time value range in the diagram is [ -500ms,500ms]The sampling interval is 2ms, the number of samples is 10 ten thousand, and the optimal data separation effect cannot be ensured by the existing method.
Based on the above, the method and the device for generating the aliasing acquisition dithering time provided by the embodiment of the invention can realize uniform and random distribution of adjacent dithering time difference values.
For the sake of understanding the present embodiment, first, a detailed description is given of a method for generating aliasing acquisition dithering time disclosed in the present embodiment.
The embodiment of the invention provides a method for generating aliasing acquisition dithering time, which is shown in a flow chart of the method for generating aliasing acquisition dithering time in fig. 1, and comprises the following steps:
step S102, acquiring a value range parameter, a value precision parameter and an initial dithering time value.
In the embodiment of the invention, the value range parameter is used for determining the duration range of the offset of the aliasing acquisition vibration time relative to the target seismic source excitation time. The value accuracy parameter is used to determine the duration interval between optional aliased acquisition dither times. The initial tremor time value may be pre-generated.
Step S104, a plurality of vibration time values to be selected are determined 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 to-be-selected tremor time values can be determined, and a plurality of to-be-selected tremor time values can be determined in the selectable range according to the value precision parameter. For example, the value range parameter includes plus or minus 500 milliseconds, the selectable range of the selected tremor time value is between minus 500 milliseconds and plus 500 milliseconds, and if the value accuracy parameter is 2 milliseconds, the selected tremor time value includes 500, 498, 496, … …, minus 498, minus 500.
And S106, determining the selection probability of each vibration time value to be selected corresponding to the target vibration source excitation moment 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 source excitation time may include a plurality of source excitation times, and for each source excitation time, a set of data is determined according to the value range parameter, the value precision parameter, and the initial tremor time value, where the set of data includes: each of the selected dithered time values and the corresponding selection probability of each of the selected dithered time values.
And S108, generating aliasing acquisition dithering time corresponding to the target source excitation moment according to the selection probability.
In the embodiment of the invention, the selection probability of the to-be-selected tremble time value is used as the selection weight, and for each target focus excitation moment, one value is selected from a plurality of to-be-selected tremble time values corresponding to the target focus excitation moment, so that the aliasing acquisition tremble time corresponding to the target focus excitation moment is obtained, and further, tremble time sequences corresponding to a plurality of target focuses are also obtained.
The embodiment of the invention provides an aliasing acquisition trembling time generation scheme, which comprises the steps of firstly acquiring a value range parameter, a value precision parameter and an initial trembling time value, determining a plurality of to-be-selected trembling time values according to the value range parameter and the value precision parameter, then determining the selection probability of each to-be-selected trembling time value corresponding to a target source excitation moment according to the value range parameter, the value precision parameter and the initial trembling time value, and generating the aliasing acquisition trembling time corresponding to the target source excitation moment according to the selection probability. The embodiment of the invention can improve the randomness of the adjacent tremor time difference value corresponding to the target seismic source excitation moment based on probability weighting, and provides support for obtaining better seismic data separation effect.
The target source excitation time includes an initial source excitation time, and the following steps may be further executed before acquiring the value range parameter, the value precision parameter and the initial dithering time value in order to improve the calculation efficiency:
randomly generating a tremble time value corresponding to the initial source excitation moment according to the value range parameter; and taking the dithering time value corresponding to the initial vibration source excitation moment as an initial dithering time value.
In the embodiment of the invention, when the tremor time sequence is acquired, the first tremor time value is acquired by adopting a conventional random method. The initial source activation time is the activation time corresponding to the source that was activated earliest among the plurality of target sources. And taking the dithering time value corresponding to the initial vibration source excitation moment as an initial dithering time value.
In order to improve randomness, the selection probability of each of the to-be-selected tremble time values corresponding to the target source excitation time is determined according to the value range parameter, the value precision parameter and the initial tremble time value, and the method can be implemented as follows.
If the target source excitation time is a second source excitation time adjacent to the initial source excitation time, determining the selection probability of each to-be-selected tremble time value corresponding to the second source excitation time according to the value range parameter, the value precision parameter and the initial tremble time value; if the target source excitation time is the target source excitation time except the initial source excitation time and the second source excitation time, determining the selection probability of each to-be-selected vibration time value corresponding to the target source excitation time according to the value range parameter, the value precision parameter and the aliasing acquisition vibration time corresponding to the last target source excitation time.
The probability of selection of each of the time values for the vibration to be selected corresponding to the target source activation time may be determined according to the following formula:
wherein p is d To select the probability, d is the value of the dithering time to be selected, d has S+1 possible values, S is the variable from 0 to S, d prev For the initial or last jitter time value, d step D for taking the value of the precision parameter max Is a value range parameter.
Referring to the general flow chart of fig. 5 when the probability weighted random method is adopted, in the embodiment of the invention, when the trembling time sequence is obtained, the first trembling time value is obtained by adopting a conventional random method, then all trembling times are obtained by adopting the probability weighted random method, and each time the last trembling time is used as a parameter to calculate the weight of each value of the current trembling time, and further the probability of each value is obtained. And finally, randomly sampling by a computer to obtain the current trembling time, and circularly reciprocating the current trembling time until the number of trembling time sequences meets the requirement.
The method may further perform the steps of: and adjusting the target source excitation time according to the aliasing acquisition dithering time so as to excite the target source according to the adjusted target source excitation time.
In the embodiment of the invention, the aliasing acquisition dithering time is a time offset, and the target vibration source excitation time is adjusted according to the aliasing acquisition dithering time, so that the adjusted target vibration source excitation time can be obtained. After obtaining the aliasing acquisition dithering time, for example, the aliasing acquisition dithering time corresponding to the target focus a, the target focus b and the target focus c is sequentially 5 seconds, 8 seconds and minus 2 seconds, the excitation time corresponding to the target focus a, the target focus b and the target focus c is sequentially 13 points, 14 points and 15 points, and the adjusted focus excitation time is sequentially 13 points, 5 seconds, 14 points, 8 seconds and 15 points with a difference of 2 seconds.
The invention can be realized as follows.
1) And (3) analysis of a value rule of the tremor time: and analyzing a dithering time value principle and a mathematical model.
2) And (3) analysis of a value rule of the vibration time difference value: and analyzing a vibration time difference value principle and a mathematical model.
3) Analysis of the probability distribution of the difference between the tremor time: and determining the probability distribution rule according to the mathematical model of the tremor time difference distribution.
4) Calculating a dithering time difference weight distribution: and adjusting weights according to the probability distribution rule of the tremor time difference values, so that the tremor time difference values realize equal probabilities in different values.
5) Calculating the weighted probability of the tremor time: and converting the weight value of the single sample point in the linear equation determined by the vibration time difference value into the integral weight and probability when the vibration time is different.
6) Cyclic acquisition of tremor time series: and circularly acquiring all the tremor time sequences according to a set flow.
7) And (3) analysis of a tremor time and tremor time difference result: the final tremor time series was analyzed and evaluated.
8) Analyzing results of time difference among cannons: and analyzing and evaluating the final inter-gun time difference distribution.
The implementation of the method is described below with reference to the drawings and formulas.
1) And (3) analysis of a value rule of the tremor time:
if the maximum absolute value of the tremor time is d max The value accuracy is d step 。
The value range is-d max To d max Is provided with
Taking the current dithering time as the vertical axis y and the next dithering time as the horizontal axis x, see the dithering time value range and the adjacent dithering time value combination schematic diagram shown in fig. 3, wherein the example is d max ＝500，d step =100. All points within a black square area represent all possible combinations of values for (x, y) that exhibit a uniform random distribution within this square area.
2) And (3) analysis of a value rule of the vibration time difference value:
let the difference between the current and last tremor times be Δd: Δd=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 delta d, and all grid points in a square area penetrated by the straight line are all possible (x, y) value combinations with a difference value equal to delta d.
The delta d is within the range of-2 d max To 2d max The total (2S+1) kinds of values are possible. Referring to the schematic diagram of the rule of values of adjacent dithering time difference values (set of linear equations) shown in fig. 4, as shown by oblique lines in the figure, each value corresponds to a unique linear equation, and the set of all equations is:
3) Analysis of the probability distribution of the difference between the tremor time:
from FIG. 4, the total number of all the dots within the black square area is (S+1) 2 The grid points are uniformly and randomly distributed.
The number of grid points traversed by the line y=x+Δd, i.e. all possible (x, y) value combinations, is equal to:the value probability is +.>Therefore:
when Δd=0, there are s+1 possibilities in total, the probability isThe highest;
when Δd=2d max Or Δd= -2d max When in use, the probability is that 1 kind of probability is onlyLowest.
4) Calculating a dithering time difference weight distribution:
because the grid points are uniformly and randomly distributed, the weight of the probability of each grid point value is equal in the default case, and the probability is set to be 1.
The number of grid points through which the straight line y=x+Δd passes varies with the Δd value, and to achieve equal probability of occurrence of all straight lines, the individual probability weights of the grid points through which they pass are set to beI.e. the opposite number of numbers, the weight of the straight line is the weight of a single grid point multiplied by the number of grid points, the weights are equal and are 1, and the equal probability of all straight lines can be realized.
5) Calculating the weighted probability of the tremor time:
when the last tremor time value x has been determined, the current tremor time value y has (s+1) possibilities, corresponding to (s+1) grid points on the two-dimensional coordinate axis. And combining the grid point single point weight formula calculated in the previous step.
The probability of value when y is any possible value is obtained as follows:
replacing x with the last jitter time d prev Y is replaced by any possible value d of the current dithering time, and the value probability of obtaining any d value is as follows:
wherein-d max ≤d≤d max D may have a common (S+1) number.
6) Cyclic acquisition of tremor time series:
when the trembling time sequence is obtained, the first trembling time value is still obtained by adopting a conventional random method, then all trembling times are obtained by adopting a probability weighted random method, the last trembling time is used as a parameter each time to calculate the weight of each value of the current trembling time, the probability of each value is further obtained, and finally the current trembling time is obtained by random sampling of a computer, and the steps are repeated until the number of the trembling time sequences meets the requirement. The flow is shown in fig. 5.
7) And (5) analyzing the tremor time and the difference value result of the tremor time.
The final tremor time distribution is regular with more sides, less middle, but the whole is symmetrical, and the average value approaches 0.
The distribution of the difference in the tremor time conforms to a uniform random characteristic, which is the final objective to be achieved by the present method. Fig. 6 and 7 are respectively a histogram of the jitter time and a histogram of the jitter time difference when a probability weighted random method is adopted (jitter time value range [ -500ms,500ms ] in the illustration, sampling interval 2ms, 10 ten thousand samples).
8) And analyzing the time difference results among the cannons.
Referring to the effect diagram without trembling time shown in fig. 8, the effect diagram of the conventional stochastic method shown in fig. 9 and the effect diagram of the probability weighted stochastic method shown in fig. 10, compared with the conventional stochastic method, the probability weighted stochastic method has a more obvious effect of improving the time difference among the cannons, and the time difference among the cannons becomes more uniform and random and meets the distribution requirement of the mixed mining seismic data separation on the time difference among the cannons.
The embodiment of the invention provides a method and a device for generating the aliasing acquisition tremble time, which can generate tremble time sequences required by the aliasing acquisition of a marine seismic exploration multi-source ship, wherein the tremble time sequences generated by the method can be distributed symmetrically and randomly, and adjacent difference values of the tremble time sequences generated by the method can be distributed symmetrically, randomly and uniformly.
The embodiment of the invention also provides an aliasing acquisition tremble time generation device, which is described in the following embodiment. Because the principle of the device for solving the problem is similar to that of the method for generating the aliasing acquisition dithering time, the implementation of the device can refer to the implementation of the method for generating the aliasing acquisition dithering time, and the repetition is not repeated. Referring to fig. 11, a block diagram of an aliasing acquisition dithering time generation apparatus is shown, the apparatus comprising:
an acquisition module 71, configured to acquire a value range parameter, a value precision parameter, and an initial jitter time value; a determining module 72, configured to determine a plurality of values of the vibration time to be selected according to the value range parameter and the value accuracy parameter; the probability module 73 is configured to determine a selection probability of each of the to-be-selected tremble time values corresponding to the target source excitation time according to the value range parameter, the value precision parameter and the initial tremble time value; the generating module 74 is configured to generate an aliasing acquisition dithering time corresponding to the target source excitation moment according to the selection probability.
In one embodiment, referring to another structural block diagram of an aliased acquisition vibroseis time generation apparatus shown in fig. 12, the target source activation time includes an initial source activation time, and the apparatus further includes a randomization module 75 for: randomly generating a tremble time value corresponding to the initial source excitation moment according to the value range parameter; and taking the dithering time value corresponding to the initial vibration source excitation moment as an initial dithering time value.
In one embodiment, the probability module is specifically configured to: if the target source excitation time is a second source excitation time adjacent to the initial source excitation time, determining the selection probability of each to-be-selected tremble time value corresponding to the second source excitation time according to the value range parameter, the value precision parameter and the initial tremble time value; if the target source excitation time is the target source excitation time except the initial source excitation time and the second source excitation time, determining the selection probability of each to-be-selected vibration time value corresponding to the target source excitation time according to the value range parameter, the value precision parameter and the aliasing acquisition vibration time corresponding to the last target source excitation time.
In one embodiment, the probability module is specifically configured to: the selection probability of each vibration time value to be selected corresponding to the target source excitation time is determined according to the following formula: wherein p is d To select the probability, d is the value of the dithering time to be selected, d has S+1 possible values, S is the variable from 0 to S, d prev For the initial or last jitter time value, d step D for taking the value of the precision parameter max Is a value range parameter.
In one embodiment, referring to another structural block diagram of the apparatus for generating aliased acquisition and dithering time shown in fig. 12, the apparatus further includes an adjustment module 76 for adjusting the target source activation time according to the aliased acquisition and dithering time so that the target source is activated according to the adjusted target source activation time.
The embodiment of the present invention further provides a computer device, referring to a schematic block diagram of a computer device structure shown in fig. 13, where the computer device includes a memory 81, a processor 82, and a computer program stored on the memory and capable of running on the processor, and when the processor executes the computer program, the processor implements the steps of any one of the aliasing acquisition and dithering time generation methods described above.
It will 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, which is not repeated herein.
The embodiment of the invention also provides a computer readable storage medium, and the computer readable storage medium stores a computer program for executing any one of the aliasing acquisition tremble time generation methods.
It will be appreciated by those skilled in the art that 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 flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations 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 examples are only specific embodiments of the present invention, and are not intended to limit the scope of the present invention, but it should be understood by those skilled in the art that the present invention is not limited thereto, and that the present invention is described in detail with reference to the foregoing examples: any person skilled in the art may modify or easily conceive of the technical solution described in the foregoing embodiments, or perform equivalent substitution of some of the technical features, while remaining within the technical scope of the present disclosure; such modifications, changes or substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention, and are intended to be included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (10)
1. A method of generating aliased acquisition vibroseis time, for use in a seismic exploration multi-source vessel, the method comprising:
acquiring a value range parameter, a value precision parameter and an initial tremor time value;
determining a plurality of to-be-selected tremor time values according to the value range parameter and the value precision parameter;
determining the selection probability of each vibration time value to be selected corresponding to the target source excitation moment according to the value range parameter, the value precision parameter and the initial vibration time value;
generating aliasing acquisition trembling time corresponding to the target seismic source excitation moment according to the selection probability;
the determining the selection probability of each to-be-selected tremble time value corresponding to the target source excitation moment according to the value range parameter, the value precision parameter and the initial tremble time value comprises the following steps:
if the target source excitation time is a second source excitation time adjacent to the initial source excitation time, determining the selection probability of each to-be-selected tremble time value corresponding to the second source excitation time according to the value range parameter, the value precision parameter and the initial tremble time value; if the target source excitation time is the target source excitation time except the initial source excitation time and the second source excitation time, determining the selection probability of each to-be-selected vibration time value corresponding to the target source excitation time according to the value range parameter, the value precision parameter and the aliasing acquisition vibration time corresponding to the last target source excitation time.
2. The method of claim 1, wherein the target source activation time comprises an initial source activation time;
before acquiring the value range parameter, the value precision parameter and the initial tremor time value, the method further comprises the following steps:
randomly generating a tremble time value corresponding to the initial source excitation moment according to the value range parameter;
and taking the tremor time value corresponding to the initial vibration source excitation moment as an initial tremor time value.
3. The method of claim 1, comprising determining a selection probability for each of the selected time-to-dither values corresponding to a target source activation time according to the following equation:
wherein p is d To select the probability, d is the value of the dithering time to be selected, d has S+1 possible values, S is the variable from 0 to S, d prev For the initial or last jitter time value, d step D for taking the value of the precision parameter max Is a value range parameter.
4. A method according to any one of claims 1-3, further comprising:
and adjusting the target focus excitation time according to the aliasing acquisition dithering time so as to excite the target focus according to the adjusted target focus excitation time.
5. An aliased acquisition vibroseis time generation device, characterized in that it is applied to a seismic exploration multi-source ship, comprising:
the acquisition module is used for acquiring the value range parameter, the value precision parameter and the initial tremor time value;
the determining module is used for determining a plurality of to-be-selected tremor 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 tremble time value corresponding to the target source excitation moment according to the value range parameter, the value precision parameter and the initial tremble time value;
the generation module is used for generating aliasing acquisition trembling time corresponding to the target seismic source excitation moment according to the selection probability;
wherein, probability module is specifically used for:
if the target source excitation time is a second source excitation time adjacent to the initial source excitation time, determining the selection probability of each to-be-selected tremble time value corresponding to the second source excitation time according to the value range parameter, the value precision parameter and the initial tremble time value; if the target source excitation time is the target source excitation time except the initial source excitation time and the second source excitation time, determining the selection probability of each to-be-selected vibration time value corresponding to the target source excitation time according to the value range parameter, the value precision parameter and the aliasing acquisition vibration time corresponding to the last target source excitation time.
6. The apparatus of claim 5, wherein the target source activation time comprises an initial source activation time; the apparatus further comprises a random module for:
randomly generating a tremble time value corresponding to the initial source excitation moment according to the value range parameter;
and taking the tremor time value corresponding to the initial vibration source excitation moment as an initial tremor time value.
7. The apparatus of claim 5, wherein the probability module is specifically configured to: determining the selection probability of each to-be-selected tremor time value corresponding to the target source excitation time according to the following formula:
wherein p is d To select the probability, d is the value of the dithering time to be selected, d has S+1 possible values, S is the variable from 0 to S, d prev For the initial or last jitter time value, d step D for taking the value of the precision parameter max Is a value range parameter.
8. The apparatus of any of claims 5-7, further comprising an adjustment module to adjust the target source activation time according to the aliased acquisition dither time to cause the target source to be activated according to the adjusted target source activation time.
9. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the processor implements the method of any of claims 1 to 4 when executing the computer program.
10. 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 4.
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 CN112051613A (en) | 2020-12-08 |
| CN112051613B true 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 (5)
| 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 |
| 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 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012097122A1 (en) * | 2011-01-12 | 2012-07-19 | Bp Corporation North America Inc. | Shot scheduling limits for seismic acquisition with simultaneous source shooting |
| US10871584B2 (en) * | 2013-06-17 | 2020-12-22 | Westerngeco L.L.C. | Seismic data processing |
| US10670753B2 (en) * | 2014-03-03 | 2020-06-02 | Saudi Arabian Oil Company | History matching of time-lapse crosswell data using ensemble kalman filtering |
| US11009618B2 (en) * | 2018-09-04 | 2021-05-18 | 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 (5)
| 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 |
| 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 |
Non-Patent Citations (3)
| Title |
|---|
| 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;Fátima Gouveia ，Rui Carrilho Gomes ， Isabel Lopes;Journal of Applied Geophysics;全文 * |
| 强环噪地区可控震源高效混叠采集应用效果;韩志雄;倪宇东;余建鹏;马涛;卢秀丽;张华;;地球物理学进展(第01期);全文 * |
| 马竹，侯昆鹏，柳兴刚等.多震源船混叠釆集导航模拟方法及应用.物探装备.2020,正文第180-183页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112051613A (en) | 2020-12-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110031896B (en) | Seismic random inversion method and device based on multi-point geostatistics prior information | |
| CN104021431B (en) | Robust Optimal methods based on average gradient value and improvement multi-objective particle swarm optimization | |
| AU2002302834A1 (en) | A method of processing data | |
| WO2002103397A1 (en) | A method of processing data | |
| CN106842306A (en) | The staggered-mesh finite difference analogy method and device of a kind of global optimization | |
| CA2969101A1 (en) | Seismic elastic wave simulation for tilted transversely isotropic media using adaptive lebedev staggered grid | |
| Valentini et al. | Relaxing segmentation on the Wasatch fault zone: Impact on seismic hazard | |
| CN112051613B (en) | Method and device for generating aliasing acquisition dithering time | |
| CN114578417A (en) | Pulse type earthquake motion fitting method based on wavelet adjustment | |
| Yu et al. | Lppn: A lightweight network for fast phase picking | |
| US20150134308A1 (en) | Method and device for acquiring optimization coefficient, and related method and device for simulating wave field | |
| CN108680968B (en) | Evaluation method and device for seismic exploration data acquisition observation system in complex structural area | |
| CN108957545B (en) | Air-gun array wavelet directionality the Method of Deconvolution and system | |
| Zhao et al. | Study on air gun wavelet characteristics under different source and towed-streamer sinking depths in the Bohai sea | |
| CN111435174B (en) | Method and device for compensating amplitude of seismic data in strong reflection area | |
| CN109633781B (en) | Geological property acquisition method and device, electronic equipment and storage medium | |
| CN109581484B (en) | Method and device for analyzing total illumination index of seismic longitudinal wave component multiples | |
| CN104217099B (en) | A method of for generating and handling the scanning signal of controlled source | |
| CN110376642B (en) | Three-dimensional seismic velocity inversion method based on conical surface waves | |
| CN113411140B (en) | Channel self-adaptive channel characteristic learning method, simulator thereof and satellite navigation system | |
| Haavik et al. | Estimation of source signatures from air guns fired at various depths: A field test of the source scaling law | |
| Richards | Spatial correlation test for chaotic dynamics in political science | |
| CN109001804A (en) | A kind of method, apparatus and system determining effective stress based on 3D seismic data | |
| CN113970789B (en) | Full waveform inversion method and device, storage medium and electronic equipment | |
| CN112539829B (en) | Arrangement method and system of self-adaptive blade tip timing sensors |
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 |