Controllable source mixed wave field separation method and device, computer equipment and storage medium

Abstract

The application provides a controllable source mixed wave field separation method, comprising the following steps: obtaining a Doppler effect correction operator and a deconvolution operator; based on the Doppler effect correction operator and the deconvolution operator, performing first separation processing on original data to obtain data domain signal data; a curved wave domain operator converts the data domain signal data to curved wave domain signal data; the sparse terms of the curved wave domain signal data are determined; based on the Doppler effect correction operator, the deconvolution operator and the curved wave domain operator, a function formula about the sparse terms of the curved wave domain signal data and the original data is constructed based on the curved wave domain; sparse inversion is performed on the target function, the sparse terms of the original data are used, and based on the Doppler effect correction operator and the deconvolution operator, second separation processing is performed on the original data to obtain mixed wave field separation data. The controllable source mixed wave field separation method simultaneously realizes Doppler moveout correction and deconvolution, and ensures the separation effect of the controllable source mixed wave field.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of seismic data processing, and particularly relates to a controllable source mixed wave field separation method and device, computer equipment and a storage medium. BACKGROUND

[0002] In marine geophysical exploration, compared with the traditional air gun source, the marine controllable source is a new type of source that is continuously excited by designing different scanning signals, which not only meets the current economic and environmental requirements of marine exploration, but also can obtain marine controllable source mixed wave field data by phase coding combination of the mixed source, so as to realize lower frequency and more efficient acquisition of offshore seismic data. However, the addition of the mixed source will also bring a series of challenges to the processing of marine seismic data. As a continuous long-time excitation source, the marine controllable source can be excited by different scanning signals at the same time through design, and the different scanning signals overlap with each other in the excitation process, producing a series of aliasing noises, so that the seismic wave field excited by the marine mixed source is more complex than that of the air gun source. In addition, due to the change of the position between the marine controllable source and the receiver during the acquisition process, there is a time error between different sources and receivers due to the Doppler effect, which also makes it difficult to separate the wave field of the mixed source.

[0003] Therefore, to realize the marine controllable source mixed wave field separation, not only the time error between the moving source and the receiver due to the Doppler effect needs to be solved, but also how to use deconvolution to convert the source scanning signal into a pulse signal needs to be considered. Therefore, a new marine controllable source mixed wave field separation method is urgently needed to solve the scanning signal excitation problem of the marine controllable source while performing Doppler correction, so as to realize the marine controllable source mixed wave field separation. SUMMARY

[0004] Therefore, it is necessary to provide a controllable source mixed wave field separation method aiming at the above technical problems, comprising:

[0005] obtaining original data, obtaining a Doppler effect correction operator and a deconvolution operator, performing first separation processing on the original data based on the Doppler effect correction operator and the deconvolution operator to obtain scanning signal data based on a data domain, and determining the data domain signal data;

[0006] obtaining a curved wave domain operator, converting the data domain signal data to scanning signal data based on a curved wave domain by using the curved wave domain operator, and determining the curved wave domain signal data;

[0007] increasing the sparsity of the curved wave domain signal data to determine the sparse terms of the curved wave domain signal data;

[0008] Based on the Doppler effect correction operator, the deconvolution operator and the curvelet domain operator, a function formula about the sparse term of the curvelet domain signal data and the original data is constructed based on the curvelet domain, and a target function is determined;

[0009] The target function is sparsely inverted to determine the sparse term of the original data, and the original data is processed for a second time by using the sparse term of the original data and based on the Doppler effect correction operator and the deconvolution operator, to obtain mixed wave field separation data.

[0010] In one embodiment, the step of obtaining original data, obtaining a Doppler effect correction operator and a deconvolution operator, processing the original data for a first time based on the Doppler effect correction operator and the deconvolution operator to obtain data domain based scanning signal data, and determining the data domain signal data comprises:

[0011] Obtaining original data;

[0012] Obtaining a Doppler effect correction operator and a deconvolution operator;

[0013] Using the Doppler effect correction operator and the deconvolution operator, a data processing calculation formula is constructed;

[0014] Processing the original data for a first time based on the data processing calculation formula to obtain data domain based scanning signal data, and determining the data domain signal data.

[0015] In one embodiment, the data processing calculation formula is:

[0016] d = L d L s d0

[0017] In the formula, d is the original data, L d is the deconvolution operator, L s is the Doppler effect correction operator, and d0 is the data domain signal data.

[0018] In one embodiment, the step of increasing the sparsity of the curvelet domain signal data to determine the sparse term of the curvelet domain signal data comprises:

[0019] Obtaining a sparse processing calculation formula;

[0020] Based on the sparse processing calculation formula, fine grid interpolation is performed on missing seismic traces in the curvelet domain signal data to increase the sparsity of the curvelet domain signal data, so as to determine the sparse term of the curvelet domain signal data.

[0021] In one embodiment, the sparse processing calculation formula is:

[0022] d0=L c m

[0023] wherein d0 is the data domain signal data, L c is the curvelet domain operator, and m is the sparse term of the curvelet domain signal data.

[0024] In one embodiment, the step of constructing a function formula based on the curvelet domain about the sparse term of the curvelet domain signal data and the original data as a target function based on the Doppler effect correction operator, the deconvolution operator and the curvelet domain operator comprises:

[0025] constructing a function formula based on the curvelet domain about the sparse term of the curvelet domain signal data and the original data as a target function based on the Doppler effect correction operator, the deconvolution operator and the curvelet domain operator;

[0026] adding a sparse constraint to the original data based on the sparsity of the original data.

[0027] In one embodiment, the target function is:

[0028] J=||d-L d L s L c m|| 2 +μ||m||1

[0029] wherein f is the original data, d is the original data, L d is the deconvolution operator, L s is the Doppler effect correction operator, μ is an inversion regularization operator, m is the sparse term of the curvelet domain signal data, and ||m||1 represents adding a sparse constraint to the original data.

