Method and device for determining multi-seismic-source aliasing acquisition parameters of land well gun

Abstract

The invention relates to a land well shot multi-seismic-source aliasing acquisition parameter determination method and device, and the method comprises the steps: dividing a target work area into a plurality of target regions, and determining an initial T-D curve and shot coordinate information corresponding to the current target region; grouping the shot coordinate information corresponding to the current target area, and determining the delay time of the secondary shot based on the initial T-D curve; performing analog aliasing processing on the three-dimensional seismic data respectively corresponding to the main gun and the auxiliary gun in each group of gun coordinate information according to the delay time of the auxiliary gun in each group of gun coordinate information to obtain analog aliasing data of the target area, separating the aliasing data, and determining the similarity between the data of each gun before aliasing and the data after aliasing separation in the aliasing data; and under the condition that the similarity is greater than or equal to a preset threshold value, updating the initial T-D curve to obtain an updated T-D curve, taking the updated T-D curve as an aliasing acquisition parameter of the target area and a mixed acquisition shot meeting the T-D curve, and carrying out effective aliasing interference separation indoors.

Description

Technical Field

[0001] The present disclosure relates to the technical field of oil exploration seismic data acquisition and processing, and particularly to a method and device for determining multi-source aliasing acquisition parameters for land well shots. Background Art

[0002] In onshore vibrator oil and gas exploration acquisition, multiple vibrators are often used to simultaneously or with short delays to excite, and long-array nodal instruments are used to receive long-record data excited by multiple shot points. This multi-source mixed acquisition method is an effective way for onshore vibrators to improve acquisition efficiency and reduce acquisition costs. The key to implementing this technology is to conduct effective aliasing acquisition parameter demonstration for the acquisition work area. These parameters can not only ensure the implementation of efficient field acquisition but also ensure high-precision aliasing separation indoors.

[0003] However, the exploration environment in the western mountainous areas is complex. In oil and gas exploration, explosive well-shot vibrators are mostly used for excitation. The differences in each shot record excited by well shots in different regions are relatively large, and the non-uniformity is strong. The parameter method for vibrator mixed acquisition demonstration cannot meet the needs of complex mountain well-shot mixed acquisition. Therefore, it is necessary to study the method and process for demonstrating well-shot multi-source aliasing acquisition parameters for the unique acquisition environment of complex mountain oil and gas exploration to meet the requirements of efficient field well-shot acquisition and effective indoor aliasing separation. Summary of the Invention

[0004] To solve the above technical problems or at least partially solve the above technical problems, embodiments of the present disclosure provide a method and device for determining multi-source aliasing acquisition parameters for land well shots.

[0005] In a first aspect, embodiments of the present disclosure provide a method for determining multi-source aliasing acquisition parameters for land well shots, including:

[0006] Dividing a target work area into multiple target regions. For each target region, determining an initial T-D curve and shot coordinate information corresponding to the current target region according to the 3D seismic data of the current target region;

[0007] Grouping the shot coordinate information corresponding to the current target region according to a preset number of aliasing shots and shot spacing, and based on the initial T-D curve, determining the delay time of the secondary shot according to the distance between the primary shot and the secondary shot in each group of shot coordinate information, where each group of shot coordinate information includes primary shot coordinate information and secondary shot coordinate information;

[0008] According to the distance between the main gun and the secondary gun in each group of gun coordinate information and the delay time of the secondary gun, perform simulated aliasing processing on the three-dimensional seismic data corresponding to the main gun coordinate information and the secondary gun coordinate information in each group of gun coordinate information to obtain the simulated aliased data of the target area, and separate the simulated aliased data of the target area to determine the similarity between the data before aliasing and the data after aliasing separation for each gun in the simulated aliased data;

[0009] In the case where the similarity is less than the preset threshold, return to the step of determining 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 based on the initial T-D curve until the similarity is greater than or equal to the preset threshold;

[0010] In the case where the similarity is greater than or equal to the preset threshold, update the initial T-D curve 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 to obtain the updated T-D curve, which is used as the multi-source aliasing acquisition parameter of the well-gun in the target area.

[0011] In a possible implementation manner, the determining the initial T-D curve corresponding to the current target area according to the three-dimensional seismic data of the current target area includes:

[0012] Select, from the historical three-dimensional seismic data of the current target area, the three-dimensional seismic data with the acquisition time length exceeding the preset duration and the maximum offset distance greater than or equal to the preset distance as the target three-dimensional seismic data of the current target area;

[0013] Select the seismic data of the nearest offset receiving line of the target three-dimensional seismic data, and divide the selected seismic data into two parts according to the positive and negative offset distances, and select the data with slower energy attenuation in the two parts as the data to be processed;

[0014] Taking 1 s as the time window size and 1 km as the distance window size, divide the data to be processed into multiple time window data according to the time window size, where the time window data is selected with an overlap of half a time window;

[0015] Calculate the root mean square amplitude value of the data in each time window, take 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, and determine the effective reflection energy and the background energy;

[0016] Pick up the energy attenuation trend curve in the energy attenuation transition zone between the effective reflection energy and the background energy, and draw the initial T-D curve according to the inflection point of the trend curve. Among them, according to the inflection point value of the initial T-D curve, set multiple T-D value ranges.

[0017] In a possible implementation manner, the determining the gun coordinate information corresponding to the current target area according to the three-dimensional seismic data of the current target area includes:

[0018] Extract shot gather data from the 3D seismic data of the current target area according to the shot file number, perform stacking processing on the data of multiple identical geophone points 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;

[0019] Sort the shot coordinate information according to the shot file number, and delete the duplicate shot coordinate information to ensure the uniqueness of the shot coordinate information for each shot.

[0020] In a possible implementation manner, the grouping of the shot coordinate information corresponding to the current target area according to the preset aliased shot number and shot spacing includes:

[0021] First step, mark each ungrouped shot coordinate information as 0;

[0022] Second step, select the shot coordinate information marked as 0 and ranked first as the main shot of the current group 1, mark it as 1, and the group number is 1;

