Partitioning method and device for joint excitation of well gun and vibroseis

Abstract

The invention relates to a partition method and device for joint excitation of a well gun and a vibroseis. The method comprises the following steps: acquiring elevation data, earth surface feature data and slope data of a target mountain work area and climbing parameters of a vibroseis vehicle; overlapping the elevation data, the surface feature data and the slope data of the target mountain work area to obtain three-dimensional image data; determining a first communication area on the three-dimensional image data according to the climbing parameters of the vibroseis vehicle, and adjusting the first communication area to obtain a second communication area; a vibroseis construction area and a well shot construction area are respectively determined according to the second communication area and the non-communication area outside the second communication area on the three-dimensional image data, the mountainous area well shot construction area and the vibroseis construction area can be accurately divided, the mountainous area construction proportion of the vibroseis is improved, and the safety of seismic exploration acquisition projects in the mountainous area is ensured. The purposes of reducing cost and increasing efficiency are achieved.

Description

Technical Field

[0001] The present disclosure relates to the technical field of onshore oil seismic exploration data acquisition, and particularly to a zoning method and device for combined excitation of well shots and vibrators. Background Art

[0002] Seismic exploration is the main method for prospecting oil and gas, and its main work includes three stages: seismic data acquisition, processing, and interpretation. Seismic data acquisition is the first process of an oil and gas seismic exploration project, mainly by artificially exciting and receiving seismic waves to obtain original seismic data containing underground geological information. Onshore seismic acquisition work mainly includes three parts: seismic wave excitation, seismic wave reception, and observation methods. Common onshore excitation sources include vibrator excitation and explosive excitation. All along, explosive excitation has been used in the excitation link of mountain seismic acquisition. When the relatively shallow strata are the target, a relatively weakened acquisition scheme and a relatively low shot density are adopted, and the data is relatively easy to obtain, while the acquisition input cost is relatively controllable.

[0003] In recent years, with the continuous deepening of exploration targets and the gradual increase in external human environment interference, it has become increasingly difficult to obtain deep and ultra-deep data. Adopting a relatively intensive seismic acquisition scheme and a multiple increase in shot density are currently the most effective methods for obtaining deep and ultra-deep data. A multiple increase in shot density means that more excitation wells need to be drilled, resulting in higher costs. To ensure data quality and reduce acquisition costs, taking into account the advantage of the low frequency of vibrators in imaging deep structures, a construction method with a relatively low-cost vibrator entering mountain construction as much as possible can be adopted to reduce the proportion of well shots.

[0004] With the development and progress of unmanned aerial vehicle (UAV) technology, UAV aerial survey has been widely used in the field of field seismic exploration. In related technologies, UAV aerial survey images and DEM (digital terrain model) data are used to guide indoor reconnaissance, geomorphic zoning, indoor physical point layout, production of driving route maps, and production management. However, this method is not fine enough, and the production of driving route maps is too simple. In related technologies, a method of using slopes, contour lines, and bulldozer-repaired flight tracks to draw acquisition route maps is also proposed, which first introduced in detail the process of designing field construction paths using aerial survey data. In related technologies, a construction process of "repairing roads first and then setting out" for path design operations and a method of designing acquisition paths based on high-definition geographical images and high-precision elevation data obtained from UAV aerial surveys are also proposed, which solves the problem of difficult uniform layout of excitation points in desert areas. However, high-definition aerial survey image data in mountainous areas has not been effectively utilized. Summary of the Invention

[0005] To solve the above technical problems or at least partially solve the above technical problems, embodiments of the present disclosure provide a zoning method and device for combined excitation of well shots and vibrators.

[0006] In a first aspect, embodiments of the present disclosure provide a zoning method for combined excitation of well shots and vibrators, including:

[0007] Obtain elevation data, surface feature data, and slope data of the target mountainous work area, as well as the climbing parameters of the vibrator vehicle;

[0008] Overlay the elevation data, surface feature data, and slope data of the target mountainous work area to obtain three-dimensional image data;

[0009] Determine a first connected area on the three-dimensional image data according to the climbing parameters of the vibrator vehicle, and adjust the first connected area to obtain a second connected area;

[0010] Determine the vibrator construction area and the well shot construction area respectively according to the second connected area on the three-dimensional image data and the unconnected areas outside the second connected area.

[0011] In a possible implementation manner, the climbing parameters of the vibrator vehicle are obtained through the following steps:

[0012] Obtain the tire radius of the vibrator vehicle;

[0013] Determine the maximum driving force of the vibrator vehicle according to the tire radius of the vibrator vehicle and the engine parameters;

[0014] Determine the climbing parameters of the vibrator vehicle according to the maximum driving force of the vibrator vehicle and the climbing resistance coefficient.

[0015] In a possible implementation manner, the tire radius of the vibrator vehicle is obtained through the following expression:

[0016]

[0017] where R is the tire radius of the vibrator vehicle, in meters, D is the tire width of the vibrator vehicle, in inches, a is the tire aspect ratio of the vibrator vehicle, and d is the wheel hub diameter of the vibrator vehicle, in inches.

[0018] In a possible implementation manner, the maximum driving force of the vibrator vehicle is determined according to the tire radius of the vibrator vehicle and the engine parameters through the following expression:

[0019]

[0020] Among them, C is the maximum driving force of the vibratory source vehicle, with the unit of Newton, N is the maximum torque of the engine, with the unit of Newton-meter, B is the gear ratio, and W is the final drive of the engine.

[0021] In a possible implementation manner, through the following expression, the climbing parameters of the vibratory source vehicle are determined according to the maximum driving force and the climbing resistance coefficient of the vibratory source vehicle:

[0022]