[0030] A controllable source hybrid wave field separation device comprises:

[0031] a first separation processing module configured to obtain original data, obtain a Doppler effect correction operator and a deconvolution operator, perform first separation processing on the original data based on the Doppler effect correction operator and the deconvolution operator, obtain data domain scanning signal data, and determine the data domain scanning signal data as data domain signal data;

[0032] a data conversion module configured to obtain a curvelet domain operator, convert the data domain signal data to curvelet domain scanning signal data by using the curvelet domain operator, and determine the curvelet domain scanning signal data as curvelet domain signal data;

[0033] a sparsity processing module, configured to increase the sparsity of the curvelet domain signal data to determine sparse terms of the curvelet domain signal data;

[0034] a target function construction module, configured to construct a function formula based on the curvelet domain and the sparse terms of the curvelet domain signal data and the original data based on the Doppler effect correction operator, the deconvolution operator and the curvelet domain operator, and determine a target function;

[0035] a second separation processing module, configured to perform sparse inversion on the target function to determine sparse terms of the original data, and perform second separation processing on the original data based on the Doppler effect correction operator and the deconvolution operator by using the sparse terms of the original data, to obtain mixed wave field separation data.

[0036] A computer device, comprising a memory and a processor, the memory stores a computer program, and the processor implements the steps of the method according to any one of the above embodiments when executing the computer program.

[0037] A computer readable storage medium, which stores a computer program, and the computer program implements the steps of the method according to any one of the above embodiments when executed by a processor.

[0038] The application provides a controllable source mixed wave field separation method.

[0039] Compared with the prior art, the controllable source mixed wave field separation method has the following advantages or beneficial effects:

[0040] The controllable source mixed wave field separation method, device, computer device and storage medium have the following advantages: the Doppler effect correction operator solves the time error between the moving source and the detector due to the Doppler effect, the deconvolution operator well suppresses aliasing noise, and the original data can be well processed for the first time, so that the data domain signal data is optimized, that is, the target function is constructed based on the optimized seismic data, so the target function can optimize the inversion of the seismic data, and the target function is constructed based on the Doppler effect correction operator and the deconvolution operator, so that the Doppler time difference correction and deconvolution are realized at the same time when the original data is processed for the second time by using the target function, and the separation effect of the controllable source mixed wave field is ensured. BRIEF DESCRIPTION OF DRAWINGS

[0041] Figure 1 It is a flowchart of the controllable source mixed wave field separation method in one embodiment;

[0042] Figure 2A flowchart of a controllable source mixed wave field separation method in another embodiment;

[0043] Figure 3 A schematic diagram of controllable source wave field data in a comparative example;

[0044] Figure 4 A schematic diagram of original controllable source mixed wave field data in a comparative example;

[0045] Figure 5 A schematic diagram of traditional first controllable source mixed wave field data;

[0046] Figure 6 A schematic diagram of second controllable source mixed wave field data in an embodiment;

[0047] Figure 7 A schematic diagram of original marine controllable source mixed wave field data in a comparative example;

[0048] Figure 8 A schematic diagram of marine controllable source mixed wave field separation data in an embodiment;

[0049] Figure 9 A structural block diagram of a controllable source mixed wave field separation device in an embodiment;

[0050] Figure 10 An internal structural diagram of a computer device in an embodiment. DETAILED DESCRIPTION

[0051] In order to make the purposes, technical solutions and advantages of the present application clearer, the present application is further described in detail below with reference to the drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and not to limit the present application.

[0052] Embodiment One

[0053] In this embodiment, as shown in Figure 1 , a controllable source mixed wave field separation method is provided, comprising:

[0054] S110, obtaining original data, obtaining a Doppler effect correction operator and a deconvolution operator, performing a first separation processing on the original data based on the Doppler effect correction operator and the deconvolution operator, obtaining data domain based scanning signal data, and determining the data domain signal data.

[0055] In this embodiment, the Doppler effect correction operator solves the time error caused by the Doppler effect between the moving source and the receiver, the deconvolution operator well suppresses the aliasing noise, and can perform a better first separation processing on the original data, so that the data domain signal data is optimized.

[0056] S120, obtaining a curvelet domain operator, converting the data domain signal data to curvelet domain signal data based on the curvelet domain operator, and determining the curvelet domain signal data.

[0057] In the embodiment, the scanning signal data is converted from the data domain to the curvelet domain, which can remove the noise and interference in the seismic signal and ensure the separation effect of the controllable source mixed wave field.

[0058] S130, increasing the sparsity of the curvelet domain signal data to determine the sparse term of the curvelet domain signal data.

[0059] In the embodiment, increasing the sparsity of the curvelet domain signal data can highlight the sparse term, making it easier to be observed and identified.

[0060] S140, based on the Doppler effect correction operator, the deconvolution operator and the curvelet domain operator, constructing a function formula of the sparse term of the curvelet domain signal data and the original data based on the curvelet domain, and determining a target function.

[0061] In the embodiment, the wave field separation problem of the controllable source is converted into a least square solution problem of the target function.

[0062] Further, the least square method is used to construct the target function, so that the sparse term of the curvelet domain signal data and the original data are well fitted on the fitting curve of the target function, and then the optimal inversion of the original data is realized in the subsequent second separation processing.

[0063] S150, performing sparse inversion on the target function to determine the sparse term of the original data, using the sparse term of the original data, and based on the Doppler effect correction operator and the deconvolution operator, performing second separation processing on the original data to obtain mixed wave field separation data.

[0064] In the embodiment, the target function is constructed by using the least square method, and the target function is constructed based on the Doppler effect correction operator and the deconvolution operator, which simultaneously realizes the Doppler moveout correction and the deconvolution, and realizes the separation and solution of the mixed wave field by using the multi-dimensional inversion process. Thus, the optimal inversion of the original data is realized by using the target function, which ensures the separation effect of the controllable source mixed wave field and improves the resolution of the marine controllable source seismic data.

