Land well multi-source mixed seismic acquisition parameter determination method and device

Abstract

The present disclosure relates to a land well gun multi-source aliasing acquisition parameter determination method and device, the method comprising: dividing a target work area into multiple target regions, determining the initial T-D curve and gun coordinate information corresponding to the current target region; grouping the gun coordinate information corresponding to the current target region, and determining the delay time of the secondary gun based on the initial T-D curve; according to the delay time of the secondary gun in each group of gun coordinate information, simulating aliasing processing on the three-dimensional seismic data corresponding to the primary gun and the secondary gun in each group of gun coordinate information, obtaining the simulated aliasing data of the target region, and separating the aliasing data to determine the similarity of the pre-aliasing data and the post-aliasing separation data of each gun in the aliasing data; in the case where the similarity is greater than or equal to a preset threshold, updating the initial T-D curve to obtain the updated T-D curve as the aliasing acquisition parameter of the target region, and the mixed shooting gun satisfying the T-D curve can effectively separate the aliasing interference in the laboratory.

Description

Technical Field

[0001] This disclosure relates to the field of seismic data acquisition and processing technology in petroleum exploration, and in particular to a method and apparatus for determining the parameters of multi-source superimposed acquisition of land well shots. Background Technology

[0002] In onshore controlled-source oil and gas exploration and acquisition, multiple seismic sources are often used to excite simultaneously or with short delays. Long arrays of nodal instruments are used to receive long-recorded data from multiple shot points. This multi-source mixed acquisition method is an effective way to improve acquisition efficiency and reduce acquisition costs for onshore controlled-source systems. The key to implementing this technology is to conduct effective verification of the mixed acquisition parameters for the acquisition area. These parameters can ensure both efficient acquisition in the field and high-precision mixed separation in the laboratory.

[0003] However, the exploration environment in the western mountainous region is complex. Oil and gas exploration mostly uses explosive well shot sources for excitation. However, the differences between each shot record from well shots in different areas are relatively large and the non-uniformity is strong. The controllable source mixed extraction parameter demonstration method cannot meet the needs of mixed well shot extraction in complex mountainous areas. Therefore, it is necessary to study the parameter demonstration method and process of multi-source mixed acquisition of well shot for oil and gas exploration in complex mountainous areas, so as to meet the requirements of efficient acquisition of well shot in the field and effective mixed separation in the laboratory. Summary of the Invention

[0004] To solve the above-mentioned technical problems, or at least partially solve them, embodiments of this disclosure provide a method and apparatus for determining parameters of multi-source superimposed acquisition in land well shots.

[0005] In a first aspect, embodiments of this disclosure provide a method for determining parameters of multi-source superimposed acquisitions from land-based well shots, including:

[0006] The target work area is divided into multiple target regions. For each target region, the initial TD curve and shot coordinate information corresponding to the current target region are determined based on the three-dimensional seismic data of the current target region.

[0007] According to the preset number of overlapping guns and gun spacing, the gun coordinate information corresponding to the current target area is grouped, and based on the initial TD curve, the delay time of the secondary gun is determined according to the distance between the main gun and the secondary gun in each group of gun coordinate information. Each group of gun coordinate information includes the main gun coordinate information and the secondary gun coordinate information.

[0008] Based on the distance between the main gun and the secondary gun and the delay time of the secondary gun in each set of gun coordinate information, simulated aliasing processing is performed on the three-dimensional seismic data corresponding to the main gun coordinate information and the secondary gun coordinate information in each set of gun coordinate information to obtain simulated aliasing data of the target area. The simulated aliasing data of the target area is then separated to determine the similarity between the data before aliasing and the data after aliasing separation for each shot in the simulated aliasing data.

[0009] If the similarity is less than a preset threshold, return to the step of determining the delay time of the secondary gun based on the distance between the main gun and the secondary gun in each set of gun coordinate information, based on the initial TD curve, until the similarity is greater than or equal to the preset threshold.

[0010] If the similarity is greater than or equal to a preset threshold, the initial TD curve is updated based on the distance between the main gun and the secondary gun corresponding to the current simulated aliasing and the delay time of the secondary gun, and the updated TD curve is obtained as the well-shot multi-source aliasing acquisition parameters for the target area.

[0011] In one possible implementation, determining the initial TD curve corresponding to the current target area based on the three-dimensional seismic data of the current target area includes:

[0012] From the historical 3D seismic data of the current target area, select 3D seismic data with an acquisition time exceeding the preset duration and a maximum offset greater than or equal to the preset distance, and use them as the target 3D seismic data of the current target area.

[0013] Select the seismic data from the nearest offset receiver line of the target 3D seismic data, and divide the selected seismic data into two parts according to the positive and negative offsets. Select the data with slower energy decay from the two parts as the data to be processed.

[0014] The data to be processed is divided into multiple time windows with a time window size of 1 second and a distance of 1 km. The time windows are selected by overlapping half of the time windows.

[0015] Calculate the root mean square amplitude value of the data within each time window, use the current time window position as the coordinate of the current root mean square amplitude value, and use the root mean square amplitude value to replace the time window data to display the root mean square amplitude value of each time window, thereby determining the effective reflected energy and background energy.

[0016] The energy attenuation trend curve is picked out in the energy attenuation transition zone between the effective reflected energy and the background energy. The initial TD curve is plotted based on the inflection point of the trend curve. Multiple TD value ranges are set based on the inflection point value of the initial TD curve.

[0017] In one possible implementation, determining the shot coordinate information corresponding to the current target area based on the three-dimensional seismic data of the current target area includes:

[0018] From the 3D seismic data of the current target area, the shot gather data is extracted according to the shot file number. The data of multiple identical receiver points in the shot gather data are overlaid to obtain the overlaid data. The file number, x and y coordinates corresponding to the shot file number in the overlaid data are extracted as shot coordinate information.

[0019] The gun coordinate information is sorted by the gun file number, and duplicate gun coordinate information is deleted to ensure the uniqueness of the gun coordinate information for each shot.

[0020] In one possible implementation, the step of grouping the gun coordinate information corresponding to the current target area according to a preset number of overlapping shots and shot spacing includes:

[0021] The first step is to mark each shot coordinate information that is not grouped as 0;

[0022] The second step is to select the gun coordinate information marked as 0 and ranked first as the main gun of the current group 1, and mark it as 1 with group number 1;

[0023] Third, among the remaining gun coordinate information marked as 0, select the gun coordinate information whose distance from the main gun is within the first TD value range, and use it as the first secondary gun of the main gun, mark it as 1, group number 1, and record the distance value from the main gun; among the remaining gun coordinate information marked as 0, select the gun coordinate information whose distance from the main gun is within the second TD value range, and use it as the second secondary gun of the main gun, mark it as 1, group number 1, and record the distance value from the main gun.

