[go: up one dir, main page]

CN107065007B - A kind of seismic data amplitude adjustment method and device - Google Patents

A kind of seismic data amplitude adjustment method and device Download PDF

Info

Publication number
CN107065007B
CN107065007B CN201710372329.4A CN201710372329A CN107065007B CN 107065007 B CN107065007 B CN 107065007B CN 201710372329 A CN201710372329 A CN 201710372329A CN 107065007 B CN107065007 B CN 107065007B
Authority
CN
China
Prior art keywords
measured
gather
seismic
signal
adjustment factor
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
Application number
CN201710372329.4A
Other languages
Chinese (zh)
Other versions
CN107065007A (en
Inventor
杜磊
郭平
高树生
柳世光
冯海涛
孙晶波
卢明德
张高
蓝阔
董兵波
桑雨
王国雷
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201710372329.4A priority Critical patent/CN107065007B/en
Publication of CN107065007A publication Critical patent/CN107065007A/en
Application granted granted Critical
Publication of CN107065007B publication Critical patent/CN107065007B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • 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. for interpretation or for event detection
    • G01V1/34—Displaying seismic recordings or visualisation of seismic data or attributes
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/50—Corrections or adjustments related to wave propagation
    • G01V2210/51—Migration
    • G01V2210/512—Pre-stack

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention provides a seismic data amplitude adjusting method and a seismic data amplitude adjusting device, wherein the method comprises the following steps: acquiring the signal-to-noise ratio of a gather to be detected; performing dynamic correction stretching excision processing on the gather to be detected, and determining excision time points of each seismic channel in the gather to be detected; determining a plurality of time periods according to the cutting time points of all seismic channels in the gather of the to-be-detected channels; counting the number of seismic channels of each time period of the removed gather to be detected falling into a plurality of time periods; calculating to obtain an energy adjustment factor of the gather to be measured in each time period according to the signal-to-noise ratio and the number of seismic channels in each time period; and adjusting the amplitude of the gather to be detected based on the energy adjustment factor to obtain an adjusted seismic signal. In the embodiment of the invention, the problem of uneven signal energy distribution caused by the change of the offset distance of the gather to be detected is considered, and exploration signals which are more accordant with geological rules can be obtained in the splicing area and the edge area of the connected block.

Description

Seismic data amplitude adjusting method and device
Technical Field
The invention relates to the technical field of geological exploration, in particular to a seismic data amplitude adjusting method and device.
Background
At present, because construction factors, excitation and receiving conditions and signal acquisition time are different among the blocks, the difference of parameters acquired among the blocks is large, for example: covering times, single shot energy intensity and other parameters. When the seismic data of each block are used for continuous imaging processing, in the splicing area of the block, because the seismic data are obtained by stacking two or more blocks of seismic data, the coverage times of the seismic data at the stacking position are obviously higher than those of other places, and if effective amplitude adjustment cannot be carried out, an arc strike phenomenon is inevitably caused; in the edge area of the block, due to acquisition reasons, the coverage times of the seismic gathers are obviously lower than those of a full coverage area, and strong energy in other areas can strike arcs to the edge of the block.
In general, when the slice block is imaged, amplitude adjustment may be performed on the slice seismic data by using methods such as geometric diffusion compensation and earth surface consistency amplitude compensation. However, when the method is adopted, the amplitude energy difference of the original seismic information in each block and among the blocks is too large, and the far and near offset energy is weak, so the adjusted seismic data has poor effect.
In view of the above problems, no effective solution has been proposed.
Disclosure of Invention
The invention provides a seismic data amplitude adjusting method and device, aiming at achieving the purpose of consistent energy among blocks under the condition of ensuring that seismic data information is not lost when a connected block is subjected to imaging processing.
The embodiment of the invention provides a seismic data amplitude adjusting method, which comprises the following steps: acquiring the signal-to-noise ratio of a gather to be detected; performing dynamic correction stretching excision processing on the gather to be detected, and determining excision time points of all seismic channels in the gather to be detected; determining a plurality of time periods according to the cutting time points of the seismic channels in the channel set to be detected; counting the number of seismic channels of the removed to-be-detected gather falling into each time period in the multiple time periods; calculating to obtain energy adjustment factors of the to-be-detected gather in each time period according to the signal-to-noise ratio and the number of seismic channels in each time period; and adjusting the amplitude of the gather to be detected based on the energy adjustment factor to obtain an adjusted seismic signal.
In one embodiment, the energy adjustment factor of the gather to be measured in each time period may be calculated according to the following formula based on the signal-to-noise ratio and the number of seismic traces in each time period:
cof(t)＝P*cov(t)-1/K+(1-P)*snr-1/L
the energy adjustment factor of the gather to be measured in the time period t is represented by cof (t), the seismic trace number of the gather to be measured in the time period t is represented by cov (t), the snr represents the signal-to-noise ratio, P is more than or equal to 0 and less than or equal to 1, K is more than 1, and L is more than 0.
In one embodiment, before acquiring the signal-to-noise ratio of the gather to be measured, the method may further include: acquiring seismic data of a work area to be detected; preprocessing the seismic data to obtain processed seismic data, wherein the preprocessing may include, but is not limited to, at least one of: pre-stack denoising, ground surface consistency amplitude compensation, geometric diffusion absorption compensation and channel sorting; and determining a plurality of to-be-detected gathers according to the processed seismic data.
In one embodiment, after calculating the energy adjustment factor of the gather to be measured in each time period according to the signal-to-noise ratio and the number of seismic traces in each time period, the method may further include: and smoothing the energy adjustment factor of the gather to be detected to obtain the energy adjustment factor after smoothing of the gather to be detected.
In one embodiment, the amplitude of the gather to be measured may be adjusted based on the energy adjustment factor according to the following formula to obtain an adjusted seismic signal:
s'(t)＝s(t)·pcof(t)
wherein s' (t) represents the adjusted seismic signal in the time period t, s (t) represents the seismic signal in the time period t before the adjustment of the to-be-measured gather, and pcof (t) represents the energy adjustment factor of the to-be-measured gather after the smoothing processing in the time period t.
In an embodiment, smoothing the energy adjustment factor of the gather to be measured to obtain the energy adjustment factor after smoothing the gather to be measured may include: sequentially calibrating the positions of all to-be-measured gathers in the to-be-measured work area according to a preset sequence to obtain all to-be-measured gathers after position calibration; acquiring energy adjustment factors after smoothing treatment of each to-be-measured gather according to the following modes: taking the position of the current gather to be measured as a center and a preset smooth radius as a radius to obtain a smooth area; and calculating to obtain the energy adjustment factor after the smoothing of the current to-be-measured gather according to the energy adjustment factor of each to-be-measured gather in the smoothing area.
In an embodiment, the energy adjustment factor after the smoothing of the current gather to be measured may be calculated according to the following formula:
wherein whenWhen wij＝0；
When in useWhen the temperature of the water is higher than the set temperature,
therein, pcofklEnergy adjustment factor after smoothing of the current gather to be measured with the representation position (k, l), cofijEnergy adjustment factor, w, before smoothing of the gather to be measured at position (i, j)ijAnd (3) representing the smoothing coefficient of the gather to be measured with the position (i, j), wherein R represents the smoothing radius, R is more than or equal to 1, i is k-R, k-R +1, …, k + R, j is l-R, l-R +1, … and l + R.
In one embodiment, the smoothing radius may be determined according to the number of to-be-measured gathers in the to-be-measured work area.
In one embodiment, the snr of the gather under test can be obtained according to the following formula:
wherein,
wherein snr represents the treatSignal-to-noise ratio, R (x), of the gather of tracesc,xd) And the cross-correlation value of the c-th seismic trace and the d-th seismic trace in the trace set to be detected is represented, wherein c is 1,2 …, N, d is 1,2 …, N and N represent the total number of sample points in the analysis time window of the c-th seismic trace.
In one embodiment, the cross-correlation value of the c-th seismic trace and the d-th seismic trace in the trace set to be tested can be determined as follows: determining the cross-correlation value of the c-th seismic channel and the d-th seismic channel in the to-be-detected channel set corresponding to each sampling point serial number in the seismic channel analysis time window according to the following formula, and obtaining a plurality of cross-correlation results corresponding to each sampling point serial number:
wherein t is 0,1, …, L represents the total number of samples in the c-th seismic trace analysis time window, τ represents the number of samples in the seismic trace analysis time window, τ is 0,1, …, L;
and selecting the maximum value in the multiple cross-correlation results as the cross-correlation value of the c-th seismic trace and the d-th seismic trace.
The embodiment of the invention also provides a seismic data amplitude adjusting device, which comprises: the signal-to-noise ratio acquisition module is used for acquiring the signal-to-noise ratio of the gather to be detected; the trace gather cutting module is used for performing dynamic correction stretching cutting processing on the trace gather to be detected and determining cutting time points of all seismic traces in the trace gather to be detected; the time period determination module is used for determining a plurality of time periods according to the cutting time points of all seismic channels in the channel set to be detected; the seismic channel number counting module is used for counting the number of seismic channels of each time period of the cut to-be-detected gather falling into the plurality of time periods; the adjusting factor calculating module is used for calculating and obtaining an energy adjusting factor of the gather to be measured in each time period according to the signal-to-noise ratio and the number of seismic channels in each time period; and the amplitude adjusting module is used for adjusting the amplitude of the to-be-detected gather based on the energy adjusting factor to obtain an adjusted seismic signal.
In the embodiment of the invention, because the dynamic correction stretching and cutting processing can effectively cut off the noise signals of the to-be-detected gather which lose the fluctuation characteristics due to the change of the offset distance, the dynamic correction stretching and cutting processing can be firstly carried out on the to-be-detected gather, the cutting time point of each seismic channel in the to-be-detected gather is determined, and the problem of the uneven signal energy distribution caused by the change of the offset distance of the to-be-detected gather is considered by counting the number of seismic channels of each time period in a plurality of time periods determined by the cutting time point of the cut-off to-be-detected gather. Furthermore, the amplitude adjustment is carried out by utilizing the energy adjustment factors of the gathers to be detected determined by the seismic channel numbers and the signal-to-noise ratios, so that exploration signals which are more in accordance with geological rules can be obtained in splicing areas and marginal areas of the connected blocks, and when the exploration signals are used for imaging the connected blocks, geological imaging results with higher quality can be obtained, so that the energy among the blocks is consistent under the condition that seismic data information is not lost.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, it is obvious that the drawings in the following description are only some embodiments described in the present application, and for those skilled in the art, other drawings can be obtained according to the drawings without any creative effort.
FIG. 1 is a flow chart of a seismic data amplitude adjustment method provided by the present application;
FIG. 2 is a schematic diagram of migration imaging of seismic signals after amplitude adjustment of a region based on a conventional amplitude adjustment method;
FIG. 3 is a schematic diagram of migration imaging of seismic signals after amplitude adjustment of a region based on the amplitude adjustment method provided herein;
fig. 4 is a block diagram of a seismic data amplitude adjustment apparatus according to the present invention.
Detailed Description
In order to make those skilled in the art better understand the technical solutions in the present application, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only a part of the embodiments of the present application, and not all of the embodiments. 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 application.
Considering the problem of uneven signal energy distribution caused by different offset distance changes in each seismic gather in the prior art, the dynamic correction stretching excision process can excise noise signals generated by offset distance changes in the to-be-measured gather, the inventor proposes a mode of combining the dynamic correction stretching excision process and multiple times of seismic trace counting, firstly counting the seismic trace number of each time period of the excised to-be-measured gather in a plurality of time periods determined by excision time points, then determining an energy adjustment factor of the to-be-measured gather by using each seismic trace number, and finally carrying out amplitude adjustment on the to-be-measured gather according to the energy adjustment factor. Based on this, a method for adjusting amplitude of seismic data is proposed, as shown in fig. 1, which may include the following steps:
s101: and acquiring the signal-to-noise ratio of the gather to be detected.
The gather to be measured can be a common central point gather, a common shot point gather, a common receiving point gather, a common reflection point gather, a common depth point gather and the like.
In an embodiment of the present application, before acquiring the signal-to-noise ratio of the gather to be measured, the method may further include: acquiring seismic data of a work area to be detected; preprocessing the seismic data to obtain processed seismic data, wherein the preprocessing may include, but is not limited to, at least one of: pre-stack denoising, ground surface consistency amplitude compensation, geometric diffusion absorption compensation and channel sorting; and determining a plurality of to-be-detected gather sets according to the processed seismic data.
After the preprocessing mode is adopted to carry out denoising processing on the to-be-detected gather, common noise signals in the to-be-detected gather can be removed, and therefore the to-be-detected gather which better accords with geological characteristics can be obtained.
In an embodiment of the present application, the signal-to-noise ratio of the gather to be measured may be obtained according to the following formula:
wherein,
wherein snr represents the signal-to-noise ratio of the gather to be measured, R (x)c,xd) And (3) representing the cross-correlation value of the c-th seismic trace and the d-th seismic trace in the trace set to be detected, wherein c is 1,2 …, N, d is 1,2 …, N and N represents the total number of sample points in the analysis time window of the c-th seismic trace.
In the process of determining the signal-to-noise ratio of the gather to be detected, the cross-correlation value of the c-th seismic channel and the d-th seismic channel in the gather to be detected can be determined in the following way:
s1-1: determining the cross-correlation value of the c-th seismic channel and the d-th seismic channel in the to-be-detected channel set corresponding to each sampling point serial number in the seismic channel analysis time window according to the following formula, and obtaining a plurality of cross-correlation results corresponding to each sampling point serial number:
wherein t is 0,1, …, L represents the total number of samples in the c-th seismic trace analysis time window, τ represents the number of samples in the seismic trace analysis time window, τ is 0,1, …, L;
s1-2: and selecting the maximum value in the multiple cross-correlation results as the cross-correlation value of the c-th seismic trace and the d-th seismic trace.
And substituting the obtained cross correlation values into a signal-to-noise ratio calculation formula, thereby calculating the signal-to-noise ratio of the gather to be measured.
S102: and performing dynamic correction stretching excision processing on the gather to be detected, and determining excision time points of all seismic channels in the gather to be detected.
S103: and determining a plurality of time periods according to the cutting time points of the seismic channels in the channel set to be detected.
S104: and counting the number of seismic channels of each time period in the removed gather to be detected in the plurality of time periods.
The seismic data dynamic correction relates to the precision and application effect of the seismic data processing result. In the dynamic correction process, the signal time can be elongated, the frequency can move to the low frequency, and the waveform variation phenomenon can be generated, when the stretching distortion is serious, the superposition effect can be damaged, the seismic signal can lose the wave dynamic characteristics, the resolution ratio is reduced, and the precision and the geological effect of seismic exploration can be influenced. Therefore, the waveform after the dynamic correction stretching can be subjected to the cutting processing, and the purpose of eliminating the above influence can be achieved.
The method is applied to the application, namely, stretching and cutting processing can be performed through dynamic correction on the to-be-detected gather, cutting time points of all seismic channels in the to-be-detected gather are determined, a plurality of time periods are determined according to the cutting time points of all seismic channels in the to-be-detected gather, and the to-be-detected gather after cutting is counted and falls into the seismic channel number of all the time periods in the plurality of time periods.
In the practical application process, the cutting time of each seismic channel can be obtained by reading the information in the head of each seismic channel, and then the statistics of the number of seismic channels is carried out corresponding to the channel set to be detected. By counting the number of seismic channels subjected to cutting processing in each time period determined by the cutting time point, the influence of offset change on the to-be-detected gather in the to-be-detected gather is considered, and the imaging processing effect of the connected block is improved.
Further, in the present application, the number of seismic channels included in the gather to be tested is the number of coverage times of the gather to be tested. Therefore, the above method for counting the number of seismic channels is to count the number of coverage times of the gather to be measured. Thus, for a single gather under test, the number of coverages counted in the manner described above varies over a period of time determined by the ablation time. For example: a certain gather to be tested is provided with 10 seismic channels, and after dynamic correction stretching excision processing is carried out, 10 seismic channels are not excised within 0-5 s; in 5s-7s, 8 seismic traces are not cut; 5 seismic traces were not cut at 7s-8 s; 3 seismic traces were not cut at 8s-9 s; in 9s-10s, 1 trace was not cut. For the gather to be measured in this example, the number of coverage is 10 in 0-5 s; in 5s-7s, the covering times are 8; in 7s-8s, the covering times are 5; in 8s-9s, the covering times are 3; the number of coverage was 1 from 9s to 10 s.
S105: and calculating the energy adjustment factor of the to-be-measured gather in each time period according to the signal-to-noise ratio and the number of seismic traces in each time period.
After the signal-to-noise ratio of the to-be-measured gather and the number of seismic channels that change with time periods are obtained by distribution by using the methods described in S101 to S105, an energy adjustment factor of the to-be-measured gather in each time period can be calculated according to the signal-to-noise ratio and the number of seismic channels in each time period by using the following formula:
cof(t)＝P*cov(t)-1/K+(1-P)*snr-1/L
the energy adjustment factor of the to-be-detected gather in the time period t is represented by cof (t), the seismic trace number of the to-be-detected gather in the time period t is represented by cov (t), the snr represents the signal to noise ratio, P is more than or equal to 0 and less than or equal to 1, K is more than 1, and L is more than 0.
For a gather to be measured, the number of coverages varies over a period of time, but the signal-to-noise ratio is constant.
In an embodiment of the application, after obtaining the energy adjustment factor, smoothing may be performed on the energy adjustment factor of the gather to be measured, so as to obtain the energy adjustment factor after smoothing of the gather to be measured. Before the smoothing processing, each gather to be measured in the work area to be measured is calculated to obtain energy adjustment factors corresponding to each gather to be measured. When the smoothing processing is performed, the smoothing operation is performed on the energy adjustment factors corresponding to the respective gathers to be measured. Specifically, the method can comprise the following steps:
s5-1: sequentially calibrating the position of each to-be-measured gather in the to-be-measured work area according to a preset sequence to obtain each to-be-measured gather after position calibration;
s5-2: acquiring energy adjustment factors after smoothing treatment of each to-be-measured gather according to the following modes: taking the position of the current gather to be measured as a center and a preset smooth radius as a radius to obtain a smooth area; and calculating to obtain the energy adjustment factor after the smoothing of the current to-be-measured gather according to the energy adjustment factor of each to-be-measured gather in the smoothing area.
That is, the energy adjustment factor after the smoothing of the current gather to be measured can be calculated according to the following formula:
wherein whenWhen wij＝0；
When in useWhen the temperature of the water is higher than the set temperature,
therein, pcofklEnergy adjustment factor after smoothing of the current gather to be measured with the position (k, l), cofijEnergy adjustment factor, w, before smoothing of the gather to be measured at position (i, j)ijAnd (3) representing the smoothing coefficient of the gather to be measured with the position (i, j), wherein R represents the smoothing radius, i is k-R, k-R +1, …, k + R, j is l-R, l-R +1, … and l + R.
The smooth radius can be determined according to the number of the gathers to be measured in the work area to be measured. That is, when there are 10 gathers to be measured in the work area to be measured, the smoothing radius may be 10.
S106: and adjusting the amplitude of the gather to be detected based on the energy adjustment factor to obtain an adjusted seismic signal.
In an embodiment of the present application, after obtaining the energy adjustment factor or the smoothed energy adjustment factor, amplitude adjustment may be performed on the gather to be measured based on the energy adjustment factor or the smoothed energy adjustment factor according to the following formula, so as to obtain an adjusted seismic signal:
s'(t)＝s(t)·pcof(t)
wherein s' (t) represents the adjusted seismic signal in the time period t, s (t) represents the seismic signal in the time period t before the adjustment of the to-be-measured gather, and pcof (t) represents the energy adjustment factor when the smoothing processing is performed on the to-be-measured gather in the time period t.
By adjusting the amplitude of the seismic data, the total energy of the surface element gather is approximately kept at the same level, and a prestack gather with uniform coverage times is obtained.
As shown in fig. 2, a schematic diagram of migration imaging of seismic signals after amplitude adjustment is performed on a certain area based on a conventional amplitude adjustment method, and as shown in fig. 3, a schematic diagram of migration imaging of seismic signals after amplitude adjustment according to an embodiment of the present invention is shown, it can be found by comparing fig. 2 and fig. 3: the 6 sections are clearly visible at the 6 white arrows in fig. 3, i.e. fig. 3 may better show the details of the formation. Therefore, the amplitude adjustment method provided by the embodiment can effectively solve the problems of pre-stack offset arc drawing and boundary effect, effectively improves the signal-to-noise ratio of the offset section, enables seismic data interpreters to more accurately judge and explain the underground geological structure condition, and improves the precision of searching favorable oil and gas trapping and determining well position.
Based on the same inventive concept, the embodiment of the present invention further provides an amplitude adjustment apparatus for seismic data, as described in the following embodiments. Because the principle of solving the problem of the seismic data amplitude adjusting device is similar to that of the seismic data amplitude adjusting method, the implementation of the seismic data amplitude adjusting device can refer to the implementation of the seismic data amplitude adjusting method, and repeated details are not repeated. As used hereinafter, the term "unit" or "module" may be a combination of software and/or hardware that implements a predetermined function. Although the means described in the embodiments below are preferably implemented in software, an implementation in hardware, or a combination of software and hardware is also possible and contemplated. Fig. 4 is a block diagram of a seismic data amplitude adjustment apparatus according to an embodiment of the present invention, and as shown in fig. 4, the apparatus may include: the system comprises a signal-to-noise ratio acquisition module 401, a gather cutting module 402, a time period determination module 403, a seismic channel number statistics module 404, an adjustment factor calculation module 405, and an amplitude adjustment module 406, which are described below.
The signal-to-noise ratio acquisition module 401 may be configured to acquire a signal-to-noise ratio of a gather to be detected;
a gather excision module 402, configured to perform dynamic correction stretching excision processing on the gather to be detected, and determine excision time points of each seismic channel in the gather to be detected;
the time period determining module 403 may be configured to determine a plurality of time periods according to the cutting time point of each seismic channel in the to-be-detected channel set;
a seismic channel number counting module 404, configured to count the number of seismic channels of the removed to-be-detected gather falling into each time period in the multiple time periods;
an adjustment factor calculation module 405, configured to calculate an energy adjustment factor of the to-be-measured gather in each time period according to the signal-to-noise ratio and the number of seismic channels in each time period;
the amplitude adjustment module 406 may be configured to perform amplitude adjustment on the gather to be detected based on the energy adjustment factor, so as to obtain an adjusted seismic signal.
In an embodiment of the application, the adjustment factor calculation module may be specifically configured to calculate an energy adjustment factor of the to-be-measured gather in each time period according to the signal-to-noise ratio and the number of seismic traces in each time period according to the following formula:
cof(t)＝P*cov(t)-1/K+(1-P)*snr-1/L
the energy adjustment factor of the gather to be measured in the time period t is represented by cof (t), the seismic trace number of the gather to be measured in the time period t is represented by cov (t), the snr represents the signal-to-noise ratio, P is more than or equal to 0 and less than or equal to 1, K is more than 1, and L is more than 0.
In an embodiment of the present application, the snr obtaining module may include: the seismic data acquisition unit can be used for acquiring seismic data of a work area to be detected before acquiring the signal-to-noise ratio of the gather to be detected; a preprocessing unit, which may be configured to preprocess the seismic data to obtain processed seismic data, wherein the preprocessing may include, but is not limited to, at least one of the following: pre-stack denoising, ground surface consistency amplitude compensation, geometric diffusion absorption compensation and channel sorting; and the gather determining unit can be used for determining a plurality of to-be-detected gathers according to the processed seismic data.
In an embodiment of the application, the adjustment factor calculation module calculates the energy adjustment factor of the to-be-measured gather in each time period according to the signal-to-noise ratio and the number of seismic channels in each time period, and then can be used for smoothing the energy adjustment factor of the to-be-measured gather to obtain the energy adjustment factor after smoothing the to-be-measured gather.
In an embodiment of the application, the amplitude adjustment module may be specifically configured to perform amplitude adjustment on the gather to be measured based on the energy adjustment factor according to the following formula to obtain an adjusted seismic signal:
s'(t)＝s(t)·pcof(t)
wherein s' (t) represents the adjusted seismic signal in the time period t, s (t) represents the seismic signal in the time period t before the adjustment of the to-be-measured gather, and pcof (t) represents the energy adjustment factor of the to-be-measured gather after the smoothing processing in the time period t.
In one embodiment of the present application, the adjustment factor calculation module may include: the position calibration unit can be used for sequentially calibrating the position of each to-be-tested gather in the to-be-tested work area according to a preset sequence to obtain each to-be-tested gather after position calibration is carried out; the factor calculating unit may be configured to obtain the energy adjustment factor after smoothing processing of each to-be-measured gather according to the following manner: taking the position of the current gather to be measured as a center and a preset smooth radius as a radius to obtain a smooth area; and calculating to obtain the energy adjustment factor after the smoothing of the current to-be-measured gather according to the energy adjustment factor of each to-be-measured gather in the smoothing area.
In an embodiment of the present application, the factor calculating unit may be specifically configured to calculate an energy adjustment factor after the smoothing processing of the current gather to be measured according to the following formula:
wherein whenWhen wij＝0；
When in useWhen the temperature of the water is higher than the set temperature,
therein, pcofklEnergy adjustment factor after smoothing of the current gather to be measured with the representation position (k, l), cofijEnergy adjustment factor, w, before smoothing of the gather to be measured at position (i, j)ijAnd (3) representing a smoothing coefficient of the gather to be measured with the position (i, j), wherein R represents the smoothing radius, i is k-R, k-R +1, …, k + R, j is l-R, l-R +1, … and l + R.
In an embodiment of the present application, the smoothing radius may be determined according to the number of to-be-measured gathers in the to-be-measured work area.
In an embodiment of the present application, the snr obtaining module may be specifically configured to obtain an snr of a gather to be measured according to the following formula:
wherein,
wherein snr represents the signal-to-noise ratio of the gather to be measured, R (x)c,xd) Representing the c-th seismic channel and the second seismic channel in the channel set to be detectedAnd d, c is 1,2 …, N, d is 1,2 …, N and N represents the total number of sample points in the analysis time window of the c-th seismic trace.
In an embodiment of the application, the snr obtaining module may be specifically configured to determine a cross-correlation value between a c-th seismic trace and a d-th seismic trace in the trace set to be detected according to the following manner: determining the cross-correlation value of the c-th seismic channel and the d-th seismic channel in the to-be-detected channel set corresponding to each sampling point serial number in the seismic channel analysis time window according to the following formula, and obtaining a plurality of cross-correlation results corresponding to each sampling point serial number:
wherein t is 0,1, …, L represents the total number of samples in the c-th seismic trace analysis time window, τ represents the number of samples in the seismic trace analysis time window, τ is 0,1, …, L;
and selecting the maximum value in the multiple cross-correlation results as the cross-correlation value of the c-th seismic trace and the d-th seismic trace.
By using the implementation mode of the seismic data amplitude adjustment device provided by each embodiment, the seismic data amplitude adjustment method can be automatically implemented, amplitude adjustment can be performed on the to-be-detected gather, specific participation of implementation personnel can be omitted, the amplitude adjustment result can be directly output, the operation is simple and rapid, and the user experience is effectively improved.
In the seismic data amplitude adjusting device, the energy adjustment factor calculation, the smoothing processing, and the signal-to-noise ratio acquisition may be extended as described in the foregoing method.
From the above description, it can be seen that the embodiments of the present invention achieve the following technical effects: the dynamic correction stretching and cutting processing can effectively cut off noise signals of the to-be-detected gather which lose fluctuation characteristics due to offset distance change, so that the dynamic correction stretching and cutting processing can be firstly carried out on the to-be-detected gather, the cutting time point of each seismic channel in the to-be-detected gather is determined, the number of seismic channels of each time period in a plurality of time periods determined by the cutting time point of the cut-off to-be-detected gather is counted, and the problem of uneven signal energy distribution caused by the offset distance change of the to-be-detected gather is considered. Furthermore, the amplitude adjustment is carried out by utilizing the energy adjustment factors of the gathers to be detected determined by the seismic channel numbers and the signal-to-noise ratios, so that exploration signals which are more in accordance with geological rules can be obtained in splicing areas and marginal areas of the connected blocks, and when the exploration signals are used for imaging the connected blocks, geological imaging results with higher quality can be obtained, so that the energy among the blocks is consistent under the condition that seismic data information is not lost.
Although the present application is described with reference to signal-to-noise ratio acquisition, energy adjustment factor calculation, smoothing processing, etc., the present application is not limited to what is necessarily described in the embodiments of the present application. Certain industry standards, or implementations modified slightly from those described using custom modes or examples, may also achieve the same, equivalent, or similar, or other, contemplated implementations of the above-described examples. Examples of data calculations, processing, etc. obtained using these modifications or variations may still fall within the scope of alternative embodiments of the present application.
Although the present application provides method steps as described in an embodiment or flowchart, more or fewer steps may be included based on conventional or non-inventive means. The order of steps recited in the embodiments is merely one manner of performing the steps in a multitude of orders and does not represent the only order of execution. When an actual apparatus or end product executes, it may execute sequentially or in parallel (e.g., parallel processors or multi-threaded environments, or even distributed data processing environments) according to the method shown in the embodiment or the figures. The terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, the presence of additional identical or equivalent elements in a process, method, article, or apparatus that comprises the recited elements is not excluded.
The units, devices, modules, etc. set forth in the above embodiments may be implemented by a computer chip or an entity, or by a product with certain functions. For convenience of description, the above devices are described as being divided into various modules by functions, and are described separately. Of course, in implementing the present application, the functions of each module may be implemented in one or more software and/or hardware, or a module implementing the same function may be implemented by a combination of a plurality of sub-modules or sub-units, and the like. The above-described embodiments of the apparatus are merely illustrative, and for example, the division of the units is only one logical division, and other divisions may be realized in practice, for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or units, and may be in an electrical, mechanical or other form.
Those skilled in the art will also appreciate that, in addition to implementing the controller as pure computer readable program code, the same functionality can be implemented by logically programming method steps such that the controller is in the form of logic gates, switches, application specific integrated circuits, programmable logic controllers, embedded microcontrollers and the like. Such a controller may therefore be considered as a hardware component, and the means included therein for performing the various functions may also be considered as a structure within the hardware component. Or even means for performing the functions may be regarded as being both a software module for performing the method and a structure within a hardware component.
The application may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, classes, etc. that perform particular tasks or implement particular abstract data types. The application may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.
From the above description of the embodiments, it is clear to those skilled in the art that the present application can be implemented by software plus necessary general hardware platform. Based on such understanding, the technical solutions of the present application may be embodied in the form of a software product, which may be stored in a storage medium, such as a ROM/RAM, a magnetic disk, an optical disk, or the like, and includes several instructions for enabling a computer device (which may be a personal computer, a mobile terminal, a server, or a network device) to execute the method according to the embodiments or some parts of the embodiments of the present application.
The embodiments in the present specification are described in a progressive manner, and the same or similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. The application is operational with numerous general purpose or special purpose computing system environments or configurations. For example: personal computers, server computers, hand-held or portable devices, tablet-type devices, multiprocessor systems, microprocessor-based systems, set top boxes, programmable electronic devices, network PCs, minicomputers, mainframe computers, distributed computing environments that include any of the above systems or devices, and the like.
While the present application has been described with examples, those of ordinary skill in the art will appreciate that there are numerous variations and permutations of the present application without departing from the spirit of the application, and it is intended that the appended claims encompass such variations and permutations without departing from the spirit of the application.

