CN114814930A - Method and device for determining distortion of vibroseis force signal - Google Patents

Method and device for determining distortion of vibroseis force signal Download PDF

Info

Publication number
CN114814930A
CN114814930A CN202110066751.3A CN202110066751A CN114814930A CN 114814930 A CN114814930 A CN 114814930A CN 202110066751 A CN202110066751 A CN 202110066751A CN 114814930 A CN114814930 A CN 114814930A
Authority
CN
China
Prior art keywords
amplitude
vibroseis
time window
value
positive
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202110066751.3A
Other languages
Chinese (zh)
Inventor
周恒�
白杰
汪长辉
蒋忠
陈永胜
徐学峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN202110066751.3A priority Critical patent/CN114814930A/en
Publication of CN114814930A publication Critical patent/CN114814930A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/02—Generating seismic energy
    • G01V1/04—Details
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/02—Generating seismic energy

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 method and a device for determining the distortion of a vibroseis force signal, wherein the method comprises the following steps: selecting discrete shot point data with the same number as the sampling points of the vibroseis force signal as a reference signal; determining the number of sampling points with different sampling point symbol values between two signals; averagely dividing the vibroseis force signal into a plurality of time windows, determining the maximum positive amplitude and the maximum negative amplitude of the vibroseis force signal in each time window, and determining the average amplitude ratio of the time windows according to the maximum positive amplitude and the maximum negative amplitude of the reference signal; and determining the distortion characteristic value of the vibroseis force signal according to the average amplitude ratio of the time window, the number of sampling points with different sampling point symbol values between the two signals and the number of the sampling points of the vibroseis force signal. Compared with the prior art, the method has the advantages that comprehensive judgment is not needed by combining the distortion attribute value and the earth surface condition within a period of time, the vibroseis force signal distortion result can be simply and effectively determined, and the accuracy of the determination result is improved.

Description