[0023] Among them, α max is the maximum climbing angle of the vibratory source vehicle, k is the climbing resistance coefficient, and Z is the sum of the weights of the vibratory source vehicle and the driver, with the unit of ton.

[0024] In a possible implementation manner, determining the first connected region on the three-dimensional image data according to the climbing parameters of the vibratory source vehicle includes:

[0025] Regarding the region where the surface dip angle is less than or equal to the maximum climbing angle of the vibratory source and the minimum width is greater than or equal to the preset width threshold as the first connected region on the three-dimensional image data.

[0026] In a possible implementation manner, through the following expression, the climbing resistance coefficient is obtained:

[0027] k = (A * cosθ + m * g * cosθ * h) / (m * g * sinθ)

[0028] Among them, A is the contact area between the tires of the vibratory source vehicle and the ground surface, θ is the surface dip angle, m is the mass of the vibratory source vehicle, g is the acceleration due to gravity, and h is the center of gravity height of the vibratory source vehicle.

[0029] In a possible implementation manner, adjusting the first connected region to obtain a second connected region includes:

[0030] For the first target unconnected region, in the case where the number of vibratory source excitation points arranged after building a road in this first target unconnected region exceeds the first preset number threshold, regarding this first target unconnected region as the region to be built, where the surface dip angle of the first target unconnected region is less than or equal to the maximum climbing angle of the vibratory source and the minimum width is less than the preset width threshold;

[0031] Regarding the first connected region surrounded by the unconnected region as the third connected region;

[0032] For the second target disconnected area, when the two first connected areas are on the same horizontal plane, the second target disconnected area is taken as the area to be repaired. When the two first connected areas are not on the same horizontal plane, the second target disconnected area with a slope less than the preset angle threshold is taken as the area to be repaired, where the surface dip angle of the second target disconnected area is greater than the maximum climbing angle of the vibrator and the minimum width is greater than or equal to the preset width threshold, and it is between two adjacent first connected areas where the vibrator excitation points both exceed the second preset quantity threshold;

[0033] Include the area to be repaired into the first connected area, and exclude the third connected area from the first connected area to obtain the second connected area.

[0034] In a possible implementation manner, for the second target disconnected area with a slope less than the preset angle threshold when the two first connected areas are not on the same horizontal plane, the road repair length of the second target disconnected area is calculated through the following steps:

[0035] Calculate the height difference between the two road repair connection points at the top and bottom of the slope of the second target disconnected area;

[0036] Through the following expression, calculate the longest horizontal distance in the horizontal direction between two points:

[0037]

[0038] where α max is the maximum climbing angle of the vibrator vehicle, i max is the actual maximum climbing slope of the vibrator considering the friction coefficients of different terrains, H is the height difference between the top and bottom of the slope, and L is the longest horizontal distance in the direction perpendicular to the slope direction,

[0039] Through the following expression, calculate the road repair length along the slope according to the calculated height difference and horizontal distance:

[0040]

[0041] where T is the shortest road repair length along the slope.

[0042] In a possible implementation manner, the determining the vibrator construction area and the downhole shot construction area according to the second connected area and the disconnected areas outside the second connected area on the three-dimensional image data includes:

[0043] Respectively take the second connected area and the disconnected areas outside the second connected area on the three-dimensional image data as the vibrator construction area and the downhole shot construction area, and form the construction boundary between the downhole shot and the vibrator according to the vibrator construction area and the downhole shot construction area;

[0044] Layout specific data points according to the designed observation system with well shots and vibroseis.

[0045] At the boundary between well shots and vibroseis, using the offset principle of the shot point, under the condition of allowing regular offset, offset the points within the well shot construction area and design them as vibroseis shot points, obtaining the optimized construction boundary between well shots and vibroseis trucks for re-dividing the vibroseis construction area and well shot construction area on the 3D image data.

[0046] In a second aspect, an embodiment of the present disclosure provides a zoning device for combined excitation of well shots and vibroseis, including:

[0047] An acquisition module for acquiring elevation data, surface feature data, slope data of the target mountain work area, and the climbing parameters of the vibroseis truck;

[0048] An overlay module for overlaying the elevation data, surface feature data, and slope data of the target mountain work area to obtain 3D image data;

[0049] An adjustment module for determining the first connected area on the 3D image data according to the climbing parameters of the vibroseis truck and adjusting the first connected area to obtain a second connected area;

[0050] A determination module for determining the vibroseis construction area and the well shot construction area respectively according to the second connected area on the 3D image data and the unconnected areas outside the second connected area.

[0051] 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 mutual communication through the communication bus;

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

[0053] The processor, when executing the program stored on the memory, implements the above-mentioned zoning method for combined excitation of well shots and vibroseis.

[0054] In a fourth aspect, an embodiment of the present disclosure provides a computer-readable storage medium, on which a computer program is stored. It is characterized in that the computer program, when executed by a processor, implements the above-mentioned zoning method for combined excitation of well shots and vibroseis.

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

[0056] The zoning method for combined excitation of well shots and vibrators according to the embodiments of the present disclosure obtains elevation data, surface feature data, slope data of a target mountain work area, and the climbing parameters of vibrator vehicles; superimposes the elevation data, surface feature data, and slope data of the target mountain work area to obtain three-dimensional image data; determines a first connected area on the three-dimensional image data according to the climbing parameters of the vibrator vehicles, and adjusts the first connected area to obtain a second connected area; determines the vibrator construction area and the well shot construction area respectively according to the second connected area on the three-dimensional image data and the unconnected areas outside the second connected area, which can accurately divide the well shot and vibrator construction areas in the mountain area, improve the proportion of vibrator construction in the mountain area and ensure the safety of seismic exploration acquisition projects in the mountain area, achieving the purpose of cost reduction and efficiency increase. Description of the Drawings

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

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

