Seismic data radon filtering method and system

Abstract

The application discloses a kind of seismic data Radon filtering method and system, system includes acquisition module, positive transformation module, filtering module and inverse transformation module.First, acquisition module obtains two-way travel time vector, slowness vector, seismic data, source three-dimensional coordinates, geophone three-dimensional coordinates.Then, positive transformation module can be converted into Radon domain value data without prior height correction with space-time domain value data.Afterwards, filtering module generates corresponding filter according to the slowness range to be retained, and filters the positive transformation result by filter.Finally, inverse transformation module converts back to space-time domain value data by the mapping relationship of each space-time domain value recorded by positive transformation module conversion to Radon domain value, converts Radon domain value data after filtering out interference part back to space-time domain value data.In this way, without the analytical formula of Radon inverse transformation, Radon filtering of seismic data under the condition of undulating surface can be realized, and the real reflected wave signal is obtained.

Description

Technical Field

[0001] This invention relates to the field of seismic data processing technology, specifically to a Radon filtering method and system for seismic data. Background Technology

[0002] Seismic reflection data contains not only reflected waves but also other interfering waves, such as surface waves, acoustic waves, and multiples. These interfering waves need to be filtered out to obtain a cleaner seismic reflection record. Radon filtering is one effective method for removing these interfering waves. By converting the spatiotemporal domain data to the Radon domain, filtering out the interfering components, and then converting it back to the spatiotemporal domain, the effective waves and interfering waves can be separated.

[0003] Traditional Radon filtering algorithms are applicable to horizontal terrain. For undulating terrain, height correction is required to align the waves to a common reference level. However, this method necessitates the wave velocity of the shallow, low-velocity zone, which is difficult to accurately obtain. Since undulating terrain significantly impacts the travel time of shallow reflected waves more than deep reflected waves, inaccurate shallow wave velocities will lead to even more severe errors after height correction in shallow reflected wave exploration. Therefore, traditional Radon filtering algorithms are unsuitable for shallow reflected wave seismic exploration on undulating terrain. Summary of the Invention

[0004] To address the shortcomings of existing technologies, this invention proposes a Radon filtering method and system for seismic data. It can obtain near-realistic reflected wave signals without the need for height correction. The specific technical solution is as follows:

[0005] Firstly, a Radon filtering method for seismic data is provided, including:

[0006] Obtain the two-way travel time vector, slowness vector, as well as seismic data and the corresponding three-dimensional coordinates of the source and the detector;

[0007] Based on the three-dimensional coordinates of the seismic source, the three-dimensional coordinates of the geophone, the two-way travel time vector, and the slowness vector, the spatiotemporal domain values ​​of each trace in the seismic data are converted into Radon domain values ​​to obtain the forward transformation result data, and the mapping relationship between the spatiotemporal domain values ​​of each trace and the corresponding Radon domain values ​​is recorded.

[0008] The forward transform result data is filtered by the slowness range that needs to be retained to obtain the filtered result data;

[0009] Based on the mapping relationship between the spatiotemporal domain values ​​and the corresponding Radon domain values, the filtered result data is inversely transformed to obtain the inverse transformation result data.

[0010] In conjunction with the first aspect, in the first possible implementation of the first aspect, converting spatiotemporal domain values ​​into Radon domain values ​​includes:

[0011] By using the three-dimensional coordinates of the seismic source, the three-dimensional coordinates of the geophone, the two-way travel time vector, and the slowness vector, the two-way travel time corresponding to each trace in the seismic data at different time nodes and different slowness nodes is calculated.

[0012] Based on the corresponding two-way travel time, the amplitude values ​​corresponding to each trace are extracted from the seismic data and summed to obtain the Radon domain values ​​of each trace at different time nodes and different slowness nodes.

[0013] In conjunction with the first possible implementation of the first aspect, in the second possible implementation of the first aspect, calculating the round-trip journey includes:

[0014]

[0015]

[0016] d = CR z -CP z ;