[0065] The controllable source mixed wave field separation method described above is based on pure data driving and does not require prior information such as known medium structure and velocity model, and is easy to implement. Compared with the traditional mixed wave field separation method based on deconvolution, the method can effectively suppress the aliasing noise generated by the mixed source.

[0066] The controllable source mixed wave field separation method can obtain the mixed wave field separated seismic data through iterative inversion, and the fidelity of the mixed wave field separation of the marine controllable source is improved in the process of iterative inversion.

[0067] The controllable source mixed wave field separation method, the Doppler effect correction operator solves the time error between the moving source and the receiver due to the Doppler effect, the deconvolution operator well suppresses the aliasing noise, and can well perform the first separation processing on the original data, so that the data domain signal data is optimized, that is, the objective function is constructed on the basis of the optimized seismic data, so that the objective function can perform optimized inversion on the seismic data, and the objective function is constructed on the basis of the Doppler effect correction operator and the deconvolution operator, so that the Doppler moveout correction and the deconvolution are simultaneously realized when the original data is processed for the second time by using the objective function, and the separation effect of the controllable source mixed wave field is ensured.

[0068] In one embodiment, the step of obtaining the original data, obtaining the Doppler effect correction operator and the deconvolution operator, performing the first separation processing on the original data on the basis of the Doppler effect correction operator and the deconvolution operator, and obtaining the data domain based scanning signal data and determining the data domain signal data comprises:

[0069] S111, obtaining the original data;

[0070] S112, obtaining the Doppler effect correction operator and the deconvolution operator;

[0071] S113, constructing a data processing calculation formula by using the Doppler effect correction operator and the deconvolution operator;

[0072] S114, performing the first separation processing on the original data on the basis of the data processing calculation formula, obtaining the data domain based scanning signal data, and determining the data domain signal data.

[0073] In this embodiment, the data processing calculation formula is constructed by using the Doppler effect correction operator and the deconvolution operator, so that the Doppler moveout correction and the deconvolution are simultaneously realized when the original data is processed for the first time by using the data processing calculation formula, not only the marine controllable source mixed wave field is well separated, but also the Doppler effect and the aliasing noise are significantly suppressed.

[0074] In one embodiment, the data processing calculation formula is:

[0075] d=L d L s d0

[0076] wherein d is the original data, L d is the deconvolution operator, L s is the Doppler effect correction operator, d0 is the data domain signal data.

[0077] In the embodiment, when the original data is processed for the first time by using the data processing calculation formula, the Doppler time difference correction and the deconvolution are simultaneously realized, the mixed wave field of the marine vibroseis is well separated, and the Doppler effect and aliasing noise are significantly suppressed.

[0078] In one embodiment, the step of increasing the sparsity of the curved wave domain signal data to determine the sparse term of the curved wave domain signal data comprises:

[0079] S131, obtaining a sparse processing calculation formula;

[0080] S132, based on the sparse processing calculation formula, performing fine grid interpolation on the missing seismic trace in the curved wave domain signal data, increasing the sparsity of the curved wave domain signal data, and determining the sparse term of the curved wave domain signal data.

[0081] In the embodiment, the fine grid interpolation on the missing seismic trace in the curved wave domain signal data can increase the sparsity of the curved wave domain signal data, so as to highlight the sparse term and easily determine the sparse term, and make the curved wave domain signal data relatively dense and complete, thereby improving the fidelity of the mixed wave field separation of the marine vibroseis in the iterative inversion process.

[0082] In one embodiment, the sparse processing calculation formula is:

[0083] d0=L c m

[0084] wherein d0 is the data domain signal data, L c is the curved wave domain operator, and m is the sparse term of the curved wave domain signal data.

[0085] In the embodiment, the curved wave domain operator is composed of curved wave domain coefficients.

[0086] In the embodiment, the fine grid interpolation between the observation points of the curved wave domain signal data can improve the effect and accuracy of the fine grid interpolation, increase the sparsity of the curved wave domain signal data, and supplement the data of the missing seismic trace in the curved wave domain signal data, so that the curved wave domain signal data is relatively dense and complete, and the objective function obtained by the curved wave domain signal data has a good fitting effect with the actual seismic data, thereby improving the fidelity of the mixed wave field separation of the marine vibroseis in the iterative inversion process.

[0087] In one embodiment, the function formula of the sparse term of the curved wave field signal data based on the curved wave field and the original data is constructed based on the Doppler effect correction operator, the deconvolution operator and the curved wave field operator, and the step of determining the target function comprises:

[0088] S141, the function formula of the sparse term of the curved wave field signal data based on the curved wave field and the original data is constructed based on the Doppler effect correction operator, the deconvolution operator and the curved wave field operator, and the target function is determined.

[0089] S142, the sparse constraint is added to the original data based on the sparsity of the original data.

[0090] In this embodiment, the sparse constraint is added to the original data, which can optimize the original data and ensure the effect of sparse inversion, and improve the fidelity of the mixed wave field separation of the marine vibroseis in the process of iterative inversion.

[0091] In one embodiment, the target function is: J = ||d-L d L s L c m|| 2 + μ||m||1, wherein d is the original data, d is the original data, L d is the deconvolution operator, L s is the Doppler effect correction operator, μ is the inversion regularization operator, m is the sparse term of the curved wave field signal data, and ||m||1 represents the sparse constraint added to the original data.

[0092] In this embodiment, the sparse term and the sparse constraint optimize the original data at the same time, ensure the effect of sparse inversion, and the inversion regularization operator is obtained by using the least square method, which realizes the optimal inversion of the original data in the second separation processing, and further improves the fidelity of the mixed wave field separation of the marine vibroseis in the process of iterative inversion.