[0059] Figure 1 Schematically shows a schematic flow chart of the zoning method for combined excitation of well shots and vibrators according to the embodiments of the present disclosure;

[0060] Figure 2 Schematically shows a schematic diagram of three-dimensional surface feature data according to the embodiments of the present disclosure;

[0061] Figure 3 Schematically shows a schematic diagram of three-dimensional elevation data according to the embodiments of the present disclosure;

[0062] Figure 4 Schematically shows a schematic diagram of three-dimensional slope data according to the embodiments of the present disclosure;

[0063] Figure 5 Schematically shows a schematic diagram of the superimposed three-dimensional image data according to the embodiments of the present disclosure;

[0064] Figure 6 Schematically shows a schematic diagram of the first connected area in the superimposed three-dimensional image data according to the embodiments of the present disclosure;

[0065] Figure 7 Schematically shows a schematic diagram of the area where the mountain gully narrows and hinders connectivity according to the embodiments of the present disclosure;

[0066] Figure 8 Schematically shows a schematic diagram of local obstacles in a large - area connected area according to an embodiment of the present disclosure;

[0067] Figure 9 Schematically shows a schematic diagram of a road - building trajectory on a high - steep slope designed by the maximum climbing gradient according to an embodiment of the present disclosure;

[0068] Figure 10 Schematically shows a partition schematic diagram of well shots and vibroseis according to an embodiment of the present disclosure;

[0069] Figure 11 Schematically shows a structural block diagram of a partition device for combined excitation of well shots and vibroseis according to an embodiment of the present disclosure;

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

[0071] To make the objectives, technical solutions, and advantages of the embodiments of the present disclosure clearer, the technical solutions in the embodiments of the present disclosure will be clearly and completely described below 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 in the present disclosure without creative efforts shall fall within the protection scope of the present disclosure.

[0072] See Figure 1 , an embodiment of the present disclosure provides a partition method for combined excitation of well shots and vibroseis, including the following steps:

[0073] S1, obtain elevation data, surface feature data, slope data of the target mountain work area, and the climbing parameters of the vibroseis vehicle.

[0074] As Figure 2 shown, the surface feature data is high - precision image data of surface feature information obtained by low - altitude aerial photography of an unmanned aerial vehicle (UAV), processed by image processing software, and registered and corrected indoors. From the figure, surface feature information can be seen, especially the surface of different land types in the mountain, ground equipment information, and the driving paths of the old survey lines in the work area. The resolution accuracy of the surface feature data is greater than or equal to 0.1m * 0.1m. As Figure 3 shown, the elevation data is high - precision elevation image data obtained by low - altitude aerial photography of an unmanned aerial vehicle (UAV) and processed by image processing software. From the figure, the undulation of the ground surface can be seen. The resolution accuracy of the elevation data is greater than or equal to 1m * 1m. As Figure 4 shown, the slope data is data reflecting the slope size of the surface undulation state calculated from the high - precision elevation image data. From the figure, the steepness of the mountain can be seen.

[0075] S2. Overlay the elevation data, surface feature data, and slope data of the target mountain work area to obtain three-dimensional image data.

[0076] In this embodiment, the overlay can be performed in the following manner: The first layer is the elevation data with the transparency set to 0%; the second layer is the surface feature data with the transparency set to 0%; the third layer is the slope data, and the slope is divided into the following slope segments: less than the maximum climbing angle of the vibrator; between the maximum climbing angle of the vibrator and 30 degrees; greater than 30 degrees. The area less than the maximum climbing angle of the vibrator is set to colorless; for the other slope segments, from small to large, the color is set from light to deep, and the color transparency of all slope segments is set to 30%. Through the above overlay method, the information of each layer can be taken into account to form a 3D high-definition image with a three-dimensional intuitive feeling.

[0077] See Figure 5 , the vehicle driving trajectory can be used as the fourth layer and overlaid with the first layer, the second layer, and the third layer together to obtain three-dimensional image data. Among them, the bottom layer of the three-dimensional image data is the elevation data layer, the surface feature data is the second layer covering the elevation data layer, the third layer is the slope layer, and the fourth layer is the vehicle driving trajectory. Among them, the vehicle driving trajectory can be formed by picking up the coordinates of the gravel road and cement road in the mountain gully and on the gentle slope of the surface feature data.

[0078] S3. Determine the first connected area on the three-dimensional image data according to the climbing parameters of the vibrator vehicle, and adjust the first connected area to obtain the second connected area.

[0079] In this embodiment, in step S3, the determining the first connected area on the three-dimensional image data according to the climbing parameters of the vibrator vehicle includes:

[0080] Regard the area where the surface dip angle is less than or equal to the maximum climbing angle of the vibrator and the minimum width is greater than or equal to the preset width threshold as the first connected area on the three-dimensional image data. Usually, the width of the vibrator vehicle is 3.5m. To ensure the safe passage of the vibrator, the preset width threshold can be 4m.

[0081] As Figure 6 shown, the area where the continuous sample point slope in the mountain gully is equal to or the width greater than or equal to the maximum climbing angle of the vibrator is greater than or equal to the preset width threshold forms the first connected area, and vibrator construction is adopted.

[0082] In this embodiment, regarding the area where the surface dip angle is less than or equal to the maximum climbing angle of the vibrator and the minimum width is greater than or equal to the preset width threshold as the first connected area on the three-dimensional image data includes:

[0083] A 0.5m×0.5m grid is established, where the coordinate system of the grid data points is the Beijing 54 coordinate system, and the coordinate data includes the stake number, east coordinate, north coordinate, and surface elevation;

