阅读说明：本技术 地震信号的数据处理方法、装置、电子设备及存储介质 (Data processing method and device for seismic signals, electronic equipment and storage medium ) 是由 黄炜霖 贾靖 赵立松 卢钢 高菲 李洁丽 刘伟杰 刘科 于 2021-08-23 设计创作，主要内容包括：本申请公开了一种地震信号的数据处理方法、装置、电子设备及可读存储介质。其中,方法包括利用奇异谱分析方法对原始地震信号数据进行滤波,得到初始滤波数据；计算原始地震信号数据和初始滤波数据之间的差异性和相似性；基于差异性、相似性、预先设置且用于预测不稳定扰动的能量衰减的加权强度控制参数和用于确定信号保留程度的局部相似性阈值计算加权矩阵。基于加权矩阵对原始地震信号数据再次进行滤波处理,得到最终地震信号的处理结果,从而有效地提高了地震信号的信噪比,可以更加高效且更加稳定地压制地震噪声。(The application discloses a data processing method and device for seismic signals, electronic equipment and a readable storage medium. The method comprises the steps of filtering original seismic signal data by using a singular spectrum analysis method to obtain initial filtering data; calculating differences and similarities between the original seismic signal data and the initial filtering data; and calculating a weighting matrix based on the difference, the similarity, a preset weighting strength control parameter for predicting the energy attenuation of the unstable disturbance and a local similarity threshold value for determining the signal retention degree. And filtering the original seismic signal data again based on the weighting matrix to obtain a final seismic signal processing result, so that the signal-to-noise ratio of the seismic signal is effectively improved, and the seismic noise can be suppressed more efficiently and more stably.)

1. A method of data processing of seismic signals, comprising:

filtering original seismic signal data by using a singular spectrum analysis method to obtain initial filtering data;

determining differences and similarities between the original seismic signal data and the initial filtered data;

calculating a weighting matrix based on the difference, the similarity, a preset weighting strength control parameter and a local similarity threshold; the weighted intensity control parameter is used for predicting the energy attenuation of the unstable disturbance, and the local similarity threshold is used for determining the signal retention degree; the weighting matrix is used for predicting the position and the energy level of the unstable disturbance;

and filtering the original seismic signal data again based on the weighting matrix.

2. The method of data processing of seismic signals of claim 1, wherein said determining differences and similarities between said original seismic signal data and said initially filtered data comprises:

computing an absolute deviation U between the raw seismic signal data and the initial filtered data by calling an absolute deviation computation relation, U ═ U [ [ U ] ]_i,j]The absolute deviation calculation relation is as follows:u_i,jfor the element at the (i, j) position in the absolute deviation correspondence matrix, d_ijFor the element at the (i, j) position in the corresponding matrix of the raw seismic signal data,the initial filtered data is mapped to the element at the (i, j) position in the matrix.

3. The method of data processing of seismic signals of claim 2, wherein said determining differences and similarities between said original seismic signal data and said initially filtered data comprises:

computing local similarity between the original seismic signal data and the initial filtered data by calling a local similarity computation relation, V ═ V_i,j]The local similarity calculation relation is as follows:

in the formula, v_i,jFor the element at the (i, j) position in the local similarity correspondence matrix, ω is the local smoothing window function, τ is the current position (i, j) shifted by τ units along i, κ is the current position (i, j) shifted by κ units along j, d_τ,κFor an element at a (τ, κ) location in the corresponding matrix of the raw seismic signal data,the initial filtered data corresponds to the element at the (τ, κ) position in the matrix.

4. The method of data processing of seismic signals according to claim 1, wherein said calculating a weighting matrix based on said differences, said similarities, a preset weighting strength control parameter and a local similarity threshold comprises:

calling a weighting element calculation relational expression to calculate the weighting matrix, wherein the weighting element calculation relational expression is as follows:

in the formula, w_i,jFor the element of the weighting matrix at position (i, j), u_i,jFor the element of the absolute deviation correspondence matrix at position (i, j) for reflecting said discrepancy, ε is the normalized absolute deviation, v_i,jη is the local similarity threshold, and P is the weighted strength control parameter, for the elements of the local similarity correspondence matrix at position (i, j) for reflecting the similarity.

