Physical data signal processing method and apparatus, computer device, and storage medium

Abstract

This invention provides a physical data signal processing method, apparatus, computer equipment, and storage medium. The method is based on a bispectral operation-based signal reconstruction approach. It calculates the harmonic coefficients, phase angle, and number of harmonic terms of the original signal's harmonic expression, as well as the harmonic coefficients, phase angle, and number of harmonic terms of the harmonic coefficients of the original signal's harmonic expression. Then, it reconstructs the harmonic coefficients of the original signal's harmonic expression and the original signal. During the reconstruction of the harmonic coefficients of the original signal's harmonic expression, the sample spacing of the harmonic coefficients is amplified, effectively increasing the dominant frequency of the physical data signal. This invention is applied to seismic exploration, effectively increasing the dominant frequency of the seismic wavelet signal, thus creating favorable conditions for improving the resolution of seismic exploration.

Description

Technical Field

[0001] This invention relates to the field of applied geophysics, and in particular to a physical data signal processing method, apparatus, computer equipment, and storage medium. Background Technology

[0002] Seismic exploration data, as one of the geophysical exploration signals, is a comprehensive reflection of various underground geological structures and the composition of various rocks and minerals within strata. To more accurately identify and describe thin-layer structures and small-sized rock masses and ore bodies, it is necessary to improve the resolution of seismic exploration data, a goal that geophysicists have been striving for. Currently, commonly used methods to improve the resolution of seismic exploration data include digital signal filtering, frequency decomposition, and deconvolution. However, these commonly used methods still cannot significantly improve the resolution of geophysical exploration data, nor can they more accurately identify and describe underground geological structures and rock properties and lithology. Summary of the Invention

[0003] Therefore, it is necessary to provide a physical data signal processing method, apparatus, computer equipment, and storage medium to address the aforementioned technical problems.

[0004] A physical data signal processing method, comprising:

[0005] Acquire physical data;

[0006] The physical data is preprocessed and bispectral calculations are performed to obtain the first harmonic signal expression of the physical data;

[0007] Calculate the first harmonic coefficient, the first phase angle, and the number of first harmonic terms in the expression of the first harmonic signal;

[0008] The first harmonic coefficient is preprocessed and bispectral calculated to obtain the second harmonic signal expression of the first harmonic coefficient;

[0009] Calculate the second harmonic coefficient, second phase angle, and second harmonic term number in the expression of the second harmonic signal;

[0010] The second harmonic coefficient, the second phase angle, the second harmonic term number, and the sample spacing of the first harmonic coefficient are amplified based on a preset amplification factor, and the amplified result is substituted into the expression of the second harmonic signal to obtain the third harmonic coefficient.

[0011] Substituting the third harmonic coefficient, the first phase angle, and the first harmonic term number into the first harmonic signal expression, the reconstructed signal data is obtained.

[0012] In one embodiment, the step of substituting the third harmonic coefficient, the first phase angle, and the first harmonic term number into the first harmonic signal expression to obtain the reconstructed physical data includes:

[0013] Substitute the third harmonic coefficient, the first phase angle, and the first harmonic term number into the first harmonic signal expression to calculate the reconstructed signal data.

[0014] The reconstructed signal data is subjected to inverse preprocessing to obtain the reconstructed physical data.

[0015] In one embodiment, the step of preprocessing the first harmonic coefficient and performing bispectral calculation to obtain the bispectral data of the first harmonic coefficient includes:

[0016] The first harmonic coefficient is preprocessed;

[0017] Perform a Fourier transform on the preprocessed first harmonic coefficient to calculate the spectrum of the first harmonic coefficient;

[0018] Based on the spectrum of the first harmonic coefficient, the bispectrum of the first harmonic coefficient in the first quadrant is calculated;

[0019] The fundamental frequency of the bispectral spectrum of the first harmonic coefficient is estimated from the bispectral spectrum of the first harmonic coefficient in the first quadrant.

[0020] Based on the fundamental frequency of the bispectral spectrum of the first harmonic coefficient, the first harmonic coefficient is phase-shifted to obtain the phase-shifted first harmonic coefficient.

[0021] Based on the phase-shifted first harmonic coefficient, the bispectral data of the phase-shifted first harmonic coefficient in the first quadrant is calculated.

[0022] Based on the bispectrum of the first harmonic coefficient in the first quadrant, the expression for the first harmonic signal of the bispectrum of the first harmonic coefficient is calculated.

[0023] In one embodiment, the step of preprocessing the physical data and performing bispectral calculations to obtain the first harmonic signal expression of the physical data includes:

[0024] The physical data is preprocessed and bispectral calculations are performed to obtain the bispectral data of the physical data;

[0025] Based on the bispectral data of the physical data, the expression for the first harmonic signal of the bispectral data of the physical data is calculated.

[0026] In one embodiment, the step of preprocessing the physical data and performing bispectral calculation to obtain the bispectral data includes:

[0027] The physical data is preprocessed;

[0028] The Fourier transform is performed on the preprocessed physical data to calculate the spectrum of the physical data;

[0029] Based on the spectrum of the physical data, the bispectrum of the physical data in the first quadrant is calculated;

[0030] The fundamental frequency of the bispectral spectrum of the physical data is estimated from the bispectral spectrum in the first quadrant of the physical data;

[0031] Based on the fundamental frequency, the physical data is phase-shifted to obtain the phase-shifted physical data;

[0032] Based on the phase-shifted physical data, the bispectral data in the first quadrant is calculated.

[0033] In one embodiment, the step of calculating the bispectrum of the physical data in the first quadrant based on the spectrum of the physical data includes:

[0034] Based on the spectrum of the physical data, the bispectrum of the physical data is calculated;

[0035] Based on the bispectral data of the physical data, the bispectral data in the first quadrant is calculated using symmetry.

[0036] In one embodiment, the step of preprocessing the physical data includes:

[0037] The physical data is subjected to delinear background processing, zero-mean processing, and period expansion processing.

[0038] A physical data signal processing device, comprising:

[0039] The physical data acquisition module is used to acquire physical data;

[0040] The first expression calculation module is used to preprocess the physical data and perform bispectral calculation to obtain the first harmonic signal expression of the physical data.