[0024] Fourth step: If the number of mixed guns in the first group is 3, complete the grouping of the first group according to the first to third steps. If the number of mixed guns in the first group is more than 3, select the gun coordinate information with a distance from the main gun within the range of the N TD value from the remaining gun coordinate information marked as 0, and use it as the Nth secondary gun of the main gun, until the grouping of the first group is completed.

[0025] Fifth step: Take the gun coordinate information that is marked as 0 and ranked first as the main gun of the current group 2, and mark it as 1 and group number 2;

[0026] Step 6: From the remaining gun coordinate information marked as 0, select the gun coordinate information whose distance from the main gun is within a preset progressive distance range after the first TD value range, and use it as the first secondary gun of the main gun, marked as 1, group number 2, and record the distance value from the main gun; From the remaining gun coordinate information marked as 0, select the gun coordinate information whose distance from the main gun is within a preset progressive distance range after the second TD value range, and use it as the second secondary gun of the main gun, marked as 1, group number 2, and record the distance value from the main gun.

[0027] Step 7: If the number of overlapping guns in the second group is 3, complete the grouping of the second group according to steps 5 and 6. If the number of overlapping guns in the second group exceeds 3, select the gun coordinate information within the range of the Nth TD value from the remaining gun coordinate information marked as 0, and use it as the Nth secondary gun of the main gun, until the grouping of the second group is completed.

[0028] Step 8: Following steps 6 and 7, complete the selection of gun coordinate information for each remaining group. The distance between each group and the main gun is increased according to the preset progressive distance until the maximum offset distance is reached. Starting from the first TD value range, if no gun coordinate information that meets the first TD value range requirement for the distance from the main gun is found among the remaining gun coordinate information marked as 0, then only the main gun coordinate information is selected for that group, and no secondary gun coordinate information is selected.

[0029] In one possible implementation, determining the delay time of the secondary gun based on the initial TD curve and the distance between the main gun and the secondary gun in each set of gun coordinate information includes:

[0030] Obtain multiple TD value ranges set according to the inflection point value of the initial TD curve. Each TD value range includes: the distance range between the main gun and the secondary gun in each set of gun coordinate information and the corresponding delay time of the secondary gun.

[0031] Based on the distance between the main gun and the secondary gun in each set of gun coordinate information, determine the range of TD values ​​for that distance;

[0032] Based on the TD value range of the distance, the delay time of the secondary gun is randomly obtained from the delay time of the secondary gun corresponding to the distance range, wherein the delay time between secondary guns in the same group is greater than 0.5s.

[0033] In one possible implementation, the step of performing simulated aliasing processing on the three-dimensional seismic data corresponding to the main gun coordinate information and the secondary gun coordinate information in each set of gun coordinate information, based on the distance between the main gun and the secondary gun and the delay time of the secondary gun in each set of gun coordinate information, to obtain simulated aliasing data of the target area, includes:

[0034] The first step is to store the trace data and trace file number that match the main gun file number in each group of shot coordinate information from the 3D seismic data into the current main gun coordinate information.

[0035] The second step is to store the trace data and trace file number that match the secondary shot file number in each group of shot coordinate information from the 3D seismic data into the current secondary shot coordinate information.

[0036] The third step is that if there is no trace data and trace file number matching the main gun file number or the secondary gun file number in the 3D seismic data, then there is no trace data and trace file number in the shot coordinate information of the main gun or the secondary gun.

[0037] Fourth, for each group of gun coordinate information, if at least one of the main gun coordinate information and the secondary gun coordinate information in the current group lacks track data and track file number, then main gun / secondary gun aliasing processing is not performed; if both the main gun coordinate information and the secondary gun coordinate information in the current group match track data and track file number, then the following main gun / secondary gun aliasing processing is performed on the current group:

[0038] Fifth step, extract all sample values ​​of the matching track data in the current group of main guns, take the delay time of the first secondary gun in the current group as the start time of the first secondary gun sample value, and add it with all sample values ​​of the current main gun in a staggered manner. Take the delay time of the second secondary gun in the current group as the start time of the second secondary gun sample value, and add it with the sum of the current main gun and the first secondary gun again in a staggered manner, until all secondary guns in the current group have been added in a staggered manner, thus completing the data aliasing of the current group;

[0039] Step 6: Perform truncation processing on the aliased data of the current group. Obtain the sample point values ​​of the main gun and secondary gun respectively according to the actual sampling time before aliasing. Replace the sample point data value before aliasing with the truncation main gun sample point data value according to the track file number matched by the main gun in the current group. For the first secondary gun, truncate sample point data of the same length according to its start time and delay time, and replace the sample point data value before aliasing with the track file number matched by the first secondary gun in the current group. For the second secondary gun, truncate sample point data of the same length according to its start time and delay time, and replace the sample point data value before aliasing with the track file number matched by the second secondary gun in the current group. Continue until all secondary guns have undergone aliased data truncation, thus completing the aliased data truncation of the current group.

[0040] Step 7: After completing the aliasing and truncation of the current group of data, perform aliasing and truncation on the next group of data according to steps 5 and 6, until all groups have completed aliasing and truncation, and obtain the simulated aliasing data of the target area.

[0041] In one possible implementation, separating the simulated aliased data of the target region includes:

[0042] For the simulated aliasing data of the target area, the non-aliased guns at adjacent positions are searched based on the gun coordinate information of the simulated aliasing data as model guns, and multi-scale curvelet transform is performed on the aliased guns and model guns in the gun domain.

[0043] At each scale, structural pattern recognition and threshold separation are performed on the aliased gun and the model gun. Then, the separated data are transformed into gun domain data through multi-scale curvelet inverse transform to complete the separation of aliased data.

[0044] In one possible implementation, the similarity between the pre-aliasing data and the post-aliasing separation data for each shot in the simulated aliasing data is determined by the following expression:

[0045]

[0046]

[0047] Where SE(x,y) represents the similarity between the pre-aliasing and post-aliasing data of each shot in the simulated aliasing data, x and y are the pre-aliasing and post-aliasing data respectively, and Ax and Ay are the Fourier amplitudes of data x and y respectively. and Let x and y represent the Fourier phases of the data respectively, and N be the number of shots.

[0048] Secondly, embodiments of this disclosure provide a device for determining parameters of multi-source superimposed acquisition data from land-based well shots, comprising:

[0049] The determination module is used to divide the target work area into multiple target regions. For each target region, the initial TD curve and shot coordinate information corresponding to the current target region are determined based on the 3D seismic data of the current target region.

[0050] The grouping module is used to group the gun coordinate information corresponding to the current target area according to the preset number of overlapping guns and gun spacing, and based on the initial TD curve, determine the delay time of the secondary gun according to the distance between the main gun and the secondary gun in each group of gun coordinate information. Each group of gun coordinate information includes the main gun coordinate information and the secondary gun coordinate information.

[0051] The aliasing module is used to simulate aliasing the 3D seismic data corresponding to the main gun coordinate information and the secondary gun coordinate information in each set of gun coordinate information based on the distance between the main gun and the secondary gun and the delay time of the secondary gun in each set of gun coordinate information. This process yields simulated aliased data of the target area, and the simulated aliased data of the target area is separated to determine the similarity between the data before aliasing and the data after aliasing separation for each shot in the simulated aliased data.

[0052] The return module is used to return the steps of determining the delay time of the secondary gun based on the distance between the main gun and the secondary gun in each group of gun coordinate information, based on the initial TD curve, when the similarity is less than a preset threshold, until the similarity is greater than or equal to the preset threshold.

[0053] The update module is used to update the initial TD curve based on the distance between the main gun and the secondary gun and the delay time of the secondary gun corresponding to the current simulated aliasing when the similarity is greater than or equal to a preset threshold. The updated TD curve is then used as the acquisition parameters for well-shot multi-source aliasing in the target area.

[0054] Thirdly, embodiments of this disclosure provide an electronic device, including a processor, a communication interface, a memory, and a communication bus, wherein the processor, the communication interface, and the memory communicate with each other via the communication bus;

[0055] Memory, used to store computer programs;

[0056] The processor, when executing the program stored in the memory, implements the above-mentioned method for determining the multi-source superimposed acquisition parameters of land well shots.

[0057] Fourthly, embodiments of this disclosure provide a computer-readable storage medium having a computer program stored thereon, characterized in that the computer program, when executed by a processor, implements the above-described method for determining the multi-source superimposed acquisition parameters of land well shots.

[0058] Compared with the prior art, the technical solutions provided in this disclosure have at least some or all of the following advantages:

[0059] The method for determining parameters of multi-source cascaded acquisition of land well shots according to the embodiments of this disclosure divides the target work area into multiple target regions. For each target region, the initial TD curve and shot coordinate information corresponding to the current target region are determined based on the 3D seismic data of the current target region. The shot coordinate information corresponding to the current target region is grouped according to a preset number of cascaded shots and shot spacing. Based on the initial TD curve, the delay time of the secondary shot is determined according to the distance between the main shot and the secondary shot in each group of shot coordinate information. Each group of shot coordinate information includes the main shot coordinate information and the secondary shot coordinate information. Based on the distance between the main shot and the secondary shot in each group of shot coordinate information and the delay time of the secondary shot, the 3D seismic data corresponding to the main shot coordinate information and the secondary shot coordinate information in each group of shot coordinate information are subjected to simulated cascading processing to obtain the simulated target region. The system analyzes and separates the simulated aliased data for the target area, determining the similarity between the pre-aliasing data and the post-aliasing data for each shot. If the similarity is less than a preset threshold, it returns to the initial TD curve, determining the delay time of the secondary gun based on the distance between the main gun and the secondary gun in each set of shot coordinate information, until the similarity is greater than or equal to the preset threshold. If the similarity is greater than or equal to the preset threshold, it updates the initial TD curve based on the distance between the main gun and the secondary gun and the delay time of the secondary gun corresponding to the current simulated aliasing, obtaining the updated TD curve as the well-shot multi-source aliasing acquisition parameter for the target area. Shots that meet this TD curve can effectively separate aliasing interference indoors, turning field aliasing data into single-shot data, achieving efficient field aliasing and high-quality indoor separation processing. Attached Figure Description

[0060] The accompanying drawings, which are incorporated in and form a part of this specification, illustrate embodiments consistent with this disclosure and, together with the description, serve to explain the principles of this disclosure.

[0061] To more clearly illustrate the technical solutions in the embodiments of this disclosure or the prior art, the accompanying drawings used in the description of the embodiments or related technologies will be briefly introduced below. Obviously, those skilled in the art can obtain other drawings based on these drawings without creative effort.

[0062] Figure 1 This illustration schematically shows a flowchart of a method for determining multi-source superimposed acquisition parameters for land well shooting according to an embodiment of the present disclosure;

[0063] Figure 2 A schematic diagram of an initial TD curve according to an embodiment of the present disclosure is shown;

[0064] Figure 3 A schematic diagram of three-dimensional seismic data before aliasing according to an embodiment of the present disclosure is shown.

[0065] Figure 4 A schematic diagram of simulated aliasing data according to an embodiment of the present disclosure is shown.

[0066] Figure 5 A schematic diagram illustrating aliasing separation data according to an embodiment of the present disclosure is shown.

[0067] Figure 6 A schematic diagram of the updated TD curve after well-shot combined production according to an embodiment of the present disclosure is shown.

[0068] Figure 7 A schematic diagram illustrates the structural block diagram of a land well shot multi-source superimposed acquisition parameter determination device according to an embodiment of the present disclosure;

[0069] Figure 8 A schematic block diagram of an electronic device according to an embodiment of the present disclosure is shown. Detailed Implementation

[0070] To make the objectives, technical solutions, and advantages of the embodiments of this disclosure clearer, the technical solutions of the embodiments of this disclosure will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of this disclosure. Based on the embodiments of this disclosure, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this disclosure.

[0071] See Figure 1The embodiments of this disclosure provide a method for determining the parameters of multi-source superimposed acquisitions from land-based well shots, including the following steps:

[0072] S1 divides the target work area into multiple target regions. For each target region, the initial TD curve and shot coordinate information corresponding to the current target region are determined based on the three-dimensional seismic data of the current target region.

[0073] In this embodiment, the target work area is divided into multiple target regions based on the ignition lithology and surface undulation height of the target work area.

[0074] S2, according to the preset number of overlapping guns and gun spacing, group the gun coordinate information corresponding to the current target area, and based on the initial TD curve, determine the delay time of the secondary gun according to the distance between the main gun and the secondary gun in each group of gun coordinate information. Each group of gun coordinate information includes the main gun coordinate information and the secondary gun coordinate information.

