A two-dimensional survey line first and last shot point automatic full coverage fast and batch calculation method

Abstract

The application discloses a kind of two-dimensional survey line first and last shot point automatic full coverage fast, batch calculation method, steps are as follows: step 1, establish all two-dimensional survey line preset parameter;Step 2, calculate a survey line azimuth, determine full end point starting point position;Step 3, determine first shot point position, calculate additional section length, first geophone point coordinate;Step 4, calculate shot number, construction length, geophone point number, calculate last shot point, geophone point coordinate;Step 5, obtain all two-dimensional survey line first and last shot point coordinate by loop.The application combines geophysical exploration industry standard and production practice, uses VB computer programming language, and innovatively realizes a kind of two-dimensional survey line first and last shot point automatic full coverage fast, batch calculation methods, so that technical personnel is freed from the inconvenience caused by manually designing each two-dimensional survey line first and last shot point position, and realizes the efficient layout of exploration design.

Description

Technical Field

[0001] This invention relates to a rapid, batch calculation method for automatic full coverage of shot detection points at the beginning and end of a two-dimensional seismic survey line, belonging to the field of seismic exploration and development technology. Background Technology

[0002] In 2D seismic exploration design, determining the coordinates of the first and last shot-receiver points is a crucial step. However, when calculating these coordinates, the process often involves extrapolating from smaller to larger coordinate directions. Given the full-scale deployment coordinates, smaller coordinates cannot always be prioritized, requiring significant effort to determine the first and last coordinates. Furthermore, human factors (skill level, work attitude, etc.) can easily lead to inaccurate calculations, necessitating repeated checks and increasing workload. A search of relevant patents reveals limited methods for rapid, batch calculation of first and last shot-receiver points for 2D seismic lines in China. For example, methods and devices for locating shot point coordinates (CN201911302675.0) and methods and devices for obtaining 2D seismic lines (CN201410778762.4) are only applicable to the later stages of 2D seismic line construction and mostly lack the capability for batch calculation of first and last shot-receiver points. Therefore, existing calculation methods are mostly cumbersome and cannot quickly and efficiently calculate the coordinates of all first and last shot-receiver points for all 2D seismic lines, thus affecting the efficiency of exploration design.

[0003] Therefore, it is necessary to develop a method that can automatically, quickly, and in batches calculate the coordinates of the first and last shot checkpoints of a two-dimensional survey line, thereby improving the efficiency and accuracy of two-dimensional survey line exploration design. Summary of the Invention

[0004] The purpose of this invention is to solve the problems existing in the prior art and to provide a method for automatic full coverage and rapid batch calculation of two-dimensional survey line start and end shot detection points.

[0005] To achieve the above objectives, the technical solution adopted by the present invention is: a method for automatic full-coverage rapid batch calculation of two-dimensional survey line start and end shot detection points, which includes the following steps:

[0006] Step 1: Establish preset parameters for all two-dimensional survey lines;

[0007] Step 2: Calculate the azimuth of a survey line and determine the starting point of the full-scale endpoint;

[0008] Step 3: Determine the location of the first shot, and calculate the length of the additional section and the coordinates of the first receiver point;

[0009] Step 4: Calculate the number of shots, construction length, and number of geophones; calculate the coordinates of the final shot point and geophones.

[0010] Step 5: Obtain the coordinates of all the first and last shot checkpoints of the two-dimensional survey line by looping.

[0011] Furthermore, in step 1, the preset parameters of the two-dimensional survey line include the maximum shot-receiver distance, minimum shot-receiver distance, full-cycle endpoint coordinates, full-cycle coverage count, shot-point distance, and receiver-point distance for all survey lines.

[0012] Furthermore, in step 2, the azimuth angle of a survey line is calculated based on the coordinates of the full-scale endpoint established in step 1, and the direction of the survey line is determined according to the range of the azimuth angle, thereby determining the starting point position of the full-scale endpoint.

[0013] Furthermore, in step 3, the firing point position of the Mth shot is determined based on the number of full coverages and the starting point position of the full coverage endpoint from steps 1 and 2. Then, the distance LPD from the first firing point to the starting point of the full coverage endpoint is calculated, thereby determining the position of the first firing point. The coordinates of the first firing point are calculated using the coordinates of the starting point of the full coverage endpoint and the azimuth angle; specifically as follows:

[0014] M = number of full coverage times / 2, where M is the first shot after the starting point of the full coverage endpoint;

[0015] LPD = 0.5 * minimum shot-receiver distance + (M-1) * shot distance;

[0016] The X-coordinate of the first shot point = the X-coordinate of the starting point of the full-scale endpoint – LPD × Sin (azimuth angle);

[0017] The Y-coordinate of the first shot point = the Y-coordinate of the starting point of the full-scale endpoint – LPD × Cos (azimuth).

[0018] Furthermore, in step 3, the length of the additional segment is the distance from the first receiver point to the starting point of the full-scale endpoint. The coordinates of the first receiver point are then calculated based on the length of the additional segment, as follows:

[0019] Additional section length = LPD + maximum shot-receiver distance;

[0020] The X-coordinate of the first detection point = X-coordinate of the starting point of the full-scale endpoint – additional segment length × Sin (azimuth); The Y-coordinate of the first detection point = Y-coordinate of the starting point of the full-scale endpoint – additional segment length × Cos (azimuth).

[0021] Furthermore, in step 4, the full-scale length is calculated based on the coordinates of the full-scale endpoint established in step 1, and then the number of shots is calculated. Finally, the coordinates of the final shot are calculated by combining the coordinates of the first shot point from step 3, as detailed below:

[0022] Number of shots = ROUNDUP(full length / shot distance + 1) + number of coverages, where ROUNDUP is the round-up function;

[0023] Shot point length = (number of shots - 1) × shot point distance;

[0024] The X-coordinate of the final shot point = the X-coordinate of the first shot point + the length of the shot point × Sin (azimuth angle);

[0025] The Y-coordinate of the final shot point = the Y-coordinate of the first shot point + the length of the shot point × Cos (azimuth).

[0026] Furthermore, in step 4, the construction length is the distance from the first receiver point to the last receiver point. The number of receiver points is calculated, and combined with the coordinates of the first receiver point in step 3, the coordinates of the last receiver point are calculated, as follows:

[0027] Construction length = 2 * maximum shot-receiver distance + (number of shots - 1) * shot point distance;

[0028] Number of detector points = Construction length / Detector point spacing + 1;

[0029] The X-coordinate of the final receiver point = the X-coordinate of the first receiver point + the construction length × Sin (azimuth angle);

[0030] The Y-coordinate of the final receiver point = the Y-coordinate of the first receiver point + the construction length × Cos (azimuth).

[0031] Furthermore, in step 5, the loop calculates the coordinates of the first and last shot checkpoints for all two-dimensional survey lines using the methods from steps 2 to 4.

[0032] Compared with the prior art, the beneficial effects of the present invention are as follows:

[0033] This invention combines geophysical exploration industry standards and production practices, and uses VB computer programming language to innovatively realize a fast and batch calculation method for automatic full coverage of the first and last shot checkpoints of two-dimensional survey lines. This frees technicians from the inconvenience of manually designing the positions of the first and last shot checkpoints of each two-dimensional survey line, and realizes efficient layout of exploration design. Attached Figure Description

[0034] Figure 1 This is a flowchart illustrating a specific embodiment of the present invention;

[0035] Figure 2 This is a schematic diagram showing the relative positions of the full-score endpoint, shot point, and receiver point in this invention.

[0036] Figure 3 This is a schematic diagram illustrating the calculation of the additional segment length in this invention;

[0037] Figure 4 This is a schematic diagram illustrating the calculation of the number of shots in this invention. Detailed Implementation

[0038] The present invention will now be described in detail with reference to the accompanying drawings and specific embodiments.

[0039] like Figure 1 A method for automatic full-coverage, rapid, and batch calculation of two-dimensional survey line start and end shot detection points includes the following steps:

[0040] Step 1: Establish preset parameters for all two-dimensional survey lines;

[0041] Step 2: Calculate the azimuth of a survey line and determine the starting point of the full-scale endpoint;

[0042] Step 3: Determine the location of the first shot, and calculate the length of the additional section and the coordinates of the first receiver point;

[0043] Step 4: Calculate the number of shots, construction length, and number of geophones, and calculate the coordinates of the final shot point and geophones.