[0084] Create data points at each grid node;

[0085] Attach the high-definition image slope value to each data point;

[0086] Extract the slope values in a fixed format (line number + point number + E coordinate + N coordinate + slope data) to form a data file (*.txt);

[0087] Extract the data points with slopes less than the maximum climbing angle of the vibroseis in the longitudinal (transverse) direction respectively. Set the areas with slopes greater than or equal to 8 consecutive points (i.e., 4m) as the areas where the vibroseis can pass through, and set them as the first connected area;

[0088] Set the positions of discontinuous points or continuous points with less than 8 excitation points (i.e., 4m) as the points to be judged for road construction. It is necessary to further judge whether road construction is needed, so as to further expand the connected area; Optimize the points to be judged for road construction and accurately divide the boundary between mountain well shots and vibroseis.

[0089] In this embodiment, the maximum climbing angle refers to the maximum slope that the vibroseis can climb when driving with the maximum driving force on a smooth road surface. It represents the ultimate climbing ability of the vibroseis vehicle. Connect the areas with slopes less than the maximum climbing angle of the vibroseis as the vibroseis construction area, and this area is called the connected area.

[0090] S4. Determine the vibroseis construction area and the well shot construction area respectively according to the second connected area and the unconnected areas outside the second connected area on the three-dimensional image data.

[0091] In this embodiment, the vibroseis construction area and the well shot construction area are used in well-seismic hybrid acquisition. Among them, well-seismic hybrid acquisition: refers to the construction method of a seismic acquisition project using the combined excitation of well shots and vibroseis, and can also be called well-seismic combined excitation.

[0092] In this embodiment, the climbing ability of the vibroseis vehicle in mountainous areas is closely related to the maximum driving force of the vibroseis vehicle, the weight of the vibroseis vehicle, the width and radius of the tires. And the maximum driving force of the vibroseis vehicle is closely related to factors such as its power system, engine power, driving device, and transmission system. The greater the power, the more powerful the driving device, and the more efficient the transmission system, the greater the maximum driving force of the vibroseis vehicle usually is, and its climbing ability in mountainous areas will also be improved accordingly. In this embodiment, in step S1, the climbing parameters of the vibroseis vehicle are obtained through the following steps:

[0093] Obtain the tire radius of the vibroseis vehicle;

[0094] Determine the maximum driving force of the vibroseis truck according to the tire radius and engine parameters of the vibroseis truck;

[0095] Determine the climbing parameters of the vibroseis truck according to the maximum driving force and climbing resistance coefficient of the vibroseis truck.

[0096] In this embodiment, the tire radius of the vibroseis truck is obtained through the following expression:

[0097]

[0098] where R is the tire radius of the vibroseis truck in meters, D is the tire width of the vibroseis truck in inches, a is the tire aspect ratio of the vibroseis truck, and d is the wheel hub diameter of the vibroseis truck in inches. Usually, the tire width of the vibroseis truck is 43 inches, the tire aspect ratio is 0.4767, and the wheel hub diameter is 25 inches. 1 inch = 25.4 mm.

[0099] In this embodiment, the maximum driving force of the vibroseis truck is determined through the following expression according to the tire radius and engine parameters of the vibroseis truck:

[0100]

[0101] where C is the maximum driving force of the vibroseis truck in Newtons, N is the maximum engine torque in Newton-meters, B is the gear ratio, and W is the engine final drive. Usually, the maximum engine torque is 2000 Newton-meters, the gear ratio is 4.5, and the engine final drive is 8.42.

[0102] In this embodiment, a climbing resistance coefficient k is set. k is related to factors such as surface lithology and surface softness. Conduct climbing tests on different terrains in the mountains, and obtain the climbing resistance coefficient of the vibroseis on different surfaces. Combine the maximum driving force of the vibroseis and the weight of the vibroseis, and determine the climbing parameters of the vibroseis truck according to the maximum driving force and climbing resistance coefficient of the vibroseis truck through the following expression:

[0103]

[0104] where α max is the maximum climbing angle of the vibroseis truck, k is the climbing resistance coefficient, and Z is the sum of the weights of the vibroseis truck and the driver in tons. Usually, Z is 28.1 tons, and the gravity of 1 kg is 9.8 Newtons.

[0105] In this embodiment, the dynamic test method is used to calculate the climbing resistance coefficient k of the vibroseis on different surfaces:

[0106] Step 1: Prepare the test equipment and surface samples. Prepare a vibroseis truck as the test vehicle and install appropriate test equipment, such as force sensors, accelerometers, etc. At the same time, prepare surface samples of different surface types, such as loess, gravel, sandstone, etc., and ensure that their area is large enough to ensure the accuracy and reliability of the test;

[0107] Step 2: Conduct on-site tests. Drive the test vehicle onto different surface samples and conduct vehicle motion tests. By controlling parameters such as the vehicle's speed, acceleration, and tire pressure, record the vehicle's motion state on different surfaces, such as acceleration, speed, and tire force. At the same time, record the climbing angle of the vibroseis truck on different surfaces;

[0108] Step 3: Process the test data. Process the data obtained from the tests, such as calculating parameters such as the friction force, speed, and acceleration of the vehicle on different surfaces. Use statistical analysis methods to obtain an estimated value of the climbing resistance coefficient of the vibroseis;

[0109] Step 4: Calculate the climbing resistance coefficient of the vibroseis. Use the calculation formula, substitute the parameters obtained from the tests into the calculation, and obtain the estimated values of the climbing resistance coefficients for different surface types. Through the following expression, obtain the climbing resistance coefficient:

[0110] k = (A * cosθ + m * g * cosθ * h) / (m * g * sinθ)

[0111] where A is the contact area between the tires of the vibroseis truck and the surface, θ is the surface inclination angle, m is the mass of the vibroseis truck, g is the acceleration due to gravity, and h is the height of the center of gravity of the vibroseis truck.

[0112] Based on multiple comparisons of the processing results and calculation results of multiple rounds of test data, obtain the climbing resistance coefficients of the vibroseis on mountain loess, grassland, gravel, and rock surfaces, as shown in Table 1 below.

[0113] Table 1

[0114]

[0115]

[0116] In this embodiment, in step S3, the adjusting the first communication area to obtain a second communication area includes:

[0117] For the first target unconnected area, in the case where the number of vibroseis excitation points arranged after building a road in the first target unconnected area exceeds the first preset number threshold, regard the first target unconnected area as the area to be built, where the surface inclination angle of the first target unconnected area is less than or equal to the maximum climbing angle of the vibroseis and the minimum width is less than the preset width threshold;

[0118] For the first connected region surrounded by a disconnected region, it is used as the third connected region;

[0119] For the second target disconnected region, when the two first connected regions are on the same horizontal plane, the second target disconnected region is used as the area to be repaired; when the two first connected regions are not on the same horizontal plane, the second target disconnected region with a slope less than the preset angle threshold is used as the area to be repaired, where the surface dip angle of the second target disconnected region is greater than the maximum climbing angle of the vibroseis and the minimum width is greater than or equal to the preset width threshold, and it is between two adjacent first connected regions where the vibroseis excitation points all exceed the second preset quantity threshold;

[0120] The area to be repaired is incorporated into the first connected region, and the third connected region is excluded from the first connected region to obtain the second connected region.

[0121] In practical applications, the first target disconnected region can be an area where the mountain gully becomes narrow and hinders connectivity. When encountering a situation where the connected region such as a gully becomes narrow by less than 4m and the vibroseis cannot pass through, it is necessary to determine whether to repair the road for connectivity, so as to determine the scope of the connected region. The specific method is as follows:

[0122] Establish a 0.5m×0.5m grid;

[0123] Create data points at each grid node;

[0124] Attach the high-definition image slope value to each data point;

[0125] Extract the slope values in a fixed format (line number + point number + E coordinate + N coordinate + slope data) to form a data file (*.txt);

[0126] Extract the data points with slopes less than the maximum climbing angle of the vibroseis in the longitudinal (transverse) direction respectively, and set the area with more than or equal to 8 consecutive points (i.e., 4m) as the area where the vibroseis can pass through, which is set as the first connected region;

[0127] Set the discontinuous points or the positions with less than 8 consecutive excitation points (i.e., 4m) as the points to judge whether to repair the road, and it is necessary to further determine whether to repair the road to further expand the connected region; preferably select the points to repair the road and accurately divide the boundary between the mountain well shots and the vibroseis. The judgment of the points to repair the road mainly depends on the cost of road repair. First, the area of the region with a slope less than the maximum climbing angle of the vibroseis blocked by the point to repair the road should be large enough to ensure that at least 10 excitation points can be arranged after the point to repair the road is repaired. If the condition is met, it is set as a point where the road can be repaired, and the blocked region is set as the connected region, regarded as the construction area of the vibroseis; if the condition is not met, it is set as a point where the road cannot be repaired, the point to repair the road is automatically closed, and the region blocked by the point to repair the road is also automatically set as the construction area of the well shots.

[0128] See Figure 7 , the base map is the layer of Figure 5 . There are two relatively large separated connected areas, separated by narrow gullies. Determine the location of the road construction point and make it connected.

[0129] In practical applications, the first connected area surrounded by unconnected areas can be a mountainous high-steep terrain, that is, there is a local connected area within a large area of unconnected areas and it cannot be connected. In a mountain area with a large slope, there are areas with local flat terrain, small slopes, and areas where vibrator sources can be constructed. Such terrain is obviously not suitable for road construction and cannot be connected to adjacent connected areas. Set such terrain as an obstacle area, draw the obstacles, and no road construction points are set within the obstacle area. Close the connected area, and well-shot excitation points are allowed to be arranged within the obstacles.

[0130] In practical applications, obstacles with a relatively large slope that hinder connection and are surrounded by connected areas are regarded as unconnected areas, that is, within the range completely surrounded by the connected area, there are local terrains with relatively large slopes, such as river channel steep slopes, dams, water channels, single high-steep mountains, tall buildings, etc. (a. Obstacles with a longitudinal or transverse width less than 250m; b. Obstacles with both longitudinal and transverse widths less than 500m). Excitation points cannot be arranged. Set such terrain as an obstacle area, draw the obstacles. The obstacle area is included in the connected area, no road construction points are set within the obstacles, and vibrator source excitation points are not allowed to be arranged.

[0131] See Figure 8 , local obstacles within a large-area connected area, the base map is the layer of Figure 5 . Set the area within the connected area where vibrator sources cannot be arranged as an obstacle area. No road construction points are set within the obstacles, and vibrator source excitation points are not allowed to be arranged.