[0017] CP x CP y CP z CR represents the three-dimensional coordinates of the earthquake source. x CR y CR z Let t be the three-dimensional coordinates of the detector, t be the two-way travel time vector node, p be the slowness vector node, and T be the two-way travel time of the target track at that time node and the slowness node.

[0018] In conjunction with the first aspect, in the third possible implementation of the first aspect, the forward transform result data is filtered by the slowness range that needs to be retained, including:

[0019] Generate a filter based on the required range of slowness.

[0020] The forward transform result data is filtered using the filter to obtain the filtered result data.

[0021] In conjunction with the third possible implementation of the first aspect, in the fourth possible implementation of the first aspect, a filter is generated according to the required retention of the slowness range, including:

[0022] Determine the vector position of the slowness range that needs to be retained within the slowness vector;

[0023] The initial filter is adjusted based on the vector position and the vector length during two-way travel to obtain the adjusted filter.

[0024] Secondly, a Radon filtering system for seismic data is provided, including:

[0025] The acquisition module is configured to acquire two-way travel time vectors, slowness vectors, as well as seismic data and the corresponding three-dimensional coordinates of the source and the detector.

[0026] The forward transformation module is configured to convert the spatiotemporal domain values ​​of each trace in the seismic data into Radon domain values ​​based on the three-dimensional coordinates of the seismic source, the three-dimensional coordinates of the detector, the two-way travel time vector, and the slowness vector, to obtain the forward transformation result data, and record the mapping relationship between the spatiotemporal domain values ​​of each trace and the corresponding Radon domain values.

[0027] The filtering module is configured to filter the forward transform result data through the slowness range that needs to be retained, and obtain the filtered result data.

[0028] The inverse transform module is configured to perform an inverse transform on the filtered result data based on the mapping relationship between the spatiotemporal domain values ​​and the corresponding Radon domain values ​​to obtain the inverse transform result data.

[0029] In conjunction with the second aspect, in a first possible implementation of the second aspect, the forward transformation module includes:

[0030] The calculation unit is configured to calculate the two-way travel time of each trace in the seismic data at different time nodes and different slowness nodes using the three-dimensional coordinates of the source, the three-dimensional coordinates of the detector, the two-way travel time vector, and the slowness vector.

[0031] The summation unit is configured to extract the amplitude values ​​corresponding to each trace from the seismic data based on the corresponding two-way travel time and sum them to obtain the Radon domain values ​​of each trace at different time nodes and different slowness nodes.

[0032] In conjunction with the first possible implementation of the second aspect, in the second possible implementation of the second aspect, the computing unit calculates the round-trip journey by including:

[0033]

[0034]

[0035] d = CR z -CP z ;

[0036] CP x CP y CP z CR represents the three-dimensional coordinates of the earthquake source. x CR y CR zLet t be the three-dimensional coordinates of the detector, t be the two-way travel time vector node, p be the slowness vector node, and T be the two-way travel time of the target track at that time node and the slowness node.

[0037] In conjunction with the second aspect, in a third possible implementation of the second aspect, the filtering module includes:

[0038] The generation unit is configured to generate filters for the slowness range that needs to be retained.

[0039] The filtering unit is configured to filter the forward transform result data through the filter to obtain filtered result data.

[0040] In conjunction with the third possible implementation of the second aspect, in the fourth possible implementation of the second aspect, the generating unit includes:

[0041] The vector position determination unit is configured to determine the vector position of the slowness range to be retained within the slowness vector;

[0042] The filter adjustment unit is configured to adjust the initial filter based on the vector position and the vector length of the two-way travel time vector to obtain the adjusted filter.

[0043] Beneficial Effects: The seismic data Radon filtering method and system of this invention, by introducing the three-dimensional coordinates of the source and detector during the Radon forward transform, can convert spatiotemporal domain data into Radon domain data without prior height correction. Furthermore, by recording the mapping relationship between each spatiotemporal domain transformation and the Radon domain transformation, the Radon domain data after filtering out interference can be converted back to spatiotemporal domain data without the need for an inverse transform formula. This achieves Radon filtering of seismic data under undulating surface conditions, obtaining the true reflected wave signal. Attached Figure Description