[0075] S3. Based on the distance between the main gun and the secondary gun and the delay time of the secondary gun in each set of gun coordinate information, simulate aliasing is performed on the three-dimensional seismic data corresponding to the main gun coordinate information and the secondary gun coordinate information in each set of gun coordinate information to obtain simulated aliasing data of the target area. The simulated aliasing data of the target area is then separated to determine the similarity between the data before aliasing and the data after aliasing separation for each shot in the simulated aliasing data.

[0076] S4, if the similarity is less than the preset threshold, return the step of determining the delay time of the secondary gun based on the distance between the main gun and the secondary gun in each group of gun coordinate information according to the initial TD curve, until the similarity is greater than or equal to the preset threshold.

[0077] S5. When the similarity is greater than or equal to the preset threshold, the initial TD curve is updated according to the distance between the main gun and the secondary gun corresponding to the current simulated aliasing and the delay time of the secondary gun, and the updated TD curve is obtained as the well-shot multi-source aliasing acquisition parameters of the target area.

[0078] In this embodiment, step S1, determining the initial TD curve corresponding to the current target area based on the three-dimensional seismic data of the current target area, includes:

[0079] Load the observation system for the current target area, and select 3D seismic data from the historical 3D seismic data of the current target area that have an acquisition time exceeding the preset duration and a maximum offset greater than or equal to the preset distance, as the target 3D seismic data for the current target area. The preset duration can be 30s and the preset distance can be 12km.

[0080] Select the seismic data from the nearest offset receiver line of the target 3D seismic data, and divide the selected seismic data into two parts according to the positive and negative offsets. Select the data with slower energy decay from the two parts as the data to be processed.

[0081] The data to be processed is divided into multiple time windows with a time window size of 1 second and a distance of 1 km. The time windows can be overlapped by half a time window.

[0082] Calculate the root mean square amplitude value of the data within each time window, use the current time window position as the coordinate of the current root mean square amplitude value, and use the root mean square amplitude value to replace the time window data. Display the root mean square amplitude value of each time window through the interactive panel tool to determine the effective reflected energy and background energy.

[0083] The energy attenuation trend curve is picked out in the energy attenuation transition zone between the effective reflected energy and the background energy. The initial TD curve is plotted based on the inflection point of the trend curve. Multiple TD value ranges are set according to the inflection point value of the initial TD curve for subsequent verification of the aliasing separation effect.

[0084] See Figure 2 In the interactive panel tools, set the color mark display range (5-20). The root mean square amplitude value of different time windows will display different colors. There is an energy decay transition zone between the effective reflected strong energy color and the background energy color. Pick the energy decay trend curve in the energy transition zone and draw the initial TD curve value according to the inflection point of the trend curve.

[0085] In this embodiment, step S1, determining the shot coordinate information corresponding to the current target area based on the three-dimensional seismic data of the current target area, includes:

[0086] Load the observation system for the current target area, extract shot gather data from the 3D seismic data of the current target area according to the shot file number, stack multiple identical receiver point data in the shot gather data to obtain stacked data, and extract the file number, x and y coordinates corresponding to the shot file number in the stacked data as shot coordinate information;

[0087] The gun coordinate information is sorted by the gun file number, and duplicate gun coordinate information is deleted to ensure the uniqueness of the gun coordinate information for each shot.

[0088] In this embodiment, the data from multiple identical receiver points in the shot gather are superimposed, including:

[0089] Multiple identical receiver data points in the shot gather data are superimposed into one receiver data point, and duplicate receiver data points are deleted.

[0090] In this embodiment, step S2, grouping the gun coordinate information corresponding to the current target area according to the preset number of overlapping shots and the shot spacing, includes:

[0091] The first step is to mark each shot coordinate information that is not grouped as 0;

[0092] The second step is to select the gun coordinate information marked as 0 and ranked first as the main gun of the current group 1, and mark it as 1 with group number 1;

[0093] Third, from the remaining gun coordinate information marked 0, select the gun coordinate information whose distance from the main gun is within the first TD value range, and designate it as the first secondary gun of the main gun, marking it as 1, group number 1, and record its distance value from the main gun; from the remaining gun coordinate information marked 0, select the gun coordinate information whose distance from the main gun is within the second TD value range, and designate it as the second secondary gun of the main gun, marking it as 1, group number 1, and record its distance value from the main gun, see [link to previous steps]. Figure 2 The first TD value range is 1 to 3 km, and the second TD value range is 2 to 4 km.

[0094] Fourth step: If the number of mixed guns in the first group is 3, complete the grouping of the first group according to the first to third steps. If the number of mixed guns in the first group is more than 3, select the gun coordinate information with a distance from the main gun within the range of the N TD value from the remaining gun coordinate information marked as 0, and use it as the Nth secondary gun of the main gun, until the grouping of the first group is completed.

[0095] Fifth step: The gun coordinates already selected in the previous group 1 will not participate in the selection of the next group. The gun coordinates with the remaining mark of 0 and the first position will be used as the main gun of the current group 2, and the mark will be 1, and the group number will be 2.

[0096] Step 6: From the remaining gun coordinate information marked as 0, select the gun coordinate information within the preset progressive distance range after the first TD value range from the main gun, and use it as the first secondary gun of the main gun, marked as 1, group number 2, and record the distance value from the main gun; From the remaining gun coordinate information marked as 0, select the gun coordinate information within the preset progressive distance range after the second TD value range from the main gun, and use it as the second secondary gun of the main gun, marked as 1, group number 2, and record the distance value from the main gun. Wherein, when the first TD value range is 1~3km and the preset progressive distance is 2km, the preset progressive distance range after the first TD value range is 3~5km; when the second TD value range is 2~4km and the preset progressive distance is 2km, the preset progressive distance range after the second TD value range is 4~6km.

[0097] Step 7: If the number of overlapping guns in the second group is 3, complete the grouping of the second group according to steps 5 and 6. If the number of overlapping guns in the second group exceeds 3, select the gun coordinate information within the range of the Nth TD value from the remaining gun coordinate information marked as 0, and use it as the Nth secondary gun of the main gun, until the grouping of the second group is completed.

[0098] Step 8: Following steps 6 and 7, complete the selection of gun coordinate information for each remaining group. The distance between each group and the main gun is increased according to the preset progressive distance until the maximum offset distance is reached. Starting from the first TD value range, if no gun coordinate information that meets the first TD value range requirement for the distance from the main gun is found among the remaining gun coordinate information marked as 0, then only the main gun coordinate information is selected for that group, and no secondary gun coordinate information is selected.

[0099] In this embodiment, step S2, determining the delay time of the secondary gun based on the distance between the main gun and the secondary gun in each set of gun coordinate information according to the initial TD curve, includes:

[0100] Obtain multiple TD value ranges set based on the inflection point value of the initial TD curve. Each TD value range includes: the distance range between the main gun and the secondary gun in each set of gun coordinate information and the corresponding delay time of the secondary gun. See [link to relevant documentation]. Figure 2 If the distance range is less than 6km, a random delay time is obtained within the time range of 6s-15s; if the distance range is greater than or equal to 6km, a random delay time is obtained within the time range of 1s-6s.

[0101] Based on the distance between the main gun and the secondary gun in each set of gun coordinate information, determine the range of TD values ​​for that distance;

[0102] Based on the TD value range of the distance, the delay time of the secondary gun is randomly obtained from the delay time of the secondary gun corresponding to the distance range, wherein the delay time between secondary guns in the same group is greater than 0.5s.

[0103] In this embodiment, in step S3, the simulation aliasing processing is performed on the three-dimensional seismic data corresponding to the main gun coordinate information and the secondary gun coordinate information in each set of gun coordinate information, based on the distance between the main gun and the secondary gun and the delay time of the secondary gun in each set of gun coordinate information, to obtain the simulated aliasing data of the target area, including:

[0104] The first step is to store the trace data and trace file number that match the main gun file number in each set of shot coordinate information from the 3D seismic data into the current main gun shot coordinate information. Specifically, the trace data and trace file number that match the main gun file number in each set of shot coordinate information are obtained by using all shot line data and trace head information on the trace set data in the 3D seismic data.

[0105] The second step is to store the trace data and trace file number that match the secondary shot file number in each group of shot coordinate information from the 3D seismic data into the current secondary shot coordinate information.

[0106] The third step is that if there is no trace data and trace file number matching the main gun file number or the secondary gun file number in the 3D seismic data, then there is no trace data and trace file number in the shot coordinate information of the main gun or the secondary gun.

[0107] Fourth, for each group of gun coordinate information, if at least one of the main gun coordinate information and the secondary gun coordinate information in the current group lacks track data and track file number, then main gun / secondary gun aliasing processing is not performed; if both the main gun coordinate information and the secondary gun coordinate information in the current group match track data and track file number, then the following main gun / secondary gun aliasing processing is performed on the current group:

[0108] Fifth step, extract all sample values ​​of the matching track data in the current group of main guns, take the delay time of the first secondary gun in the current group as the start time of the first secondary gun sample value, and add it with all sample values ​​of the current main gun in a staggered manner. Take the delay time of the second secondary gun in the current group as the start time of the second secondary gun sample value, and add it with the sum of the current main gun and the first secondary gun again in a staggered manner, until all secondary guns in the current group have been added in a staggered manner, thus completing the data aliasing of the current group;

[0109] Step 6: Perform truncation processing on the aliased data of the current group. Obtain the sample point values ​​of the main gun and secondary gun respectively according to the actual sampling time before aliasing. Replace the sample point data value before aliasing with the truncation main gun sample point data value according to the track file number matched by the main gun in the current group. For the first secondary gun, truncate sample point data of the same length according to its start time and delay time, and replace the sample point data value before aliasing with the track file number matched by the first secondary gun in the current group. For the second secondary gun, truncate sample point data of the same length according to its start time and delay time, and replace the sample point data value before aliasing with the track file number matched by the second secondary gun in the current group. Continue until all secondary guns have undergone aliased data truncation, thus completing the aliased data truncation of the current group.

[0110] Step 7: After completing the aliasing and truncation of the current group of data, perform aliasing and truncation on the next group of data according to steps 5 and 6, until all groups have completed aliasing and truncation, and obtain the simulated aliasing data of the target area.

[0111] In this embodiment, step S3, separating the simulated aliased data of the target region, includes:

[0112] Since most land-based well shots are preprocessed in the form of shot collections, for the simulated aliasing data of the target area, the non-aliased shots at adjacent positions are searched based on the shot coordinate information of the simulated aliasing data as model shots, and multi-scale curvelet transform is performed on the aliased shots and model shots in the shot domain.

[0113] At each scale, structural pattern recognition and threshold separation are performed on the aliased gun and the model gun. Then, the separated data are transformed into gun domain data through multi-scale curvelet inverse transform to complete the separation of aliased data.

[0114] In this embodiment, in step S3, the quality control of the separation effect is based on calculating the aliased data acquisition points as the basis for the actual TD curve values. Similarity analysis is performed on the data before and after aliasing for quality control. The similarity between the data before and after aliasing separation for each shot in the simulated aliasing data is determined by the following expression:

[0115]

[0116] Where SE(x,y) represents the similarity between the pre-aliasing and post-aliasing data of each shot in the simulated aliasing data, x and y are the pre-aliasing and post-aliasing data respectively, and Ax and Ay are the Fourier amplitudes of data x and y respectively. and Let x and y represent the Fourier phases of the data respectively, and N be the number of shots.

[0117] In this embodiment, using Figure 3 The receiver point gathers of the three-dimensional data are generated according to the method for determining the parameters of multi-source superimposed acquisitions from land wells and shots disclosed in this invention, as follows: Figure 4 The simulated aliasing data at different distances and times shown are subjected to aliasing separation to obtain the following results: Figure 5 The data shown.

[0118] Through similarity analysis Figure 3 Data and Figure 5 The data similarity was 98.3%, meeting the standard of aliasing separation greater than 95%. Using the separated data that met the standard, different distance and time pairs were obtained and plotted as follows: Figure 6 The TD parameter curves for well-blast combined mining are shown.

[0119] This disclosed method for determining parameters of multi-source cascading acquisition in land-based well-shot seismic fields involves collecting previously acquired single-shot data from the area, preprocessing it to extract shot coordinate information, and grouping the data according to the number of cascading shots and the distance between shots. For each grouped shot coordinate, the delay time of the secondary shot is obtained according to the excitation time range. Data matching and simulated cascading processing are performed on the grouped shot coordinate information and previously acquired old data. The simulated cascading data undergoes cascading interference separation processing, and the separated data is analyzed for similarity. Based on the time and distance information showing good separation results, a TD curve for the area is plotted as the basis for field cascading construction parameters. Only cascading shots that meet this TD curve can achieve effective cascading interference separation indoors. Only then can the field cascading data be transformed into a more complete cascading parameter demonstration for complex well-shot areas in mountainous regions. The generated TD cascading curve meets the requirements of efficient field cascading construction and high-quality indoor separation processing.