[0093] In one embodiment, the function formula of the sparse term of the curved wave field signal data based on the curved wave field and the original data is constructed based on the Doppler effect correction operator, the deconvolution operator and the curved wave field operator, and the target function is determined; the sparse inversion is performed on the target function to determine the sparse term of the original data, the sparse term of the original data is used, and the second separation processing is performed on the original data based on the Doppler effect correction operator and the deconvolution operator to obtain the mixed wave field separation data, which is determined as iterative inversion. The number of iterative inversions is at least twice, and the mixed wave field separation data obtained in the previous iterative inversion is taken as the original data for each iterative inversion from the second time.

[0094] In the embodiment, the at least two times of iterative inversion can constantly optimize the interpretation result of the seismic data, and ensure the fidelity of the vibroseis mixed wave field separation method to the mixed wave field separation.

[0095] It should be understood that, although Figure 1 The steps in the flowchart of the method can be displayed in sequence according to the arrows, but these steps are not necessarily executed in the order indicated by the arrows. Unless otherwise specified herein, the execution of these steps is not strictly limited in sequence, and these steps can be executed in other orders. Moreover, Figure 1 At least part of the steps in the method can include multiple sub-steps or multiple stages, which are not necessarily executed at the same time, but can be executed at different times, and the execution order of these sub-steps or stages is not necessarily sequential, but can be executed in rotation or alternation with other steps or sub-steps or stages of other steps.

[0096] Embodiment two

[0097] In the embodiment, as shown in Figure 2 A vibroseis mixed wave field separation method is provided, comprising:

[0098] Step one, input data;

[0099] Step two, based on the data, construct a target function;

[0100] Step three, use the target function to solve the optimal solution of the data by least squares;

[0101] Step four, perform iterative inversion through the optimal solution of the data to obtain optimized seismic data;

[0102] Step five, the number of times of iterative inversion is multiple, and each iteration uses the optimized seismic data obtained by the previous iteration as data.

[0103] In the embodiment, multiple times of iterative inversion can constantly optimize the interpretation result of the seismic data, and ensure the fidelity of the vibroseis mixed wave field separation method to the mixed wave field separation.

[0104] Embodiment three

[0105] In the embodiment, a method for verifying the fidelity of vibroseis mixed wave field separation is provided, comprising:

[0106] Set the data of three source ships;

[0107] The controllable source wave field data obtained without considering the mixed source and the Doppler effect is shown in Figure 3 .

[0108] The original controllable source mixed wave field data obtained by considering the mixed source among the three source ships is shown in Figure 4 .

[0109] The first controllable source mixed wave field data is obtained by using the traditional deconvolution method to perform the controllable source mixed wave field separation processing on the original controllable source mixed wave field data, as shown in Figure 5 , and the first controllable source mixed wave field data is compared with the controllable source wave field data and the original controllable source mixed wave field data respectively.

[0110] The second controllable source mixed wave field data is obtained by using the above controllable source mixed wave field separation method to perform the controllable source mixed wave field separation processing on the original controllable source mixed wave field data, as shown in Figure 6 , and the second controllable source mixed wave field data is compared with the controllable source wave field data, the original controllable source mixed wave field data and the first controllable source mixed wave field data respectively.

[0111] In the embodiment, the ship speed of each source ship is 4.8 m / s, the shot point is set to 80, the number of geophones is set to 15, the distance between the source and the geophone is set to 24 m, and the maximum frequency of the source and the geophone is set to 25 Hz.

[0112] As shown in Figure 3 , Sx(m) represents the amplitude value of the seismic wave, Rx(m) represents the amplitude value of the anti-phase axis of the seismic wave, Figures 4 to 8 , and the same applies hereinafter.

[0113] In the embodiment, the original controllable source mixed wave field data is different from the conventional seismic trace data and is formed by linear scanning signal correlation, and needs to be subjected to Doppler correction and deconvolution processing in the processing process. It can be seen that, compared with the first controllable source mixed wave field data obtained by the traditional deconvolution method, the second controllable source mixed wave field data obtained by using the above controllable source mixed wave field separation method can well suppress the aliasing noise.

[0114] In combination with the controllable source wave field data and the original controllable source mixed wave field data, it can be known that, without considering the mixed source and the Doppler effect, the source has a great influence on the wave field separation effect, the traditional deconvolution method has a poor elimination effect on the interference of the influence, and the second controllable source mixed wave field data obtained by using the above controllable source mixed wave field separation method can well eliminate the interference of the influence, that is, can well suppress the aliasing noise.

[0115] Embodiment Four

[0116] The marine actual data of a work area is selected to further test the controllable source mixed wave field separation method:

[0117] As shown in Figure 7 and Figure 8 , 80 mixed shots are set, the number of geophones is 150, the distance between the shot points and the geophones is 12.5 m, and the speed of the source boat is 1.8 m / s.

[0118] As can be seen from Figure 7 , the aliasing noise of the marine controllable source mixed wave field original data is large, and as shown in Figure 8 , the aliasing noise of the marine controllable source mixed wave field separation data obtained after the above controllable source mixed wave field separation method is small. It can be seen that after the above controllable source mixed wave field separation method is separated and processed, not only the marine controllable source mixed wave field is well separated, but also the Doppler effect and aliasing noise are significantly suppressed, proving the effectiveness of the above controllable source mixed wave field separation method.

[0119] Embodiment Five

[0120] In this embodiment, as shown in Figure 9 , a controllable source mixed wave field separation device is provided, comprising:

[0121] The first separation processing module 210 is configured to obtain original data, obtain a Doppler effect correction operator and a deconvolution operator, perform first separation processing on the original data based on the Doppler effect correction operator and the deconvolution operator, obtain data domain-based scanning signal data, and determine the data domain-based scanning signal data as data domain signal data;

[0122] The data conversion module 220 is configured to obtain a curved wave domain operator, convert the data domain signal data to curved wave domain-based scanning signal data by using the curved wave domain operator, and determine the curved wave domain-based scanning signal data as curved wave domain signal data;

[0123] The sparse processing module 230 is configured to increase the sparsity of the curved wave domain signal data to determine sparse terms of the curved wave domain signal data;