[0044] To more clearly illustrate the specific embodiments of the present invention, the accompanying drawings used in the specific embodiments will be briefly described below. In all the drawings, the elements or parts are not necessarily drawn to scale.

[0045] Figure 1 A flowchart of a Radon filtering method for seismic data provided in an embodiment of the present invention;

[0046] Figure 2 A flowchart of the forward transformation provided in an embodiment of the present invention;

[0047] Figure 3 This is a flowchart illustrating the computer's calculation of seismic data filtering according to an embodiment of the present invention.

[0048] Figure 4This is a system block diagram of a Radon filtering system for seismic data provided in an embodiment of the present invention;

[0049] Figure 5 This is a system block diagram of a Radon filtering system for seismic data provided in an embodiment of the present invention;

[0050] Figure 6 This is a system block diagram of a Radon filtering system for seismic data provided in an embodiment of the present invention;

[0051] Figure 7 This is a system block diagram of a Radon filtering system for seismic data provided in an embodiment of the present invention. Detailed Implementation

[0052] The embodiments of the technical solution of the present invention will now be described in detail with reference to the accompanying drawings. These embodiments are merely illustrative of the technical solution of the present invention and are therefore intended to limit the scope of protection of the present invention.

[0053] like Figure 1 The flowchart shown illustrates the Radon filtering method for seismic data. This filtering method includes:

[0054] Step 1: Obtain the two-way travel time vector, slowness vector, seismic data, and corresponding source 3D coordinates and detector 3D coordinates;

[0055] Step 2: Based on the three-dimensional coordinates of the seismic source, the three-dimensional coordinates of the detector, the two-way travel time vector, and the slowness vector, convert the spatiotemporal domain values ​​of each trace in the seismic data into Radon domain values ​​to obtain the forward transformation result data, and record the mapping relationship between the spatiotemporal domain values ​​of each trace and the corresponding Radon domain values.

[0056] Step 3: Filter the forward transform result data according to the slowness range that needs to be retained to obtain the filtered result data;

[0057] Step 4: Based on the mapping relationship between the spatiotemporal domain values ​​and the corresponding Radon domain values, perform an inverse transform on the filtered result data to obtain the inverse transform result data.

[0058] Specifically, firstly, the spatiotemporal domain data included in the seismic data can be converted into Radon threshold data using the Radon forward transform. During the forward transform, the three-dimensional coordinates of the source and detector can be introduced, allowing the spatiotemporal domain data to be converted into Radon threshold data without prior height correction. The mapping relationship between each spatiotemporal domain data point and the converted Radon threshold data is recorded during the transformation process. Then, a corresponding filter can be generated based on the desired slowness range, and the forward transform result is filtered using this filter. Because three-dimensional coordinates are introduced, the analytical formula for the inverse Radon transform cannot be derived. Therefore, the Radon threshold data after filtering out interference can be converted back to spatiotemporal domain data by recording the mapping relationship between each spatiotemporal domain value and the Radon threshold value. Thus, Radon filtering of seismic data under undulating surface conditions can be achieved without the need for the analytical formula for the inverse Radon transform, obtaining the true reflected wave signal.

[0059] In this embodiment, optionally, in step 2, converting the spatiotemporal domain values ​​into Radon domain values ​​includes:

[0060] Step 2-1: Calculate the two-way travel time of each trace in the seismic data at different time nodes and different slowness nodes using the three-dimensional coordinates of the source, the three-dimensional coordinates of the detector, the two-way travel time vector, and the slowness vector.

[0061] Step 2-2: Extract the amplitude values ​​corresponding to each trace from the seismic data according to the corresponding two-way travel time and sum them to obtain the Radon threshold values ​​of each trace at different time nodes and different slowness nodes.

[0062] Specifically, firstly, based on the three-dimensional coordinates of the seismic source, the three-dimensional coordinates of the geophone, the two-way travel time vector, and the slowness vector, the two-way travel time T of each trace in the seismic data at different time nodes and different slowness nodes can be calculated. The specific formula for calculating the two-way travel time T is as follows:

[0063]

[0064]

[0065] d = CR z -CP z ;

[0066] CP x CP y CP z CR represents the three-dimensional coordinates of the earthquake source. x CR y CR zLet t be the three-dimensional coordinates of the detector, t be the two-way travel time vector node, p be the slowness vector node, and T be the two-way travel time of the target track at that time node and the slowness node.

[0067] Then, the amplitude values ​​corresponding to each trace at time T can be extracted from the seismic data, and the amplitude values ​​corresponding to each trace can be added together to obtain the sum of the amplitude values ​​corresponding to that time node and the slowness node, which is the Radon threshold value corresponding to that time node and the slowness node. Combining the Radon threshold values ​​corresponding to different time nodes and different slowness nodes yields the forward transform result data.

[0068] The following will combine Figure 3 The specific calculation process of the Radon forward transform is explained in detail.

[0069] In this embodiment, the computer can perform a Radon forward transform on the seismic data according to the following calculation process:

[0070] S1, Input earthquake data to the computer: Data, three-dimensional coordinates of the epicenter, three-dimensional coordinates of the detector, and two-way travel time vector t(t1, t2, t3...t...). n And the slowness vector p(p1, p2, p3 ... p) n );

[0071] S2. Let the initial cycle number j = 1 for the two-way time vector t;

[0072] S3. Let the initial cycle number i = 1 for the slowness vector p;

[0073] S4. Set the initial cycle number k of the seismic data trace number to 1;

[0074] S5. Set the initial amplitude value Amp = 0, and calculate the round-trip travel time T of the kth track according to the above formula;

[0075] S6. Retrieve the index T of T in the two-way time vector t. num And extract the amplitude value Data(T) of the k-th channel at time T. num ,k);

[0076] S7, Transfer Data(T) nyY Add (k) to Amp and replace the value of Amp, then save Y. num The value is entered into the corresponding position (i, k, j) of the three-dimensional array trecord;

[0077] S8. Determine if k is less than the maximum number of traces in the seismic data Data. If not, proceed to the next step. If yes, set k = k + 1 and repeat S5-S8.

[0078] S9. Save the Amp value to the corresponding position (j,i) of the positive transformation result Radon(j,i);

[0079] S10. Determine if i is less than the maximum index of the slowness vector p. If not, proceed to the next step. If yes, set i = i + 1 and repeat S4-S10.

[0080] S11. Determine if j is less than the maximum index of the two-way time vector t. If not, end the process and obtain the forward transform result data. It should be understood that the obtained forward transform result data is a set of Radon threshold arrays. If yes, let j = j + 1 and repeat S3-S11.

[0081] In this embodiment, optionally, in step 3, the forward transform result data is filtered according to the slowness range that needs to be retained, including:

[0082] Generate a filter based on the required range of slowness.

[0083] The forward transform result data is filtered using the filter to obtain the filtered result data.

[0084] Specifically, firstly, a filter can be generated based on the desired range of slowness. Then, the forward transform result data is filtered using the generated filter to obtain the filtered result data after removing interference. When generating the filter, the vector position of the desired range of slowness within the slowness vector can be determined first. Then, the initial filter is adjusted based on the vector position and the length of the two-way travel time vector to obtain the adjusted filter.

[0085] For example, the slowness ranges that need to be retained are divided into p1, p2, p3, and p4. Here, p1 and p2 represent the slowness ranges corresponding to the minimum value of the two-way travel time vector t (corresponding to the first row of data in the forward transformation result), physically representing the surface slowness intervals. p3 and p4 represent the slowness ranges corresponding to the maximum value of the two-way travel time vector t (corresponding to the last row of data in the forward transformation result), physically representing the deepest slowness intervals reached by seismic waves.