5. The method of processing seismic signal data according to any of claims 1 to 4, wherein said performing a filtering process again on said original seismic signal data based on said weighting matrix comprises:

calculating the Hadamard product of the weighting matrix to the original seismic signal data to obtain primary seismic signal modification data;

carrying out forward one-dimensional Fourier transform on the primary seismic signal modification data along the time direction to obtain a frequency slice;

sequentially performing Hankel matrix embedding operation, rank reduction operation and averaging operation on the frequency slices to obtain secondary seismic signal modification data;

and carrying out inverse one-dimensional Fourier transform on the secondary seismic signal modification data along the time direction to obtain a data processing result of the original seismic signal data.

6. A seismic signal data processing apparatus, comprising:

the initial filtering module is used for filtering the original seismic signal data by using a singular spectrum analysis method to obtain initial filtering data;

a parameter calculation module for determining differences and similarities between the original seismic signal data and the initial filtered data;

the weight calculation module is used for calculating a weighting matrix based on the difference, the similarity, a preset weighting strength control parameter and a local similarity threshold; the weighted intensity control parameter is used for predicting the energy attenuation of the unstable disturbance, and the local similarity threshold is used for determining the signal retention degree; the weighting matrix is used for predicting the position and the energy level of the unstable disturbance;

and the denoising module is used for filtering the original seismic signal data again based on the weighting matrix.

7. The seismic signal data processing apparatus of claim 6, wherein the parameter calculation module is further configured to:

computing an absolute deviation U between the raw seismic signal data and the initial filtered data by calling an absolute deviation computation relation, U ═ U [ [ U ] ]_i,j]The absolute deviation calculation relation can be expressed as:u_i,jfor the element at the (i, j) position in the absolute deviation correspondence matrix, d_ijFor the element at the (i, j) position in the corresponding matrix of the raw seismic signal data,mapping the (i, j) position in the matrix for the initial filtered dataThe elements of (1).

8. The seismic signal data processing apparatus of claim 6 or 7, wherein the denoising module is further configured to:

calculating the Hadamard product of the weighting matrix to the original seismic signal data to obtain primary seismic signal modification data;

carrying out forward one-dimensional Fourier transform on the primary seismic signal modification data along the time direction to obtain a frequency slice;

sequentially performing Hankel matrix embedding operation, rank reduction operation and averaging operation on the frequency slices to obtain secondary seismic signal modification data;

and carrying out inverse one-dimensional Fourier transform on the secondary seismic signal modification data along the time direction to obtain a data processing result of the original seismic signal data.

9. An electronic device comprising a processor and a memory, the processor being adapted to carry out the steps of the method of data processing of seismic signals according to any of claims 1 to 5 when executing a computer program stored in the memory.

10. A readable storage medium, characterized in that the readable storage medium has stored thereon a computer program which, when being executed by a processor, carries out the steps of the method for data processing of seismic signals according to any of claims 1 to 5.

Technical Field

The present application relates to the field of signal processing technologies, and in particular, to a method and an apparatus for processing seismic signal data, an electronic device, and a readable storage medium.

Background

Noise in seismic data brings errors to subsequent attribute analysis, wave impedance inversion, crack prediction and the like, and accordingly works such as underground structure interpretation, geophysical inversion and the like are seriously affected. Since the seismic data collected in the field are inevitably interfered by noise, the signal-to-noise ratio, i.e., the signal-to-noise ratio, in the signal system needs to be increased.

To improve the signal-to-noise ratio of the seismic data, incoherent noise in the acquired seismic signal data may be attenuated. The singular spectrum analysis method can detect the low-rank structure of the data, and can effectively inhibit the seismic random noise based on different expressions of the seismic signal and the random noise on the singular spectrum. However, since the singular spectrum analysis method adopts a quadratic fitting curve which is very sensitive to non-gaussian interference, when seismic data is interfered by unstable noise, the performance of the singular spectrum analysis method is unstable, abnormal values cannot be identified and suppressed based on a weighting criterion of absolute differences, a large amount of noise energy is easily left or considered to be interfered, the denoising effect of the seismic data is not obvious, the signal-to-noise ratio is not high, and if a better denoising effect is obtained in a mode of n (n is more than 2) iterations, the data processing efficiency is low due to the larger calculation cost required by the iteration step, especially when a function with robustness is difficult to approximate.