[0023] Third step, among the remaining shot coordinate information marked as 0, select the shot coordinate information whose distance from the main shot is within the first T-D value range as the first sub-shot of the main shot, mark it as 1, and the group number is 1, and record the distance value from the main shot; among the remaining shot coordinate information marked as 0, select the shot coordinate information whose distance from the main shot is within the second T-D value range as the second sub-shot of the main shot, mark it as 1, and the group number is 1, and record the distance value from the main shot;

[0024] Fourth step, if the aliased shot number of the first group is 3 shots, complete the grouping of the first group according to the first to third steps. If the aliased shot number of the first group exceeds 3 shots, execute the step of selecting the shot coordinate information whose distance from the main shot is within the Nth T-D value range as the Nth sub-shot of the main shot among the remaining shot coordinate information marked as 0 until the grouping of the first group is completed;

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

[0026] Sixth step, among the remaining shot coordinate information marked as 0, select the shot coordinate information whose distance from the main shot is within the range after shifting the preset progressive distance within the first T-D value range as the first sub-shot of the main shot, mark it as 1, and the group number is 2, and record the distance value from the main shot; among the remaining shot coordinate information marked as 0, select the shot coordinate information whose distance from the main shot is within the range after shifting the preset progressive distance within the second T-D value range as the second sub-shot of the main shot, mark it as 1, and the group number is 2, and record the distance value from the main shot;

[0027] Step 7: If the aliasing gun count of the second group is 3 guns, complete the grouping of the second group according to Steps 5 to 6. If the aliasing gun count of the second group exceeds 3 guns, perform the step of selecting the gun coordinate information within the N T-D value range from the remaining gun coordinate information marked as 0 and having a distance from the main gun as the Nth secondary gun of the main gun until the grouping of the second group is completed;

[0028] Step 8: Complete the selection of the gun coordinate information for each remaining group according to Steps 6 and 7. The distance between each group and the main gun progresses by a preset progressive distance until the maximum offset is reached. Starting from the first T-D value range, if no gun coordinate information that meets the requirement of the first T-D value range of the distance from the main gun is found in the remaining gun coordinate information marked as 0, then only the gun coordinate information of the main gun is selected for this group, and there is no gun coordinate information of the secondary gun.

[0029] In a possible implementation manner, determining the delay time of the secondary gun based on the initial T-D curve according to the distance between the main gun and the secondary gun in each group of gun coordinate information includes:

[0030] Obtain a plurality of T-D value ranges set according to the inflection point values of the initial T-D curve, where each T-D value range includes: the distance range between the main gun and the secondary gun in each group of gun coordinate information and the corresponding delay time of the secondary gun;

[0031] Determine the T-D value range where the distance is located according to the distance between the main gun and the secondary gun in each group of gun coordinate information;

[0032] Based on the T-D value range where the distance is located, randomly obtain the delay time of the secondary gun from the delay times of the secondary guns corresponding to the distance range, where the delay time between the secondary guns in the same group is greater than 0.5 s.

[0033] In a possible implementation manner, performing simulated aliasing processing on the 3D seismic data corresponding to the main gun coordinate information and the secondary gun coordinate information in each group of gun coordinate information according to the distance between the main gun and the secondary gun in each group of gun coordinate information and the delay time of the secondary gun to obtain the simulated aliased data of the target area includes:

[0034] Step 1: Store the trace data and trace file numbers in the 3D seismic data that match the main gun file number in each group of gun coordinate information into the current main gun coordinate information;

[0035] Step 2: Store the trace data and trace file numbers in the 3D seismic data that match the secondary gun file number in each group of gun coordinate information into the current secondary gun coordinate information;

[0036] Step 3: If there is no trace data and trace file number in the 3D seismic data that match the main gun file number or the secondary gun file number, then there is no trace data and trace file number in the gun coordinate information of this main gun or secondary gun;

[0037] Step 4: For each group of gun coordinate information, when at least one of the main gun coordinate information and the secondary gun coordinate information in the current group of gun coordinate information has no trace data and trace file number, the main gun and secondary gun aliasing processing is not performed; when both the main gun coordinate information and the secondary gun coordinate information in the current group of gun coordinate information match the trace data and trace file number, the following main gun and secondary gun aliasing processing is performed on the current group:

[0038] Step 5: Take out all the sample point values that match the trace data in the main gun of the current group, use the delay time of the first secondary gun in the current group as the start time of the first secondary gun sample point values, perform misaligned addition with all the sample point values of the current main gun, use the delay time of the second secondary gun in the current group as the start time of the second secondary gun sample point values, and perform misaligned addition with the sum of the current main gun and the first secondary gun again until all the secondary guns in the current group have been misaligned and added, completing the data aliasing of the current group;

[0039] Step 6: Perform truncation processing on the aliased data of the current group, obtain the main and secondary gun sample point values according to the actual sampling time before aliasing respectively, and use the truncated main gun sample data values to replace the sample data values before aliasing according to the trace file number matched by the main gun of the current group; the truncation start time of the first secondary gun is based on its delay time, truncate the sample data of the same length, and use the trace file number matched by the first secondary gun of the current group to replace the sample data values before aliasing; the truncation start time of the second secondary gun is based on its delay time, truncate the sample data of the same length, and use the trace file number matched by the second secondary gun of the current group to replace the sample data values before aliasing until all the secondary guns have completed the truncation of the aliased data, completing the truncation of the aliased data of the current group;

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

[0041] In a possible implementation manner, separating the simulated aliased data of the target area includes:

[0042] For the simulated aliased data of the target area, search for the non-aliased guns at adjacent positions according to the gun coordinate information of the simulated aliased data as the model guns, and perform multi-scale curvelet transform on the aliased guns and the model guns in the gun domain;

[0043] Perform structure pattern recognition and threshold separation processing on the aliased guns and the model guns at each scale, and then perform multi-scale curvelet inverse transform on the separated data to the gun domain data, completing the separation of the aliased data.

[0044] In a possible implementation manner, the similarity between the data before aliasing and the data after aliasing separation of each gun in the simulated aliased data is determined through the following expression:

[0045]

[0046]

[0047] Among them, SE(x, y) is the similarity between the data before aliasing and the data after aliasing separation for each shot in the simulated aliased data. x and y are the data before aliasing and the data after separation respectively, Ax and Ay are the Fourier amplitudes of the data x and y respectively, and represent the Fourier phases of the data x and y respectively, and N is the number of shots.