[0086] When generating the filter, firstly, the positions of p1, p2, p3, and p4 in the slowness vector p are sequentially retrieved, resulting in Cp1, Cp2, Cp3, and Cp4 respectively. Then, the initial filter (filter) and the forward transform result data (i.e., the Radon threshold array) are set to the same size, with all values ​​set to 1. Next, the left node (Cleft) and right node (Cright) of the retention range for each row of the initial filter (filter) are calculated using the following formula, ensuring that the data within the left and right nodes of each row remain unchanged, while the data outside these nodes are set to 0, thus obtaining the adjusted filter. The specific calculation formula is as follows:

[0087]

[0088]

[0089] Where r is the row number of the initial filter and N is the length of the time vector t. Finally, performing a matrix dot product between the adjusted filter and the Radon threshold array yields the filtered Radon array, which is the filtered result data.

[0090] In this embodiment, after obtaining the filtered result data, the Radon domain value data after filtering out the interference can be converted back to the spatiotemporal domain value data by recording the mapping relationship between each spatiotemporal domain value and the Radon domain value.

[0091] The following will combine Figure 3 The specific calculation process of the Radon inverse transform is explained in detail.

[0092] In this embodiment, the computer can perform an inverse Radon transform on the filtered result data according to the following calculation process:

[0093] SS1. Set the initial loop number I of the original data channel number to 1;

[0094] SS2. Set the initial cycle number J of the time vector t to 1;

[0095] SS3. Set the initial page number K of the three-dimensional array trecord to 1;

[0096] SS4. Set the initial amplitude value Amp = 0;

[0097] SS5. Set the initial row number kk of the three-dimensional array trecord to 1;

[0098] SS6. Determine if trecord(kk,I,K) equals J. If yes, proceed to the next step. If no, jump to SS8.

[0099] SS7. Determine if Radon(K,kk) is equal to 0. If yes, add Radon(K,kk) to Amp and replace the value of Amp. If no, proceed to the next step.

[0100] SS8. Determine if kk is less than the maximum row number of trecord. If not, proceed to the next step. If yes, then kk = kk + 1, and repeat SS6-SS8.

[0101] SS9. Determine if K is less than the maximum page number of trecord. If not, proceed to the next step. If yes, then...

[0102] K = K + 1, repeat SS5-S9;

[0103] SS10. Save the Amp value to the corresponding position (J,I) of the inverse transform result data;

[0104] SS11. Determine if J is less than the maximum index of t. If not, proceed to the next step. If yes, then J = J + 1.

[0105] Repeat SS3-SS11;

[0106] SS12. Determine if I is less than the maximum number of channels in the original data Data. If not, proceed to the next step. If yes, then I = I + 1, and repeat SS2-SS12.

[0107] SS13. Save the inverse transform result data.

[0108] like Figure 4 The diagram shown is a system block diagram of the Radon filtering system for seismic data. The filtering system includes:

[0109] The acquisition module is configured to acquire two-way travel time vectors, slowness vectors, as well as seismic data and the corresponding three-dimensional coordinates of the source and the detector.

[0110] The forward transformation module is configured to convert the spatiotemporal domain values ​​of each trace in the seismic data into Radon domain values ​​based on the three-dimensional coordinates of the seismic source, the three-dimensional coordinates of the detector, the two-way travel time vector, and the slowness vector, to obtain the forward transformation result data, and record the mapping relationship between the spatiotemporal domain values ​​of each trace and the corresponding Radon domain values.

[0111] The filtering module is configured to filter the forward transform result data through the slowness range that needs to be retained, and obtain the filtered result data.

[0112] The inverse transform module is configured to perform an inverse transform on the filtered result data based on the mapping relationship between the spatiotemporal domain values ​​and the corresponding Radon domain values ​​to obtain the inverse transform result data.

[0113] Specifically, the filtering system consists of an acquisition module, a forward transform module, a filtering module, and an inverse transform module. The acquisition module acquires the two-way travel time vector, the slowness vector, and the seismic data along with the corresponding three-dimensional coordinates of the source and the detector. The forward transform module uses Radon forward transform to convert the spatiotemporal domain data included in the seismic data into Radon domain data. During the forward transform, the three-dimensional coordinates of the source and detector can be incorporated, allowing the conversion of spatiotemporal domain data into Radon domain data without prior height correction. The module also records the mapping relationship between each spatiotemporal domain data point and the converted Radon domain data during the transformation process.