[0124] The objective function construction module 240 is configured to construct a function formula about the sparse terms of the curved wave domain signal data and the original data based on the Doppler effect correction operator, the deconvolution operator, and the curved wave domain operator, and determine the function formula as an objective function;

[0125] The second separation processing module 250 is configured to perform sparse inversion on the target function, determine sparse terms of the original data, perform second separation processing on the original data based on the Doppler effect correction operator and the deconvolution operator by using the sparse terms of the original data, and obtain mixed wave field separation data.

[0126] In the embodiment, the first separation processing module 210 not only solves the time error between the moving source and the receiver due to the Doppler effect, but also suppresses the aliasing noise well, and can perform better first separation processing on the original data, so that the data domain signal data is optimized.

[0127] In the embodiment, the data conversion module 220 converts the scanning signal data from the data domain to the curved wave domain, which can remove the noise and interference in the seismic signal and ensure the separation effect of the controllable source mixed wave field separation device on the controllable source mixed wave field.

[0128] In the embodiment, the sparsity of the curved wave domain signal data is increased, which can highlight the sparse terms and make the sparse terms more easily observed and identified.

[0129] In the embodiment, the wave field separation problem of the controllable source is converted into a least square solution problem of the target function.

[0130] Further, the least square method is used to construct the target function, so that the sparse terms of the curved wave domain signal data are well fitted on the fitting curve of the target function with the original data, and then the optimal inversion of the original data by the subsequent second separation processing is realized.

[0131] In the embodiment, the target function is constructed by using the least square method, and the first separation processing module simultaneously realizes the Doppler time difference correction and the deconvolution, and the separation and solution of the mixed wave field are realized by using the multi-dimensional inversion process. Therefore, the optimal inversion of the original data is realized by using the target function, the separation effect of the controllable source mixed wave field is ensured, and the resolution of the marine controllable source seismic data is improved.

[0132] The controllable source mixed wave field separation device has the advantages that the first separation processing module not only solves the time error between the moving source and the receiver due to the Doppler effect, but also suppresses the aliasing noise well, and can perform better first separation processing on the original data, so that the data domain signal data is optimized. The target function is constructed based on the optimized seismic data, so that the target function can perform optimal inversion on the seismic data. The target function is constructed based on the Doppler effect correction operator and the deconvolution operator, so that the Doppler time difference correction and the deconvolution are simultaneously realized when the original data is processed by using the target function for the second separation processing, and the separation effect of the controllable source mixed wave field is ensured.

[0133] In one embodiment, the first separation processing module comprises:

[0134] an original data acquisition unit configured to acquire original data;

[0135] an operator acquisition unit configured to acquire a Doppler effect correction operator and a deconvolution operator;

[0136] a data construction unit configured to construct a data processing formula by using the Doppler effect correction operator and the deconvolution operator;

[0137] a data processing unit configured to perform first separation processing on the original data based on the data processing formula to obtain data domain scanning signal data, which is determined as data domain signal data.

[0138] In this embodiment, the data processing formula is constructed by using the Doppler effect correction operator and the deconvolution operator, so that Doppler time difference correction and deconvolution are simultaneously achieved when the data processing unit performs first separation processing on the original data, and not only the marine vibrator mixed wave field is well separated, but also the Doppler effect and aliasing noise are significantly suppressed.

[0139] In one embodiment, the data processing formula is:

[0140] d = L d L s d0

[0141] In the formula, d is the original data, L d is the deconvolution operator, L s is the Doppler effect correction operator, and d0 is the data domain signal data.

[0142] In this embodiment, Doppler time difference correction and deconvolution are simultaneously achieved when the data processing formula is used to perform first separation processing on the original data, and not only the marine vibrator mixed wave field is well separated, but also the Doppler effect and aliasing noise are significantly suppressed.

[0143] In one embodiment, the vibrator mixed wave field separation device further comprises a sparsity processing module, and the sparsity processing module comprises:

[0144] a sparsity formula acquisition unit configured to acquire a sparsity processing formula;

[0145] a fine grid interpolation unit configured to perform fine grid interpolation on missing seismic traces in the curved wave domain signal data based on the sparsity processing formula.

[0146] In the embodiment, after the sparse calculation formula acquisition unit acquires the sparse processing calculation formula, the fine grid interpolation unit performs fine grid interpolation on the missing seismic trace in the curvelet domain signal data, which can increase the sparsity of the curvelet domain signal data, so as to facilitate subsequent highlighting of the sparse term, i.e., easy determination of the sparse term, and make the curvelet domain signal data more dense and complete, thereby improving the fidelity of the marine vibrator mixed wave field separation in the iterative inversion process.

[0147] In one embodiment, the sparse processing calculation formula is:

[0148] d0=L c m

[0149] In the formula, d0 is the data domain signal data, L c is the curvelet domain operator, and m is the sparse term of the curvelet domain signal data.

[0150] In the embodiment, the curvelet domain operator is composed of curvelet domain coefficients.

[0151] In the embodiment, fine grid interpolation is performed between the observation points of the curvelet domain signal data, which can improve the effect and accuracy of fine grid interpolation, increase the sparsity of the curvelet domain signal data, and supplement the data of the missing seismic trace in the curvelet domain signal data, so that the curvelet domain signal data is more dense and complete, thereby improving the fidelity of the marine vibrator mixed wave field separation in the iterative inversion process.

[0152] In one embodiment, the vibrator mixed wave field separation device further comprises:

[0153] A sparse constraint module is configured to add a sparse constraint to the original data based on the sparsity of the original data.

[0154] In the embodiment, the sparse constraint module adds a sparse constraint to the original data, which can optimize the original data and ensure the effect of sparse inversion, thereby improving the fidelity of the marine vibrator mixed wave field separation in the iterative inversion process.

[0155] In one embodiment, the target function is:

[0156] J=||d-L d L s L c m|| 2 +μ||m||1