[0048] In a second aspect, an embodiment of the present disclosure provides a device for determining multi-source aliasing acquisition parameters for land well shots, including:

[0049] A determination module, configured to divide a target work area into multiple target areas. For each target area, determine the initial T-D curve and shot coordinate information corresponding to the current target area according to the three-dimensional seismic data of the current target area;

[0050] A grouping module, configured to group the shot coordinate information corresponding to the current target area according to a preset number of aliased shots and shot spacing, and based on the initial T-D curve, determine the delay time of the secondary shot according to the distance between the primary shot and the secondary shot in each group of shot coordinate information. Among them, each group of shot coordinate information includes the primary shot coordinate information and the secondary shot coordinate information;

[0051] An aliasing module, configured to perform simulated aliasing processing on the three-dimensional seismic data corresponding to the primary shot coordinate information and the secondary shot coordinate information in each group of shot coordinate information according to the distance between the primary shot and the secondary shot in each group of shot coordinate information and the delay time of the secondary shot, obtain the simulated aliased data of the target area, and separate the simulated aliased data of the target area to determine the similarity between the data before aliasing and the data after aliasing separation for each shot in the simulated aliased data;

[0052] A return module, configured to, when the similarity is less than a preset threshold, return the step of determining the delay time of the secondary shot according to the distance between the primary shot and the secondary shot in each group of shot coordinate information based on the initial T-D curve until the similarity is greater than or equal to the preset threshold;

[0053] An update module, configured to, when the similarity is greater than or equal to the preset threshold, update the initial T-D curve according to the distance between the primary shot and the secondary shot corresponding to the current simulated aliasing and the delay time of the secondary shot to obtain an updated T-D curve, which is used as the multi-source aliasing acquisition parameters for well shots in the target area.

[0054] In a third aspect, an embodiment of the present disclosure provides an electronic device, including a processor, a communication interface, a memory, and a communication bus. Among them, the processor, the communication interface, and the memory complete communication with each other through the communication bus;

[0055] The memory is used to store a computer program;

[0056] The processor is configured to implement the above-mentioned method for determining the multi-source aliasing acquisition parameters of land well guns when executing the program stored on the memory.

[0057] In a fourth aspect, an embodiment of the present disclosure provides a computer-readable storage medium, on which a computer program is stored. The computer program, when executed by a processor, implements the above-mentioned method for determining the multi-source aliasing acquisition parameters of land well guns.

[0058] The above technical solutions provided by the embodiments of the present disclosure have at least some or all of the following advantages compared with the prior art:

[0059] In the method for determining the multi-source aliasing acquisition parameters of land well guns according to the embodiment of the present disclosure, the target work area is divided into multiple target areas. For each target area, the initial T-D curve and gun coordinate information corresponding to the current target area are determined according to the three-dimensional seismic data of the current target area; according to the preset aliasing gun number and gun spacing, the gun coordinate information corresponding to the current target area is grouped, and based on the initial T-D curve, according to the distance between the main gun and the secondary gun in each group of gun coordinate information, the delay time of the secondary gun is determined, where each group of gun coordinate information includes the main gun coordinate information and the secondary gun coordinate information; according to the distance between the main gun and the secondary gun in each group of gun coordinate information and the delay time of the secondary gun, the three-dimensional seismic data corresponding to the main gun coordinate information and the secondary gun coordinate information in each group of gun coordinate information are respectively subjected to simulated aliasing processing to obtain the simulated aliasing data of the target area, and the simulated aliasing data of the target area is separated to determine the similarity between the data before aliasing and the data after aliasing separation for each gun in the simulated aliasing data; when the similarity is less than the preset threshold, return to the step of determining 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 based on the initial T-D curve until the similarity is greater than or equal to the preset threshold; when the similarity is greater than or equal to the preset threshold, update the initial T-D curve 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 to obtain the updated T-D curve, which is used as the multi-source aliasing acquisition parameters of the well guns in the target area. The mixed acquisition guns that meet this T-D curve can effectively separate the aliasing interference indoors, and the field mixed acquisition data becomes single-shot data, realizing efficient field mixed acquisition and high-quality indoor separation processing. Description of the Drawings

[0060] The accompanying drawings here are incorporated into the specification and form a part of this specification, showing embodiments in accordance with the present disclosure, and are used together with the specification to explain the principles of the present disclosure.

[0061] To more clearly illustrate the technical solutions in the embodiments of the present disclosure or the prior art, the following will briefly introduce the accompanying drawings required for use in the description of the embodiments or related technologies. Obviously, for those of ordinary skill in the art, without creative efforts, other drawings can also be obtained based on these drawings.

[0062] Figure 1 Schematically shows a schematic flow diagram of a method for determining multi-source aliasing acquisition parameters for land well shots according to an embodiment of the present disclosure;

[0063] Figure 2 Schematically shows a schematic diagram of an initial T-D curve according to an embodiment of the present disclosure;

[0064] Figure 3 Schematically shows a schematic diagram of 3D seismic data before aliasing according to an embodiment of the present disclosure;

[0065] Figure 4 Schematically shows a schematic diagram of simulated aliased data according to an embodiment of the present disclosure;

[0066] Figure 5 Schematically shows a schematic diagram of aliasing separation data according to an embodiment of the present disclosure;

[0067] Figure 6 Schematically shows a schematic diagram of the T-D curve after well shot multi-source acquisition update according to an embodiment of the present disclosure;

[0068] Figure 7 Schematically shows a structural block diagram of a device for determining multi-source aliasing acquisition parameters for land well shots according to an embodiment of the present disclosure;

[0069] Figure 8 Schematically shows a structural block diagram of an electronic device according to an embodiment of the present disclosure. Detailed implementation manners

[0070] To make the objectives, technical solutions, and advantages of the embodiments of the present disclosure clearer, the following will clearly and completely describe the technical solutions in the embodiments of the present disclosure with reference to the accompanying drawings in the embodiments of the present disclosure. Obviously, the described embodiments are some, but not all, of the embodiments of the present disclosure. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present disclosure without creative efforts fall within the scope of protection of the present disclosure.