[0132] In practical applications, the second target non-connected area can be an area used to partition adjacent large connected areas. Two large connected areas (with at least 20 shots at the excitation points of each connected area) are partitioned by an obstacle with a large slope (greater than the maximum climbing angle of the vibroseis). This results in a situation where they cannot be connected. ① If the two connected areas are on the same horizontal plane, they are connected according to the driving track. If they cannot be connected, a road can be built to connect them and establish a connected area. ② If the two connected areas are not on the same horizontal plane, one connected area is at the bottom of the slope (lower position), and the other connected area is at the top of the slope (higher position). The top of the slope is relatively flat, resulting in a situation where they cannot be connected. It is necessary to determine whether they can be connected. First, determine the slope range of the steep slope that partitions the two connected areas: a. If the partition slope is greater than 30 degrees, it is regarded as a non-road-building point, and the connected area in this area is automatically partitioned and judged as not connectable. The connected area on the top of the mountain is automatically closed and set as a mountain well-shot construction area, where well-shot excitation points can be arranged. b. If the partition slope is less than 30 degrees, a road-building route can be designed for this steep slope. Among them, the slope-direction connection design for the connected area partitioned by the mountain steep slope.

[0133] In this embodiment, for the terrain where two connected areas not on the same horizontal plane are partitioned by a mountain steep slope with a partition slope less than 30 degrees and need to be connected by building a road, a targeted slope-direction connection design is required. A lateral route can be designed in the direction at a certain angle (β) to the slope for building a road. In this way, along the road-building direction, the climbing angle of the vibroseis vehicle can be reduced. By selecting the shortest designed route, it is possible to ensure that the climbing angle of the vibroseis vehicle along the road-building direction is lower than the maximum climbing angle of the vibroseis vehicle, ensuring that the vibroseis can pass. Use the maximum climbing slope (i max ) of the vibroseis vehicle to calculate the road-building length. First, calculate the height difference (H) between the two road-building connection points at the top and bottom of the slope. Secondly, through the following expression, calculate the longest horizontal distance (L) in the horizontal direction between the two points; through the following expression, calculate the road-building length T along the slope based on the calculated height difference and horizontal distance:

[0134]

[0135]

[0136] Among them, α max is the maximum climbing angle of the vibroseis vehicle, i max is the actual maximum climbing slope of the vibroseis considering the friction coefficient of different terrains, H is the height difference between the top and bottom of the slope, L is the longest horizontal distance in the direction perpendicular to the slope direction, and T is the shortest road-building length along the slope.

[0137] Therefore, as long as it can be ensured that it is possible to design the longest T, ensuring that after the road construction, the climbing angle of the vibroseis along the road construction direction is less than the maximum climbing angle of the vibroseis, so that the vibroseis can pass through, thereby connecting the two connected areas and further accurately defining the connected area range. On the contrary, if the designed road construction route cannot meet the requirements even when it is the longest this condition, it can be judged that the vibroseis cannot pass through after the road construction at this partition point. Cancel this road construction point, close the connected area on the mountain top, and set it as the construction area for downhole guns in the mountain, where downhole gun firing points can be arranged.

[0138] See Figure 9 , design the road construction trajectory of the high-steep slope through the maximum climbing gradient. The base map is Figure 5 the layer of. Design the gentlest road construction trajectory on the high-steep mountain slope between two large connected areas to ensure the connection of the two connected areas and ensure the normal construction of the vibroseis.

[0139] In this embodiment, in step S4, the determination of the vibroseis construction area and the downhole gun construction area according to the second connected area and the unconnected area outside the second connected area on the three-dimensional image data includes:

[0140] Respectively take the second connected area and the unconnected area outside the second connected area on the three-dimensional image data as the vibroseis construction area and the downhole gun construction area, and form the construction boundary between the downhole guns and the vibroseis according to the vibroseis construction area and the downhole gun construction area;

[0141] Layout specific data points for the downhole guns and the vibroseis according to the designed observation system;

[0142] At the boundary between the downhole guns and the vibroseis, using the offset principle of the firing points, under the condition of allowing regular offset, offset the points in the downhole gun construction area and design them as vibroseis firing points, increase the vibroseis construction area and the number of firing points, and narrow the construction boundary of the downhole guns in the mountain to obtain the optimized construction boundary between the downhole guns and the vibroseis vehicles, realizing the fine division of the construction areas of the downhole guns and the vibroseis.

[0143] See Figure 10 , the partition of the downhole guns and the vibroseis, where the black area is the construction area for downhole guns in the mountain, and downhole guns are used for excitation; the remaining area is the vibroseis construction area, and vibroseis are used for excitation.

[0144] The zoning method for combined shooting of well shots and vibrators in the present disclosure proposes a mountain well-seismic design method based on 3D stereo superimposed high-definition images, a calculation method for the "climbing ability of vibrators in mountainous areas", a design method for vibrator connectivity areas, and the concept of designing road construction points indoors. It also proposes a brand-new accurate well-seismic zoning design method and process for mountainous areas, replacing the traditional process of "humans first" for field reconnaissance and then design operations. The well-seismic zoning is accurate, and the design accuracy of excitation points is high. In addition, for the problem of accurate well-seismic design in the mixed well-seismic operation in mountainous areas on land, a brand-new well-seismic zoning design method and design process are proposed. By using surface feature data (accuracy greater than or equal to 0.1 m), high-precision elevation data (accuracy greater than or equal to 1 m), and slope data to form 3D high-definition images, combined with the vibrator connectivity area in the work area, the bulldozer road construction track, and the reconnaissance track, the precise division of the vibrator construction area in mountainous areas is carried out to determine the well-seismic demarcation.

[0145] Aiming at the problems of high construction cost, low efficiency, and poor data quality in the single well-shot excitation construction for onshore mountain seismic acquisition operations, the zoning method for combined shooting of well shots and vibrators in the present disclosure uses surface feature data, elevation data, slope data, and vehicle passage paths, combined with the mountain climbing ability of vibrators and the slope connectivity where large connected areas in the piedmont zone are blocked by mountain steep slopes, to conduct zoning for combined shooting of well shots and vibrators in mountainous areas. It can effectively support the economy, safety, and environmental protection of complex mountain seismic acquisition projects, significantly reduce costs and increase efficiency, and improve the economic feasibility of seismic acquisition construction projects in mountainous areas. The method proposed in the present invention has been verified for a certain mountain seismic acquisition project in western China.