[0120] The aliasing separation step disclosed herein can record the sequence number of each shot trace using aliased shot gather data acquired by TD curves, extract and retain the outlier and empty traces of each shot data, perform multi-scale transformation on the aliased data and model data, calculate their scale coefficients respectively, use the scale coefficients of the aliased data and model data to calculate the model characteristic coefficients at different scales between the two, set a threshold value for filtering, perform multi-scale inverse transformation on the filtered model characteristic coefficients, and merge the obtained aliased data with the extracted outlier and empty traces to obtain the final separated main shot data and aliased data, which is convenient for subsequent high-precision processing of seismic data.

[0121] See Figure 7 The present disclosure provides an apparatus for determining parameters of multi-source superimposed acquisition data from land-based well-shot seismic sources, comprising:

[0122] Module 11 is used to divide the target work area into multiple target regions. For each target region, the initial TD curve and shot coordinate information corresponding to the current target region are determined based on the three-dimensional seismic data of the current target region.

[0123] Grouping module 12 is used to group the gun coordinate information corresponding to the current target area according to the preset number of overlapping guns and gun spacing, and based on the initial TD curve, determine the delay time of the secondary gun according to the distance between the main gun and the secondary gun in each group of gun coordinate information. Each group of gun coordinate information includes the main gun coordinate information and the secondary gun coordinate information.

[0124] The aliasing module 13 is used to simulate aliasing of the three-dimensional seismic data corresponding to the main gun coordinate information and the secondary gun coordinate information in each set of gun coordinate information according to the distance between the main gun and the secondary gun and the delay time of the secondary gun in each set of gun coordinate information, to obtain the simulated aliasing data of the target area, and to separate the simulated aliasing data of the target area to determine the similarity between the data before aliasing and the data after aliasing separation in the simulated aliasing data.

[0125] Return module 14 is used to return the steps of determining the delay time of the secondary gun based on the distance between the main gun and the secondary gun in each group of gun coordinate information, based on the initial TD curve, when the similarity is less than the preset threshold, until the similarity is greater than or equal to the preset threshold.

[0126] The update module 15 is used to update the initial TD curve based on the distance between the main gun and the secondary gun and the delay time of the secondary gun corresponding to the current simulated aliasing when the similarity is greater than or equal to a preset threshold, so as to obtain the updated TD curve as the well-shot multi-source aliasing acquisition parameters of the target area.

[0127] The specific implementation process of the functions and roles of each unit in the above device can be found in the implementation process of the corresponding steps in the above method, and will not be repeated here.

[0128] For the device embodiments, since they basically correspond to the method embodiments, the relevant parts can be referred to in the description of the method embodiments. The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of the present invention according to actual needs. Those skilled in the art can understand and implement this without creative effort.

[0129] In the above embodiments, any and more of the determination module 11, grouping module 12, aliasing module 13, return module 14, and update module 15 can be combined into one module, or any one of these modules can be split into multiple modules. Alternatively, at least part of the functionality of one or more of these modules can be combined with at least part of the functionality of other modules and implemented in one module. At least one of the determination module 11, grouping module 12, aliasing module 13, return module 14, and update module 15 can be at least partially implemented as hardware circuitry, such as a field-programmable gate array (FPGA), programmable logic array (PLA), system-on-a-chip, system-on-a-substrate, system-on-package, application-specific integrated circuit (ASIC), or any other reasonable means of integrating or packaging circuitry, or implemented in software, hardware, or firmware, or in any appropriate combination of any of these three implementation methods. Alternatively, at least one of the determination module 11, grouping module 12, aliasing module 13, return module 14, and update module 15 can be at least partially implemented as a computer program module, which can perform corresponding functions when the computer program module is run.

[0130] Reference Figure 8 As shown, the electronic device provided in the embodiments of this disclosure includes a processor 1110, a communication interface 1120, a memory 1130, and a communication bus 1140, wherein the processor 1110, the communication interface 1120, and the memory 1130 communicate with each other through the communication bus 1140.

[0131] Memory 1130 is used to store computer programs;

[0132] When processor 1110 executes the program stored in memory 1130, it implements the following method for determining the parameters of multi-source cascading acquisition in land well shots:

[0133] The target work area is divided into multiple target regions. For each target region, the initial TD curve and shot coordinate information corresponding to the current target region are determined based on the three-dimensional seismic data of the current target region.

[0134] According to the preset number of overlapping guns and gun spacing, the gun coordinate information corresponding to the current target area is grouped, and based on the initial TD curve, the delay time of the secondary gun is determined according to the distance between the main gun and the secondary gun in each group of gun coordinate information. Each group of gun coordinate information includes the main gun coordinate information and the secondary gun coordinate information.

[0135] Based on the distance between the main gun and the secondary gun and the delay time of the secondary gun in each set of gun coordinate information, simulated aliasing processing is performed on the three-dimensional seismic data corresponding to the main gun coordinate information and the secondary gun coordinate information in each set of gun coordinate information to obtain simulated aliasing data of the target area. The simulated aliasing data of the target area is then separated to determine the similarity between the data before aliasing and the data after aliasing separation for each shot in the simulated aliasing data.

[0136] If the similarity is less than a preset threshold, return to the step of determining the delay time of the secondary gun based on the distance between the main gun and the secondary gun in each set of gun coordinate information, based on the initial TD curve, until the similarity is greater than or equal to the preset threshold.

[0137] If the similarity is greater than or equal to a preset threshold, the initial TD curve is updated based on the distance between the main gun and the secondary gun corresponding to the current simulated aliasing and the delay time of the secondary gun, and the updated TD curve is obtained as the well-shot multi-source aliasing acquisition parameters for the target area.

[0138] The aforementioned communication bus 1140 can be a Peripheral Component Interconnect (PCI) bus or an Extended Industry Standard Architecture (EISA) bus, etc. This communication bus 1140 can be divided into an address bus, a data bus, a control bus, etc. For ease of illustration, it is represented by only one thick line in the figure, but this does not indicate that there is only one bus or one type of bus.

[0139] The communication interface 1120 is used for communication between the above-mentioned electronic device and other devices.

[0140] The memory 1130 may include random access memory (RAM) or non-volatile memory, such as at least one disk storage device. Optionally, the memory 1130 may also be at least one storage device located remotely from the aforementioned processor 1110.