[0041] The first harmonic calculation module is used to calculate the first harmonic coefficient, the first phase angle, and the first harmonic term number of the first harmonic signal expression.

[0042] The second expression calculation module is used to preprocess the first harmonic coefficient and perform bispectral calculation to obtain the second harmonic signal expression of the first harmonic coefficient.

[0043] The second harmonic calculation module is used to calculate the second harmonic coefficient, the second phase angle, and the number of second harmonic terms in the expression of the second harmonic signal.

[0044] The amplification module is used to amplify the second harmonic coefficient, the second phase angle, the second harmonic term number and the sample spacing of the first harmonic coefficient based on a preset amplification factor, and substitute the amplified result into the expression of the second harmonic signal to obtain the third harmonic coefficient;

[0045] The reconstructed signal data acquisition module is used to substitute the third harmonic coefficient, the first phase angle, and the first harmonic term number into the first harmonic signal expression to calculate the reconstructed signal data.

[0046] A computer device includes a memory and a processor, the memory storing a computer program, characterized in that the processor executes the computer program to perform the following steps:

[0047] Acquire physical data;

[0048] The physical data is preprocessed and bispectral calculations are performed to obtain the first harmonic signal expression of the physical data;

[0049] Calculate the first harmonic coefficient, the first phase angle, and the number of first harmonic terms in the expression of the first harmonic signal;

[0050] The first harmonic coefficient is preprocessed and bispectral calculated to obtain the second harmonic signal expression of the first harmonic coefficient;

[0051] Calculate the second harmonic coefficient, second phase angle, and second harmonic term number in the expression of the second harmonic signal;

[0052] The second harmonic coefficient, the second phase angle, the second harmonic term number, and the sample spacing of the first harmonic coefficient are amplified based on a preset amplification factor, and the amplified result is substituted into the expression of the second harmonic signal to obtain the third harmonic coefficient.

[0053] Substituting the third harmonic coefficient, the first phase angle, and the first harmonic term number into the first harmonic signal expression, the reconstructed signal data is obtained.

[0054] A computer-readable storage medium having a computer program stored thereon, the computer program performing the following steps when executed by a processor:

[0055] Acquire physical data;

[0056] The physical data is preprocessed and bispectral calculations are performed to obtain the first harmonic signal expression of the physical data;

[0057] Calculate the first harmonic coefficient, the first phase angle, and the number of first harmonic terms in the expression of the first harmonic signal;

[0058] The first harmonic coefficient is preprocessed and bispectral calculated to obtain the second harmonic signal expression of the first harmonic coefficient;

[0059] Calculate the second harmonic coefficient, second phase angle, and second harmonic term number in the expression of the second harmonic signal;

[0060] The second harmonic coefficient, the second phase angle, the second harmonic term number, and the sample spacing of the first harmonic coefficient are amplified based on a preset amplification factor, and the amplified result is substituted into the expression of the second harmonic signal to obtain the third harmonic coefficient.

[0061] Substituting the third harmonic coefficient, the first phase angle, and the first harmonic term number into the first harmonic signal expression, the reconstructed signal data is obtained.

[0062] The aforementioned physical data signal processing method, apparatus, computer equipment, and storage medium, based on a bispectral operation signal reconstruction method, calculate the harmonic coefficients, phase angle, and number of harmonic terms of the original signal harmonic expression, as well as the harmonic coefficients, phase angle, and number of harmonic terms of the harmonic coefficients of the original signal harmonic expression. Then, the harmonic coefficients of the original signal harmonic expression and the original signal are reconstructed. During the reconstruction of the harmonic coefficients of the original signal harmonic expression, the sample spacing of the harmonic coefficients is amplified, effectively increasing the dominant frequency of the seismic wavelet signal, thus creating favorable conditions for improving the resolution of seismic exploration. Attached Figure Description

[0063] Figure 1 This is a flowchart illustrating a physical data signal processing method in one embodiment;

[0064] Figure 2 This is a schematic diagram of the structure of a computer device in one embodiment;

[0065] Figure 3 This is a flowchart illustrating the physical data signal processing method in another embodiment;

[0066] Figure 4 This is a schematic diagram comparing the signal of theoretical seismic wavelet data with the reconstructed signal after the seismic wavelet dominant frequency is boosted by 1.5 times in one embodiment;

[0067] Figure 5 This is a schematic diagram comparing the signal of theoretical seismic wavelet data with the reconstructed signal after the seismic wavelet dominant frequency is boosted by 9 times in one embodiment;

[0068] Figure 6A This is a schematic diagram of seismic fluctuations from the original seismic data.

[0069] Figure 6B for Figure 6A The diagram shows the seismic fluctuations in the reconstructed seismic data after the main frequency of the original seismic data was increased by 1.5 times through reconstruction. Detailed Implementation

[0070] To make the objectives, technical solutions, and advantages of this application clearer, the following detailed description is provided in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the scope of this application.

[0071] Example 1

[0072] In this embodiment, as Figure 1 As shown, a physical data signal processing method is provided, which includes:

[0073] Step 110: Obtain physical data.

[0074] Specifically, the physical data is geophysical data. In this embodiment, the data file of geophysical data x(n) is read to obtain the physical data x(n).

[0075] Step 120: Preprocess the physical data and perform bispectral calculations to obtain the first harmonic signal expression of the physical data.

[0076] In this step, the physical data x(n) is preprocessed and bispectral calculated sequentially to obtain the harmonic signal expression of the physical data, namely the first harmonic signal expression x(t).

[0077] In one embodiment, preprocessing includes delinearizing the background, zero-meaning, and periodic expansion of the physical data. Bispectral calculation involves calculating the spectrum of the physical data using FFT, and then calculating the bispectral data using a triple correlation method.

[0078] Step 130: Calculate the first harmonic coefficient, first phase angle, and first harmonic term number of the first harmonic signal expression.

[0079] In this step, the harmonic coefficient a of the first harmonic signal expression x(t) is obtained. i With phase angle The number of harmonic terms M, and the harmonic coefficient a i That is, the first harmonic coefficient, phase angle That is, the first phase angle, and the number of harmonic terms M is the number of the first harmonic terms.

[0080] Step 140: Preprocess the first harmonic coefficient and perform bispectral calculation to obtain the second harmonic signal expression of the first harmonic coefficient.