[0071] See Figure 1, embodiments of the present disclosure provide a method for determining multi-source aliasing acquisition parameters for land well shots, including the following steps:

[0072] S1. Divide the target work area into multiple target areas. For each target area, determine the initial T-D curve and shot coordinate information corresponding to the current target area according to the 3D seismic data of the current target area.

[0073] In this embodiment, the target work area is divided into multiple target areas according to the excitation lithology and surface undulation height of the target work area.

[0074] S2. Group the shot coordinate information corresponding to the current target area according to the preset number of aliasing shots and shot spacing, and based on the initial T-D curve, determine the delay time of the secondary shot according to the distance between the main shot and the secondary shot in each group of shot coordinate information, where each group of shot coordinate information includes the main shot coordinate information and the secondary shot coordinate information.

[0075] S3. According to 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, perform simulated aliasing processing on the 3D seismic data corresponding to the main shot coordinate information and the secondary shot coordinate information in each group of shot coordinate information respectively to obtain the simulated aliasing data of the target area, and separate the simulated aliasing data of the target area 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. In the case where the similarity is less than the preset threshold, return to the step of determining the delay time of the secondary shot according to the distance between the main shot and the secondary shot in each group of shot coordinate information based on the initial T-D curve until the similarity is greater than or equal to the preset threshold.

[0077] S5. In the case where the similarity is greater than or equal to the preset threshold, update the initial T-D curve according to the distance between the main shot and the secondary shot corresponding to the current simulated aliasing and the delay time of the secondary shot to obtain the updated T-D curve, which is used as the multi-source aliasing acquisition parameters for well shots in the target area.

[0078] In this embodiment, in step S1, the determining the initial T-D curve corresponding to the current target area according to the 3D seismic data of the current target area includes:

[0079] Load the observation system of the current target area, and select the 3D seismic data with a collection time length exceeding the preset duration and a maximum offset distance greater than or equal to the preset distance from the historical 3D seismic data of the current target area as the target 3D seismic data of the current target area, where the preset duration can be 30s and the preset distance can be 12km;

[0080] Select the seismic data of the nearest offset receiving line of the target 3D seismic data, divide the selected seismic data into two parts according to positive and negative offsets, and select the data with slower energy attenuation in the two parts as the data to be processed;

[0081] Taking 1 s as the time window size and 1 km as the distance window size, divide the data to be processed into multiple time window data according to the time window size, where the time window data can be selected overlapping by half a time window;

[0082] Calculate the root mean square amplitude value of the data in each time window, take 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 reflection energy and background energy;

[0083] Pick up the energy attenuation trend curve in the energy attenuation transition zone between the effective reflection energy and the background energy, and draw the initial T-D curve according to the inflection point of the trend curve. Among them, according to the inflection point value of the initial T-D curve, set multiple T-D value ranges for verifying the subsequent aliasing separation effect.

[0084] See Figure 2 , in the interactive panel tool, set the color scale display range (5 - 20), and the root mean square amplitude values of different time windows are displayed in different colors. The energy attenuation transition zone is between the color of the strong effective reflection energy and the color of the background energy. Pick up the energy attenuation trend curve in the energy transition zone and draw the initial T-D curve value according to the inflection point of the trend curve.

[0085] In this embodiment, in step S1, the determining the gun coordinate information corresponding to the current target area according to the 3D seismic data of the current target area includes:

[0086] Load the observation system of the current target area, extract the gun gather data from the 3D seismic data of the current target area according to the gun file number, perform stacking processing on the data of multiple same geophone points in the gun gather data to obtain the stacked data, and extract the file number, x and y coordinates corresponding to the gun file number in the stacked data as the gun coordinate information;

[0087] Sort the gun coordinate information according to the gun file number, and delete the duplicate gun coordinate information to ensure the uniqueness of the gun coordinate information for each gun.

[0088] In this embodiment, the performing stacking processing on the data of multiple same geophone points in the gun gather data includes:

[0089] The data of multiple same geophone points in the gun gather data are stacked into 1 geophone point data, and the duplicate geophone point information is deleted.

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

[0091] First, each gun coordinate information that has not been grouped is marked as 0;

[0092] Second, 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, and the group number is 1;

[0093] 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 T-D value range as the first secondary gun of the main gun, and mark it as 1, and the group number is 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 T-D value range as the second secondary gun of the main gun, and mark it as 1, and the group number is 1, and record the distance value from the main gun. See Figure 2 , the distance of the first T-D value range is 1 - 3 km, and the distance of the second T-D value range is 2 - 4 km;

[0094] Fourth, if the aliasing gun number of group 1 is 3 guns, complete the grouping of group 1 according to the first to third steps. If the aliasing gun number of group 1 exceeds 3 guns, execute the step of selecting the gun coordinate information whose distance from the main gun is within the Nth T-D value range among the remaining gun coordinate information marked as 0 as the Nth secondary gun of the main gun until the grouping of group 1 is completed;

[0095] Fifth, the selected gun coordinate information in the previous group 1 does not participate in the selection of the next group. Select the gun coordinate information marked as 0 and ranked first as the main gun of the current group 2, and mark it as 1, and the group number is 2;

[0096] Sixth, among the remaining gun coordinate information marked as 0, select the gun coordinate information whose distance from the main gun is within the range of the first T-D value range shifted by a preset progressive distance as the first secondary gun of the main gun, and mark it as 1, and the group number is 2, 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 range of the second T-D value range shifted by a preset progressive distance as the second secondary gun of the main gun, and mark it as 1, and the group number is 2, and record the distance value from the main gun. Among them, when the first T-D value range is 1 - 3 km and the preset progressive distance is 2 km, the range of the first T-D value range shifted by the preset progressive distance is 3 - 5 km. When the second T-D value range is 2 - 4 km and the preset progressive distance is 2 km, the range of the second T-D value range shifted by the preset progressive distance is 4 - 6 km;

[0097] Step 7: If the aliased shot count of the second group is 3 shots, complete the grouping of the second group according to Steps 5 to 6. If the aliased shot count of the second group exceeds 3 shots, perform the step of selecting shot coordinate information within the range of the Nth T-D value from the main shot among the remaining shot coordinate information marked as 0 as the Nth secondary shot of the main shot until the grouping of the second group is completed;