Claims (11)

1.一种地震资料振幅调整方法，其特征在于，包括：1. a seismic data amplitude adjustment method, is characterized in that, comprises: 获取待测道集的信噪比；Obtain the signal-to-noise ratio of the gather to be measured; 对所述待测道集进行动校正拉伸切除处理，并确定所述待测道集中各个地震道的切除时间点；Performing dynamic correction stretch cut processing on the to-be-measured gathers, and determining the cut-off time points of each seismic trace in the to-be-measured gathers; 根据所述待测道集中各个地震道的切除时间点，确定出多个时间段；According to the cut-off time point of each seismic trace in the to-be-measured trace set, a plurality of time periods are determined; 统计切除后的待测道集落入所述多个时间段中各个时间段的地震道数；Counting the number of seismic traces that the excised trace set to be measured falls into each of the multiple time periods; 根据所述信噪比和所述各个时间段的地震道数，计算得到所述待测道集在所述各个时间段的能量调整因子；According to the signal-to-noise ratio and the number of seismic traces in each time period, the energy adjustment factor of the to-be-measured gather in each time period is calculated and obtained; 基于所述能量调整因子，对所述待测道集进行振幅调整，得到调整后的地震信号；Based on the energy adjustment factor, performing amplitude adjustment on the gather to be measured to obtain an adjusted seismic signal; 其中，按照以下公式根据所述信噪比和所述各个时间段的地震道数，计算得到所述待测道集在所述各个时间段的能量调整因子：Wherein, according to the following formula, according to the signal-to-noise ratio and the number of seismic traces in each time period, the energy adjustment factor of the to-be-measured gather in each time period is calculated and obtained: cof(t)＝P*cov(t)-1/K+(1-P)*snr-1/L cof(t)=P*cov(t) -1/K +(1-P)*snr -1/L 其中，cof(t)表示所述待测道集在时间段t的能量调整因子，cov(t)表示所述待测道集在时间段t的地震道数，snr表示所述信噪比，0≤P≤1，K＞1，L＞0；where cof(t) represents the energy adjustment factor of the gather to be measured in time period t, cov(t) represents the number of seismic traces of the gather to be measured in time period t, snr represents the signal-to-noise ratio, 0≤P≤1, K>1, L>0; 其中，P表示在能量调整因子cof(t)中，覆盖次数cov(t)和信噪比snr所占权重，K表示覆盖次数cov(t)的平均程度，L表示信噪比snr的调整因子。Among them, P represents the weight of the coverage times cov(t) and the signal-to-noise ratio snr in the energy adjustment factor cof(t), K represents the average degree of the coverage times cov(t), and L represents the adjustment factor of the signal-to-noise ratio snr . 2.如权利要求1所述的方法，其特征在于，在获取待测道集的信噪比之前，所述方法还包括：2. The method of claim 1, wherein before acquiring the signal-to-noise ratio of the gather to be measured, the method further comprises: 获取待测工区的地震数据；Obtain the seismic data of the work area to be tested; 对所述地震数据进行预处理，得到处理后的地震数据，其中，所述预处理包括以下至少之一：叠前去噪、地表一致性振幅补偿、几何扩散吸收补偿和道分选；Preprocessing the seismic data to obtain processed seismic data, wherein the preprocessing includes at least one of the following: pre-stack denoising, surface consistent amplitude compensation, geometric diffusion and absorption compensation, and trace sorting; 根据所述处理后的地震数据，确定多个所述待测道集。According to the processed seismic data, a plurality of the gathers to be measured are determined. 3.如权利要求2所述的方法，其特征在于，在根据所述信噪比和所述各个时间段的地震道数，计算得到所述待测道集在所述各个时间段的能量调整因子之后，所述方法还包括：3 . The method according to claim 2 , wherein, according to the signal-to-noise ratio and the number of seismic traces in each time period, the energy adjustment of the gather to be measured in each time period is calculated and obtained. 4 . After the factor, the method further includes: 对所述待测道集的能量调整因子进行平滑处理，得到所述待测道集平滑处理后的能量调整因子。Smoothing the energy adjustment factors of the gathers to be measured is performed to obtain the smoothed energy adjustment factors of the gathers to be measured. 4.如权利要求3所述的方法，其特征在于，按照以下公式基于所述能量调整因子，对所述待测道集进行振幅调整，得到调整后的地震信号：4. The method according to claim 3, characterized in that, based on the energy adjustment factor, amplitude adjustment is performed on the gather to be measured according to the following formula to obtain an adjusted seismic signal: s'(t)＝s(t)·pcof(t)s'(t)=s(t)·pcof(t) 其中，s'(t)表示调整后的在时间段t的地震信号，s(t)表示所述待测道集调整前在时间段t的地震信号，pcof(t)表示所述待测道集在时间段t进行平滑处理后的能量调整因子。Among them, s'(t) represents the adjusted seismic signal in the time period t, s(t) represents the seismic signal in the time period t before the adjustment of the to-be-measured gather, and pcof(t) represents the to-be-measured trace Set the energy adjustment factor after smoothing at time period t. 5.如权利要求3所述的方法，其特征在于，对所述待测道集的能量调整因子进行平滑处理，得到所述待测道集平滑处理后的能量调整因子，包括：5. The method of claim 3, wherein the energy adjustment factor of the gather to be measured is smoothed to obtain the smoothed energy adjustment factor of the gather to be measured, comprising: 按照预定的顺序依次对所述待测工区中各个待测道集的位置进行标定，得到进行位置标定后的各个待测道集；According to a predetermined order, the positions of each gather to be measured in the work area to be measured are calibrated in turn, so as to obtain each gather to be measured after the position calibration is carried out; 按照以下方式获取各个待测道集平滑处理后的能量调整因子：以当前待测道集的位置为中心，预定的平滑半径为半径，得到平滑区域；根据所述平滑区域中各个待测道集的能量调整因子，计算得到所述当前待测道集平滑处理后的能量调整因子。Obtain the smoothed energy adjustment factor of each gather to be measured in the following manner: take the position of the current gather to be measured as the center and the predetermined smoothing radius as the radius to obtain a smooth area; The energy adjustment factor is calculated to obtain the energy adjustment factor after the smoothing of the current gather to be measured. 6.如权利要求5所述的方法，其特征在于，按照以下公式计算得到所述当前待测道集平滑处理后的能量调整因子：6. The method of claim 5, wherein the energy adjustment factor after the smoothing of the current gather to be measured is obtained by calculating according to the following formula: 其中，当时，wij＝0；Among them, when , w ij = 0; 当时， when hour, 其中，pcofkl表示位置为(k,l)的所述当前待测道集进行平滑处理后的能量调整因子，cofij表示位置为(i,j)的待测道集进行平滑处理前的能量调整因子，wij表示位置为(i,j)的待测道集的平滑系数，R表示所述平滑半径，R≥1，i＝k-R,k-R+1,…,k+R，j＝l-R,l-R+1,…,l+R。Wherein, pcof kl represents the energy adjustment factor of the current gather to be measured whose position is (k, l) after smoothing, and cof ij represents the energy of the gather to be tested whose position is (i, j) before smoothing. Adjustment factor, w ij represents the smoothing coefficient of the gather to be measured whose position is (i, j), R represents the smoothing radius, R≥1, i=kR, k-R+1,...,k+R,j =lR,l-R+1,...,l+R. 7.如权利要求6所述的方法，其特征在于，所述平滑半径是根据所述待测工区中待测道集的个数确定的。7 . The method of claim 6 , wherein the smoothing radius is determined according to the number of gathers to be measured in the work area to be measured. 8 . 8.如权利要求1所述的方法，其特征在于，按照以下公式获取待测道集的信噪比：8. The method of claim 1, wherein the signal-to-noise ratio of the gather to be measured is obtained according to the following formula: 其中， in, 其中，snr表示所述待测道集的信噪比，R(xc,xd)表示所述待测道集中第c个地震道与第d个地震道的互相关值，c＝1,2…,N，d＝1,2…,N，N表示所述第c个地震道分析时窗内的总样点数。Wherein, snr represents the signal-to-noise ratio of the gather to be measured, R(x c , x d ) represents the cross-correlation value of the c-th seismic trace and the d-th seismic trace in the to-be-measured gather, c=1, 2..., N, d=1, 2..., N, N represents the total number of samples in the analysis time window of the c-th seismic trace. 9.如权利要求8所述的方法，其特征在于，按照以下方式确定所述待测道集中第c个地震道与第d个地震道的互相关值：9. The method of claim 8, wherein the cross-correlation value of the c-th seismic trace and the d-th seismic trace in the gathers to be measured is determined in the following manner: 按照以下公式确定所述待测道集中第c个地震道与第d个地震道对应于地震道分析时窗内的各个样点序号的互相关值，得到与所述各个样点序号相对应的多个互相关结果：Determine the cross-correlation values of the c-th seismic trace and the d-th seismic trace in the to-be-measured trace set corresponding to each sample point number in the seismic trace analysis time window according to the following formula, and obtain the corresponding value of each sample point number Multiple cross-correlation results: 其中，t＝0,1,…,L，L表示所述第c个地震道分析时窗内的总样点数，τ表示地震道分析时窗内的样点序号，τ＝0,1,…,L；Among them, t=0,1,...,L, L represents the total number of sample points in the cth seismic trace analysis time window, τ represents the sample point number in the seismic trace analysis time window, τ=0,1,... ,L; 选取所述多个互相关结果中的最大值作为所述第c个地震道与第d个地震道的互相关值。The maximum value among the multiple cross-correlation results is selected as the cross-correlation value between the c-th seismic trace and the d-th seismic trace. 10.一种地震资料振幅调整装置，其特征在于，包括：10. A seismic data amplitude adjustment device, comprising: 信噪比获取模块，用于获取待测道集的信噪比；The signal-to-noise ratio acquisition module is used to obtain the signal-to-noise ratio of the gathers to be measured; 道集切除模块，用于对所述待测道集进行动校正拉伸切除处理，并确定所述待测道集中各个地震道的切除时间点；a gather excision module, configured to perform dynamic correction stretch excision processing on the gather to be measured, and determine the excision time point of each seismic trace in the gather to be measured; 时间段确定模块，用于根据所述待测道集中各个地震道的切除时间点，确定出多个时间段；a time period determination module, configured to determine a plurality of time periods according to the excision time points of each seismic trace in the to-be-measured trace set; 地震道数统计模块，用于统计切除后的待测道集落入所述多个时间段中各个时间段的地震道数；A seismic trace number statistics module, configured to count the seismic trace numbers of the excised trace sets to be measured that fall into each of the multiple time periods; 调整因子计算模块，用于根据所述信噪比和所述各个时间段的地震道数，计算得到所述待测道集在所述各个时间段的能量调整因子；an adjustment factor calculation module, configured to calculate the energy adjustment factor of the gathers to be measured in the respective time periods according to the signal-to-noise ratio and the number of seismic traces in the respective time periods; 振幅调整模块，用于基于所述能量调整因子，对所述待测道集进行振幅调整，得到调整后的地震信号；an amplitude adjustment module, configured to perform amplitude adjustment on the gathers to be measured based on the energy adjustment factor to obtain an adjusted seismic signal; 其中，所述调整因子计算模块按照以下公式根据所述信噪比和所述各个时间段的地震道数，计算得到所述待测道集在所述各个时间段的能量调整因子：Wherein, the adjustment factor calculation module calculates the energy adjustment factors of the gathers to be measured in the respective time periods according to the signal-to-noise ratio and the number of seismic traces in the respective time periods according to the following formula: cof(t)＝P*cov(t)-1/K+(1-P)*snr-1/L cof(t)=P*cov(t)-1/K+(1-P)*snr -1/L 其中，cof(t)表示所述待测道集在时间段t的能量调整因子，cov(t)表示所述待测道集在时间段t的地震道数，snr表示所述信噪比，0≤P≤1，K＞1，L＞0；where cof(t) represents the energy adjustment factor of the gather to be measured in time period t, cov(t) represents the number of seismic traces of the gather to be measured in time period t, snr represents the signal-to-noise ratio, 0≤P≤1, K>1, L>0; 其中，P表示在能量调整因子cof(t)中，覆盖次数cov(t)和信噪比snr所占权重，K表示覆盖次数cov(t)的平均程度，L表示信噪比snr的调整因子。Among them, P represents the weight of the coverage times cov(t) and the signal-to-noise ratio snr in the energy adjustment factor cof(t), K represents the average degree of the coverage times cov(t), and L represents the adjustment factor of the signal-to-noise ratio snr . 11.一种地震资料振幅调整装置，包括处理器以及用于存储处理器可执行指令的存储器，所述处理器执行所述指令时实现：11. An apparatus for adjusting the amplitude of seismic data, comprising a processor and a memory for storing instructions executable by the processor, the processor implementing the instructions to achieve: 获取待测道集的信噪比；对所述待测道集进行动校正拉伸切除处理，并确定所述待测道集中各个地震道的切除时间点；根据所述待测道集中各个地震道的切除时间点，确定出多个时间段；统计切除后的待测道集落入所述多个时间段中各个时间段的地震道数；根据所述信噪比和所述各个时间段的地震道数，计算得到所述待测道集在所述各个时间段的能量调整因子；基于所述能量调整因子，对所述待测道集进行振幅调整，得到调整后的地震信号；Acquire the signal-to-noise ratio of the gathers to be measured; perform dynamic correction and stretch excision processing on the gathers to be measured, and determine the excision time point of each seismic trace in the gathers to be measured; According to the cut-off time points of the traces, multiple time periods are determined; the number of seismic traces in each of the multiple time periods after the excision of the traces to be measured is counted; according to the signal-to-noise ratio and the The number of seismic traces, the energy adjustment factors of the gathers to be measured in the respective time periods are obtained by calculation; based on the energy adjustment factors, the amplitudes of the gathers to be measured are adjusted to obtain the adjusted seismic signals; 其中，按照以下公式根据所述信噪比和所述各个时间段的地震道数，计算得到所述待测道集在所述各个时间段的能量调整因子：Wherein, according to the following formula, according to the signal-to-noise ratio and the number of seismic traces in each time period, the energy adjustment factor of the to-be-measured gather in each time period is calculated and obtained: cof(t)＝P*cov(t)-1/K+(1-P)*snr-1/L cof(t)=P*cov(t) -1/K +(1-P)*snr -1/L 其中，cof(t)表示所述待测道集在时间段t的能量调整因子，cov(t)表示所述待测道集在时间段t的地震道数，snr表示所述信噪比，0≤P≤1，K＞1，L＞0；where cof(t) represents the energy adjustment factor of the gather to be measured in time period t, cov(t) represents the number of seismic traces of the gather to be measured in time period t, snr represents the signal-to-noise ratio, 0≤P≤1, K>1, L>0; 其中，P表示在能量调整因子cof(t)中，覆盖次数cov(t)和信噪比snr所占权重，K表示覆盖次数cov(t)的平均程度，L表示信噪比snr的调整因子。Among them, P represents the weight of the coverage times cov(t) and the signal-to-noise ratio snr in the energy adjustment factor cof(t), K represents the average degree of the coverage times cov(t), and L represents the adjustment factor of the signal-to-noise ratio snr .
CN201710372329.4A 2017-05-24 2017-05-24 A kind of seismic data amplitude adjustment method and device Active CN107065007B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710372329.4A CN107065007B (en) 2017-05-24 2017-05-24 A kind of seismic data amplitude adjustment method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710372329.4A CN107065007B (en) 2017-05-24 2017-05-24 A kind of seismic data amplitude adjustment method and device