[0081] In this step, the first harmonic coefficient a is processed in the same way as in step 120. i Preprocessing and bispectral calculations were performed to obtain the first harmonic coefficient a. i The expression for the second harmonic signal.

[0082] In one embodiment, preprocessing includes delinearization, zero-mean normalization, and periodic expansion. Bispectral calculation involves calculating the spectrum of the first harmonic coefficient using FFT, and then calculating the bispectral density of the first harmonic coefficient using triple correlation.

[0083] Step 150: Calculate the second harmonic coefficient, second phase angle, and second harmonic term number of the second harmonic signal expression.

[0084] In this step, the harmonic coefficient aa in the second harmonic expression is obtained. i Phase angle With the number of harmonic terms M′, the harmonic coefficient aa i That is, the second harmonic coefficient, phase angle This is the second phase angle, and the number of harmonic terms M′ is the number of the second harmonic terms.

[0085] Step 160: Amplify the second harmonic coefficient, the second phase angle, the second harmonic term number, and the sample spacing of the first harmonic coefficient based on a preset amplification factor, and substitute the amplified result into the expression of the second harmonic signal to obtain the third harmonic coefficient.

[0086] In this embodiment, the preset amplification factor is coefficient α, and the preset amplification factor is greater than or equal to 1, so that the second harmonic coefficient aa i The second phase angle The second harmonic term M′ and the first harmonic coefficient a i Multiplying the sample spacing by a preset amplification factor α yields the third harmonic coefficient a′. i The third harmonic coefficient a′ i It can also be called the reconstructed harmonic coefficient.

[0087] Step 170: Substitute the third harmonic coefficient, the first phase angle, and the first harmonic term number into the first harmonic signal expression to calculate the reconstructed signal data.

[0088] In this step, the third harmonic coefficient a′ i The first phase angle Substituting the first harmonic term number M into the first harmonic signal expression, the reconstructed signal data x′(t) is obtained.

[0089] In this embodiment, the main frequency of the reconstructed signal x′(t) is increased by α (α≥1) times compared to the main frequency of x(t), thereby realizing the increase of the main frequency of the geophysical data signal. The magnitude of the increase is controlled by the input parameter α (α≥1).

[0090] In the above embodiments, the signal reconstruction method based on bispectral operation calculates the harmonic coefficients, phase angle, and number of harmonic terms of the original signal harmonic expression, as well as the harmonic coefficients, phase angle, and number of harmonic terms of the harmonic coefficients of the original signal harmonic expression. Then, the harmonic coefficients of the original signal harmonic expression and the original signal are reconstructed respectively. When reconstructing the harmonic coefficients of the original signal harmonic expression, the sample spacing of the harmonic coefficients of the original signal harmonic expression is enlarged, which effectively improves the main frequency of the seismic wavelet signal, thereby creating favorable conditions for improving the resolution of seismic exploration.

[0091] In one embodiment, the step of substituting the third harmonic coefficient, the first phase angle, and the first harmonic term number into the first harmonic signal expression to obtain the reconstructed physical data includes: substituting the third harmonic coefficient, the first phase angle, and the first harmonic term number into the first harmonic signal expression to calculate and obtain reconstructed signal data; and performing inverse preprocessing on the reconstructed signal data to obtain reconstructed physical data.

[0092] In this embodiment, the inverse preprocessing includes linear background compensation and mean compensation processing on the reconstructed signal data to obtain the reconstructed physical data and output the calculation results.

[0093] In one embodiment, the step of preprocessing the physical data and calculating the bispectral density to obtain the bispectral density of the physical data includes: preprocessing the physical data; performing a Fourier transform on the preprocessed physical data to calculate the spectrum of the physical data; calculating the bispectral density of the physical data in the first quadrant based on the spectrum of the physical data; estimating the fundamental frequency of the bispectral density of the physical data from the bispectral density of the physical data in the first quadrant; performing a phase shift on the physical data according to the fundamental frequency to obtain the phase-shifted physical data; and calculating the bispectral density of the phase-shifted physical data in the first quadrant based on the phase-shifted physical data.

[0094] In one embodiment, the step of calculating the bispectrum of the physical data in the first quadrant based on the spectrum of the physical data includes: calculating the bispectrum of the physical data based on the spectrum of the physical data; and calculating the bispectrum of the physical data in the first quadrant based on the bispectrum of the physical data using symmetry.

[0095] In the above embodiments, the process of preprocessing the physical data and performing bispectral calculations is as follows:

[0096] Step 1: Read the data file containing geophysical data x(n).

[0097] Step 2: Perform preprocessing on the geophysical data x(n) including delinearization, zero-meaning, and periodic expansion.

[0098] Step 3: Calculate the spectrum X(f) of the data x(n) using Fourier transform.

[0099] Step 4: Calculate the bispectral B(f1,f2) of the data x(n) through triple correlation.

[0100] Step 5: Calculate the bispectral B(f1,f2) of the data x(n) in the entire first quadrant using symmetry.

[0101] Step 6: Take the first peak value of |B(f,f)| as the estimated fundamental frequency f0.

[0102] Step 7: Add a half-fundamental frequency cosine signal with a phase of 0 to the original data signal x(n), i.e., y(n) = x(n) + cos(πf0n). In this step, the physical data is phase-shifted according to the fundamental frequency f0 to obtain the phase-shifted physical data y(n) = x(n) + cos(πf0n).

[0103] Step 8: Calculate the spectrum Y(f) of the data signal y(n) using Fourier transform.

[0104] Step nine: Calculate the bispectral B of the data signal y(n) through triple correlation. y (f1,f2).

[0105] Step 10: Calculate the bispectral B of the data signal y(n) across the entire first quadrant using symmetry. y (f1,f2).

[0106] In one embodiment, the step of preprocessing the first harmonic coefficient and calculating its bispectral density to obtain the bispectral density of the first harmonic coefficient includes: preprocessing the first harmonic coefficient; performing a Fourier transform on the preprocessed first harmonic coefficient to calculate its spectrum; calculating the bispectral density of the first harmonic coefficient in the first quadrant based on its spectrum; estimating the fundamental frequency of the bispectral density of the first harmonic coefficient from its bispectral density in the first quadrant; performing a phase shift on the first harmonic coefficient based on its fundamental frequency to obtain the phase-shifted first harmonic coefficient; calculating the bispectral density of the phase-shifted first harmonic coefficient in the first quadrant based on its bispectral density; and calculating the first harmonic signal expression of the bispectral density of the first harmonic coefficient based on its bispectral density in the first quadrant.