[0098] Step 8: Complete the selection of shot coordinate information for each remaining group according to Steps 6 and 7. The distance between each group and the main shot progresses according to a preset progressive distance until the maximum offset is reached. Starting from the first T-D value range, if no shot coordinate information that meets the requirement of the first T-D value range of the distance from the main shot is found among the remaining shot coordinate information marked as 0, then only the shot coordinate information of the main shot is selected for this group, and there is no shot coordinate information of the secondary shot.

[0099] In this embodiment, in step S2, the determining the delay time of the secondary shot based on the initial T-D curve according to the distance between the main shot and the secondary shot in each group of shot coordinate information includes:

[0100] Obtain multiple T-D value ranges set according to the inflection point values of the initial T-D curve, where each T-D value range includes: the distance range between the main shot and the secondary shot in each group of shot coordinate information and the corresponding delay time of the secondary shot. Refer to Figure 2 If the distance range is less than 6 km, obtain a random delay time within the time range of 6 s - 15 s. If the distance range is greater than or equal to 6 km, obtain a random delay time within the time range of 1 s - 6 s;

[0101] Determine the T-D value range where the distance is located according to the distance between the main shot and the secondary shot in each group of shot coordinate information;

[0102] Based on the T-D value range where the distance is located, randomly obtain the delay time of the secondary shot from the delay times of the secondary shot corresponding to the distance range, where the delay time between the secondary shots in the same group is greater than 0.5 s.

[0103] In this embodiment, in step S3, the performing simulated aliasing processing on the three-dimensional seismic data corresponding to the main shot coordinate information and the secondary shot coordinate information in each group of shot coordinate information according to 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 to obtain the simulated aliased data of the target area includes:

[0104] First step, store the trace data and trace file number in the three-dimensional seismic data that match the main shot file number in each group of shot coordinate information into the current main shot coordinate information, where the trace data and trace file number that match the main shot file number in each group of shot coordinate information are obtained through all shot line data and trace header information on the trace gather data in the three-dimensional seismic data;

[0105] Step 2: Store the trace data and trace file numbers in the 3D seismic data that match the secondary shot file numbers in each group of shot coordinate information into the current secondary shot coordinate information;

[0106] Step 3: If there are no trace data and trace file numbers in the 3D seismic data that match the primary shot file number or the secondary shot file number, then there are no trace data and trace file numbers in the shot coordinate information of this primary shot or secondary shot;

[0107] Step 4: For each group of shot coordinate information, when there are no trace data and trace file numbers in at least one of the primary shot coordinate information and the secondary shot coordinate information in the current group of shot coordinate information, do not perform primary-secondary shot aliasing processing; when there are trace data and trace file numbers in both the primary shot coordinate information and the secondary shot coordinate information in the current group of shot coordinate information, perform the following primary-secondary shot aliasing processing on the current group:

[0108] Step 5: Take out all the sample values of the traces that match in the current group of primary shots, use the delay time of the first secondary shot in the current group as the start time of the sample values of the first secondary shot, perform staggered addition with all the sample values of the current primary shot, use the delay time of the second secondary shot in the current group as the start time of the sample values of the second secondary shot, and perform staggered addition again with the added value of the current primary shot and the first secondary shot until all the secondary shots in the current group have been subjected to staggered addition, 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 values of the primary and secondary shots respectively according to the actual sampling time before aliasing. According to the trace file number matched by the primary shot in the current group, use the truncated sample data values of the primary shot to replace the sample data values before aliasing; the truncation start time of the first secondary shot is based on its delay time, truncate the sample data of the same length, and use the trace file number matched by the first secondary shot in the current group to replace the sample data values before aliasing; the truncation start time of the second secondary shot is based on its delay time, truncate the sample data of the same length, and use the trace file number matched by the second secondary shot in the current group to replace the sample data values before aliasing until all the secondary shots have completed the truncation of the aliased data, thus completing the truncation of the aliased data of the current group;

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

[0111] In this embodiment, in step S3, the separation of the simulated aliased data of the target area includes:

[0112] Since most of the onshore well shots are preprocessed in the form of shot gathers, for the simulated aliased data of the target area, search for the unaliased shots at adjacent positions as model shots according to the shot coordinate information of the simulated aliased data, and perform multi-scale curvelet transform on the aliased shots and the model shots in the shot domain;

[0113] In each scale, structural pattern recognition and threshold separation processing are performed on aliased shots and model shots, and then the separated data is multi-scale curvelet inverse transformed to shot domain data to complete the separation of aliased data.

[0114] In this embodiment, in step S3, the quality control of the separation effect is to calculate the aliased data acquisition points as the basis for the actual T-D curve values, and perform similarity analysis on the data before and after aliasing to perform quality control. Through the following expression, the similarity between the data before aliasing and the data after aliasing separation for each shot in the simulated aliased data is determined:

[0115]

[0116] where SE(x,y) is the similarity between the data before aliasing and the data after aliasing separation for each shot in the simulated aliased data, x and y are the data before and after aliasing respectively, Ax and Ay are the Fourier amplitudes of the data x and y, and represent the Fourier phases of the data x and y respectively, and N is the number of shots.

[0117] In this embodiment, using Figure 3 the receiver gather of 3D data, simulated aliased data with different distances and times as shown in Figure 4 is generated according to the method for determining the multi-source aliased acquisition parameters of land well shots disclosed in this disclosure, and aliasing separation is performed to obtain the data as shown in Figure 5 shown.

[0118] Through similarity analysis, Figure 3 the similarity between the data and Figure 5 the data is 98.3%, meeting the standard of aliasing separation greater than 95%. Through the separated data that meets the standard, the obtained pairs of different distances and times are used to plot the well shot multi-source acquisition T-D parameter curve as shown in Figure 6 shown.