In view of this, how to more efficiently and more stably compress the seismic noise is a technical problem to be solved by those skilled in the art.

Disclosure of Invention

The application provides a data processing method and device for seismic signals, electronic equipment and a readable storage medium, which can more efficiently and more stably suppress seismic noise and effectively improve the signal-to-noise ratio of seismic signal data.

In order to solve the above technical problems, embodiments of the present invention provide the following technical solutions:

one aspect of the embodiments of the present invention provides a data processing method for seismic signals, including:

filtering original seismic signal data by using a singular spectrum analysis method to obtain initial filtering data;

determining differences and similarities between the original seismic signal data and the initial filtered data;

calculating a weighting matrix based on the difference, the similarity, a preset weighting strength control parameter and a local similarity threshold; the weighted intensity control parameter is used for predicting the energy attenuation of the unstable disturbance, and the local similarity threshold is used for determining the signal retention degree;

and filtering the original seismic signal data again based on the weighting matrix.

Optionally, the determining the differences and similarities between the original seismic signal data and the initial filtered data includes:

computing an absolute deviation U between the raw seismic signal data and the initial filtered data by calling an absolute deviation computation relation, U ═ U [ [ U ] ]_i,j]The absolute deviation calculation relation is as follows:u_i,jfor the element at the (i, j) position in the absolute deviation correspondence matrix, d_ijFor the element at the (i, j) position in the corresponding matrix of the raw seismic signal data,the initial filtered data is mapped to the element at the (i, j) position in the matrix.

Optionally, the determining the differences and similarities between the original seismic signal data and the initial filtered data includes:

computing local similarity between the original seismic signal data and the initial filtered data by calling a local similarity computation relationshipV, V ═ V_i,j]The local similarity calculation relation is as follows:

in the formula, v_i,jFor the element at the (i, j) position in the local similarity correspondence matrix, ω is the local smoothing window function, τ is the current position (i, j) shifted by τ units along i, κ is the current position (i, j) shifted by κ units along j, d_τ,κFor an element at a (τ, κ) location in the corresponding matrix of the raw seismic signal data,the initial filtered data corresponds to the element at the (τ, κ) position in the matrix.

Optionally, the calculating a weighting matrix based on the difference, the similarity, a preset weighting strength control parameter, and a local similarity threshold includes:

calling a weighting element calculation relational expression to calculate the weighting matrix, wherein the weighting element calculation relational expression is as follows:

in the formula, w_i,jFor the element of the weighting matrix at position (i, j), u_i,jFor the element of the absolute deviation correspondence matrix at position (i, j) for reflecting said discrepancy, ε is the normalized absolute deviation, v_i,jη is the local similarity threshold, and P is the weighted strength control parameter, for the elements of the local similarity correspondence matrix at position (i, j) for reflecting the similarity.

Optionally, the filtering the original seismic signal data based on the weighting matrix again includes:

calculating the Hadamard product of the weighting matrix to the original seismic signal data to obtain primary seismic signal modification data;

carrying out forward one-dimensional Fourier transform on the primary seismic signal modification data along the time direction to obtain a frequency slice;

sequentially performing Hankel matrix embedding operation, rank reduction operation and averaging operation on the frequency slices to obtain secondary seismic signal modification data;

and carrying out inverse one-dimensional Fourier transform on the secondary seismic signal modification data along the time direction to obtain a data processing result of the original seismic signal data.

Another aspect of the embodiments of the present invention provides a data processing apparatus for seismic signals, including:

the initial filtering module is used for filtering the original seismic signal data by using a singular spectrum analysis method to obtain initial filtering data;

a parameter calculation module for determining differences and similarities between the original seismic signal data and the initial filtered data;

the weight calculation module is used for calculating a weighting matrix based on the difference, the similarity, a preset weighting strength control parameter and a local similarity threshold; the weighted intensity control parameter is used for predicting the energy attenuation of the unstable disturbance, and the local similarity threshold is used for determining the signal retention degree;