Publications (2)

Publication Number Publication Date
CN107065007A CN107065007A (en) 2017-08-18
CN107065007B true CN107065007B (en) 2019-02-15

Family

ID=59610602

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710372329.4A Active CN107065007B (en) 2017-05-24 2017-05-24 A kind of seismic data amplitude adjustment method and device

Country Status (1)

Country Link
CN (1) CN107065007B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109581506B (en) * 2018-12-24 2021-01-05 核工业北京地质研究院 A Data-Driven Stretch-Free Motion Correction Method
CN112666601B (en) * 2019-10-15 2024-03-22 中国石油化工股份有限公司 Method and system for fitting effective time window to seismic data amplitude
CN113109872B (en) * 2021-04-07 2022-06-21 中国煤炭地质总局物测队 Pre-stack amplitude compensation processing method based on coverage frequency weighting
CN114325826A (en) * 2021-11-30 2022-04-12 中交公路规划设计院有限公司 Automatic equalization method and device for engineering seismic exploration energy
CN120233428A (en) * 2023-12-28 2025-07-01 中国石油天然气集团有限公司 A method, device, equipment and storage medium for adjusting energy of pre-stacked joints

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102636813A (en) * 2012-05-04 2012-08-15 郭平 Normal moveout stretch cutting method for processing geophysical exploitation seismic data
CN104536048A (en) * 2014-07-08 2015-04-22 中国石油集团东方地球物理勘探有限责任公司 Method for automatically cutting seismic wave arrival

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6748329B2 (en) * 2000-12-08 2004-06-08 Halliburton Energy Services, Inc. Acoustic signal processing method using array coherency
WO2010059036A1 (en) * 2008-11-21 2010-05-27 Magali Hadar Rebound-effector

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102636813A (en) * 2012-05-04 2012-08-15 郭平 Normal moveout stretch cutting method for processing geophysical exploitation seismic data
CN104536048A (en) * 2014-07-08 2015-04-22 中国石油集团东方地球物理勘探有限责任公司 Method for automatically cutting seismic wave arrival

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
泌阳凹陷陡坡带三维地震资料叠前时间偏移连片处理实践;朱景修 等;《石油地质与工程》;20111130;第25卷(第6期);第42-45页