[0107] In one embodiment, the step of preprocessing the physical data and performing bispectral calculation to obtain the first harmonic signal expression of the physical data includes: preprocessing the physical data and performing bispectral calculation to obtain the bispectrum of the physical data; and calculating the first harmonic signal expression of the bispectrum of the physical data based on the bispectrum of the physical data.

[0108] In the above embodiments, the process of preprocessing the first harmonic data and performing bispectral calculations can be referred to the process of preprocessing the physical data and performing bispectral calculations. The preprocessing and bispectral calculations of the first harmonic coefficients are performed using the same method as the preprocessing and bispectral calculations of the physical data, employing steps two to ten of the above embodiments. This embodiment will not elaborate on these steps.

[0109] Example 2

[0110] (1) Principle of physical data signal reconstruction method based on bispectral operation:

[0111] Let x(t) be a continuous real signal with a mean of 0, and its Fourier transform be:

[0112] X(f)=A(f)·exp[jφ(f)] (1)

[0113] A(f) and φ(f) are the amplitude and phase of the spectrum X(f), respectively.

[0114] The expression for the third-order cumulant function of x(t) is:

[0115]

[0116] t1 and t2 represent two different time variables.

[0117] The two-dimensional Fourier transform of C(t1,t2) is the bispectral B(f1,f2).

[0118]

[0119] X * (f1+f2) is the conjugate of the spectrum X(f1+f2) (Note: the same applies below). Where f1 and f2 are the frequency variables corresponding to the time variables t1 and t2.

[0120] Furthermore, the bispectral expression for x(t) can be obtained:

[0121]

[0122] Let the seismic exploration data x(n) be the discrete form of the continuous real signal x(t), which can be described by the expression for harmonic signals, i.e.

[0123]

[0124] In the formula: M is the number of harmonic terms, a i and These are the harmonic coefficients and phase angles, respectively, and f0 is the fundamental frequency.

[0125] Based on the properties of the Fourier transform, the relationships between the bispectral amplitude spectrum, phase spectrum, harmonic coefficients, and phase angle of x(n) can be obtained from equations (4) and (5):

[0126]

[0127]

[0128] The unknowns in equations (6a) and (6b) are also the unknowns in equation (5), including the fundamental frequency f0 and the harmonic coefficient a. i With phase angle Obviously, the number of unknowns in (6a) and (6b) is greater than the number of equations. In practical applications, the fundamental frequency f0 must be accurately estimated first. This is the key to the successful implementation of this signal reconstruction method. Moreover, as can be seen from equation (5), the accuracy of the fundamental frequency f0 estimation determines the fidelity of the reconstructed signal.

[0129] Based on the estimated fundamental frequency f0, a half-fundamental frequency cosine signal z(t) = cos(πf0t) with zero phase is added to the signal x(t), i.e.

[0130] y(t)=x(t)+z(t)=x(t)+cos(πf0t) (7)

[0131] Since the Fourier transform of y(t) can be expressed as:

[0132] Y(f)=X(f)+Z(f) (8)

[0133] Therefore, the bispectral representation of y(t) can be expressed as:

[0134]

[0135] For equation (9), if f1 = f2 = f0 / 2, then according to the properties of the Fourier transform, we have

[0136]

[0137] Therefore,

[0138]

[0139]

[0140] If neither f1 nor f2 is f0 / 2, and f1 + f2 ≠ f0 / 2, then we have

[0141] Z(f1)=Z(f2)=Z * (f1+f2)=0 (12)

[0142] Therefore, from equations (12) and (9), we can obtain

[0143] B y (f1,f2)=X(f1)·X(f2)·X * (f1+f2)=B x (f1,f2) (13)

[0144] That is, the bispectral density of signal y(t) is now the same as the bispectral density of signal x(t). Therefore, the harmonic coefficients a in the desired harmonic signal expression x(n) can be obtained. i With phase angle The recurrence relation for solving the problem is:

[0145]

[0146]

[0147] The harmonic coefficient a obtained from equation (14) i and phase angle Substituting into equation (5), the signal x(n) can be reconstructed.

[0148] Before performing bispectral operations on the original data signal, three preprocessing steps are performed on the original signal: delinearization, zero-mean normalization, and period expansion (see the invention patent "An Optimization Method for Signal Recovery", patent number: 201010508001.9).

[0149] The core content of this invention:

[0150] The bispectral density of the physical data signal after the above preprocessing is calculated, and then the harmonic coefficients a of its harmonic signal expression x(t) are obtained according to equation (5). i With phase angle

[0151] For harmonic coefficient a i Preprocessing and bispectral calculation are performed, and then a is obtained according to equation (5). i The harmonic coefficients aa in the expression of the harmonic signal i With phase angle

[0152] aa i , and a i Multiply the sample spacing by a coefficient α (α≥1) and substitute it into a i The expression for the harmonic signal is a′. i ;

[0153] a′ i , Substituting the harmonic signal expression of x(t) into the expression, we obtain the new physical data harmonic signal expression x′(t);

[0154] The main frequency of signal x′(t) is increased by α (α≥1) times compared to the main frequency of x(t). Thus, the main frequency of the physical data signal is increased, and the magnitude of the increase is controlled by the input parameter α (α≥1).

[0155] This application is applied to the reconstruction of seismic wavelets in geophysical exploration, effectively realizing the enhancement of the dominant frequency of the seismic wavelet signal. The enhancement magnitude can be controlled by the input parameter α (α≥1), thereby creating favorable conditions for improving the resolution of seismic exploration.

[0156] This invention patent application effectively increases the main frequency of physical data digital signals, providing a very useful tool for improving the resolution and precision of complex signal interpretation.

[0157] The specific calculation process is explained below. Please refer to the following text. Figure 3 Physical data signal processing methods include:

[0158] Step 1: Read the data file containing the physical data x(n);

[0159] Step 2: Perform preprocessing on the physical data x(n), including delinearization, zero-mean normalization, and periodic expansion.