Method and device for determining distortion of vibroseis force signal
Technical Field
The invention relates to the technical field of geophysical exploration, in particular to a method and a device for determining the distortion of a vibroseis force signal.
Background
In land oil seismic exploration, a controllable seismic source is a safe and environment-friendly excitation device, the controllable seismic source vibrates by taking a regular reference signal as a reference, a ground force signal is generated through the interaction of a heavy hammer and a flat plate and the response of the ground, and the generated force signal is inevitably distorted due to the nonlinear characteristic of a mechanical driving device of the controllable seismic source and the coupling relation between a vibration device and the ground, wherein the distortion is reflected in two aspects, namely phase and amplitude. Currently, in field seismic exploration, the evaluation of vibroseis distortion is generally realized by 4 variables, namely: average phase, maximum phase, average amplitude distortion and maximum amplitude distortion, wherein the variables are in the range of 0-100, and the larger the value is, the larger the distortion is, the worse the quality of the signal is.
However, in practical applications, it has been found that there are some problems in determining the distortion of the vibroseis force signal by using the above 4 variables. In the prior art, threshold values are generally set for the 4 variables, and once the set threshold values are exceeded, the distortion degree of the vibroseis force signal is considered to be too large. However, under some surface conditions, these 4 variables easily exceed the threshold, and the actual vibroseis force signal is not distorted to a high degree. Therefore, when the distortion of the vibroseis force signal is determined, the distortion attribute value and the earth surface condition within a period of time are combined for comprehensive judgment, so that the difficulty of determining the distortion of the vibroseis force signal is improved, and the determination result is not accurate enough.
Disclosure of Invention
The embodiment of the invention provides a method for determining the distortion of a vibroseis force signal, which is used for simply and effectively determining the distortion result of the vibroseis force signal and improving the accuracy of the determination result, and comprises the following steps:
selecting discrete shot point data with the same number as the sampling points of the vibroseis force signal as a reference signal;
determining the number of sampling points with different sampling point symbol values between two signals according to the sampling point symbol values of the controllable seismic source force signal and the reference signal;
averagely dividing the vibroseis force signal and the reference signal into a plurality of time windows, and determining the maximum positive amplitude value and the maximum negative amplitude value of the vibroseis force signal and the maximum positive amplitude value and the maximum negative amplitude value of the reference signal in each time window;
determining a time window average amplitude ratio according to the positive amplitude maximum value and the negative amplitude maximum value of the vibroseis force signals in the multiple time windows and the positive amplitude maximum value and the negative amplitude maximum value of the reference signal;
and determining the distortion characteristic value of the vibroseis force signal according to the average amplitude ratio of the time window, the number of sampling points with different sampling point symbol values between the two signals and the number of the sampling points of the vibroseis force signal.
The embodiment of the invention also provides a device for determining the distortion of the vibroseis force signal, which is used for simply and effectively determining the distortion result of the vibroseis force signal and improving the accuracy of the determination result, and comprises the following steps:
the reference signal selecting module is used for selecting discrete shot point data with the same number as the sampling points of the vibroseis force signal as a reference signal;
the sampling point symbol counting module is used for determining the number of sampling points with different sampling point symbol values between two signals according to the sampling point symbol values of the vibroseis force signal and the reference signal;
the amplitude maximum value determining module is used for averagely dividing the vibroseis force signal and the reference signal into a plurality of time windows, and determining the positive amplitude maximum value and the negative amplitude maximum value of the vibroseis force signal, and the positive amplitude maximum value and the negative amplitude maximum value of the reference signal in each time window;
the time window average amplitude ratio determining module is used for determining the time window average amplitude ratio according to the maximum positive amplitude value and the maximum negative amplitude value of the vibroseis force signal in the time windows and the maximum positive amplitude value and the maximum negative amplitude value of the reference signal;
and the distortion characteristic value determining module is used for determining the distortion characteristic value of the vibroseis force signal according to the time window average amplitude ratio, the number of sampling points with different sampling point symbol values between the two signals and the number of the sampling points of the vibroseis force signal.
The embodiment of the invention also provides computer equipment which comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor executes the computer program to realize the method for determining the distortion of the vibroseis force signal.
Embodiments of the present invention also provide a computer-readable storage medium storing a computer program for executing the method for determining distortion of a vibroseis force signal.
In the embodiment of the invention, discrete shot point data with the same number as the sampling points of the vibroseis force signal is selected as a reference signal; determining the number of sampling points with different sampling point symbol values between two signals according to the sampling point symbol values of the vibroseis force signal and the reference signal; averagely dividing the vibroseis force signal and the reference signal into a plurality of time windows, and determining the maximum positive amplitude value and the maximum negative amplitude value of the vibroseis force signal and the maximum positive amplitude value and the maximum negative amplitude value of the reference signal in each time window; determining a time window average amplitude ratio according to the positive amplitude maximum value and the negative amplitude maximum value of the vibroseis force signals in the multiple time windows and the positive amplitude maximum value and the negative amplitude maximum value of the reference signal; and determining the distortion characteristic value of the vibroseis force signal according to the average amplitude ratio of the time window, the number of sampling points with different sampling point symbol values between the two signals and the number of the sampling points of the vibroseis force signal. The reference signal is closest to the theoretical value of the vibroseis force signal, and the time window average amplitude ratio is determined by taking the reference signal as a reference, so that the distortion degree of the vibroseis force signal can be accurately reflected; compared with the prior art, the method has the advantages that comprehensive judgment is not needed by combining the distortion attribute value and the earth surface condition within a period of time, the vibroseis force signal distortion result can be simply and effectively determined, and the accuracy of the determination result is improved.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
Fig. 1 is a schematic diagram of a method for determining distortion of a vibroseis force signal in an embodiment of the invention.
Fig. 2 is a schematic diagram of a specific implementation method of step 103 in an embodiment of the present invention.
Fig. 3 is a schematic diagram of a specific implementation method of step 104 in an embodiment of the present invention.
Fig. 4 is a schematic diagram of a vibroseis reference signal in an implementation of the present invention.
Fig. 5 is a schematic diagram of a vibroseis force signal in an implementation of an embodiment of the invention.
Fig. 6 is a schematic diagram of a vibroseis force signal with a distortion eigenvalue of 8.01 in an implementation of one embodiment of the present invention.
Fig. 7 is a schematic diagram of a vibroseis force signal with a distortion eigenvalue of 157.1 in an implementation of one embodiment of the present invention.
Fig. 8 is a schematic representation of a vibroseis force signal with a distortion signature value of 315.8 in an implementation of one embodiment of the present invention.
Fig. 9 is a schematic diagram of a vibroseis force signal distortion determining apparatus in an embodiment of the invention.
Fig. 10 is a schematic structural diagram of the maximum amplitude determination module 903 according to an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, 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 invention.
The embodiment of the invention provides a method for determining the distortion of a vibroseis force signal, which is used for simply and effectively determining the distortion result of the vibroseis force signal and improving the accuracy of the determination result, and as shown in figure 1, the method comprises the following steps:
step 101: selecting discrete shot point data with the same number as the sampling points of the vibroseis force signal as a reference signal;
step 102: determining the number of sampling points with different sampling point symbol values between two signals according to the sampling point symbol values of the vibroseis force signal and the reference signal;
step 103: averagely dividing the vibroseis force signal and the reference signal into a plurality of time windows, and determining the maximum positive amplitude value and the maximum negative amplitude value of the vibroseis force signal and the maximum positive amplitude value and the maximum negative amplitude value of the reference signal in each time window;
step 104: determining a time window average amplitude ratio according to the positive amplitude maximum value and the negative amplitude maximum value of the vibroseis force signals in the multiple time windows and the positive amplitude maximum value and the negative amplitude maximum value of the reference signal;
step 105: and determining the distortion characteristic value of the vibroseis force signal according to the average amplitude ratio of the time window, the number of sampling points with different sampling point symbol values between the two signals and the number of the sampling points of the vibroseis force signal.
As can be known from the process shown in fig. 1, in the embodiment of the present invention, the discrete shot data with the same number of sampling points as the vibroseis force signal is selected as the reference signal; determining the number of sampling points with different sampling point symbol values between two signals according to the sampling point symbol values of the vibroseis force signal and the reference signal; averagely dividing the vibroseis force signal and the reference signal into a plurality of time windows, and determining the maximum positive amplitude value and the maximum negative amplitude value of the vibroseis force signal and the maximum positive amplitude value and the maximum negative amplitude value of the reference signal in each time window; determining a time window average amplitude ratio according to the positive amplitude maximum value and the negative amplitude maximum value of the vibroseis force signals in the multiple time windows and the positive amplitude maximum value and the negative amplitude maximum value of the reference signal; and determining the distortion characteristic value of the vibroseis force signal according to the average amplitude ratio of the time window, the number of sampling points with different sampling point symbol values between the two signals and the number of the sampling points of the vibroseis force signal. The reference signal is closest to the theoretical value of the vibroseis force signal, and the time window average amplitude ratio is determined by taking the reference signal as a reference, so that the distortion degree of the vibroseis force signal can be accurately reflected; compared with the prior art, the method has the advantages that comprehensive judgment is not needed by combining the distortion attribute value and the earth surface condition within a period of time, the vibroseis force signal distortion result can be simply and effectively determined, and the accuracy of the determination result is improved.
In specific implementation, firstly, the discrete shot point data with the same number of sampling points as the vibroseis force signal is selected as a reference signal. In practice, the reference signal is generally the signal closest to the ideal state of the vibroseis force signal in the actual measurement.
And after the reference signal is selected, determining the number of sampling points with different sampling point symbol values between the two signals according to the sampling point symbol values of the vibroseis force signal and the reference signal. The time is used as an abscissa of the vibroseis force signal and the time is used as a reference signal, the signal amplitude is used as an ordinate, the sign values of the sampling points of the vibroseis force signal and the reference signal at the same moment are positive or negative, and the number n of the sampling points of the vibroseis force signal and the reference signal at the same moment are determined.
Then, the vibroseis force signal and the reference signal are averagely divided into a plurality of time windows, and the maximum positive amplitude value and the maximum negative amplitude value of the vibroseis force signal and the maximum positive amplitude value and the maximum negative amplitude value of the reference signal are determined in each time window. The specific implementation process, as shown in fig. 2, includes:
step 201: dividing the vibroseis force signal and the reference signal into a plurality of time windows, and respectively carrying out curve fitting of discrete sampling points on the vibroseis force signal and the reference signal in each time window to obtain a vibration curve of the vibroseis force signal and the reference signal in each time window;
step 202: according to the positive and negative of the sign value of the amplitude on the vibration curve of the reference signal, respectively selecting the amplitude value with the largest numerical value from the positive sign value and the negative sign value, and determining the amplitude value as the maximum positive amplitude value and the maximum negative amplitude value of the vibroseis force signal in each time window;
step 203: and according to the positive and negative of the sign value of the amplitude on the vibration curve of the vibroseis force signal, respectively selecting the amplitude value with the largest numerical value from the positive sign value and the negative sign value, and determining the amplitude value as the maximum positive amplitude value and the maximum negative amplitude value of the reference signal in each time window.
And after determining the maximum positive amplitude value and the maximum negative amplitude value of the vibroseis force signal and the maximum positive amplitude value and the maximum negative amplitude value of the reference signal in each time window, determining the average amplitude ratio of the time windows according to the maximum positive amplitude value and the maximum negative amplitude value of the vibroseis force signal and the maximum positive amplitude value and the maximum negative amplitude value of the reference signal in the plurality of time windows. In specific implementation, as shown in fig. 3, the method includes:
step 301: determining a positive amplitude ratio in a time window according to the maximum positive amplitude of the vibroseis force signal in the time window and the maximum positive amplitude of the reference signal;
step 302: determining a negative amplitude ratio in a time window according to the maximum value of the negative amplitude of the vibroseis force signal in the time window and the maximum value of the negative amplitude of the reference signal;
step 303: determining a positive-negative amplitude ratio difference value corresponding to the time window according to the positive amplitude ratio in the time window and the negative amplitude ratio in the time window;
step 304: determining the positive amplitude ratio difference value of the adjacent time window corresponding to the time window according to the positive amplitude ratio in the time window and the positive amplitude ratio of the next time window;
step 305: determining the negative amplitude ratio difference value of the adjacent time window corresponding to the time window according to the negative amplitude ratio in the time window and the negative amplitude ratio of the next time window;
step 306: and determining the average amplitude ratio of the time windows according to the positive and negative amplitude ratio difference values corresponding to the plurality of time windows, the positive amplitude ratio difference value of the adjacent time windows, the negative amplitude ratio difference value of the adjacent time windows and the number of the time windows.
In step 306, when the method is implemented, the average amplitude ratio of the time windows is determined according to the positive-negative amplitude ratio difference corresponding to the plurality of time windows, the positive amplitude ratio difference of the adjacent time windows, the negative amplitude ratio difference of the adjacent time windows, and the number of the time windows according to the following formula:
Figure BDA0002904483960000061
Figure BDA0002904483960000062
wherein V represents the average amplitude ratio of the time windows; w represents the number of time windows; v i Representing the amplitude ratio of each time window; i represents the number of time windows; f + (max) i Represents the maximum value of the positive amplitude of the vibroseis force signal in the ith time window; r + (max) i Represents the maximum value of the positive amplitude of the reference signal in the ith time window; f - (max) i Represents the maximum value of the negative amplitude of the vibroseis force signal in the ith time window; r - (max) i Represents the maximum value of the negative amplitude of the reference signal in the ith time window; f + (max) i+1 Represents the maximum value of the positive amplitude of the vibroseis force signal in the (i + 1) th time window; r + (max) i+1 The maximum value of the positive amplitude of the reference signal in the (i + 1) th time window; f - (max) i+1 The maximum value of the negative amplitude of the vibroseis force signal in the (i + 1) th time window; r - (max) i+1 The negative amplitude maximum of the reference signal in the (i + 1) th time window.
And after the time window average amplitude ratio is determined, determining the distortion characteristic value of the vibroseis force signal according to the time window average amplitude ratio, the number of sampling points with different sampling point symbol values between the two signals and the number of the sampling points of the vibroseis force signal. In the specific implementation process, the distortion characteristic value Q of the vibroseis force signal is determined according to the following formula:
Figure BDA0002904483960000063
wherein V represents the average amplitude ratio of the time window; n represents the number of sampling points with different sampling point symbol values between two signals; n represents the number of samples of the vibroseis force signal.
The larger the distortion characteristic value of the vibroseis force signal is, the larger the distortion degree of the vibroseis force signal is, namely, the worse the quality of the vibroseis force signal is.
A specific example is given below to illustrate how embodiments of the present invention determine the degree of distortion of the vibroseis force signal. The working steps of the example are as follows:
given a discrete reference signal and a vibroseis force signal, as shown in fig. 4 and 5, the number of samples of both signals is N, and a time window length is set such that both signals are equally divided into w time windows.
And counting the number N of sampling points with different sign positions between the vibroseis force signal and the reference signal, and solving the proportion N/N of the sampling points in the total sampling points.
Calculating the maximum amplitude value in each time window by respectively using the reference signal and the vibroseis force signal: in the ith time window, the maximum value of numerical values with positive and negative numerical values in the vibroseis force signal is respectively marked as F + (max) i And F - (max) i The maximum values of the positive and negative values in the reference signal are denoted as R + (max) i And R - (max) i Calculating the ratio of the maximum value in the time window of the reference signal and the force signal in the ith time window, and recording the ratio with positive value as
Figure BDA0002904483960000071
The ratio of negative values is recorded
Figure BDA0002904483960000072
The absolute value of the difference between the two ratios is recorded as
Figure BDA0002904483960000073
abs is a function for calculating the absolute value, and the absolute value of the ratio difference between the value of the ith time window and the value of the (i + 1) th time window is calculated and recorded as
Figure BDA0002904483960000074
And calculating the absolute value of the ratio difference between the value of the ith time window and the value of the (i + 1) th time window, and recording the absolute value as
Figure BDA0002904483960000075
Summing these ratios, we note:
Figure BDA0002904483960000076
then averaging all the time windows to obtain the average amplitude ratio of the time windows
Figure BDA0002904483960000077
Multiplying the average amplitude ratio of the time window and the result of N/N to obtain a statistic value
Figure BDA0002904483960000078
The larger this statistic is, the worse the vibroseis force signal quality is, the smaller the value is, the smaller the signal distortion is, and the better the quality is.
Fig. 6, 7 and 8 represent the distortion of three different force signals, with distortion characteristic values Q of 8.01, 157.1 and 315.8, respectively. Meanwhile, the distortion degree of the vibroseis force signal can be visually observed from the graph, the distortion degree of the force signal shown in fig. 6 is minimum, so that compared with a reference signal, the coincidence ratio is higher, the distortion of the force signal shown in fig. 8 is the most serious, the distortion characteristic value is also the largest, and the distortion condition of the force signal shown in fig. 7 is between the two. Thus proving the feasibility of the method for quantifying and representing the distortion magnitude of the vibroseis force signal by using a characteristic value.
The implementation of the above specific application is only an example, and the rest of the embodiments are not described in detail.
Based on the same inventive concept, embodiments of the present invention further provide a device for determining distortion of a vibroseis force signal, where the principle of the problem solved by the device for determining distortion of a vibroseis force signal is similar to that of the method for determining distortion of a vibroseis force signal, so that the implementation of the device for determining distortion of a vibroseis force signal may refer to the implementation of the method for determining distortion of a vibroseis force signal, and repeated details are omitted, and the specific structure is shown in fig. 9:
a reference signal selecting module 901, configured to select discrete shot point data with the same number of sampling points as the vibroseis force signal as a reference signal;
a sampling point symbol counting module 902, configured to determine, according to sampling point symbol values of the vibroseis force signal and the reference signal, the number of sampling points with different sampling point symbol values between the two signals;
an amplitude maximum determination module 903, configured to averagely divide the vibroseis force signal and the reference signal into a plurality of time windows, and determine a positive amplitude maximum and a negative amplitude maximum of the vibroseis force signal, and a positive amplitude maximum and a negative amplitude maximum of the reference signal in each time window;
a time window average amplitude ratio determining module 904 for determining a time window average amplitude ratio based on the positive amplitude maximum and the negative amplitude maximum of the vibroseis force signal and the positive amplitude maximum and the negative amplitude maximum of the reference signal within the plurality of time windows;
and a distortion characteristic value determining module 905, configured to determine a distortion characteristic value of the vibroseis force signal according to the time window average amplitude ratio, the number of sampling points with different sampling point symbol values between the two signals, and the number of sampling points of the vibroseis force signal.
In a specific embodiment, the structure of the maximum amplitude value determining module 903 is shown in fig. 10, and includes:
a curve fitting unit 1001, configured to divide the vibroseis force signal and the reference signal into multiple time windows, and perform curve fitting of discrete sampling points on the vibroseis force signal and the reference signal in each time window, to obtain a vibration curve of the vibroseis force signal and the reference signal in each time window;
a reference signal amplitude maximum value determining unit 1002, configured to select, according to the positive and negative of the sign value of the amplitude on the vibration curve of the reference signal, an amplitude value with the largest magnitude from the positive sign value and the negative sign value, and determine the amplitude value as a positive amplitude maximum value and a negative amplitude maximum value of the vibroseis force signal in each time window;
the maximum force signal amplitude determination unit 1003 is configured to select, according to the positive and negative signs of the sign value of the amplitude on the vibration curve of the vibroseis force signal, an amplitude value with the largest magnitude from the positive sign value and the negative sign value, and determine the amplitude value as the maximum positive amplitude value and the maximum negative amplitude value of the reference signal in each time window.
In a specific embodiment, the time window average amplitude ratio determining module 904 is specifically configured to:
determining a positive amplitude ratio in a time window according to the maximum positive amplitude of the vibroseis force signal in the time window and the maximum positive amplitude of the reference signal;
determining a negative amplitude ratio in a time window according to the maximum value of the negative amplitude of the vibroseis force signal in the time window and the maximum value of the negative amplitude of the reference signal;
determining a positive-negative amplitude ratio difference value corresponding to the time window according to the positive amplitude ratio in the time window and the negative amplitude ratio in the time window;
determining the positive amplitude ratio difference value of the adjacent time window corresponding to the time window according to the positive amplitude ratio in the time window and the positive amplitude ratio of the next time window;
determining the negative amplitude ratio difference value of the adjacent time window corresponding to the time window according to the negative amplitude ratio in the time window and the negative amplitude ratio of the next time window;
and determining the average amplitude ratio of the time windows according to the positive and negative amplitude ratio difference values corresponding to the plurality of time windows, the positive amplitude ratio difference value of the adjacent time windows, the negative amplitude ratio difference value of the adjacent time windows and the number of the time windows.
In a specific implementation, the time window average amplitude ratio determining module 904 is specifically configured to:
determining the average amplitude ratio of the time windows according to the positive and negative amplitude ratio difference values corresponding to the time windows, the positive amplitude ratio difference value of the adjacent time windows, the negative amplitude ratio difference value of the adjacent time windows and the number of the time windows according to the following formula:
Figure BDA0002904483960000091
Figure BDA0002904483960000092
wherein V represents the average amplitude ratio of the time window; w represents the number of time windows; v i Representing the amplitude ratio of each time window; i represents the number of time windows; f + (max) i Represents the maximum value of the positive amplitude of the vibroseis force signal in the ith time window; r + (max) i Represents the maximum value of the positive amplitude of the reference signal in the ith time window; f - (max) i Represents the maximum value of the negative amplitude of the vibroseis force signal in the ith time window; r - (max) i Represents the maximum value of the negative amplitude of the reference signal in the ith time window; f + (max) i+1 Represents the maximum value of the positive amplitude of the vibroseis force signal in the (i + 1) th time window; r + (max) i+1 The maximum value of the positive amplitude of the reference signal in the (i + 1) th time window; f - (max) i+1 The maximum value of the negative amplitude of the vibroseis force signal in the (i + 1) th time window; r - (max) i+1 The negative amplitude maximum of the reference signal in the (i + 1) th time window.
In a specific embodiment, the distortion characteristic value determining module 905 is specifically configured to:
determining the distortion characteristic value of the vibroseis force signal according to the average amplitude ratio of the time window, the number of sampling points with different sampling point symbol values between two signals and the number of the sampling points of the vibroseis force signal according to the following formula:
Figure BDA0002904483960000093
wherein V represents the average amplitude ratio of the time window; n represents the number of sampling points with different sampling point symbol values between two signals; n represents the number of samples of the vibroseis force signal.
The embodiment of the invention also provides computer equipment which comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor executes the computer program to realize the method for determining the distortion of the vibroseis force signal.
Embodiments of the present invention further provide a computer-readable storage medium storing a computer program for executing the method for determining distortion of a vibroseis force signal.
In summary, the method and the device for determining the distortion of the vibroseis force signal provided by the embodiment of the invention have the following advantages:
selecting discrete shot point data with the same number as the sampling points of the vibroseis force signal as a reference signal; determining the number of sampling points with different sampling point symbol values between two signals according to the sampling point symbol values of the vibroseis force signal and the reference signal; averagely dividing the vibroseis force signal and the reference signal into a plurality of time windows, and determining the maximum positive amplitude value and the maximum negative amplitude value of the vibroseis force signal and the maximum positive amplitude value and the maximum negative amplitude value of the reference signal in each time window; determining a time window average amplitude ratio according to the positive amplitude maximum value and the negative amplitude maximum value of the vibroseis force signals in the multiple time windows and the positive amplitude maximum value and the negative amplitude maximum value of the reference signal; and determining the distortion characteristic value of the vibroseis force signal according to the average amplitude ratio of the time window, the number of sampling points with different sampling point symbol values between the two signals and the number of the sampling points of the vibroseis force signal. The reference signal is closest to the theoretical value of the vibroseis force signal, and the time window average amplitude ratio is determined by taking the reference signal as a reference, so that the distortion degree of the vibroseis force signal can be accurately reflected; compared with the prior art, the method has the advantages that comprehensive judgment is not needed by combining the distortion attribute value and the earth surface condition within a period of time, the vibroseis force signal distortion result can be simply and effectively determined, and the accuracy of the determination result is improved.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, apparatus, 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 and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above description is only a preferred embodiment of the present invention, and is not intended to limit the present invention, and various modifications and changes may be made to the embodiment of the present invention by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (12)