[0119] The method for determining the multi-source aliasing acquisition parameters of land well guns in the present disclosure collects the single-shot data acquired in the past in this area, extracts the gun coordinate information after preprocessing, groups them according to the number of aliasing guns and the distance between guns, and obtains the delay time of each secondary gun for the grouped gun coordinate information within the excitation time range; performs data matching and simulated aliasing processing on the grouped gun coordinate information and the old data acquired in the past; performs aliasing interference separation processing on the simulated aliasing data, and conducts similarity analysis on the separated data; based on the time and distance information with better separation effect, draws the T-D curve of this area as the basis for field mixed acquisition construction parameters. Only the mixed acquisition guns that meet this T-D curve can be effectively separated from aliasing interference indoors, and the field mixed acquisition data can become raw data. It can conduct a relatively complete argumentation on the mixed acquisition parameters for the complex well gun work area in the mountains. The generated T-D mixed acquisition curve can meet the requirements of high-efficiency field mixed acquisition construction and high-quality separation processing indoors.

[0120] The aliasing separation step in the present disclosure can collect the aliasing gun set data through the T-D curve, record the sequence number of each gun trace, extract and retain the abnormal traces and empty traces of each gun data, perform multi-scale transformation on the aliasing data and the model data, respectively obtain their scale coefficients, use the scale coefficients of the aliasing data and the model data to obtain the model feature coefficients of different scales between the two, set a threshold for filtering processing, perform multi-scale inverse transformation on the model feature coefficients after filtering processing, and merge the obtained aliasing data with the extracted abnormal traces and empty traces to obtain the finally separated main gun data and aliasing data, which is convenient for subsequent high-precision processing of seismic data.

[0121] See Figure 7 , the embodiment of the present disclosure provides a device for determining multi-source aliasing acquisition parameters of land well guns, including:

[0122] A determination module 11, configured to divide the target work area into multiple target regions. For each target region, determine the initial T-D curve and gun coordinate information corresponding to the current target region according to the three-dimensional seismic data of the current target region;

[0123] A grouping module 12, configured to group the gun coordinate information corresponding to the current target region according to the preset number of aliasing guns and the gun spacing, and based on the initial T-D 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, where each group of gun coordinate information includes the main gun coordinate information and the secondary gun coordinate information;

[0124] An aliasing module 13 is configured to perform simulated aliasing processing on the 3D seismic data corresponding to the main gun coordinate information and the secondary gun coordinate information in each group 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 group of gun coordinate information, obtain the simulated aliased data of the target area, and separate the simulated aliased data of the target area to determine the similarity between the data before aliasing and the data after aliasing separation for each gun in the simulated aliased data;

[0125] A return module 14 is configured to, when the similarity is less than a preset threshold, return the step of determining the delay time of the secondary gun based on the initial T-D curve according to the distance between the main gun and the secondary gun in each group of gun coordinate information until the similarity is greater than or equal to the preset threshold;

[0126] An update module 15 is configured to, when the similarity is greater than or equal to the preset threshold, update the initial T-D curve according to the distance between the main gun and the secondary gun and the delay time of the secondary gun corresponding to the current simulated aliasing, obtain the updated T-D curve, and use it as the multi-source aliasing acquisition parameter of the well-gun in the target area.

[0127] The implementation processes of the functions and effects of each unit in the above device are specifically described in detail in the implementation processes of the corresponding steps in the above method, and will not be elaborated here.

[0128] For the device embodiment, since it basically corresponds to the method embodiment, the relevant parts can refer to the partial description of the method embodiment. The device embodiments described above are only illustrative. The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed to multiple network units. Some or all of the modules can be selected according to actual needs to achieve the purpose of the solution of the present invention. Those of ordinary skill in the art can understand and implement it without creative efforts.

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

[0130] Referring to Figure 8 As shown, the electronic device provided by the embodiments of the present disclosure includes a processor 1110, a communication interface 1120, a memory 1130, and a communication bus 1140. Among them, the processor 1110, the communication interface 1120, and the memory 1130 complete communication with each other through the communication bus 1140;

[0131] The memory 1130 is used to store computer programs;

[0132] When the processor 1110 is used to execute the program stored on the memory 1130, the method for determining the multi-source aliasing acquisition parameters of land well guns as shown below is implemented:

[0133] The target work area is divided into multiple target areas. For each target area, the initial T-D curve and gun coordinate information corresponding to the current target area are determined according to the three-dimensional seismic data of the current target area;

[0134] According to the preset number of aliasing guns and gun spacing, the gun coordinate information corresponding to the current target area is grouped, and based on the initial T-D 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, where each group of gun coordinate information includes the main gun coordinate information and the secondary gun coordinate information;

[0135] According to the distance between the main gun and the secondary gun in each group of gun coordinate information and the delay time of the secondary gun, perform simulated aliasing processing on the three-dimensional seismic data corresponding to the main gun coordinate information and the secondary gun coordinate information in each group of gun coordinate information to obtain the simulated aliasing data of the target area, and separate the simulated aliasing data of the target area to determine the similarity between the data before aliasing and the data after aliasing separation for each gun in the simulated aliasing data;

[0136] In the case where the similarity is less than the preset threshold, return to the step of determining 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 based on the initial T-D curve until the similarity is greater than or equal to the preset threshold;

[0137] In the case where the similarity is greater than or equal to the preset threshold, update the initial T-D curve 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 to obtain the updated T-D curve, which is used as the multi-source aliasing acquisition parameter of the well-gun in the target area.

[0138] The above communication bus 1140 can be a Peripheral Component Interconnect (PCI) bus or an Extended Industry Standard Architecture (EISA) bus, etc. The communication bus 1140 can be divided into an address bus, a data bus, a control bus, etc. For the sake of simplicity, only a thick line is shown in the figure, but it does not mean that there is only one bus or one type of bus.

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

[0140] The memory 1130 can include a Random Access Memory (RAM), and can also include a non-volatile memory, such as at least one disk memory. Optionally, the memory 1130 can also be at least one storage device located far from the aforementioned processor 1110.

[0141] The above-mentioned processor 1110 may be a general-purpose processor, including a Central Processing Unit (CPU), a Network Processor (NP), etc.; it may 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, discrete hardware components.

[0142] Embodiments of the present disclosure also provide a computer-readable storage medium. A computer program is stored on the above-mentioned computer-readable storage medium, and when the computer program is executed by a processor, the method for determining the multi-source aliasing acquisition parameters of a land well gun as described above is implemented.