[0160] Step 3: Calculate the spectrum X(f) of the data x(n) using Fourier transform;

[0161] Step 4: According to equation (4), calculate the bispectral B(f1,f2) of the data x(n) through triple correlation;

[0162] Step 5: Calculate the bispectral B(f1,f2) of the data x(n) in the entire first quadrant using symmetry.

[0163] Step 6: Take the first peak value of |B(f,f)| as the estimated fundamental frequency f0;

[0164] Step 7: Add a half-fundamental frequency cosine signal with a phase of 0 to the original data signal x(n), i.e., y(n) = x(n) + cos(πf0n);

[0165] Step 8: Calculate the spectrum Y(f) of the data signal y(n) using Fourier transform;

[0166] Step nine, according to equation (4), calculate the bispectral B of the data signal y(n) through triple correlation. y (f1,f2);

[0167] Step 10: Calculate the bispectral B of the data signal y(n) across the entire first quadrant using symmetry. y (f1,f2);

[0168] Step 11: Calculate the harmonic coefficient a according to equation (14). i Phase angle And the number of harmonic terms M;

[0169] Step 12: The harmonic coefficients a obtained in Step 11 are... i Repeat steps two through ten for the calculation;

[0170] Step 13, calculate a according to equation (14) i The harmonic coefficients aa in the harmonic expression i Phase angle With harmonic term number M′;

[0171] Step fourteen, input parameter α (α≥1);

[0172] Step 15, use the harmonic coefficients aa obtained in Step 13. i Phase angle With harmonic term number M′, and a i Multiply the sample spacing by the parameter α input in step fourteen, and substitute it into equation (5) to obtain the reconstructed harmonic coefficient a′. i ;

[0173] Step sixteen, the harmonic coefficients a′ obtained in step fifteen i The phase angle obtained in step eleven Substituting the number of harmonic terms M into equation (5), we obtain the reconstructed signal data.

[0174] Step 17: Perform linear background compensation and mean compensation on the reconstructed signal data to obtain the reconstructed physical data.

[0175] Step 18: Output the calculation results.