[0146] See Figure 11 , the embodiments of the present disclosure provide a zoning device for combined shooting of well shots and vibrators, including:

[0147] An acquisition module 11, configured to acquire elevation data, surface feature data, slope data of a target mountain work area, and the climbing parameters of a vibrator vehicle;

[0148] A superimposing module 12, configured to superimpose the elevation data, surface feature data, and slope data of the target mountain work area to obtain three-dimensional image data;

[0149] An adjustment module 13, configured to determine a first connected area on the three-dimensional image data according to the climbing parameters of the vibrator vehicle, and adjust the first connected area to obtain a second connected area;

[0150] A determination module 14, configured to determine a vibrator construction area and a well-shot construction area respectively according to the second connected area on the three-dimensional image data and the unconnected area outside the second connected area.

[0151] The partition device for combined excitation of well shots and vibrators in the present disclosure precisely divides the construction areas of well shots and vibrator trucks in the well-vibrator mixed acquisition project, improves the proportion of vibrator operations in mountainous areas, enhances the acquisition efficiency, and reduces the acquisition cost.

[0152] The partition device for combined excitation of well shots and vibrators in the present disclosure creates 3D stereoscopic superimposed high-definition images, experimentally quantifies the climbing ability of vibrators, establishes and delineates the connected areas, designs the slope connection where the connected areas are blocked by mountain slopes, and precisely divides the construction areas of well shots and vibrators. It focuses on improving quality and efficiency in the field and conducts research on how to implement the method for dividing the excitation positions of well-vibrator in mountainous areas and the method for increasing the proportion of vibrators in mountainous areas, etc.

[0153] For the realization processes of the functions and roles of each unit in the above device, please refer to the realization processes of the corresponding steps in the above method for details, and they will not be elaborated here.

[0154] For the device embodiments, since they basically correspond to the method embodiments, the relevant parts can refer to the partial descriptions of the method embodiments. The device embodiments described above are merely 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.

[0155] In the above embodiments, any combination of the acquisition module 11, the superimposing module 12, the adjustment module 13, and the determination module 14 can be combined and implemented in one module, or any one of the modules can be split into multiple modules. Or, 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 acquisition module 11, the superimposing module 12, the adjustment module 13, and the determination module 14 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 a chip, a system on a substrate, a system on a package, an application specific integrated circuit (ASIC), or can be implemented by any other reasonable way of integrating or packaging circuits, etc., in hardware or firmware, or implemented in any one of the three implementation methods of software, hardware, and firmware, or in an appropriate combination of any several of them. Or, at least one of the acquisition module 11, the superimposing module 12, the adjustment module 13, and the determination module 14 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.

[0156] Refer to Figure 12As shown in the figure, 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;

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

[0158] When the processor 1110 is used to execute the program stored on the memory 1130, the following partition method for the combined excitation of well shots and vibrators is implemented:

[0159] Obtain the elevation data, surface feature data, slope data of the target mountain work area, and the climbing parameters of the vibrator vehicle;

[0160] Overlay the elevation data, surface feature data, and slope data of the target mountain work area to obtain three-dimensional image data;

[0161] Determine the first connected area on the three-dimensional image data according to the climbing parameters of the vibrator vehicle, and adjust the first connected area to obtain a second connected area;

[0162] Determine the vibrator construction area and the well shot construction area according to the second connected area on the three-dimensional image data and the unconnected areas outside the second connected area respectively.

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

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

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

[0166] 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.

[0167] 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 above-mentioned zoning method for combined excitation of well shots and vibrators is implemented.

[0168] 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 above-mentioned one or more programs are executed, the zoning method for combined excitation of well shots and vibrators according to the embodiments of the present disclosure is implemented.

[0169] According to an embodiment 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 can be used by or combined with an instruction execution system, device, or apparatus.

[0170] It should be noted that, in this document, 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 terms "comprising", "including" or any other variant thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device comprising a series of elements not only includes those elements, but also includes other elements not expressly listed, or further includes elements inherent to such process, method, article or device. Without further limitation, an element defined by the statement "comprising an..." does not exclude the existence of additional identical elements in the process, method, article or device comprising the said element.

[0171] 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 should be accorded the widest scope consistent with the principles and novel features claimed herein.

Claims

1. A zoning method for combined excitation of well shots and vibroseis, characterized in that, the method includes: Obtain the elevation data, surface feature data and slope data of the target mountain work area, as well as the climbing parameters of the vibroseis vehicle; Overlay the elevation data, surface feature data and slope data of the target mountain work area to obtain three-dimensional image data; Determine the first connected area on the three-dimensional image data according to the climbing parameters of the vibroseis vehicle, and adjust the first connected area to obtain the second connected area; Determine the vibroseis construction area and the well shot construction area according to the second connected area on the three-dimensional image data and the unconnected areas outside the second connected area respectively.

2. The method according to claim 1, characterized in that, the climbing parameters of the vibroseis vehicle are obtained through the following steps: Obtain the tire radius of the vibroseis vehicle; Determine the maximum driving force of the vibroseis vehicle according to the tire radius and engine parameters of the vibroseis vehicle; Determine the climbing parameters of the vibroseis vehicle according to the maximum driving force and climbing resistance coefficient of the vibroseis vehicle.

3. The method according to claim 2, characterized in that, obtain the tire radius of the vibroseis vehicle through the following expression: wherein, R is the tire radius of the vibroseis vehicle, in meters, D is the tire width of the vibroseis vehicle, in inches, a is the tire aspect ratio of the vibroseis vehicle, and d is the wheel hub diameter of the vibroseis vehicle, in inches.