[0143] The computer-readable storage medium may be included in the device / apparatus described in the above embodiments; it may also exist alone without being assembled into the device / apparatus. The above-mentioned computer-readable storage medium carries one or more programs, and when the one or more programs are executed, the method for determining the multi-source aliasing acquisition parameters of a land well gun according to the embodiments of the present disclosure is implemented.

[0144] According to the embodiments of the present disclosure, the computer-readable storage medium may be a non-volatile computer-readable storage medium, for example, it may include but is not limited to: portable computer disks, hard disks, random access memories (RAMs), read-only memories (ROMs), erasable programmable read-only memories (EPROMs or flash memories), portable compact disk read-only memories (CD-ROMs), optical storage devices, magnetic storage devices, or any suitable combination of the above. In the present disclosure, the computer-readable storage medium may be any tangible medium that contains or stores a program, and the program may be used by or in combination with an instruction execution system, device, or device.

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

[0146] The above are only specific embodiments of the present disclosure, enabling those skilled in the art to understand or implement the present disclosure. Various modifications to these embodiments will be obvious to those skilled in the art, and the general principles defined herein can be implemented in other embodiments without departing from the spirit or scope of the present disclosure. Therefore, the present disclosure will not be limited to these embodiments shown herein, but rather will be accorded the widest scope consistent with the principles and novel features claimed herein.

Claims

1. A method for determining multi-source aliasing acquisition parameters of land well guns, characterized in that The method includes: Dividing the target work area into multiple target regions. For each target region, determining the initial T-D curve and gun coordinate information corresponding to the current target region according to the 3D seismic data of the current target region; Grouping the gun coordinate information corresponding to the current target region according to the preset aliased gun number and gun spacing, and based on the initial T-D curve, determining 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, where each group of gun coordinate information includes the main gun coordinate information and the secondary gun coordinate information; According to the distance between the main gun and the secondary gun in each group of gun coordinate information and the delay time of the secondary gun, performing simulated aliasing processing on the 3D seismic data corresponding to the main gun coordinate information and the secondary gun coordinate information in each group of gun coordinate information respectively to obtain the simulated aliased data of the target region, and separating the simulated aliased data of the target region to determine the similarity between the data before aliasing and the data after aliasing separation for each gun in the simulated aliased data; In the case where the similarity is less than the preset threshold, return to the step of determining 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 based on the initial T-D curve until the similarity is greater than or equal to the preset threshold; In the case where the similarity is greater than or equal to the preset threshold, update the initial T-D curve 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 to obtain the updated T-D curve, which is used as the multi-source aliased acquisition parameter of the well gun in the target region.

2. The method according to claim 1, characterized in that The determining the initial T-D curve corresponding to the current target region according to the 3D seismic data of the current target region includes: Selecting, from the historical 3D seismic data of the current target region, the 3D seismic data with the acquisition time length exceeding the preset time length and the maximum offset distance greater than or equal to the preset distance as the target 3D seismic data of the current target region; Selecting the seismic data of the nearest offset receiving line of the target 3D seismic data, dividing the selected seismic data into two parts according to the positive and negative offsets, and selecting the data with slower energy attenuation in the two parts as the data to be processed; Taking 1 s for time and 1 km for distance as the time window size, dividing the data to be processed into multiple time window data according to the time window size, where the time window data is selected with an overlap of half a time window; Calculating the root mean square amplitude value of the data in each time window, taking the current time window position as the coordinate of the current root mean square amplitude value, and replacing the time window data with the root mean square amplitude value, displaying the root mean square amplitude value of each time window, and determining the effective reflection energy and the background energy; Picking up the energy attenuation trend curve in the energy attenuation transition zone between the effective reflection energy and the background energy, and drawing the initial T-D curve according to the inflection point of the trend curve, where multiple T-D value ranges are set according to the inflection point value of the initial T-D curve.

3. The method according to claim 1, wherein The determining the gun coordinate information corresponding to the current target region according to the 3D seismic data of the current target region includes: Extract shot gather data from the 3D seismic data of the current target area according to the shot file number, perform stacking processing on the data of multiple identical geophone points 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; Sort the shot coordinate information according to the shot file number, and delete the duplicate shot coordinate information to ensure the uniqueness of the shot coordinate information for each shot.

4. The method according to claim 1, wherein Group the shot coordinate information corresponding to the current target area according to the preset aliased shot number and shot spacing, including: First step, mark each ungrouped shot coordinate information as 0; Second step, select the shot coordinate information marked as 0 and ranked first as the main shot of the current group 1, mark it as 1, and the group number is 1; Third step, among the remaining shot coordinate information marked as 0, select the shot coordinate information whose distance from the main shot is within the first T-D value range as the first sub-shot of the main shot, mark it as 1, and the group number is 1, and record the distance value from the main shot; among the remaining shot coordinate information marked as 0, select the shot coordinate information whose distance from the main shot is within the second T-D value range as the second sub-shot of the main shot, mark it as 1, and the group number is 1, and record the distance value from the main shot; Fourth step, if the aliased shot number of the first group is 3 shots, complete the grouping of the first group according to the first to third steps. If the aliased shot number of the first group exceeds 3 shots, perform the step of selecting the shot coordinate information whose distance from the main shot is within the Nth T-D value range as the Nth sub-shot of the main shot among the remaining shot coordinate information marked as 0 until the grouping of the first group is completed; Fifth step, select the shot coordinate information marked as 0 and ranked first as the main shot of the current group 2, mark it as 1, and the group number is 2; Sixth step, among the remaining shot coordinate information marked as 0, select the shot coordinate information whose distance from the main shot is within the range after shifting the preset progressive distance of the first T-D value range as the first sub-shot of the main shot, mark it as 1, and the group number is 2, and record the distance value from the main shot; among the remaining shot coordinate information marked as 0, select the shot coordinate information whose distance from the main shot is within the range after shifting the preset progressive distance of the second T-D value range as the second sub-shot of the main shot, mark it as 1, and the group number is 2, and record the distance value from the main shot; Seventh step, if the aliased shot number of the second group is 3 shots, complete the grouping of the second group according to the fifth to sixth steps. If the aliased shot number of the second group exceeds 3 shots, perform the step of selecting the shot coordinate information whose distance from the main shot is within the Nth T-D value range as the Nth sub-shot of the main shot among the remaining shot coordinate information marked as 0 until the grouping of the second group is completed; Eighth step, complete the selection of the shot coordinate information for each remaining group according to the sixth and seventh steps. The distance between each group and the main shot progresses according to the preset progressive distance. When the maximum offset is reached, starting from the first T-D value range, if no shot coordinate information that meets the requirement of the first T-D value range of the distance from the main shot is found among the remaining shot coordinate information marked as 0, then only the shot coordinate information of the main shot is selected for this group, and there is no shot coordinate information of the sub-shot.