[0176] For theoretical seismic wavelet data ( Figure 4 , Figure 5 The theoretical model and model signal in the model were used to perform seismic wavelet reconstruction calculations using the bispectral signal reconstruction method developed in this invention. Figure 4 , Figure 5 (Reconstructed signals in the data). For example, Figure 4 As shown, the reconstructed seismic wavelet dominant frequency is increased by 1.5 times compared to the original signal, that is, the seismic wavelet dominant frequency is increased from 30Hz to 45Hz. Figure 5 As shown, the dominant frequency of the seismic wavelet increased from 30Hz to 270Hz, a ninefold increase. Figure 6A As shown, the reconstructed frequency of actual seismic data using the bispectral signal reconstruction method developed in this invention was increased by 1.5 times. Figure 6B As shown in the figure, the continuity of the reconstructed data reflection layers has been significantly improved.

[0177] This invention applies a signal reconstruction method based on bispectral operations. It calculates the harmonic coefficients, phase angle, and number of harmonic terms of the original signal harmonic expression, as well as the harmonic coefficients, phase angle, and number of harmonic terms of the harmonic coefficients of the original signal harmonic expression. Then, it reconstructs the harmonic coefficients of the original signal harmonic expression and the original signal. When reconstructing the harmonic coefficients of the original signal harmonic expression, the sample spacing of the harmonic coefficients of the original signal harmonic expression is multiplied by an input parameter α (α≥1). The value of the input parameter α is the multiple by which the main frequency of the original signal is boosted.

[0178] The method of this invention is applied to the reconstruction of seismic wavelets in geophysical exploration, which effectively improves the dominant frequency of the wavelets. The application of actual data has significantly improved the continuity of the reflection horizon.

[0179] Example 3

[0180] In this embodiment, a physical data signal processing device is provided, comprising:

[0181] The physical data acquisition module is used to acquire physical data;

[0182] The first expression calculation module is used to preprocess the physical data and perform bispectral calculation to obtain the first harmonic signal expression of the physical data.

[0183] The first harmonic calculation module is used to calculate the first harmonic coefficient, the first phase angle, and the first harmonic term number of the first harmonic signal expression.

[0184] The second expression calculation module is used to preprocess the first harmonic coefficient and perform bispectral calculation to obtain the second harmonic signal expression of the first harmonic coefficient.

[0185] The second harmonic calculation module is used to calculate the second harmonic coefficient, the second phase angle, and the number of second harmonic terms in the expression of the second harmonic signal.

[0186] The amplification module is used to amplify the second harmonic coefficient, the second phase angle, the second harmonic term number and the sample spacing of the first harmonic coefficient based on a preset amplification factor, and substitute the amplified result into the expression of the second harmonic signal to obtain the third harmonic coefficient;

[0187] The reconstructed signal data acquisition module is used to substitute the third harmonic coefficient, the first phase angle, and the first harmonic term number into the first harmonic signal expression to calculate the reconstructed signal data.

[0188] Specific limitations regarding the physical data signal processing device can be found in the limitations of the physical data signal processing method described above, and will not be repeated here. Each unit in the aforementioned physical data signal processing device can be implemented entirely or partially through software, hardware, or a combination thereof. These units can be embedded in or independent of the processor in a computer device in hardware form, or stored in the memory of a computer device in software form, so that the processor can call and execute the operations corresponding to each unit.

[0189] Example 4

[0190] In this embodiment, a computer device is provided. Its internal structure diagram can be shown as follows: Figure 2As shown, the computer device includes a processor, memory, network interface, display screen, and input devices connected via a system bus. The processor provides computing and control capabilities. The memory includes non-volatile storage media and internal memory. The non-volatile storage media stores the operating system and computer programs, and also deploys a database for storing physical data. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage media. The network interface is used to communicate with other computer devices that have deployed application software. When the computer program is executed by the processor, it implements a physical data signal processing method. The display screen can be a liquid crystal display (LCD) or an e-ink display. The input devices can be a touch layer covering the display screen, buttons, a trackball, or a touchpad mounted on the computer device's casing, or an external keyboard, touchpad, or mouse.

[0191] Those skilled in the art will understand that Figure 2 The structure shown is merely a block diagram of a portion of the structure related to the present application and does not constitute a limitation on the computer device to which the present application is applied. Specific computer devices may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.

[0192] In one embodiment, a computer device is provided, including a memory and a processor, the memory storing a computer program, the processor executing the computer program to perform the following steps:

[0193] Step 110: Obtain physical data.

[0194] Specifically, the physical data is geophysical data. In this embodiment, the data file of geophysical data x(n) is read to obtain the physical data x(n).

[0195] Step 120: Preprocess the physical data and perform bispectral calculations to obtain the first harmonic signal expression of the physical data.

[0196] In this step, the physical data x(n) is preprocessed and bispectral calculated sequentially to obtain the harmonic signal expression of the physical data, namely the first harmonic signal expression x(t).

[0197] In one embodiment, preprocessing includes delinearizing the background, zero-meaning, and periodic expansion of the physical data. Bispectral calculation involves calculating the spectrum of the physical data using FFT, and then calculating the bispectral data using a triple correlation method.

[0198] Step 130: Calculate the first harmonic coefficient, first phase angle, and first harmonic term number of the first harmonic signal expression.

[0199] In this step, the harmonic coefficient a of the first harmonic signal expression x(t) is obtained. i With phase angle The number of harmonic terms M, and the harmonic coefficient a i That is, the first harmonic coefficient, phase angle That is, the first phase angle, and the number of harmonic terms M is the number of the first harmonic terms.

[0200] Step 140: Preprocess the first harmonic coefficient and perform bispectral calculation to obtain the second harmonic signal expression of the first harmonic coefficient.

[0201] In this step, the first harmonic coefficient a is processed in the same way as in step 120. i Preprocessing and bispectral calculations were performed to obtain the first harmonic coefficient a. i The expression for the second harmonic signal.

[0202] In one embodiment, preprocessing includes delinearization, zero-mean normalization, and periodic expansion. Bispectral calculation involves calculating the spectrum of the first harmonic coefficient using FFT, and then calculating the bispectral density of the first harmonic coefficient using triple correlation.

[0203] Step 150: Calculate the second harmonic coefficient, second phase angle, and second harmonic term number of the second harmonic signal expression.

[0204] In this step, the harmonic coefficient aa in the second harmonic expression is obtained. i Phase angle With the number of harmonic terms M′, the harmonic coefficient aa i That is, the second harmonic coefficient, phase angle This is the second phase angle, and the number of harmonic terms M′ is the number of the second harmonic terms.

[0205] Step 160: Amplify the second harmonic coefficient, the second phase angle, the second harmonic term number, and the sample spacing of the first harmonic coefficient based on a preset amplification factor, and substitute the amplified result into the expression of the second harmonic signal to obtain the third harmonic coefficient.

[0206] In this embodiment, the preset amplification factor is coefficient α, and the preset amplification factor is greater than or equal to 1, so that the second harmonic coefficient aa i The second phase angle The second harmonic term M′ and the first harmonic coefficient a i Multiplying the sample spacing by a preset amplification factor α yields the third harmonic coefficient a′. iThe third harmonic coefficient a′ i It can also be called the reconstructed harmonic coefficient.

[0207] Step 170: Substitute the third harmonic coefficient, the first phase angle, and the first harmonic term number into the first harmonic signal expression to calculate the reconstructed signal data.

[0208] In this step, the third harmonic coefficient a′ i The first phase angle Substituting the first harmonic term number M into the first harmonic signal expression, the reconstructed signal data x′(t) is obtained.

[0209] In this embodiment, the main frequency of the reconstructed signal x′(t) is increased by α (α≥1) times compared to the main frequency of x(t), thereby realizing the increase of the main frequency of the geophysical data signal. The magnitude of the increase is controlled by the input parameter α (α≥1).

[0210] In one embodiment, the processor, when executing a computer program, also performs the following steps:

[0211] Substituting the third harmonic coefficient, the first phase angle, and the first harmonic term number into the first harmonic signal expression, the reconstructed signal data is obtained.

[0212] The reconstructed signal data is subjected to inverse preprocessing to obtain the reconstructed physical data.

[0213] In one embodiment, the processor, when executing a computer program, also performs the following steps:

[0214] The first harmonic coefficient is preprocessed;

[0215] Perform a Fourier transform on the preprocessed first harmonic coefficient to calculate the spectrum of the first harmonic coefficient;

[0216] Based on the spectrum of the first harmonic coefficient, the bispectrum of the first harmonic coefficient in the first quadrant is calculated;

[0217] The fundamental frequency of the bispectral spectrum of the first harmonic coefficient is estimated from the bispectral spectrum of the first harmonic coefficient in the first quadrant.

[0218] Based on the fundamental frequency of the bispectral spectrum of the first harmonic coefficient, the first harmonic coefficient is phase-shifted to obtain the phase-shifted first harmonic coefficient.

[0219] Based on the phase-shifted first harmonic coefficient, the bispectral data of the phase-shifted first harmonic coefficient in the first quadrant is calculated.

[0220] Based on the bispectrum of the first harmonic coefficient in the first quadrant, the expression for the first harmonic signal of the bispectrum of the first harmonic coefficient is calculated.

[0221] In one embodiment, the processor, when executing a computer program, also performs the following steps:

[0222] The physical data is preprocessed and bispectral calculations are performed to obtain the bispectral data of the physical data;

[0223] Based on the bispectral data of the physical data, the expression for the first harmonic signal of the bispectral data of the physical data is calculated.

[0224] In one embodiment, the processor, when executing a computer program, also performs the following steps:

[0225] The physical data is preprocessed;

[0226] The Fourier transform is performed on the preprocessed physical data to calculate the spectrum of the physical data;

[0227] Based on the spectrum of the physical data, the bispectrum of the physical data in the first quadrant is calculated;

[0228] The fundamental frequency of the bispectral spectrum of the physical data is estimated from the bispectral spectrum in the first quadrant of the physical data;

[0229] Based on the fundamental frequency, the physical data is phase-shifted to obtain the phase-shifted physical data;

[0230] Based on the phase-shifted physical data, the bispectral data in the first quadrant is calculated.

[0231] In one embodiment, the processor, when executing a computer program, also performs the following steps:

[0232] Based on the spectrum of the physical data, the bispectrum of the physical data is calculated;

[0233] Based on the bispectral data of the physical data, the bispectral data in the first quadrant is calculated using symmetry.

[0234] In one embodiment, the processor, when executing a computer program, also performs the following steps:

[0235] The physical data is subjected to delinear background processing, zero-mean processing, and period expansion processing.

[0236] Example 5

[0237] In this embodiment, a computer-readable storage medium is provided, on which a computer program is stored. When the computer program is executed by a processor, it performs the following steps:

[0238] Step 110: Obtain physical data.

[0239] Specifically, the physical data is geophysical data. In this embodiment, the data file of geophysical data x(n) is read to obtain the physical data x(n).

[0240] Step 120: Preprocess the physical data and perform bispectral calculations to obtain the first harmonic signal expression of the physical data.

[0241] In this step, the physical data x(n) is preprocessed and bispectral calculated sequentially to obtain the harmonic signal expression of the physical data, namely the first harmonic signal expression x(t).

[0242] In one embodiment, preprocessing includes delinearizing the background, zero-meaning, and periodic expansion of the physical data. Bispectral calculation involves calculating the spectrum of the physical data using FFT, and then calculating the bispectral data using a triple correlation method.

[0243] Step 130: Calculate the first harmonic coefficient, first phase angle, and first harmonic term number of the first harmonic signal expression.

[0244] In this step, the harmonic coefficient a of the first harmonic signal expression x(t) is obtained. i With phase angle The number of harmonic terms M, and the harmonic coefficient a i That is, the first harmonic coefficient, phase angle That is, the first phase angle, and the number of harmonic terms M is the number of the first harmonic terms.

[0245] Step 140: Preprocess the first harmonic coefficient and perform bispectral calculation to obtain the second harmonic signal expression of the first harmonic coefficient.

[0246] In this step, the first harmonic coefficient a is processed in the same way as in step 120. i Preprocessing and bispectral calculations were performed to obtain the first harmonic coefficient a. i The expression for the second harmonic signal.

[0247] In one embodiment, preprocessing includes delinearization, zero-mean normalization, and periodic expansion. Bispectral calculation involves calculating the spectrum of the first harmonic coefficient using FFT, and then calculating the bispectral density of the first harmonic coefficient using triple correlation.

[0248] Step 150: Calculate the second harmonic coefficient, second phase angle, and second harmonic term number of the second harmonic signal expression.

[0249] In this step, the harmonic coefficient aa in the second harmonic expression is obtained. i Phase angle With the number of harmonic terms M′, the harmonic coefficient aa i That is, the second harmonic coefficient, phase angle This is the second phase angle, and the number of harmonic terms M′ is the number of the second harmonic terms.

[0250] Step 160: Amplify the second harmonic coefficient, the second phase angle, the second harmonic term number, and the sample spacing of the first harmonic coefficient based on a preset amplification factor, and substitute the amplified result into the expression of the second harmonic signal to obtain the third harmonic coefficient.

[0251] In this embodiment, the preset amplification factor is coefficient α, and the preset amplification factor is greater than or equal to 1, so that the second harmonic coefficient aa i The second phase angle The second harmonic term M′ and the first harmonic coefficient a i Multiplying the sample spacing by a preset amplification factor α yields the third harmonic coefficient a′. i The third harmonic coefficient a′ i It can also be called the reconstructed harmonic coefficient.

[0252] Step 170: Substitute the third harmonic coefficient, the first phase angle, and the first harmonic term number into the first harmonic signal expression to calculate the reconstructed signal data.

[0253] In this step, the third harmonic coefficient a i ′、First phase angle Substituting the first harmonic term number M into the first harmonic signal expression, the reconstructed signal data x′(t) is obtained.

[0254] In this embodiment, the main frequency of the reconstructed signal x′(t) is increased by α (α≥1) times compared to the main frequency of x(t), thereby realizing the increase of the main frequency of the geophysical data signal. The magnitude of the increase is controlled by the input parameter α (α≥1).

[0255] In one embodiment, when the computer program is executed by a processor, it also performs the following steps:

[0256] Substituting the third harmonic coefficient, the first phase angle, and the first harmonic term number into the first harmonic signal expression, the reconstructed signal data is obtained.

[0257] The reconstructed signal data is subjected to inverse preprocessing to obtain the reconstructed physical data.

[0258] In one embodiment, when the computer program is executed by a processor, it also performs the following steps:

[0259] The first harmonic coefficient is preprocessed;

[0260] Perform a Fourier transform on the preprocessed first harmonic coefficient to calculate the spectrum of the first harmonic coefficient;

[0261] Based on the spectrum of the first harmonic coefficient, the bispectrum of the first harmonic coefficient in the first quadrant is calculated;

[0262] The fundamental frequency of the bispectral spectrum of the first harmonic coefficient is estimated from the bispectral spectrum of the first harmonic coefficient in the first quadrant.

[0263] Based on the fundamental frequency of the bispectral spectrum of the first harmonic coefficient, the first harmonic coefficient is phase-shifted to obtain the phase-shifted first harmonic coefficient.

[0264] Based on the phase-shifted first harmonic coefficient, the bispectral data of the phase-shifted first harmonic coefficient in the first quadrant is calculated.

[0265] Based on the bispectrum of the first harmonic coefficient in the first quadrant, the expression for the first harmonic signal of the bispectrum of the first harmonic coefficient is calculated.

[0266] In one embodiment, when the computer program is executed by a processor, it also performs the following steps:

[0267] The physical data is preprocessed and bispectral calculations are performed to obtain the bispectral data of the physical data;

[0268] Based on the bispectral data of the physical data, the expression for the first harmonic signal of the bispectral data of the physical data is calculated.

[0269] In one embodiment, when the computer program is executed by a processor, it also performs the following steps:

[0270] The physical data is preprocessed;

[0271] The Fourier transform is performed on the preprocessed physical data to calculate the spectrum of the physical data;

[0272] Based on the spectrum of the physical data, the bispectrum of the physical data in the first quadrant is calculated;

[0273] The fundamental frequency of the bispectral spectrum of the physical data is estimated from the bispectral spectrum in the first quadrant of the physical data;

[0274] Based on the fundamental frequency, the physical data is phase-shifted to obtain the phase-shifted physical data;

[0275] Based on the phase-shifted physical data, the bispectral data in the first quadrant is calculated.

[0276] In one embodiment, when the computer program is executed by a processor, it also performs the following steps:

[0277] Based on the spectrum of the physical data, the bispectrum of the physical data is calculated;

[0278] Based on the bispectral data of the physical data, the bispectral data in the first quadrant is calculated using symmetry.

[0279] In one embodiment, when the computer program is executed by a processor, it also performs the following steps:

[0280] The physical data is subjected to delinear background processing, zero-mean processing, and period expansion processing.

[0281] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium, and when executed, it can include the processes of the embodiments of the above methods. Any references to memory, storage, databases, or other media used in the embodiments provided in this application can include non-volatile and / or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in various forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), dual data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), Rambus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.