and the denoising module is used for filtering the original seismic signal data again based on the weighting matrix.

Optionally, the parameter calculating module is further configured to: computing an absolute deviation U between the raw seismic signal data and the initial filtered data by calling an absolute deviation computation relation, U ═ U [ [ U ] ]_i,j]The absolute deviation calculation relation can be expressed as:u_i,jfor the element at the (i, j) position in the absolute deviation correspondence matrix, d_ijFor the element at the (i, j) position in the corresponding matrix of the raw seismic signal data,the initial filtered data is mapped to the element at the (i, j) position in the matrix.

An embodiment of the present invention further provides an electronic device, which includes a processor, and the processor is configured to implement the steps of the data processing method for seismic signals according to any one of the preceding claims when executing the computer program stored in the memory.

Finally, an embodiment of the present invention provides a readable storage medium, on which a computer program is stored, where the computer program, when executed by a processor, implements the steps of the data processing method for seismic signals according to any of the previous claims.

The technical scheme provided by the application has the advantages that the position and the energy level of unstable disturbance can be accurately predicted by utilizing the weighting matrix obtained by calculating the difference and the similarity between the primary filtering data and the original seismic signal data, the original seismic signal data are subjected to secondary filtering based on the weighting matrix, unstable and random noise in the seismic signal data can be eliminated, the signal-to-noise ratio of the seismic signal data is effectively improved, the whole process does not need to iterate for many times, large calculation cost is not needed, and the seismic noise can be suppressed more efficiently and more stably.

In addition, the embodiment of the invention also provides a corresponding implementation device, electronic equipment and a readable storage medium for the data processing method of the seismic signals, so that the method has higher practicability, and the device, the electronic equipment and the readable storage medium have corresponding advantages.

It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the disclosure.

Drawings

In order to more clearly illustrate the embodiments of the present invention or the technical solutions of the related art, the drawings required to be used in the description of the embodiments or the related art will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.

FIG. 1 is a schematic flow chart of a seismic signal data processing method according to an embodiment of the present invention;

FIG. 2 is a block diagram of an embodiment of a seismic signal data processing apparatus according to the present invention;

fig. 3 is a block diagram of an embodiment of an electronic device according to an embodiment of the present invention.

Detailed Description

In order that those skilled in the art will better understand the disclosure, the invention will be described in further detail with reference to the accompanying drawings and specific embodiments. It is to be understood that the described embodiments are merely exemplary of the invention, and not restrictive of the full scope of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.

The terms "first," "second," "third," "fourth," and the like in the description and claims of this application and in the above-described drawings are used for distinguishing between different objects and not for describing a particular order. Furthermore, the terms "comprising" and "having," as well as any variations thereof, are intended to cover non-exclusive inclusions. For example, a process, method, system, article, or apparatus that comprises a list of steps or elements is not limited to only those steps or elements but may include other steps or elements not expressly listed.

Having described the technical solutions of the embodiments of the present invention, various non-limiting embodiments of the present application are described in detail below.

Referring to fig. 1, fig. 1 is a schematic flow chart of a seismic signal data processing method according to an embodiment of the present invention, where the embodiment of the present invention includes the following steps:

s101: and filtering the original seismic signal data by using a singular spectrum analysis method to obtain initial filtering data.

In this step, the singular spectrum analysis method is any one of the existing methods for studying nonlinear time series data, and can construct a trajectory matrix according to the observed time series, and decompose and reconstruct the trajectory matrix, thereby extracting signals representing different components of the original time series, such as noise signals, periodic signals, and the like. The filtering process of the original address signal data based on the singular spectrum analysis method may include:

the singular spectrum analysis method comprises Hankel matrix embedding operation, rank reduction operation and anti-diagonal averaging operation. If D (x, ω) is a frequency slice of the original seismic signal data, where x ═ 1,2,3_x，ω＝1,2,3,...，N_xThe number of sampling points in the spatial and frequency directions. The singular spectrum analysis method firstly embeds the frequency slice D (x, omega) into the Hankel matrix to obtain a frequency slice with the following size (N)_xHankel matrix H of-m +1) × m^ωThe value of each element of the Hankel matrix is equal to the value of the element in the minor diagonal direction of the matrix. m is a predefined integer, and the value of m is such that a Hankel matrix H^ωApproximating a square matrix, e.g. m ═ N_x-N_xA/2 whereinRepresenting the integer part of the parameter. Hankel matrix H^ωCan be expressed as:

the singular spectrum analysis method assumes that the seismic data can be regarded as superposition of K plane waves, so the Hankel matrix H^ωIs K. The addition of noise will result in H^ωThe rank of (c) is increased. Therefore, the singular spectral analysis method reduces H by^ωRank to eliminate noise, i.e. Singular Value Decomposition (SVD):

wherein in ()^HRepresenting the Hermitian transpose of the matrix.Andrepresentation matrix H^ωThe first K maximum singular values and the associated K singular vectors.Represents H^ωA low rank approximation of. Finally, the singular spectrum analysis method is applied to the matrixIs based on the relationAveraging is performed to recover the filtered data. In the formula (I), the compound is shown in the specification,an averaging operator is indicated. The filtering method of the singular spectrum analysis method may be based onThe realization method is realized in the way that,respectively representing an averaging operator, a rank reduction operator and a Hankel operator.

S102: the differences and similarities between the raw seismic signal data and the initial filtered data are determined.

In this step, the difference refers to a deviation between the original seismic signal data and the initial filtered data, and the similarity refers to a degree of similarity between the original seismic signal data and the initial filtered data. Any parameter reflecting the deviation and similarity between data may be used for this purpose, and the present application is not limited thereto.

S103: and calculating a weighting matrix based on the difference, the similarity, a preset weighting strength control parameter and a local similarity threshold.

In order to accurately predict the position and energy level of the unstable disturbance, the present application may preset a plurality of parameters, including but not limited to a weighted intensity control parameter p that can be used to predict the energy attenuation of the unstable disturbance, a local similarity threshold η for determining the degree of signal retention, a constant K that depends on the number of plane waves included in the processing window, and a constant value λ for determining the intensity of the unstable disturbance. The denoising strength of the present embodiment increases as K, η, and λ decrease, and increases as p increases. The reference values for the parameters p, η and λ may be, for example, 5, 0.6 and 4, respectively. In the actual processing, the value can be floated up and down according to specific problems, which does not influence the realization of the application. According to the similarity and difference between the initial filtering data and the original seismic signals, a weighting matrix capable of predicting the position and the energy level of the unstable disturbance can be obtained by combining parameters determining the energy attenuation and the signal retention degree of the unstable disturbance.

S104: and filtering the original seismic signal data again based on the weighting matrix.

In the step, the original seismic signal data is processed again based on the weighting matrix calculated in the last step, the Hadamard product of the weighting matrix to the original seismic signal data can be calculated firstly, namely, the weighting matrix is multiplied by the corresponding position elements of the matrix corresponding to the original seismic signal data to obtain primary seismic signal modification data. And in the time direction, carrying out forward one-dimensional Fourier transform on the primary seismic signal modification data to obtain a frequency slice. And performing Hankel matrix embedding operation, rank reduction operation and averaging operation on the frequency slices in sequence according to the S101 method to obtain secondary seismic signal modification data. And in the time direction, performing inverse one-dimensional Fourier transform on the secondary seismic signal modification data to obtain final seismic signal data, and outputting the final seismic signal data serving as a data processing result of the original seismic signal data.

In the technical scheme provided by the embodiment of the invention, the position and the energy level of unstable disturbance can be accurately predicted by using the weighting matrix obtained by calculating the difference and the similarity between the primary filtering data and the original seismic signal data, the original seismic signal data is subjected to secondary filtering based on the weighting matrix, unstable and random noise in the seismic signal data can be eliminated, the signal-to-noise ratio of the seismic signal data is effectively improved, and the whole process does not need multiple iterations and large calculation cost, so that the seismic noise can be more efficiently and more stably suppressed.

It should be noted that, in the present application, there is no strict sequential execution order among the steps, and as long as the logical order is met, the steps may be executed simultaneously or according to a certain preset order, and fig. 1 is only an exemplary manner, and does not represent that only the execution order is the order.