1. A method for determining distortion of a vibroseis force signal, comprising:
selecting discrete shot point data with the same number as the sampling points of the vibroseis force signal as a reference signal;
determining the number of sampling points with different sampling point symbol values between two signals according to the sampling point symbol values of the vibroseis force signal and the reference signal;
averagely dividing the vibroseis force signal and the reference signal into a plurality of time windows, and determining the maximum positive amplitude value and the maximum negative amplitude value of the vibroseis force signal and the maximum positive amplitude value and the maximum negative amplitude value of the reference signal in each time window;
determining a time window average amplitude ratio according to the positive amplitude maximum value and the negative amplitude maximum value of the vibroseis force signals in the multiple time windows and the positive amplitude maximum value and the negative amplitude maximum value of the reference signal;
and determining the distortion characteristic value of the vibroseis force signal according to the average amplitude ratio of the time window, the number of sampling points with different sampling point symbol values between the two signals and the number of the sampling points of the vibroseis force signal.
2. The method of claim 1, wherein the averaging of the vibroseis force signal and the reference signal into a plurality of time windows, determining a positive amplitude maximum and a negative amplitude maximum of the vibroseis force signal and a positive amplitude maximum and a negative amplitude maximum of the reference signal within each time window, comprises:
dividing the vibroseis force signal and the reference signal into a plurality of time windows, and respectively carrying out curve fitting of discrete sampling points on the vibroseis force signal and the reference signal in each time window to obtain a vibration curve of the vibroseis force signal and the reference signal in each time window;
according to the positive and negative of the sign value of the amplitude on the vibration curve of the reference signal, respectively selecting the amplitude value with the largest numerical value from the positive sign value and the negative sign value, and determining the amplitude value as the maximum positive amplitude value and the maximum negative amplitude value of the vibroseis force signal in each time window;
and according to the positive and negative of the sign value of the amplitude on the vibration curve of the vibroseis force signal, respectively selecting the amplitude value with the largest numerical value from the positive sign value and the negative sign value, and determining the amplitude value as the maximum positive amplitude value and the maximum negative amplitude value of the reference signal in each time window.
3. The method of claim 1, wherein determining the time-window average amplitude ratio based on the positive and negative amplitude maxima of the vibroseis force signal and the positive and negative amplitude maxima of the reference signal over the plurality of time windows comprises:
determining a positive amplitude ratio in a time window according to the maximum positive amplitude of the vibroseis force signal in the time window and the maximum positive amplitude of the reference signal;
determining a negative amplitude ratio in a time window according to the maximum value of the negative amplitude of the vibroseis force signal in the time window and the maximum value of the negative amplitude of the reference signal;
determining a positive-negative amplitude ratio difference value corresponding to the time window according to the positive amplitude ratio in the time window and the negative amplitude ratio in the time window;
determining the positive amplitude ratio difference value of the adjacent time window corresponding to the time window according to the positive amplitude ratio in the time window and the positive amplitude ratio of the next time window;
determining the negative amplitude ratio difference value of the adjacent time window corresponding to the time window according to the negative amplitude ratio in the time window and the negative amplitude ratio of the next time window;
and determining the average amplitude ratio of the time windows according to the positive and negative amplitude ratio difference values corresponding to the plurality of time windows, the positive amplitude ratio difference value of the adjacent time windows, the negative amplitude ratio difference value of the adjacent time windows and the number of the time windows.
4. The method of claim 3, wherein the window average amplitude ratio is determined based on the positive and negative amplitude ratio differences for a plurality of time windows, the positive amplitude ratio difference for adjacent time windows, the negative amplitude ratio difference for adjacent time windows, and the number of time windows according to the following equation:
Figure FDA0002904483950000021
Figure FDA0002904483950000022
wherein V represents the average amplitude ratio of the time window; w represents the number of time windows; v i Representing the amplitude ratio of each time window; i represents the number of time windows; f + (max) i Represents the maximum value of the positive amplitude of the vibroseis force signal in the ith time window; r + (max) i Represents the maximum value of the positive amplitude of the reference signal in the ith time window; f - (max) i Represents the maximum value of the negative amplitude of the vibroseis force signal in the ith time window; r - (max) i Represents the maximum value of the negative amplitude of the reference signal in the ith time window; f + (max) i+1 Represents the maximum value of the positive amplitude of the vibroseis force signal in the (i + 1) th time window; r + (max) i+1 The maximum value of the positive amplitude of the reference signal in the (i + 1) th time window; f - (max) i+1 The maximum value of the negative amplitude of the vibroseis force signal in the (i + 1) th time window; r - (max) i+1 The negative amplitude maximum of the reference signal in the (i + 1) th time window.
5. The method of claim 1, wherein the vibroseis force signal distortion characteristic is determined based on a time window average amplitude ratio, a number of samples having different sample symbol values between the two signals, and a number of samples of the vibroseis force signal according to the following formula:
Figure FDA0002904483950000023
wherein V represents the average amplitude ratio of the time window; n represents the number of sampling points with different sampling point symbol values between two signals; n represents the number of samples of the vibroseis force signal.
6. A vibroseis force signal distortion determination apparatus, comprising:
the reference signal selecting module is used for selecting discrete shot point data with the same number as the sampling points of the vibroseis force signal as a reference signal;
the sampling point symbol counting module is used for determining the number of sampling points with different sampling point symbol values between two signals according to the sampling point symbol values of the vibroseis force signal and the reference signal;
the amplitude maximum value determining module is used for averagely dividing the vibroseis force signal and the reference signal into a plurality of time windows, and determining the positive amplitude maximum value and the negative amplitude maximum value of the vibroseis force signal, and the positive amplitude maximum value and the negative amplitude maximum value of the reference signal in each time window;
the time window average amplitude ratio determining module is used for determining the time window average amplitude ratio according to the maximum positive amplitude value and the maximum negative amplitude value of the vibroseis force signals in the time windows and the maximum positive amplitude value and the maximum negative amplitude value of the reference signal;
and the distortion characteristic value determining module is used for determining the distortion characteristic value of the vibroseis force signal according to the time window average amplitude ratio, the number of sampling points with different sampling point symbol values between the two signals and the number of the sampling points of the vibroseis force signal.
7. The apparatus of claim 6, wherein the maximum amplitude determination module comprises:
the curve fitting unit is used for dividing the vibroseis force signal and the reference signal into a plurality of time windows, and respectively performing curve fitting of discrete sampling points on the vibroseis force signal and the reference signal in each time window to obtain vibration curves of the vibroseis force signal and the reference signal in each time window;
the maximum amplitude determination unit of the reference signal is used for selecting the amplitude value with the largest numerical value from the positive symbol value and the negative symbol value respectively according to the positive and negative of the symbol value of the amplitude on the vibration curve of the reference signal, and determining the maximum amplitude value and the maximum negative amplitude value of the vibroseis force signal in each time window;
and the force signal amplitude maximum value determining unit is used for selecting an amplitude value with the largest numerical value from the positive symbol value and the negative symbol value respectively according to the positive and negative of the symbol value of the amplitude on the vibration curve of the vibroseis force signal, and determining the amplitude value as the positive amplitude maximum value and the negative amplitude maximum value of the reference signal in each time window.
8. The apparatus of claim 6, wherein the time-window average magnitude ratio determination module is specifically configured to:
determining a positive amplitude ratio in a time window according to the maximum positive amplitude of the vibroseis force signal in the time window and the maximum positive amplitude of the reference signal;
determining a negative amplitude ratio in a time window according to the maximum value of the negative amplitude of the vibroseis force signal in the time window and the maximum value of the negative amplitude of the reference signal;
determining a positive-negative amplitude ratio difference value corresponding to the time window according to the positive amplitude ratio in the time window and the negative amplitude ratio in the time window;
determining the positive amplitude ratio difference value of the adjacent time window corresponding to the time window according to the positive amplitude ratio in the time window and the positive amplitude ratio of the next time window;
determining the negative amplitude ratio difference value of the adjacent time window corresponding to the time window according to the negative amplitude ratio in the time window and the negative amplitude ratio of the next time window;
and determining the average amplitude ratio of the time windows according to the positive and negative amplitude ratio difference values corresponding to the plurality of time windows, the positive amplitude ratio difference value of the adjacent time windows, the negative amplitude ratio difference value of the adjacent time windows and the number of the time windows.
9. The apparatus of claim 8, wherein the time-window average magnitude ratio determination module is specifically configured to:
determining the average amplitude ratio of the time windows according to the positive and negative amplitude ratio difference values corresponding to the time windows, the positive amplitude ratio difference value of the adjacent time windows, the negative amplitude ratio difference value of the adjacent time windows and the number of the time windows according to the following formula:
Figure FDA0002904483950000041
Figure FDA0002904483950000042
wherein V represents the average amplitude ratio of the time window; w represents the number of time windows; v i Representing the amplitude ratio of each time window; i represents the number of time windows; f + (max) i Represents the maximum value of the positive amplitude of the vibroseis force signal in the ith time window; r is + (max) i Represents the maximum value of the positive amplitude of the reference signal in the ith time window; f - (max) i Represents the maximum value of the negative amplitude of the vibroseis force signal in the ith time window; r - (max) i Represents the maximum value of the negative amplitude of the reference signal in the ith time window; f + (max) i+1 Represents the maximum value of the positive amplitude of the vibroseis force signal in the (i + 1) th time window; r + (max) i+1 The maximum value of the positive amplitude of the reference signal in the (i + 1) th time window; f - (max) i+1 The maximum value of the negative amplitude of the vibroseis force signal in the (i + 1) th time window; r - (max) i+1 The negative amplitude maximum of the reference signal in the (i + 1) th time window.
10. The apparatus of claim 6, wherein the distortion characteristic value determination module is specifically configured to:
determining the distortion characteristic value of the vibroseis force signal according to the average amplitude ratio of the time window, the number of sampling points with different sampling point symbol values between two signals and the number of the sampling points of the vibroseis force signal according to the following formula:
Figure FDA0002904483950000043
wherein V represents the average amplitude ratio of the time windows; n represents the number of sampling points with different sampling point symbol values between two signals; n represents the number of samples of the vibroseis force signal.
11. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method of any of claims 1 to 5 when executing the computer program.
12. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program for executing the method of any one of claims 1 to 5.
CN202110066751.3A 2021-01-19 2021-01-19 Method and device for determining distortion of vibroseis force signal Pending CN114814930A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110066751.3A CN114814930A (en) 2021-01-19 2021-01-19 Method and device for determining distortion of vibroseis force signal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110066751.3A CN114814930A (en) 2021-01-19 2021-01-19 Method and device for determining distortion of vibroseis force signal