[0044] Step 5: Obtain the coordinates of all the first and last shot checkpoints of the two-dimensional survey line by looping.

[0045] The following are several specific embodiments of the application of the present invention.

[0046] Example 1:

[0047] In a specific embodiment 1 of the present invention, such as Figure 1 The diagram shown is a flowchart of a method for automatic full-coverage rapid batch calculation of two-dimensional survey line start and end shot detection points according to the present invention.

[0048] In step 1, the preset parameters for the two-dimensional survey lines include the maximum shot-receiver distance, minimum shot-receiver distance, full-cycle endpoint coordinates, full-cycle coverage count, shot-point distance, and receiver-point distance for all survey lines.

[0049] In step 2, the azimuth angle of a survey line is calculated based on the coordinates of the full-scale endpoint established in step 1, and the direction of the survey line is determined according to the range of the azimuth angle, thereby determining the starting point position of the full-scale endpoint.

[0050] In step 3, the position of the Mth shot is determined based on the number of full coverages and the starting point of the full coverage endpoint from steps 1 and 2. Then, the distance LPD from the first shot to the starting point of the full coverage endpoint is calculated, thus determining the position of the first shot. The coordinates of the first shot are calculated using the coordinates of the starting point of the full coverage endpoint and the azimuth angle, as detailed below:

[0051] M = number of full coverage times / 2, where M is the first shot after the starting point of the full coverage endpoint, and the shot is the shot with the first shot count as the first half of the full coverage times; for example, if the full coverage times are 160 times, the corresponding M shot is the 80th shot.

[0052] LPD = 0.5 * minimum shot-receiver distance + (M-1) * shot distance;

[0053] The X-coordinate of the first shot point = the X-coordinate of the starting point of the full-scale endpoint – LPD × Sin (azimuth angle);

[0054] The Y-coordinate of the first shot point = the Y-coordinate of the starting point of the full-scale endpoint – LPD × Cos (azimuth).

[0055] In step 3, the additional segment length refers to the distance from the first receiver point to the starting point of the full-scale endpoint. The coordinates of the first receiver point are then calculated based on the additional segment length, as follows:

[0056] Additional section length = LPD + maximum shot-receiver distance;

[0057] X-coordinate of the first detection point = X-coordinate of the starting point of the full-scale endpoint – Length of the additional segment × Sin(azimuth angle);

[0058] Y-coordinate of the first detection point = Y-coordinate of the starting point of the full-scale endpoint – Length of the additional segment × Cos(azimuth).

[0059] In step 4, the full-scale length is calculated based on the coordinates of the endpoint established in step 1, and then the number of shots is calculated. Combined with the coordinates of the first shot point in step 3, the coordinates of the last shot point are calculated, as follows:

[0060] Number of shots = ROUNDUP(full length / shot distance + 1) + number of coverages, where ROUNDUP is the round-up function;

[0061] Shot point length = (number of shots - 1) × shot point distance;

[0062] The X-coordinate of the final shot point = the X-coordinate of the first shot point + the length of the shot point × Sin (azimuth angle);

[0063] The Y-coordinate of the final shot point = the Y-coordinate of the first shot point + the length of the shot point × Cos (azimuth).

[0064] In step 4, the construction length refers to the distance from the first receiver point to the last receiver point. The number of receiver points is then calculated, and combined with the coordinates of the first receiver point from step 3, the coordinates of the last receiver point are calculated, as follows:

[0065] Construction length = 2 * maximum shot-receiver distance + (number of shots - 1) * shot point distance;

[0066] Number of detector points = Construction length / Detector point spacing + 1;

[0067] The X-coordinate of the final receiver point = the X-coordinate of the first receiver point + the construction length × Sin (azimuth angle);

[0068] The Y-coordinate of the final receiver point = the Y-coordinate of the first receiver point + the construction length × Cos (azimuth).

[0069] In step 5, the loop is to calculate the coordinates of the first and last shot checkpoints for all two-dimensional survey lines using the methods from steps 2 to 4.

[0070] Example 2:

[0071] In a specific embodiment 2 of the present invention, such as Figure 2The diagram shown illustrates the relative positions of the full-scale endpoint, shot point, and receiver point; as shown... Figure 3 The diagram shown illustrates the calculation of the additional segment length. Figure 4 The diagram shown illustrates the calculation of the number of shots. The specific implementation includes the following steps:

[0072] Step 1: Establish all preset parameters for the two-dimensional survey lines, including the maximum shot-receiver distance of 6390m, the minimum shot-receiver distance of 10m, the coordinates of the endpoint of the full-cycle survey line 1 (0, 9555) and (0, 11185), the coordinates of the full-cycle survey line 2 (9555, 0) and (111855, 0), the number of full-cycle coverages of 160, the shot-point distance of 40m, and the receiver-point distance of 20m, etc.

[0073] Step 2: Calculate the azimuth angle of survey line 1 to 0 degrees. According to the range of the azimuth angle, survey line 1 is determined to be a north-south survey line. Therefore, the south direction of the survey line is the starting direction. The starting point of the full-time endpoint is (0, 9555).

[0074] Step 3: Based on the full coverage count of 160 and the starting point of the full coverage endpoint (0, 9555), determine the firing point position (0, 9390) for the M=80th shot. Then, calculate the distance LPD from the first firing point to the starting point of the full coverage endpoint as 3165m. Using the coordinates of the starting point of the full coverage endpoint (0, 9555) and the azimuth angle of 0 degrees, calculate the position of the first firing point (0, 6390), as follows:

[0075] M = 160 / 2 = 80;

[0076] LPD=0.5*10+(80-1)*40=3165;

[0077] The X coordinate of the first shot point = the X coordinate of the starting point of the full-scale endpoint – LPD × Sin (azimuth angle) = 0 - 3165 × Sin (0°) = 0;

[0078] The Y-coordinate of the first shot point = Y-coordinate of the starting point of the full-scale endpoint - LPD × Cos (azimuth angle) = 9555 - 3165 × Cos (0°) = 6390.

[0079] The additional segment length refers to the distance from the first receiver point to the starting point of the full-pass endpoint, which is 9555m. The coordinates of the first receiver point (0, 0) are then calculated based on the additional segment length, as follows:

[0080] Additional segment length = 3165 + 6390 = 9555;

[0081] X-coordinate of the first detection point = X-coordinate of the starting point of the full-scale endpoint – Length of the additional segment × Sin(azimuth) = 0 - 9555 × Sin(0°) = 0;

[0082] Y-coordinate of the first detection point = Y-coordinate of the starting point of the full-scale endpoint - length of the additional segment × Cos(azimuth angle) = 9555 - 9555 × Cos(0°) = 0.

[0083] Step 4: Based on the full-scale endpoint coordinates (0, 9555) and (0, 11185) established in Step 1, the full-scale length is calculated to be 1630m, and the number of shots is calculated to be 202. Combining the first shot coordinates (0, 6390) from Step 3, the last shot coordinates (0, 14430) are calculated, as follows:

[0084] Number of shots = ROUNDUP(1630 / 40+1) + 160 = 202, where ROUNDUP is the round-up function;

[0085] Length of firing point = (number of shots - 1) × distance between firing points = (202 - 1) × 40 = 8040;

[0086] Final shot X coordinate = First shot X coordinate + Shot length × Sin(azimuth) = 0 + 8040 × Sin(0°) = 0;

[0087] The Y-coordinate of the shot point = the Y-coordinate of the first shot point + the length of the shot point × Cos(azimuth) = 6390 + 8040 × Cos(0°) = 14430.

[0088] The construction length refers to the distance from the first receiver point to the last receiver point, which is 20820m. Therefore, the number of receiver points is calculated to be 1042. Combining the coordinates of the first receiver point (0, 0) from step 3, the coordinates of the last receiver point (0, 20820) are calculated as follows:

[0089] Construction length = 2 * 6390 + (202 - 1) * 40 = 20820;

[0090] Number of detector points = 20820 / 20 + 1 = 1042;

[0091] The X coordinate of the last receiver point = the X coordinate of the first receiver point + the construction length × Sin (azimuth angle) = 0 + 20820 × Sin (0°) = 0;

[0092] The Y-coordinate of the final receiver point = the Y-coordinate of the first receiver point + the construction length × Cos(azimuth) = 0 + 20820 × Cos(0°) = 20820.