[0157] In the formula, d is the original data, d is the original data, Ld L is the deconvolution operator, s μ is the inversion regularization operator, m is the sparse term of the curvelet domain signal data, and ||m||1 represents adding a sparse constraint to the original data.

[0158] In this embodiment, the sparse term and the sparse constraint optimize the original data at the same time, ensure the effect of sparse inversion, and the inversion regularization operator is obtained by using the least square method, realizes the optimal inversion of the original data by the second separation processing, and further improves the fidelity of the marine vibroseis mixed wave field separation in the process of iterative inversion.

[0159] The specific limitations of the marine vibroseis mixed wave field separation device can be referred to the limitations of the marine vibroseis mixed wave field separation method in the above, which will not be repeated here. Each unit in the above marine vibroseis mixed wave field separation device can be realized by software, hardware and combinations thereof, in whole or in part. The above each unit can be embedded in or independent of the processor in the computer device in hardware form, or can be stored in the memory in the computer device in software form, so as to be called and executed by the processor to perform the operations corresponding to each unit.

[0160] Embodiment six

[0161] In this embodiment, as shown in Figure 10 , a computer device is provided. Its internal structure diagram can be as shown in Figure 10 . The computer device includes a processor, a memory, a network interface, a display screen and an input device connected by a system bus. Among them, the processor of the computer device is used to provide computing and control capability. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program, and the non-volatile storage medium is deployed with a database, which is used to store original data. The internal memory provides an environment for the operating system and the computer program in the non-volatile storage medium to run. The network interface of the computer device is used to communicate with other computer devices deployed with application software. The computer program is executed by the processor to implement a marine vibroseis mixed wave field separation method. The display screen of the computer device can be a liquid crystal display screen or an electronic ink display screen, and the input device of the computer device can be a touch layer overlaid on the display screen, or a key, trackball or touchpad arranged on the shell of the computer device, or an external keyboard, touchpad or mouse, etc.

[0162] Those skilled in the art can understand, Figure 10The structure shown in the figure is only a block diagram of part of the structure related to the scheme of the present application, and does not constitute a limitation on the computer device to which the scheme of the present application is applied. The specific computer device can include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components.

[0163] In one embodiment, a computer device is provided, comprising a memory and a processor, the memory storing a computer program, and the processor implementing the following steps when executing the computer program:

[0164] S110, obtaining original data, obtaining a Doppler effect correction operator and a deconvolution operator, performing first separation processing on the original data based on the Doppler effect correction operator and the deconvolution operator to obtain data domain-based scanning signal data, and determining the data domain signal data.

[0165] In this embodiment, the Doppler effect correction operator solves the time error between the moving seismic source and the geophone due to the Doppler effect, and the deconvolution operator well suppresses the aliasing noise, so that the original data can be well separated for the first time, and the data domain signal data is optimized.

[0166] S120, obtaining a curvelet domain operator, converting the data domain signal data to curvelet domain-based scanning signal data using the curvelet domain operator, and determining the curvelet domain signal data.

[0167] In this embodiment, the scanning signal data is converted from the data domain to the curvelet domain, which can remove noise and interference in the seismic signal and ensure the separation effect of the computer device on the controllable seismic source mixed wave field.

[0168] S130, increasing the sparsity of the curvelet domain signal data to determine the sparse terms of the curvelet domain signal data.

[0169] In this embodiment, increasing the sparsity of the curvelet domain signal data can highlight the sparse terms, so that the sparse terms can be more easily observed and identified.

[0170] S140, based on the Doppler effect correction operator, the deconvolution operator and the curvelet domain operator, constructing a function formula of the sparse terms of the curvelet domain signal data and the original data based on the curvelet domain, and determining the target function.

[0171] In this embodiment, the wave field separation problem of the controllable seismic source is converted into a least square solution problem of the target function.

[0172] Further, a target function is constructed by using a least square method, so that the sparse term of the curvelet domain signal data is well fitted on a fitting curve of the target function with the original data, and then the optimal inversion of the original data is realized by subsequent second separation processing.

[0173] S150, sparse inversion is performed on the target function to determine the sparse term of the original data, and the original data is subjected to second separation processing based on the Doppler effect correction operator and the deconvolution operator by using the sparse term of the original data, to obtain mixed wave field separation data.

[0174] In the embodiment, the target function is constructed by using a least square method, and the target function is constructed based on the Doppler effect correction operator and the deconvolution operator, so that the Doppler moveout correction and the deconvolution are simultaneously realized, the separation and solution of the mixed wave field are realized by using a multi-dimensional inversion process, thus the optimal inversion of the original data is realized by using the target function, and the separation effect of the mixed wave field of the marine vibroseis is ensured, and the resolution of the marine vibroseis seismic data is improved.

[0175] In one embodiment, the processor, when executing the computer program, further implements the following steps:

[0176] S111, obtaining original data;

[0177] S112, obtaining a Doppler effect correction operator and a deconvolution operator;

[0178] S113, constructing a data processing calculation formula by using the Doppler effect correction operator and the deconvolution operator;

[0179] S114, performing first separation processing on the original data based on the data processing calculation formula to obtain scanning signal data based on a data domain, and determining the scanning signal data as data domain signal data.

[0180] In one embodiment, the data processing calculation formula is:

[0181] d = L d L s d0

[0182] In the formula, d is the original data, L d is the deconvolution operator, L s is the Doppler effect correction operator, and d0 is the data domain signal data.

[0183] In one embodiment, the processor, when executing the computer program, further implements the following steps:

[0184] S131, obtaining a sparse processing calculation formula;

[0185] S132, performing fine grid interpolation on the missing seismic trace in the curvelet domain signal data based on the sparse processing calculation formula, increasing the sparsity of the curvelet domain signal data, to determine the sparse item of the curvelet domain signal data.

[0186] In one embodiment, the sparse processing calculation formula is:

[0187] d0=L c m

[0188] In the formula, k0 is the data domain signal data, L c is the curvelet domain operator, and m is the sparse item of the curvelet domain signal data.