Publications (1)

Publication Number Publication Date
CN114814930A true CN114814930A (en) 2022-07-29

Family

ID=82524725

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110066751.3A Pending CN114814930A (en) 2021-01-19 2021-01-19 Method and device for determining distortion of vibroseis force signal

Country Status (1)

Country Link
CN (1) CN114814930A (en)

Similar Documents

Publication Publication Date Title
US10705238B2 (en) Method and apparatus for processing seismic data
CN101241192B (en) Method for eliminating pneumatic gun near-field wavelet imaginary reaction
CN105204065B (en) A kind of method and apparatus of first break picking ripple
CN108693555A (en) Intelligent time-varying blind deconvolution wideband processing method and processing device
US20180156933A1 (en) Seismic acquisition geometry full-waveform inversion
CN110879412A (en) Underground transverse wave velocity inversion method, device, computing equipment and storage medium
CA2567993A1 (en) Method and apparatus of performing amplitude variation with offset analysis
CN113341455A (en) Viscous anisotropic medium seismic wave numerical simulation method, device and equipment
CN111736222A (en) Single-shot data signal-to-noise ratio determining method and device
CN114814930A (en) Method and device for determining distortion of vibroseis force signal
CN113568044B (en) Method and device for determining head arrival time of array acoustic logging
CN110688785B (en) Krauklis wave numerical simulation method and device based on plane wave seismic source
CN113900141B (en) Oil-gas distribution prediction method and device
CN113411140B (en) Channel self-adaptive channel characteristic learning method, simulator thereof and satellite navigation system
WO2021155754A1 (en) Method and apparatus for removing tube wave interference from optical fiber acoustic wave sensing seismic data
CN115270579A (en) Second-order acoustic wave equation finite difference numerical simulation parameter selection method
Zhang et al. First break of the seismic signals in oil exploration based on information theory
CN113239636A (en) Performance evaluation method, system, equipment and storage medium of weapon equipment
CN111352158B (en) Seismic signal enhancement method and device
CN109298451A (en) A kind of automatic Picking S wave seismic phase method for improving degree of skewness
CN106054245A (en) Earthquake data high-resolution processing quantitative evaluation method
CN112415601A (en) Method and device for determining surface quality factor Q value
CN110426736B (en) Offset gather acquisition method and device
CN114002742B (en) Euler Gaussian beam offset imaging method and device
CN112051613B (en) Method and device for generating aliasing acquisition dithering 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