In the above embodiment, how to perform steps S102 and S103 is not limited, and a calculation manner of the weighting matrix in this embodiment may include the following steps:

first, an absolute deviation calculation relationship may be used to calculate the absolute deviation U between the original seismic signal data and the original filtered data, U [ [ U ] ]_i,j]The absolute deviation calculation relationship can be expressed as:u_i,jfor the element at the (i, j) position in the absolute deviation correspondence matrix, d_ijFor the element at the (i, j) position in the corresponding matrix of the raw seismic signal data,the initial filtered data corresponds to the element at the (i, j) position in the matrix.

Second, a local similarity calculation relationship may be invoked to calculate a local similarity between the original seismic signal data and the initial filtered data, V ═ V_i,j]The local similarity calculation relationship can be expressed as:

in the formula, v_i,jFor the element at the (i, j) position in the local similarity correspondence matrix, ω is the local smoothing window function, τ is the current position (i, j) shifted by τ units along i, κ is the current position (i, j) shifted by κ units along j, d_τ,κFor an element at a (τ, κ) location in the corresponding matrix of the raw seismic signal data,the initial filtered data corresponds to the element at the (τ, κ) position in the matrix. ω may be selected from bell functions, trigonometric functions or rectangular functions. In this embodiment, the local similarity may also be calculated by using a shaping regularization method, wherein the smoothness may be determined by a shaping operator.

Finally, a weighting matrix representing the position and magnitude of the unstable disturbance may be calculated using the weighting element calculation relationship, with the grid points (i, j) of the raw seismic signal data at w_i,jContamination by unstable disturbances, < 1, w_i,jThe smaller the value of (c), the stronger the disturbance. The weighted element calculation relation can be expressed as:

in the formula, w_i,jIs the element of the weighting matrix at position (i, j), u_i,jFor the absolute deviation U used to reflect the difference to correspond to the element at the (i, j) position in the matrix, ε is the normalized absolute deviation, which may be calculated asv_i,jFor the element at the (i, j) position in the local similarity correspondence matrix for reflecting the similarity, η is the local similarity threshold, and P is the weighted intensity control parameter.

In the present embodiment, w is obtained by the final calculation_i,jThe position and the amplitude of the unstable disturbance can be predicted, and the position and the amplitude of the stable disturbance are obtained based on the predictionAnd the initial seismic signal data is filtered again, so that the noise can be removed to the maximum extent, and the signal-to-noise ratio of the initial seismic signal data is improved.

The embodiment of the invention also provides a corresponding device for the data processing method of the seismic signals, so that the method has higher practicability. Wherein the means can be described separately from the functional module point of view and the hardware point of view. The following describes a seismic signal data processing apparatus according to an embodiment of the present invention, and the seismic signal data processing apparatus described below and the seismic signal data processing method described above may be referred to in correspondence with each other.

Referring to fig. 2, fig. 2 is a block diagram of a data processing apparatus for seismic signals according to an embodiment of the present invention, in an embodiment, the apparatus may include:

an initial filtering module 201, configured to filter the original seismic signal data by using a singular spectrum analysis method to obtain initial filtering data;

a parameter calculation module 202 for determining differences and similarities between the original seismic signal data and the initial filtered data;

the weight calculation module 203 is configured to calculate a weighting matrix based on the difference, the similarity, a preset weighting strength control parameter, and a local similarity threshold; the weighted intensity control parameter is used for predicting the energy attenuation of the unstable disturbance, and the local similarity threshold is used for determining the signal retention degree;

and the denoising module 204 is configured to perform filtering processing on the original seismic signal data again based on the weighting matrix.

Optionally, in some embodiments of the present embodiment, the parameter calculating module 202 may include an absolute deviation calculating unit and a similarity calculating unit;

wherein the absolute deviation calculation unit is operable to: and calling an absolute deviation calculation relation to calculate the absolute deviation U of the original seismic signal data and the initial filtering data, wherein U is [ U ═ U_i,j]The absolute deviation calculation relationship can be expressed as:u_i,jfor the element at the (i, j) position in the absolute deviation correspondence matrix, d_ijFor the element at the (i, j) position in the corresponding matrix of the raw seismic signal data,the initial filtered data corresponds to the element at the (i, j) position in the matrix.