[0189] In one embodiment, when the processor executes the computer program, the following steps are also implemented:

[0190] S141, constructing a function formula based on the curvelet domain about the sparse item of the curvelet domain signal data and the original data based on the Doppler effect correction operator, the deconvolution operator and the curvelet domain operator, to determine a target function;

[0191] S142, adding a sparse constraint to the original data based on the sparsity of the original data.

[0192] In one embodiment, the target function is:

[0193] J=||d-L d L s L c m|| 2 +μ||m||1

[0194] In the formula, d is the original data, d is the original data, L d is the deconvolution operator, L s is the Doppler effect correction operator, μ is an inversion regularization operator, m is the sparse item of the curvelet domain signal data, and ||m||1 represents adding a sparse constraint to the original data.

[0195] Embodiment seven

[0196] In this embodiment, a computer readable storage medium is provided, which stores a computer program. When the computer program is executed by a processor, the following steps are implemented:

[0197] S110, obtaining original data, obtaining a Doppler effect correction operator and a deconvolution operator, performing a first separation processing on the original data based on the Doppler effect correction operator and the deconvolution operator, to obtain data domain scanning signal data, and determining the data domain signal data.

[0198] In the embodiment, the Doppler effect correction operator solves the time error caused by the Doppler effect between the moving source and the detector, the deconvolution operator well suppresses the aliasing noise, and the original data can be well separated for the first time, so that the data domain signal data is optimized.

[0199] In S120, a curvelet domain operator is obtained, the data domain signal data is converted into curvelet domain scanning signal data based on the curvelet domain operator, and the curvelet domain signal data is determined.

[0200] In the embodiment, the scanning signal data is converted from the data domain to the curvelet domain, the noise and interference in the seismic signal can be removed, and the separation effect of the computer program on the controllable source mixed wave field is ensured.

[0201] In S130, the sparsity of the curvelet domain signal data is increased to determine the sparse items of the curvelet domain signal data.

[0202] In the embodiment, the sparsity of the curvelet domain signal data is increased, and the sparse items are highlighted, so that the sparse items are more easily observed and identified.

[0203] In S140, based on the Doppler effect correction operator, the deconvolution operator and the curvelet domain operator, a function formula about the sparse items of the curvelet domain signal data and the original data based on the curvelet domain is constructed, and a target function is determined.

[0204] In the embodiment, the wave field separation problem of the controllable source is converted into a least square solution problem of the target function.

[0205] Further, the least square method is used to construct the target function, so that the sparse items of the curvelet domain signal data and the original data are well fitted on the fitting curve of the target function, and the subsequent second separation processing optimizes the inversion of the original data.

[0206] In S150, the sparse inversion of the target function is performed to determine the sparse items of the original data, the sparse items of the original data are used, and based on the Doppler effect correction operator and the deconvolution operator, the original data is subjected to the second separation processing to obtain the mixed wave field separation data.

[0207] In the embodiment, the target function is constructed by using the least square method, and the target function is constructed based on the Doppler effect correction operator and the deconvolution operator, so that the Doppler time difference correction and the deconvolution are realized, the separation and solution of the mixed wave field are realized by using the multi-dimensional inversion, the target function is used to realize the optimized inversion of the original data, the separation effect of the controllable source mixed wave field is ensured, and the resolution of the marine controllable source seismic data is improved.

[0208] In one embodiment, the computer program, when executed by the processor, further implements the following steps:

[0209] S111, obtaining original data;

[0210] S112, obtaining a Doppler effect correction operator and a deconvolution operator;

[0211] S113, constructing a data processing calculation formula by using the Doppler effect correction operator and the deconvolution operator;

[0212] S114, performing first separation processing on the original data based on the data processing calculation formula to obtain data domain-based scanning signal data, which is determined as data domain signal data.

[0213] In this embodiment, the data processing calculation formula is constructed by using the Doppler effect correction operator and the deconvolution operator. Thus, when the original data is processed by using the data processing calculation formula, Doppler time difference correction and deconvolution are simultaneously implemented. Not only the ocean bottom seismic mixed wave field is well separated, but also the Doppler effect and aliasing noise are significantly suppressed.

[0214] In one embodiment, the data processing calculation formula is:

[0215] d = L d L s d0

[0216] In the formula, d is the original data, L d is the deconvolution operator, L s is the Doppler effect correction operator, and d0 is the data domain signal data.

[0217] In one embodiment, the computer program, when executed by the processor, further implements the following steps:

[0218] S131, obtaining a sparse processing calculation formula;

[0219] S132, performing fine grid interpolation on missing seismic traces in the curved wave domain signal data based on the sparse processing calculation formula to increase the sparsity of the curved wave domain signal data, so as to determine sparse terms of the curved wave domain signal data.

[0220] In this embodiment, fine grid interpolation is performed on the missing seismic traces in the curved wave domain signal data, which can increase the sparsity of the curved wave domain signal data, so as to highlight the sparse terms in the subsequent process, that is, the sparse terms are easily determined, and the curved wave domain signal data is relatively dense and complete, thereby improving the fitting effect of the objective function obtained by the curved wave domain signal data on the actual seismic data, and further improving the fidelity of the ocean bottom seismic mixed wave field separation in the iterative inversion process.

[0221] In one embodiment, the sparse processing calculation formula is:

[0222] d0=L c m

[0223] In the formula, d0 is the data domain signal data, L c is the curvelet domain operator, and m is the sparse term of the curvelet domain signal data.

[0224] In one embodiment, the computer program, when executed by the processor, further implements the following steps:

[0225] S141, based on the Doppler effect correction operator, the deconvolution operator and the curvelet domain operator, constructing a function formula based on the curvelet domain about the sparse term of the curvelet domain signal data and the original data, determining as the objective function;

[0226] S142, based on the sparsity of the original data, adding a sparse constraint to the original data.

[0227] In this embodiment, adding a sparse constraint to the original data can optimize the original data and ensure the effect of sparse inversion, thereby improving the fidelity of mixed wave field separation of the marine vibroseis in the process of iterative inversion.