4. The method according to claim 2, characterized in that, determine the maximum driving force of the vibroseis vehicle according to the tire radius and engine parameters of the vibroseis vehicle through the following expression: wherein, C is the maximum driving force of the vibroseis vehicle, in Newtons, N is the maximum engine torque, in Newton-meters, B is the gear ratio, and W is the engine final drive.

5. The method according to claim 2, characterized in that, determine the climbing parameters of the vibroseis vehicle according to the maximum driving force and climbing resistance coefficient of the vibroseis vehicle through the following expression: Among them, α max is the maximum climbing angle of the vibroseis truck, k is the climbing resistance coefficient, and Z is the sum of the weights of the vibroseis truck and the driver, with the unit of ton.

6. The method according to claim 5, characterized in that, the determination of the first connected area on the three-dimensional image data according to the climbing parameters of the vibroseis vehicle includes: Regard the area where the surface dip angle is less than or equal to the maximum climbing angle of the vibroseis and the minimum width is greater than or equal to the preset width threshold as the first connected area on the three-dimensional image data.

7. The method according to claim 2, characterized in that, obtain the climbing resistance coefficient through the following expression: k = (A * cosθ + m * g * cosθ * h) / (m * g * sinθ) wherein, A is the contact area between the tire of the vibroseis vehicle and the ground surface, θ is the surface dip angle, m is the mass of the vibroseis vehicle, g is the acceleration due to gravity, and h is the center of gravity height of the vibroseis vehicle.

8. The method according to claim 1, characterized in that, the adjustment of the first connected area to obtain the second connected area includes: For the first target disconnected area, when the number of vibroseis excitation points arranged after road construction in the first target disconnected area exceeds the first preset quantity threshold, the first target disconnected area is taken as the area to be road constructed, where the surface dip angle of the first target disconnected area is less than or equal to the maximum climbing angle of the vibroseis and the minimum width is less than the preset width threshold; For the first connected area surrounded by disconnected areas, it is taken as the third connected area; For the second target disconnected area, when the two first connected areas are on the same horizontal plane, the second target disconnected area is taken as the area to be road constructed, and when the two first connected areas are not on the same horizontal plane, the second target disconnected area with a slope less than the preset angle threshold is taken as the area to be road constructed, where the surface dip angle of the second target disconnected area is greater than the maximum climbing angle of the vibroseis and the minimum width is greater than or equal to the preset width threshold, and it is between two adjacent first connected areas where the number of vibroseis excitation points both exceed the second preset quantity threshold; The area to be road constructed is incorporated into the first connected area, and the third connected area is excluded from the first connected area to obtain the second connected area.

9. The method according to claim 8, wherein, For the second target disconnected area with a slope less than the preset angle threshold when the two first connected areas are not on the same horizontal plane, the road construction length of the second target disconnected area is calculated through the following steps: Calculate the height difference between the two road construction connection points at the top and bottom of the slope of the second target disconnected area; Through the following expression, calculate the longest horizontal distance in the horizontal direction between the two points: Among them, α max is the maximum climbing angle of the vibroseis truck, and i max is the actual maximum climbing gradient of the vibroseis after considering the friction coefficients of different terrains. H is the height difference between the top and bottom of the slope, and L is the longest horizontal distance in the direction perpendicular to the slope direction. Through the following expression, calculate the road construction length along the slope according to the calculated height difference and horizontal distance: where T is the shortest road construction length along the slope.

10. The method according to claim 1, wherein, The determining the vibroseis construction area and the downhole shot construction area according to the second connected area and the disconnected areas outside the second connected area on the three-dimensional image data respectively includes: Respectively taking the second connected area and the disconnected areas outside the second connected area on the three-dimensional image data as the vibroseis construction area and the downhole shot construction area, and forming the construction demarcation between the downhole shot and the vibroseis according to the vibroseis construction area and the downhole shot construction area; Layout the specific data points for the downhole shot and the vibroseis according to the designed acquisition system; At the demarcation between the downhole shot and the vibroseis, using the offset principle of the excitation points, under the condition that the regular offset is allowed, offset the points in the downhole shot construction area and design them as vibroseis excitation points to obtain the optimized construction demarcation between the downhole shot and the vibroseis vehicle, which is used to re-divide the vibroseis construction area and the downhole shot construction area on the three-dimensional image data.

11. A zoning device for combined excitation of downhole shots and vibroseis, wherein, comprising: An acquisition module, configured to acquire the elevation data, surface ground object data, slope data of the target mountain work area and the climbing parameters of the vibroseis vehicle; An overlay module, configured to overlay the elevation data, surface ground object data and slope data of the target mountain work area to obtain three-dimensional image data; An adjustment module, configured to determine a first connected region on three-dimensional image data according to the climbing parameters of a vibrator truck, and adjust the first connected region to obtain a second connected region; A determination module, configured to determine a vibrator construction area and a well shot construction area according to the second connected region on the three-dimensional image data and the unconnected regions outside the second connected region, respectively.

12. An electronic device, characterized in that it includes a processor, a communication interface, a memory, and a communication bus, wherein the processor, the communication interface, and the memory complete communication with each other through the communication bus; The memory is used for storing a computer program; The processor is configured to implement the zoning method for combined excitation of well shots and vibrators according to any one of claims 1-10 when executing the program stored on the memory.

13. A computer-readable storage medium, on which a computer program is stored, characterized in that the computer program, when executed by a processor, implements the zoning method for combined excitation of well shots and vibrators according to any one of claims 1-10.

Images