Also Published As

Publication number Publication date
CN107065007A (en) 2017-08-18

Similar Documents

Publication Publication Date Title
CN107065007B (en) A kind of seismic data amplitude adjustment method and device
Küperkoch et al. Automated determination of P-phase arrival times at regional and local distances using higher order statistics
CN103235179B (en) Frequency mask trigger with uneven bandwidth
Franzen et al. Source counts and confusion at 72–231 MHz in the MWA GLEAM survey
Du et al. Earthquake relocation using cross-correlation time delay estimates verified with the bispectrum method
RU2015143896A (en) DEVICE AND METHOD FOR DETERMINING BASIC PHYSIOLOGICAL INDICATORS OF A SUBJECT
CN105372506B (en) A kind of mutative scale gridding method calculated for region electromagnetic environment and system
CN109425894A (en) A kind of seismic anomaly road detection method and device
Persico et al. A microwave tomographic imaging approach for multibistatic configuration: The choice of the frequency step
RU2014135771A (en) METHOD AND DEVICE FOR PROCESSING SEISMIC DATA
Oelze et al. Roughness measurements of soil surfaces by acoustic backscatter
CN104122586A (en) Wireless node data monitoring method of seismic exploration
CN116577829B (en) Automatic extraction method based on background noise dispersion curve
CN108562936A (en) A kind of crack prediction method, system, storage medium and terminal
CN114924328B (en) A method and system for urban artificial source electromagnetic exploration with vertical magnetic field reference track
Abraham et al. Use of a handheld broadband EM induction system for deriving resistivity depth images
WO2015138485A2 (en) Hydrophone response compensation filter derivation, design and application
CN118444312A (en) Evaluation method of grouting effect of asphalt pavement cracks based on ground penetrating radar signal characteristics
CN108873077A (en) A kind of new natural electric field exploitation method
Lehujeur et al. Vertical seismic profiling using double-beamforming processing of nonuniform anthropogenic seismic noise: The case study of Rittershoffen, Upper Rhine Graben, France
CN103760599B (en) A kind of miniature fault detection method and fault detection device
Stoneback et al. In situ irregularity identification and scintillation estimation using wavelets and CINDI on C/NOFS
RU2685048C1 (en) Method of determining places of inhomogeneities and damage of power transmission lines
CN111142160B (en) Analysis method and device for time lapse seismic observation data
JP2010529475A (en) Method and apparatus for monitoring storm activity on the surface of the earth in real time

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