[0228] In one embodiment, the objective function is:

[0229] J=||d-L d L s L c m|| 2 +μ||m||1

[0230] In the formula, d is the original data, d is the original data, L d is the deconvolution operator, L s is the Doppler effect correction operator, μ is the inversion regularization operator, m is the sparse term of the curvelet domain signal data, and ||m||1 represents adding a sparse constraint to the original data.

[0231] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program. The computer program can be stored in a non-volatile computer readable storage medium, and when the computer program is executed, the processes of the above-mentioned embodiments of the methods can be included. Any reference to memory, storage, databases, or other media in the embodiments provided by the present 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. As an illustration but not limitation, RAM is available in many forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDR SDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.

[0232] The technical features of the above embodiments can be combined in any way. In order to make the description simple, not all possible combinations of the technical features in the above embodiments are described, but as long as the combinations of the technical features do not exist, they should be considered as the scope of the present application.

[0233] The above embodiments only express several implementation manners of the present application, and the description is more specific and detailed, but it should not be understood as a limitation on the scope of the patent. It should be pointed out that for ordinary skilled in the art, without departing from the concept of the present application, some modifications and improvements can be made, which are all within the scope of the present application. Therefore, the scope of the patent of the present application should be subject to the appended claims.

Claims

1. A vibroseis hybrid wavefield separation method, characterized in that: include: Acquiring raw data, acquiring a Doppler effect correction operator and a deconvolution operator, and performing a first separation process on the raw data based on the Doppler effect correction operator and the deconvolution operator to obtain scanning signal data based on a data domain, which is determined as data domain signal data; Obtaining a curvelet domain operator, and using the curvelet domain operator to convert the data domain signal data into scan signal data based on the curvelet domain, and determining the data as curvelet domain signal data; Increasing the sparsity of the curvelet domain signal data to determine a sparse term of the curvelet domain signal data; Based on the Doppler effect correction operator, the deconvolution operator, and the curvelet domain operator, constructing a curvelet domain-based functional expression of a sparse term of the curvelet domain signal data and the original data, and determining the functional expression as a target function; The objective function is subjected to sparse inversion to determine the sparse terms of the original data. The sparse terms of the original data are used to perform a second separation process on the original data based on the Doppler effect correction operator and the deconvolution operator to obtain mixed wavefield separation data.

2. The method according to claim 1, characterized in that The steps of obtaining raw data, obtaining a Doppler effect correction operator and a deconvolution operator, performing a first separation process on the raw data based on the Doppler effect correction operator and the deconvolution operator to obtain scanning signal data based on the data domain, and determining the data domain signal data include: Get the original data; Obtaining Doppler effect correction operators and deconvolution operators; Using the Doppler effect correction operator and the deconvolution operator, a data processing calculation formula is constructed; Based on the data processing calculation formula, the original data is subjected to a first separation process to obtain scanning signal data based on the data domain, which is determined as data domain signal data.

3. The method according to claim 2, characterized in that The data processing calculation formula is: d=L d L s d0 Wherein, d is the original data, L d is the deconvolution operator, L s is the Doppler effect correction operator, and d0 is the data domain signal data.

4. The method according to claim 1, wherein The step of increasing the sparsity of the curvelet domain signal data to determine a sparse term of the curvelet domain signal data comprises: Get the sparse processing calculation formula; Based on the sparse processing calculation formula, fine grid interpolation is performed on the missing seismic traces in the curvelet domain signal data to increase the sparsity of the curvelet domain signal data and determine the sparse items of the curvelet domain signal data.

5. The method according to claim 4, characterized in that The sparse processing calculation formula is: d0=L c m Wherein, d0 is the data domain signal data, L c is the curvelet domain operator, and m is the sparse term of the curvelet domain signal data.

6. The method according to claim 1, characterized in that The step of constructing a function formula of the sparse terms of the curvelet domain signal data and the original data based on the curvelet domain based on the Doppler effect correction operator, the deconvolution operator and the curvelet domain operator, and determining the function as the objective function includes: Based on the Doppler effect correction operator, the deconvolution operator, and the curvelet domain operator, constructing a curvelet domain-based functional expression of a sparse term of the curvelet domain signal data and the original data, and determining the functional expression as a target function; Based on the sparsity of the original data, a sparse constraint is added to the original data.

7. The method according to claim 6, characterized in that The objective function is: J=||dL d L s L c m|| 2 +μ||m||1 Wherein, d is the original data, d is the original data, L d is the deconvolution operator, L s is the Doppler effect correction operator, μ is the inversion regularization operator, m is the sparse term of the curvelet domain signal data, and ||m||1 represents adding a sparse constraint to the original data.

8. A vibroseis hybrid wavefield separation device, characterized in that: include: a first separation processing module, configured to obtain raw data, obtain a Doppler effect correction operator and a deconvolution operator, and perform a first separation processing on the raw data based on the Doppler effect correction operator and the deconvolution operator to obtain scanning signal data based on the data domain, and determine the data domain signal data; a data conversion module, configured to obtain a curvelet domain operator, and convert the data domain signal data into curvelet domain-based scanning signal data using the curvelet domain operator, and determine the data domain signal data as curvelet domain signal data; a sparse processing module, configured to increase the sparsity of the curvelet domain signal data to determine a sparse term of the curvelet domain signal data; an objective function construction module, configured to construct, based on the Doppler effect correction operator, the deconvolution operator, and the curvelet domain operator, a function formula of the sparse terms of the curvelet domain signal data and the original data in the curvelet domain, and determine the function as an objective function; The second separation processing module is used to perform sparse inversion on the objective function to determine the sparse terms of the original data, and perform a second separation processing on the original data based on the sparse terms of the original data and the Doppler effect correction operator and the deconvolution operator to obtain mixed wavefield separation data.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of the method according to any one of claims 1 to 7 are implemented.

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

Images