The similarity calculation unit may be operable to: local similarity calculation relation is called to calculate local similarity V between original seismic signal data and initial filtering data, and V is [ V ═ V [ [ V ]_i,j]The local similarity calculation relationship is:

in the formula, v_i,jFor the element at the (i, j) position in the local similarity correspondence matrix, ω is the local smoothing window function, τ is the current position (i, j) shifted by τ units along i, κ is the current position (i, j) shifted by κ units along j, d_τ,κFor an element at a (τ, κ) location in the corresponding matrix of the raw seismic signal data,the initial filtered data corresponds to the element at the (τ, κ) position in the matrix.

Optionally, in other embodiments of this embodiment, the weight calculating module 203 may be further configured to: calling a weighting element to calculate a relational calculation weighting matrix, wherein the weighted element calculation relational expression is as follows:

in the formula, w_i,jIs the element of the weighting matrix at position (i, j), u_i,jIs the absolute deviation, ε is the normalized absolute deviation, v_i,jFor local similarity, η is the local similarity threshold, and P is the weighted intensity controlAnd (4) parameters.

Optionally, in some other embodiments of this embodiment, the denoising module 204 may be further configured to:

calculating the Hadamard product of the weighting matrix to the original seismic signal data to obtain primary seismic signal modification data;

in the time direction, carrying out forward one-dimensional Fourier transform on the primary seismic signal modification data to obtain a frequency slice;

respectively performing Hankel matrix embedding, rank reduction and averaging operation on the frequency slices to obtain secondary seismic signal modification data;

and performing inverse one-dimensional Fourier transform on the secondary seismic signal modification data along the time direction to obtain final seismic signal data.

The functions of the functional modules of the data processing device of seismic signals in the embodiment of the present invention may be specifically implemented according to the method in the embodiment of the method, and the specific implementation process may refer to the related description of the embodiment of the method, which is not described herein again.

Therefore, the embodiment of the invention effectively improves the signal-to-noise ratio of the seismic signal, and can more efficiently and more stably compress the seismic noise.

The above mentioned data processing device for seismic signals is described from the perspective of functional modules, and further, the present application also provides an electronic device described from the perspective of hardware. Fig. 3 is a schematic structural diagram of an electronic device according to an embodiment of the present application. As shown in fig. 3, the electronic device comprises a memory 30 for storing a computer program; a processor 31 for implementing the steps of the method for data processing of seismic signals as mentioned in any of the embodiments above when executing a computer program.

The processor 31 may include one or more processing cores, such as a 4-core processor, an 8-core processor, and the processor 31 may also be a controller, a microcontroller, a microprocessor or other data processing chip, and the like. The processor 31 may be implemented in at least one hardware form of a DSP (Digital Signal Processing), an FPGA (Field-Programmable Gate Array), and a PLA (Programmable Logic Array). The processor 31 may also include a main processor and a coprocessor, where the main processor is a processor for Processing data in an awake state, and is also called a Central Processing Unit (CPU); a coprocessor is a low power processor for processing data in a standby state. In some embodiments, the processor 31 may be integrated with a GPU (Graphics Processing Unit), which is responsible for rendering and drawing the content required to be displayed by the display screen. In some embodiments, the processor 31 may further include an AI (Artificial Intelligence) processor for processing a calculation operation related to machine learning.

Memory 30 may include one or more computer-readable storage media, which may be non-transitory. Memory 30 may also include high speed random access memory as well as non-volatile memory, such as one or more magnetic disk storage devices, flash memory storage devices. The memory 30 may in some embodiments be an internal storage unit of the electronic device, for example a hard disk of a server. The memory 30 may also be an external storage device of the electronic device in other embodiments, such as a plug-in hard disk provided on a server, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card), and the like. Further, the memory 30 may also include both an internal storage unit and an external storage device of the electronic device. The memory 30 can be used for storing various data and application software installed in the electronic device, such as: the code of the program that executes the vulnerability handling method, etc. may also be used to temporarily store data that has been output or is to be output. In this embodiment, the memory 30 is at least used for storing the following computer program 301, wherein after being loaded and executed by the processor 31, the computer program can realize the relevant steps of the data processing method of seismic signals disclosed in any one of the foregoing embodiments. In addition, the resources stored by the memory 30 may also include an operating system 302, data 303, and the like, and the storage may be transient storage or permanent storage. Operating system 302 may include Windows, Unix, Linux, etc. Data 303 may include, but is not limited to, data corresponding to data processing results of seismic signals, and the like.