[0114] The filtering module can generate corresponding filters based on the desired slowness range, and then filter the forward transform results using these filters. Because three-dimensional coordinates are introduced, the analytical formula for the inverse Radon transform cannot be derived. Therefore, the inverse transform module can use the mapping relationship between each spatiotemporal domain value recorded by the forward transform module and the Radon domain value to convert the Radon domain value data (after filtering out interference) back to spatiotemporal domain value data. In this way, Radon filtering of seismic data under undulating surface conditions can be achieved without the need for the analytical formula for the inverse Radon transform, obtaining the true reflected wave signal.

[0115] In this embodiment, optionally, such as Figure 5 As shown, the forward transformation module includes:

[0116] The calculation unit is configured to calculate the two-way travel time of each trace in the seismic data at different time nodes and different slowness nodes using the three-dimensional coordinates of the source, the three-dimensional coordinates of the detector, the two-way travel time vector, and the slowness vector.

[0117] The summation unit is configured to extract the amplitude values ​​corresponding to each trace from the seismic data based on the corresponding two-way travel time and sum them to obtain the Radon domain values ​​of each trace at different time nodes and different slowness nodes.

[0118] Specifically, the forward transform module includes a calculation unit and a summation unit. The calculation unit can calculate the two-way travel time (T) of each trace in the seismic data at different time nodes and different slowness nodes based on the three-dimensional coordinates of the seismic source, the three-dimensional coordinates of the detector, the two-way travel time vector, and the slowness vector. The specific formula for calculating the two-way travel time T is as follows:

[0119]

[0120]

[0121] d = CR z -CP z ;

[0122] CP x CP y CP z CR represents the three-dimensional coordinates of the earthquake source. x CR y CR z Let t be the three-dimensional coordinates of the detector, t be the two-way travel time vector node, p be the slowness vector node, and T be the two-way travel time of the target track at that time node and the slowness node.

[0123] The summation unit can extract the amplitude values ​​corresponding to each trace at time T from the seismic data, and sum the amplitude values ​​of each trace to obtain the sum of the amplitude values ​​corresponding to that time node and the slowness node, which is the Radon threshold value corresponding to that time node and the slowness node. Combining the Radon threshold values ​​corresponding to different time nodes and slowness nodes yields the forward transform result data.

[0124] In this embodiment, optionally, such as Figure 6 As shown, the filtering module includes:

[0125] The generation unit is configured to generate filters for the slowness range that needs to be retained.

[0126] The filtering unit is configured to filter the forward transform result data through the filter to obtain filtered result data.

[0127] Specifically, the filtering module includes a generation unit and a filtering unit. The generation unit generates a corresponding filter based on the desired slowness range. The filtering unit uses the filter generated by the generation unit to filter the forward transform result data.

[0128] In this embodiment, as Figure 7 As shown, the generation unit includes a vector position determination unit and a filter adjustment unit. The vector position determination unit can determine the vector position of the slowness range to be retained within the slowness vector. The filter adjustment unit can adjust the initial filter based on the vector position and the vector length of the two-way travel time vector to obtain the adjusted filter.

[0129] For example, the slowness ranges that need to be retained are divided into p1, p2, p3, and p4. Here, p1 and p2 represent the slowness ranges corresponding to the minimum value of the two-way travel time vector t (corresponding to the first row of data in the forward transformation result), physically representing the surface slowness intervals. p3 and p4 represent the slowness ranges corresponding to the maximum value of the two-way travel time vector t (corresponding to the last row of data in the forward transformation result), physically representing the deepest slowness intervals reached by seismic waves.