[0282] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0283] The embodiments described above are merely illustrative of several implementation methods of this application, and while the descriptions are relatively specific and detailed, they should not be construed as limiting the scope of the invention patent. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of this application, and these all fall within the protection scope of this application. Therefore, the protection scope of this patent application should be determined by the appended claims.

Claims

1. A physical data signal processing method, characterized in that, include: Acquire physical data; The physical data is preprocessed and bispectral calculations are performed to obtain the first harmonic signal expression of the physical data; Calculate the first harmonic coefficient, the first phase angle, and the number of first harmonic terms in the expression of the first harmonic signal; The first harmonic coefficient is preprocessed and bispectral calculated to obtain the second harmonic signal expression of the first harmonic coefficient; Calculate the second harmonic coefficient, second phase angle, and second harmonic term number in the expression of the second harmonic signal; The second harmonic coefficient, the second phase angle, the second harmonic term number, and the sample spacing of the first harmonic coefficient are amplified based on a preset amplification factor, and the amplified result is substituted into the expression of the second harmonic signal to obtain the third harmonic coefficient. Substituting the third harmonic coefficient, the first phase angle, and the first harmonic term number into the first harmonic signal expression, the reconstructed signal data is obtained.

2. The method according to claim 1, characterized in that, The step of substituting the third harmonic coefficient, the first phase angle, and the first harmonic term number into the first harmonic signal expression to obtain the reconstructed physical data includes: Substitute the third harmonic coefficient, the first phase angle, and the first harmonic term number into the first harmonic signal expression to calculate the reconstructed signal data. The reconstructed signal data is subjected to inverse preprocessing to obtain the reconstructed physical data.