[0141] The processor 1110 mentioned above can be a general-purpose processor, including a central processing unit (CPU), a network processor (NP), etc.; it can also be a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field-programmable gate array (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, or discrete hardware components.

[0142] Embodiments of this disclosure also provide a computer-readable storage medium. The computer-readable storage medium stores a computer program that, when executed by a processor, implements the method for determining parameters of multi-source cascaded acquisition in land-based well-shot seismic systems as described above.

[0143] The computer-readable storage medium may be included in the device / apparatus described in the above embodiments; or it may exist independently and not assembled into the device / apparatus. The computer-readable storage medium carries one or more programs that, when executed, implement the method for determining multi-source superimposed acquisition parameters for land well shots according to embodiments of this disclosure.

[0144] According to embodiments of this disclosure, the computer-readable storage medium can be a non-volatile computer-readable storage medium, such as including, but not limited to: portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination thereof. In this disclosure, the computer-readable storage medium can be any tangible medium that contains or stores a program that can be used by or in conjunction with an instruction execution system, apparatus, or device.

[0145] It should be noted that, in this document, relational terms such as "first" and "second" are used merely to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes said element.

[0146] The above description is merely a specific embodiment of this disclosure, enabling those skilled in the art to understand or implement it. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of this disclosure. Therefore, this disclosure is not to be limited to the embodiments shown herein, but is to be accorded the widest scope consistent with the principles and novel features claimed herein.

Claims

1. A method for determining parameters of multi-source superimposed acquisition data from land-based seismic wells, characterized in that, The method includes: The target work area is divided into multiple target regions. For each target region, the initial TD curve and shot coordinate information corresponding to the current target region are determined based on the three-dimensional seismic data of the current target region. According to the preset number of overlapping guns and gun spacing, the gun coordinate information corresponding to the current target area is grouped, and based on the initial TD curve, the delay time of the secondary gun is determined according to the distance between the main gun and the secondary gun in each group of gun coordinate information. Each group of gun coordinate information includes the main gun coordinate information and the secondary gun coordinate information. Based on the distance between the main gun and the secondary gun and the delay time of the secondary gun in each set of gun coordinate information, simulated aliasing processing is performed on the three-dimensional seismic data corresponding to the main gun coordinate information and the secondary gun coordinate information in each set of gun coordinate information to obtain simulated aliasing data of the target area. The simulated aliasing data of the target area is then separated to determine the similarity between the data before aliasing and the data after aliasing separation for each shot in the simulated aliasing data. If the similarity is less than a preset threshold, return to the step of determining the delay time of the secondary gun based on the distance between the main gun and the secondary gun in each set of gun coordinate information, based on the initial TD curve, until the similarity is greater than or equal to the preset threshold. If the similarity is greater than or equal to a preset threshold, the initial TD curve is updated based on the distance between the main gun and the secondary gun corresponding to the current simulated aliasing and the delay time of the secondary gun, and the updated TD curve is obtained as the well-shot multi-source aliasing acquisition parameters for the target area.

2. The method according to claim 1, characterized in that, The step of determining the initial TD curve corresponding to the current target area based on the three-dimensional seismic data of the current target area includes: From the historical 3D seismic data of the current target area, select 3D seismic data with an acquisition time exceeding the preset duration and a maximum offset greater than or equal to the preset distance, and use them as the target 3D seismic data of the current target area. Select the seismic data from the nearest offset receiver line of the target 3D seismic data, and divide the selected seismic data into two parts according to the positive and negative offsets. Select the data with slower energy decay from the two parts as the data to be processed. The data to be processed is divided into multiple time windows with a time window size of 1 second and a distance of 1 km. The time windows are selected by overlapping half of the time windows. Calculate the root mean square amplitude value of the data within each time window, use the current time window position as the coordinate of the current root mean square amplitude value, and use the root mean square amplitude value to replace the time window data to display the root mean square amplitude value of each time window, thereby determining the effective reflected energy and background energy. The energy attenuation trend curve is picked out in the energy attenuation transition zone between the effective reflected energy and the background energy. The initial TD curve is plotted based on the inflection point of the trend curve. Multiple TD value ranges are set based on the inflection point value of the initial TD curve.

3. The method according to claim 1, characterized in that, The step of determining the shot coordinate information corresponding to the current target area based on the three-dimensional seismic data of the current target area includes: From the 3D seismic data of the current target area, the shot gather data is extracted according to the shot file number. The data of multiple identical receiver points in the shot gather data are overlaid to obtain the overlaid data. The file number, x and y coordinates corresponding to the shot file number in the overlaid data are extracted as shot coordinate information. The gun coordinate information is sorted by the gun file number, and duplicate gun coordinate information is deleted to ensure the uniqueness of the gun coordinate information for each shot.

4. The method according to claim 1, characterized in that, The step of grouping the gun coordinate information corresponding to the current target area according to the preset number of overlapping shots and the shot spacing includes: The first step is to mark each shot coordinate information that is not grouped as 0; The second step is to select the gun coordinate information marked as 0 and ranked first as the main gun of the current group 1, and mark it as 1 with group number 1; Third, among the remaining gun coordinate information marked as 0, select the gun coordinate information whose distance from the main gun is within the first TD value range, and use it as the first secondary gun of the main gun, mark it as 1, group number 1, and record the distance value from the main gun; among the remaining gun coordinate information marked as 0, select the gun coordinate information whose distance from the main gun is within the second TD value range, and use it as the second secondary gun of the main gun, mark it as 1, group number 1, and record the distance value from the main gun. Fourth step: If the number of mixed guns in the first group is 3, complete the grouping of the first group according to the first to third steps. If the number of mixed guns in the first group is more than 3, select the gun coordinate information with a distance from the main gun within the range of the N TD value from the remaining gun coordinate information marked as 0, and use it as the Nth secondary gun of the main gun, until the grouping of the first group is completed. Fifth step: Take the gun coordinate information that is marked as 0 and ranked first as the main gun of the current group 2, and mark it as 1 and group number 2; Step 6: From the remaining gun coordinate information marked as 0, select the gun coordinate information whose distance from the main gun is within a preset progressive distance range after the first TD value range, and use it as the first secondary gun of the main gun, marked as 1, group number 2, and record the distance value from the main gun; From the remaining gun coordinate information marked as 0, select the gun coordinate information whose distance from the main gun is within a preset progressive distance range after the second TD value range, and use it as the second secondary gun of the main gun, marked as 1, group number 2, and record the distance value from the main gun. Step 7: If the number of overlapping guns in the second group is 3, complete the grouping of the second group according to steps 5 and 6. If the number of overlapping guns in the second group exceeds 3, select the gun coordinate information within the range of the Nth TD value from the remaining gun coordinate information marked as 0, and use it as the Nth secondary gun of the main gun, until the grouping of the second group is completed. Step 8: Following steps 6 and 7, complete the selection of gun coordinate information for each remaining group. The distance between each group and the main gun is increased according to the preset progressive distance until the maximum offset distance is reached. Starting from the first TD value range, if no gun coordinate information that meets the first TD value range requirement for the distance from the main gun is found among the remaining gun coordinate information marked as 0, then only the main gun coordinate information is selected for that group, and no secondary gun coordinate information is selected.

5. The method according to claim 1, characterized in that, The determination of the secondary gun's delay time based on the initial TD curve and the distance between the main gun and secondary gun in each set of gun coordinate information includes: Obtain multiple TD value ranges set according to the inflection point value of the initial TD curve. Each TD value range includes: the distance range between the main gun and the secondary gun in each set of gun coordinate information and the corresponding delay time of the secondary gun. Based on the distance between the main gun and the secondary gun in each set of gun coordinate information, determine the range of TD values ​​for that distance; Based on the TD value range of the distance, the delay time of the secondary gun is randomly obtained from the delay time of the secondary gun corresponding to the distance range, wherein the delay time between secondary guns in the same group is greater than 0.5s.

6. The method according to claim 1, characterized in that, The method involves performing simulated aliasing processing on the three-dimensional seismic data corresponding to the main gun and secondary gun coordinates in each set of gun coordinate information, based on the distance between the main gun and secondary gun and the delay time of the secondary gun in each set of gun coordinate information, to obtain simulated aliasing data of the target area, including: The first step is to store the trace data and trace file number that match the main gun file number in each group of shot coordinate information from the 3D seismic data into the current main gun coordinate information. The second step is to store the trace data and trace file number that match the secondary shot file number in each group of shot coordinate information from the 3D seismic data into the current secondary shot coordinate information. The third step is that if there is no trace data and trace file number matching the main gun file number or the secondary gun file number in the 3D seismic data, then there is no trace data and trace file number in the shot coordinate information of the main gun or the secondary gun. Fourth, for each group of gun coordinate information, if at least one of the main gun coordinate information and the secondary gun coordinate information in the current group lacks track data and track file number, then main gun / secondary gun aliasing processing is not performed; if both the main gun coordinate information and the secondary gun coordinate information in the current group match track data and track file number, then the following main gun / secondary gun aliasing processing is performed on the current group: Fifth step, extract all sample values ​​of the matching track data in the current group of main guns, take the delay time of the first secondary gun in the current group as the start time of the first secondary gun sample value, and add it with all sample values ​​of the current main gun in a staggered manner. Take the delay time of the second secondary gun in the current group as the start time of the second secondary gun sample value, and add it with the sum of the current main gun and the first secondary gun again in a staggered manner, until all secondary guns in the current group have been added in a staggered manner, thus completing the data aliasing of the current group; Step 6: Perform truncation processing on the aliased data of the current group. Obtain the sample point values ​​of the main gun and secondary gun respectively according to the actual sampling time before aliasing. Replace the sample point data value before aliasing with the truncation main gun sample point data value according to the track file number matched by the main gun in the current group. For the first secondary gun, truncate sample point data of the same length according to its start time and delay time, and replace the sample point data value before aliasing with the track file number matched by the first secondary gun in the current group. For the second secondary gun, truncate sample point data of the same length according to its start time and delay time, and replace the sample point data value before aliasing with the track file number matched by the second secondary gun in the current group. Continue until all secondary guns have undergone aliased data truncation, thus completing the aliased data truncation of the current group. Step 7: After completing the aliasing and truncation of the current group of data, perform aliasing and truncation on the next group of data according to steps 5 and 6, until all groups have completed aliasing and truncation, and obtain the simulated aliasing data of the target area.

7. The method according to claim 1, characterized in that, The separation of simulated aliased data in the target region includes: For the simulated aliasing data of the target area, the non-aliased guns at adjacent positions are searched based on the gun coordinate information of the simulated aliasing data as model guns, and multi-scale curvelet transform is performed on the aliased guns and model guns in the gun domain. At each scale, structural pattern recognition and threshold separation are performed on the aliased gun and the model gun. Then, the separated data are transformed into gun domain data through multi-scale curvelet inverse transform to complete the separation of aliased data.

8. The method according to claim 1, characterized in that, The similarity between the pre-aliasing data and the post-aliasing data for each shot in the simulated aliasing data is determined using the following expression: in, To simulate the similarity between the pre-aliasing and post-aliasing data for each shot in the aliasing data, let x and y be the pre-aliasing and post-aliasing data, respectively, and let Ax and Ay be the Fourier amplitudes of data x and y, respectively. and These represent the Fourier phases of data x and y, respectively. and Let x and y represent the mean of the Fourier phases of the data x and y, respectively, and N be the number of shots.

9. A device for determining parameters of multi-source superimposed acquisition data from land-based seismic wells, characterized in that, include: The determination module is used to divide the target work area into multiple target regions. For each target region, the initial TD curve and shot coordinate information corresponding to the current target region are determined based on the 3D seismic data of the current target region. The grouping module is used to group the gun coordinate information corresponding to the current target area according to the preset number of overlapping guns and gun spacing, and based on the initial TD curve, determine the delay time of the secondary gun according to the distance between the main gun and the secondary gun in each group of gun coordinate information. Each group of gun coordinate information includes the main gun coordinate information and the secondary gun coordinate information. The aliasing module is used to simulate aliasing the 3D seismic data corresponding to the main gun coordinate information and the secondary gun coordinate information in each set of gun coordinate information based on the distance between the main gun and the secondary gun and the delay time of the secondary gun in each set of gun coordinate information. This process yields simulated aliased data of the target area, and the simulated aliased data of the target area is separated to determine the similarity between the data before aliasing and the data after aliasing separation for each shot in the simulated aliased data. The return module is used to return the steps of determining the delay time of the secondary gun based on the distance between the main gun and the secondary gun in each group of gun coordinate information, based on the initial TD curve, when the similarity is less than a preset threshold, until the similarity is greater than or equal to the preset threshold. The update module is used to update the initial TD curve based on the distance between the main gun and the secondary gun and the delay time of the secondary gun corresponding to the current simulated aliasing when the similarity is greater than or equal to a preset threshold. The updated TD curve is then used as the acquisition parameters for well-shot multi-source aliasing in the target area.

10. An electronic device, characterized in that, It includes a processor, a communication interface, a memory, and a communication bus, wherein the processor, the communication interface, and the memory communicate with each other through the communication bus; Memory, used to store computer programs; The processor, when executing a program stored in memory, implements the method for determining the multi-source superimposed acquisition parameters of land well shots as described in any one of claims 1-8.

11. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the method for determining the multi-source superimposed acquisition parameters of land well shot as described in any one of claims 1-8.

Images