[0093] Step 5: Calculate the coordinates of the first and last shot checkpoints for the two-dimensional survey line 2 using the methods from Step 2 to Step 4. Obtain the coordinates of the first shot checkpoint (6390, 0), the last shot checkpoint (14430, 0), the first receiver checkpoint (0, 0), and the last receiver checkpoint (20820, 0) for the survey line 2. This allows you to calculate the coordinates of the first and last shot checkpoints for all two-dimensional survey lines.

[0094] This invention combines geophysical exploration industry standards and production practices, and uses VB computer programming language to innovatively realize a rapid, batch calculation method for automatic full coverage of the first and last shot checkpoints of two-dimensional survey lines. This frees technicians from the inconvenience of manually designing the locations of the first and last shot checkpoints for each two-dimensional survey line, and achieves efficient layout of exploration design.

[0095] Working principle and usage process of this invention:

[0096] In use, this invention involves the following steps: Step 1 is establishing preset parameters for all two-dimensional survey lines; Step 2 is calculating the azimuth of a survey line and determining the starting point of the full-coverage endpoint; Step 3 is determining the position of the first shot point and calculating the length of the additional section and the coordinates of the first receiver point; Step 4 is calculating the number of shots, the construction length, and the number of receiver points, and calculating the coordinates of the final shot point and receiver points; Step 5 is obtaining the coordinates of the first and final shot-receiver points of all two-dimensional survey lines through iteration. In Step 1, the preset parameters for the two-dimensional survey lines include the maximum and minimum shot-receiver distances, the coordinates of the full-coverage endpoint, the number of full-coverage cycles, the shot-point distance, and the receiver-point distance for all survey lines. In Step 2, the azimuth of a survey line is calculated based on the full-coverage endpoint coordinates established in Step 1, and the survey line direction is determined according to the range of the azimuth, thereby determining the starting point of the full-coverage endpoint. In step 3, the position of the Mth shot is determined based on the number of full coverages and the starting point of the full coverage endpoint from steps 1 and 2. The distance LPD from the first shot to the starting point of the full coverage endpoint is then calculated, thus determining the position of the first shot. The coordinates of the first shot are calculated using the coordinates of the starting point of the full coverage endpoint and the azimuth, specifically as follows: M = number of full coverages / 2; LPD = 0.5 * minimum shot-receiver distance + (M-1) * shot-point distance. In step 3, the additional segment length refers to the distance from the first receiver point to the starting point of the full coverage endpoint. The coordinates of the first receiver point are calculated based on the additional segment length, specifically as follows: Additional segment length = LPD + maximum shot-receiver distance. In step 4, the full coverage length is calculated based on the full coverage endpoint coordinates established in step 1, and the number of shots is calculated. Combined with the first shot coordinates from step 3, the coordinates of the last shot are calculated, specifically as follows: Number of shots = ROUNDUP(full coverage length / shot-point distance + 1) + number of coverages, where ROUNDUP is the round-up function. In step 4, the construction length refers to the distance from the first receiver point to the last receiver point. The number of receiver points is then calculated, and combined with the coordinates of the first receiver point from step 3, the coordinates of the last receiver point are calculated as follows: Construction length = 2 * maximum shot-receiver distance + (number of shots - 1) * shot-point distance; Number of receiver points = Construction length / receiver-point distance + 1. In step 5, the loop uses the method from steps 2 to 4 to calculate the coordinates of the first and last shot-receiver points for all two-dimensional survey lines.

[0097] The foregoing has shown and described the basic principles, main features, and advantages of the present invention. Those skilled in the art should understand that the above embodiments do not limit the scope of protection of the present invention in any way, and all technical solutions obtained by equivalent substitution or other means fall within the scope of protection of the present invention.

[0098] All parts not covered in this invention are the same as or can be implemented using existing technologies.

Claims

1. A method for automatic full-coverage, rapid, and batch calculation of first and last shot detection points on a two-dimensional survey line, characterized in that, Includes the following steps: Step 1: Establish preset parameters for all two-dimensional survey lines; Step 2: Calculate the azimuth of a survey line and determine the starting point of the full-scale endpoint; Step 3: Determine the location of the first shot, and calculate the length of the additional section and the coordinates of the first receiver point; Step 4: Calculate the number of shots, construction length, and number of geophones; calculate the coordinates of the final shot point and geophones. Step 5: Obtain the coordinates of all the first and last shot checkpoints of the two-dimensional survey line by looping.