5. The method according to claim 1, wherein Based on the initial T-D curve, determining the delay time of the secondary gun according to the distance between the main gun and the secondary gun in each set of gun coordinate information, includes: Obtaining a plurality of T-D value ranges set according to the inflection point values of the initial T-D curve, where each T-D 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; Determining the T-D value range where the distance is located according to the distance between the main gun and the secondary gun in each set of gun coordinate information; Based on the T-D value range where the distance is located, randomly obtaining the delay time of the secondary gun from the delay times of the secondary gun corresponding to the distance range, where the delay time between secondary guns in the same group is greater than 0.5 s.

6. The method according to claim 1, wherein According to the distance between the main gun and the secondary gun in each set of gun coordinate information and the delay time of the secondary gun, 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 to obtain the simulated aliased data of the target area, includes: First step, storing the trace data and trace file numbers in the three-dimensional seismic data that match the main gun file number in each set of gun coordinate information into the current main gun coordinate information; Second step, storing the trace data and trace file numbers in the three-dimensional seismic data that match the secondary gun file number in each set of gun coordinate information into the current secondary gun coordinate information; Third step, if there is no trace data and trace file number in the three-dimensional seismic data that match the main gun file number or the secondary gun file number, then there is no trace data and trace file number in the gun coordinate information of the main gun or the secondary gun; Fourth step, for each set of gun coordinate information, when there is no trace data and trace file number in at least one of the main gun coordinate information and the secondary gun coordinate information in the current set of gun coordinate information, the main gun-secondary gun aliasing processing is not performed; when there are trace data and trace file numbers matching in both the main gun coordinate information and the secondary gun coordinate information in the current set of gun coordinate information, the following main gun-secondary gun aliasing processing is performed on the current set: Fifth step, taking out all the sample point values in the current set of main guns that match the trace data, using the delay time of the first secondary gun in the current set as the start time of the first secondary gun sample point values, adding them out of phase with all the sample point values of the current main gun, using the delay time of the second secondary gun in the current set as the start time of the second secondary gun sample point values, and adding them out of phase with the added value of the current main gun and the first secondary gun again until all the secondary guns in the current set are added out of phase, completing the data aliasing of the current set; Sixth step, performing truncation processing on the aliased data of the current set, obtaining the main and secondary gun sample point values respectively according to the actual sampling time before aliasing, and using the truncated main gun sample data values to replace the sample data values before aliasing according to the trace file number matched by the main gun in the current set; the truncation start time of the first secondary gun intercepts the sample data of the same length according to its delay time, and uses the trace file number matched by the first secondary gun in the current set to replace the sample data values before aliasing; the truncation start time of the second secondary gun intercepts the sample data of the same length according to its delay time, and uses the trace file number matched by the second secondary gun in the current set to replace the sample data values before aliasing until all the secondary guns have completed the truncation of the aliased data, completing the truncation of the aliased data of the current set; In the seventh step, after aliasing and truncating the data of the current group, perform aliasing and truncating on the next group of data according to the fifth and sixth steps until aliasing and truncating are completed for all groups, and obtain the simulated aliased data of the target area.

7. The method according to claim 1, characterized in that The separation of the simulated aliased data of the target area includes: For the simulated aliased data of the target area, search for the unaliased guns at adjacent positions according to the gun coordinate information of the simulated aliased data as the model guns, and perform multi-scale curvelet transform on the aliased guns and the model guns in the gun domain; In each scale, perform structural pattern recognition and threshold separation processing on the aliased guns and the model guns, and then perform multi-scale curvelet inverse transform on the separated data to the gun domain data to complete the separation of the aliased data.

8. The method according to claim 1, wherein Through the following expression, determine the similarity between the data before aliasing and the data after aliasing separation for each gun in the simulated aliased data: Among them, SE(x, y) is the similarity between the data before aliasing and the data after alias separation for each shot in the simulated aliased data. x and y are the data before aliasing and the data after separation respectively. Ax and Ay are the Fourier amplitudes of data x and y respectively, and represent the Fourier phases of data x and y respectively, and N is the number of shots.

9. A device for determining acquisition parameters of multi-source aliasing in land well shooting, characterized in that Including: A determination module for dividing the target work area into multiple target areas. For each target area, determine the initial T-D curve and gun coordinate information corresponding to the current target area according to the 3D seismic data of the current target area; A grouping module for grouping the gun coordinate information corresponding to the current target area according to the preset number of aliased guns and gun spacing, and based on the initial T-D 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, where each group of gun coordinate information includes the main gun coordinate information and the secondary gun coordinate information; An aliasing module for performing simulated aliasing processing on the 3D seismic data corresponding to the main gun coordinate information and the secondary gun coordinate information in each group of gun coordinate information according to the distance between the main gun and the secondary gun in each group of gun coordinate information and the delay time of the secondary gun, obtain the simulated aliased data of the target area, and separate the simulated aliased data of the target area to determine the similarity between the data before aliasing and the data after aliasing separation for each gun in the simulated aliased data; A return module for, when the similarity is less than the preset threshold, returning to the step of determining 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 based on the initial T-D curve until the similarity is greater than or equal to the preset threshold; An update module for, when the similarity is greater than or equal to the preset threshold, updating the initial T-D curve 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 to obtain the updated T-D curve as the well-gun multi-source aliasing acquisition parameter of the target area.

10. An electronic device, characterized in that, Including a processor, a communication interface, a memory, and a communication bus. Among them, the processor, the communication interface, and the memory complete mutual communication through the communication bus; The memory is used to store computer programs; The processor, when executing the programs stored on the memory, implements the method for determining the well-gun multi-source aliasing acquisition parameter 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 well-gun multi-source aliasing acquisition parameter described in any one of claims 1-8.

Images