3. The method according to claim 1, characterized in that, The steps of preprocessing and bispectral calculation of the first harmonic coefficient to obtain the bispectral data of the first harmonic coefficient include: The first harmonic coefficient is preprocessed; Perform a Fourier transform on the preprocessed first harmonic coefficient to calculate the spectrum of the first harmonic coefficient; Based on the spectrum of the first harmonic coefficient, the bispectrum of the first harmonic coefficient in the first quadrant is calculated; The fundamental frequency of the bispectral spectrum of the first harmonic coefficient is estimated from the bispectral spectrum of the first harmonic coefficient in the first quadrant. Based on the fundamental frequency of the bispectral spectrum of the first harmonic coefficient, the first harmonic coefficient is phase-shifted to obtain the phase-shifted first harmonic coefficient. Based on the phase-shifted first harmonic coefficient, the bispectral data of the phase-shifted first harmonic coefficient in the first quadrant is calculated. Based on the bispectrum of the first harmonic coefficient in the first quadrant, the expression for the first harmonic signal of the bispectrum of the first harmonic coefficient is calculated.

4. The method according to claim 1, characterized in that, The steps of preprocessing the physical data and performing bispectral calculations to obtain the first harmonic signal expression of the physical data include: The physical data is preprocessed and bispectral calculations are performed to obtain the bispectral data of the physical data; Based on the bispectral data of the physical data, the expression for the first harmonic signal of the bispectral data of the physical data is calculated.

5. The method according to claim 4, characterized in that, The steps of preprocessing the physical data and performing bispectral calculation to obtain the bispectral data include: The physical data is preprocessed; The Fourier transform is performed on the preprocessed physical data to calculate the spectrum of the physical data; Based on the spectrum of the physical data, the bispectrum of the physical data in the first quadrant is calculated; The fundamental frequency of the bispectral spectrum of the physical data is estimated from the bispectral spectrum in the first quadrant of the physical data; Based on the fundamental frequency, the physical data is phase-shifted to obtain the phase-shifted physical data; Based on the phase-shifted physical data, the bispectral data in the first quadrant is calculated.

6. The method according to claim 5, characterized in that, The step of calculating the bispectrum of the physical data in the first quadrant based on the spectrum of the physical data includes: Based on the spectrum of the physical data, the bispectrum of the physical data is calculated; Based on the bispectral data of the physical data, the bispectral data in the first quadrant is calculated using symmetry.

7. The method according to claim 5, characterized in that, The preprocessing steps for the physical data include: The physical data is subjected to delinear background processing, zero-mean processing, and period expansion processing.

8. A physical data signal processing device, characterized in that, include: The physical data acquisition module is used to acquire physical data; The first expression calculation module is used to preprocess the physical data and perform bispectral calculation to obtain the first harmonic signal expression of the physical data. The first harmonic calculation module is used to calculate the first harmonic coefficient, the first phase angle, and the first harmonic term number of the first harmonic signal expression. The second expression calculation module is used to preprocess the first harmonic coefficient and perform bispectral calculation to obtain the second harmonic signal expression of the first harmonic coefficient. The second harmonic calculation module is used to calculate the second harmonic coefficient, the second phase angle, and the number of second harmonic terms in the expression of the second harmonic signal. The amplification module is used to amplify the second harmonic coefficient, the second phase angle, the second harmonic term number and the sample spacing of the first harmonic coefficient based on a preset amplification factor, and substitute the amplified result into the expression of the second harmonic signal to obtain the third harmonic coefficient; The reconstructed signal data acquisition module is used to substitute the third harmonic coefficient, the first phase angle, and the first harmonic term number into the first harmonic signal expression to calculate the reconstructed signal data.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that, When the processor executes the computer program, it implements the steps of the method according to any one of claims 1 to 7.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 7.

Images