2. The method for automatic full-coverage rapid batch calculation of two-dimensional survey line start and end shot detection points according to claim 1, characterized in that, In step 1, the preset parameters of the two-dimensional survey line include the maximum shot-receiver distance, minimum shot-receiver distance, full-cycle endpoint coordinates, full-cycle coverage count, shot-point distance, and receiver-point distance for all survey lines.

3. The method for automatic full-coverage rapid batch calculation of two-dimensional survey line start and end shot detection points according to claim 1, characterized in that, In step 2, the azimuth angle of a survey line is calculated based on the coordinates of the full-scale endpoint established in step 1, and the direction of the survey line is determined according to the range of the azimuth angle, thereby determining the starting point position of the full-scale endpoint.

4. The method for automatic full-coverage rapid batch calculation of two-dimensional survey line start and end shot detection points according to claim 1, characterized in that, In step 3, the position of the Mth shot is determined based on the number of full coverages and the starting point of the full coverage endpoint from steps 1 and 2. Then, the distance LPD from the first shot to the starting point of the full coverage endpoint is calculated, thus determining the position of the first shot. The coordinates of the first shot are calculated using the coordinates of the starting point of the full coverage endpoint and the azimuth angle. Specifically, as follows: M = number of full coverage times / 2, where M is the first shot after the starting point of the full coverage endpoint; LPD = 0.5 * minimum shot-receiver distance + (M-1) * shot distance; The X-coordinate of the first shot point = the X-coordinate of the starting point of the full-scale endpoint – LPD × Sin (azimuth angle); The Y-coordinate of the first shot point = the Y-coordinate of the starting point of the full-scale endpoint – LPD × Cos (azimuth).

5. The method for automatic full-coverage rapid batch calculation of two-dimensional survey line start and end shot detection points according to claim 1, characterized in that, In step 3, the length of the additional segment is the distance from the first receiver point to the starting point of the full-scale endpoint. The coordinates of the first receiver point are then calculated based on the length of the additional segment, as follows: Additional section length = LPD + maximum shot-receiver distance; X-coordinate of the first detection point = X-coordinate of the starting point of the full-scale endpoint – Length of the additional segment × Sin(azimuth angle); Y-coordinate of the first detection point = Y-coordinate of the starting point of the full-scale endpoint – Length of the additional segment × Cos(azimuth).

6. The method for automatic full-coverage rapid batch calculation of two-dimensional survey line start and end shot detection points according to claim 1, characterized in that, In step 4, the full-scale length is calculated based on the full-scale endpoint coordinates established in step 1, and then the number of shots is calculated. Finally, the coordinates of the first shot point are combined with the coordinates of the last shot point from step 3, as detailed below: Number of shots = ROUNDUP(full length / shot distance + 1) + number of coverages, where ROUNDUP is the round-up function; Shot point length = (number of shots - 1) × shot point distance; The X-coordinate of the final shot point = the X-coordinate of the first shot point + the length of the shot point × Sin (azimuth angle); The Y-coordinate of the final shot point = the Y-coordinate of the first shot point + the length of the shot point × Cos (azimuth).

7. The method for automatic full-coverage rapid batch calculation of two-dimensional survey line start and end shot detection points according to claim 1, characterized in that, In step 4, the construction length is the distance from the first receiver point to the last receiver point. The number of receiver points is calculated, and combined with the coordinates of the first receiver point from step 3, the coordinates of the last receiver point are calculated as follows: Construction length = 2 * maximum shot-receiver distance + (number of shots - 1) * shot point distance; Number of detector points = Construction length / Detector point spacing + 1; The X-coordinate of the final receiver point = the X-coordinate of the first receiver point + the construction length × Sin (azimuth angle); The Y-coordinate of the final receiver point = the Y-coordinate of the first receiver point + the construction length × Cos (azimuth).

8. The method for automatic full-coverage rapid batch calculation of two-dimensional survey line start and end shot detection points according to claim 1, characterized in that, In step 5, the loop calculates the coordinates of the first and last shot checkpoints for all two-dimensional survey lines using the methods from steps 2 to 4.

Images