In some embodiments, the electronic device may further include a display 32, an input/output interface 33, a communication interface 34 or network interface, a power source 35, and a communication bus 36. The display 32 and the input/output interface 33, such as a Keyboard (Keyboard), belong to a user interface, and the optional user interface may also include a standard wired interface, a wireless interface, and the like. Alternatively, in some embodiments, the display may be an LED display, a liquid crystal display, a touch-sensitive liquid crystal display, an OLED (Organic Light-Emitting Diode) touch device, or the like. The display, which may also be referred to as a display screen or display unit, as appropriate, is used for displaying information processed in the electronic device and for displaying a visualized user interface. The communication interface 34 may optionally include a wired interface and/or a wireless interface, such as a WI-FI interface, a bluetooth interface, etc., typically used to establish a communication connection between an electronic device and other electronic devices. The communication bus 36 may be a Peripheral Component Interconnect (PCI) bus, an Extended Industry Standard Architecture (EISA) bus, or the like. The bus may be divided into an address bus, a data bus, a control bus, etc. For ease of illustration, only one thick line is shown in FIG. 3, but this does not mean only one bus or one type of bus.

Those skilled in the art will appreciate that the configuration shown in fig. 3 is not intended to be limiting of the electronic device and may include more or fewer components than those shown, such as a sensor 37 to perform various functions.

The functions of the functional modules of the electronic device according to the embodiments of the present invention may be specifically implemented according to the method in the above method embodiments, and the specific implementation process may refer to the description related to the above method embodiments, which is not described herein again.

Therefore, the embodiment of the invention effectively improves the signal-to-noise ratio of the seismic signal, and can more efficiently and more stably compress the seismic noise.

It is to be understood that, if the data processing method of the seismic signal in the above-described embodiment is implemented in the form of a software functional unit and sold or used as a separate product, it may be stored in a computer-readable storage medium. Based on such understanding, the technical solutions of the present application may be substantially or partially implemented in the form of a software product, which is stored in a storage medium and executes all or part of the steps of the methods of the embodiments of the present application, or all or part of the technical solutions. And the aforementioned storage medium includes: a U disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), an electrically erasable programmable ROM, a register, a hard disk, a multimedia card, a card type Memory (e.g., SD or DX Memory, etc.), a magnetic Memory, a removable magnetic disk, a CD-ROM, a magnetic or optical disk, and other various media capable of storing program codes.

Based on this, the embodiment of the present invention further provides a readable storage medium, which stores a computer program, and the computer program is executed by a processor, and the steps of the data processing method of the seismic signal according to any one of the above embodiments are provided.

The functions of the functional modules of the readable storage medium according to the embodiment of the present invention may be specifically implemented according to the method in the foregoing method embodiment, and the specific implementation process may refer to the description related to the foregoing method embodiment, which is not described herein again.

The embodiments are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same or similar parts among the embodiments are referred to each other. For hardware including devices and electronic equipment disclosed by the embodiment, the description is relatively simple because the hardware includes the devices and the electronic equipment correspond to the method disclosed by the embodiment, and the relevant points can be obtained by referring to the description of the method.

Those of skill would further appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, computer software, or combinations of both, and that the various illustrative components and steps have been described above generally in terms of their functionality in order to clearly illustrate this interchangeability of hardware and software. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the implementation. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.

The present application provides a method, an apparatus, an electronic device and a readable storage medium for processing seismic signal data. The principles and embodiments of the present invention are explained herein using specific examples, which are presented only to assist in understanding the method and its core concepts. It should be noted that, for those skilled in the art, it is possible to make various improvements and modifications to the present invention without departing from the principle of the present invention, and those improvements and modifications also fall within the scope of the claims of the present application.