[0130] When generating the filter, the vector position determination unit sequentially retrieves the positions of p1, p2, p3, and p4 in the slowness vector p, which are Cp1, Cp2, Cp3, and Cp4, respectively. Then, the initial filter (filter) and the forward transform result data (i.e., the Radon threshold array) are made to have the same array size and all values ​​are set to 1. The filter adjustment unit calculates the left node (Cleft) and right node (Cright) of the retention range for each row of the initial filter (filter) using the following formula, keeping the data within the left and right nodes unchanged and setting the data outside to 0, thus obtaining the adjusted filter. The filtering unit performs a matrix multiplication between the adjusted filter (filter) and the Radon threshold array to obtain the filtered Radon array, i.e., the filtered result data.

[0131] The above embodiments are only used to illustrate the technical solutions of the present invention, and are not intended to limit it. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some or all of the technical features therein. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present invention, and they should all be covered within the scope of the claims and specification of the present invention.

Claims

1. A Radon filtering method for seismic data, characterized in that, include: Obtain the two-way travel time vector, slowness vector, as well as seismic data and the corresponding three-dimensional coordinates of the source and the detector; Based on the three-dimensional coordinates of the seismic source, the three-dimensional coordinates of the geophone, the two-way travel time vector, and the slowness vector, the spatiotemporal domain values ​​of each trace in the seismic data are converted into Radon domain values ​​to obtain the forward transformation result data, and the mapping relationship between the spatiotemporal domain values ​​of each trace and the corresponding Radon domain values ​​is recorded. The forward transform result data is filtered by the slowness range that needs to be retained to obtain the filtered result data; Based on the mapping relationship between the spatiotemporal domain values ​​and the corresponding Radon domain values, the filtered result data is inversely transformed to obtain the inverse transform result data, including: SS1. Set the initial cycle number I of the original data channel number to 1; SS2, Let the time vector The initial loop number J = 1; SS3, initialize the page number of the three-dimensional array trecord. =1; SS4. Set the initial amplitude value Amp = 0; SS5. Initialize the row numbers of the three-dimensional array trecord. =1; SS6, Determine if trecord( Is it equal to J? If yes, proceed to the next step; otherwise, jump to SS8. SS7. Determine the result of the orthogonal transformation (Radon) Is it equal to 0? If so, then set Radon( Add to Amp and replace the value of Amp; otherwise, proceed to the next step. SS8, Judgment Is it less than the maximum row number in trecord? If not, proceed to the next step; if yes, then... = +1, repeat SS6-SS8; SS9. Judgment Is it less than the maximum page number of trecord? If not, proceed to the next step; if yes, then... = +1, repeat SS5-S9; SS10. Save the Amp value to the corresponding location in the inverse transform result. middle; SS11. Determine if J is less than The largest sequence number; if not, proceed to the next step; if yes, then J = J + 1, and repeat SS3-SS11. SS12. Determine if I is less than the maximum number of channels in the original data. If not, proceed to the next step. If yes, then I = I + 1, and repeat SS2-SS12. SS13. Save the inverse transform result data.

2. The Radon filtering method for seismic data according to claim 1, characterized in that, Converting spatiotemporal domain values ​​to Radon domain values ​​includes: By using the three-dimensional coordinates of the seismic source, the three-dimensional coordinates of the geophone, the two-way travel time vector, and the slowness vector, the two-way travel time corresponding to each trace in the seismic data at different time nodes and different slowness nodes is calculated. Based on the corresponding two-way travel time, the amplitude values ​​corresponding to each trace are extracted from the seismic data and summed to obtain the Radon domain values ​​of each trace at different time nodes and different slowness nodes.

3. The Radon filtering method for seismic data according to claim 2, characterized in that, Calculating the round trip includes: ； ； ； , , The three-dimensional coordinates of the earthquake source are: , , Let t be the three-dimensional coordinates of the detector, t be the two-way travel time vector node, p be the slowness vector node, and T be the two-way travel time of the target track at that time node and the slowness node.

4. The Radon filtering method for seismic data according to claim 1, characterized in that, The forward transform result data is filtered by retaining a specific slowness range, including: Generate a filter based on the required range of slowness. The forward transform result data is filtered using the filter to obtain the filtered result data.

5. The Radon filtering method for seismic data according to claim 4, characterized in that, Generate filters based on the required retention of the slowness range, including: Determine the vector position of the slowness range that needs to be retained within the slowness vector; The initial filter is adjusted based on the vector position and the vector length during two-way travel to obtain the adjusted filter.

6. A Radon filtering system for seismic data, characterized in that, include: The acquisition module is configured to acquire two-way travel time vectors, slowness vectors, as well as seismic data and the corresponding three-dimensional coordinates of the source and the detector. The forward transformation module is configured to convert the spatiotemporal domain values ​​of each trace in the seismic data into Radon domain values ​​based on the three-dimensional coordinates of the seismic source, the three-dimensional coordinates of the detector, the two-way travel time vector, and the slowness vector, to obtain the forward transformation result data, and record the mapping relationship between the spatiotemporal domain values ​​of each trace and the corresponding Radon domain values. The filtering module is configured to filter the forward transform result data through the slowness range that needs to be retained, and obtain the filtered result data. The inverse transform module is configured to perform an inverse transform on the filtered result data based on the mapping relationship between the spatiotemporal domain values ​​and the corresponding Radon domain values ​​to obtain inverse transform result data, including: SS1. Set the initial cycle number I of the original data channel number to 1; SS2, Let the time vector The initial loop number J = 1; SS3, initialize the page number of the three-dimensional array trecord. =1; SS4. Set the initial amplitude value Amp = 0; SS5. Initialize the row numbers of the three-dimensional array trecord. =1; SS6, Determine if trecord( Is it equal to J? If yes, proceed to the next step; otherwise, jump to SS8. SS7, Determine the filtering result data Radon( Is it equal to 0? If so, then set Radon( Add to Amp and replace the value of Amp; otherwise, proceed to the next step. SS8, Judgment Is it less than the maximum row number in trecord? If not, proceed to the next step; if yes, then... = +1, repeat SS6-SS8; SS9. Judgment Is it less than the maximum page number of trecord? If not, proceed to the next step; if yes, then... = +1, repeat SS5-S9; SS10. Save the Amp value to the corresponding location in the inverse transform result. middle; SS11. Determine if J is less than The largest sequence number; if not, proceed to the next step; if yes, then J = J + 1, and repeat SS3-SS11. SS12. Determine if I is less than the maximum number of channels in the original data. If not, proceed to the next step. If yes, then I = I + 1, and repeat SS2-SS12. SS13. Save the inverse transform result data.

7. The Radon filtering system for seismic data according to claim 6, characterized in that, The forward transformation module includes: The calculation unit is configured to calculate the two-way travel time of each trace in the seismic data at different time nodes and different slowness nodes using the three-dimensional coordinates of the source, the three-dimensional coordinates of the detector, the two-way travel time vector, and the slowness vector. The summation unit is configured to extract the amplitude values ​​corresponding to each trace from the seismic data based on the corresponding two-way travel time and sum them to obtain the Radon domain values ​​of each trace at different time nodes and different slowness nodes.

8. The Radon filtering system for seismic data according to claim 7, characterized in that, The calculation unit calculates the round-trip journey by including: ； ； ； , , The three-dimensional coordinates of the earthquake source are: , , Let t be the three-dimensional coordinates of the detector, t be the two-way travel time vector node, p be the slowness vector node, and T be the two-way travel time of the target track at that time node and the slowness node.

9. The Radon filtering system for seismic data according to claim 6, characterized in that, The filtering module includes: The generation unit is configured to generate filters for the slowness range that needs to be retained. The filtering unit is configured to filter the forward transform result data through the filter to obtain filtered result data.

10. The Radon filtering system for seismic data according to claim 9, characterized in that, The generation unit includes: The vector position determination unit is configured to determine the vector position of the slowness range to be retained within the slowness vector; The filter adjustment unit is configured to adjust the initial filter based on the vector position and the vector length of the two-way travel time vector to obtain the adjusted filter.

Images