Partition method and device for joint shooting of well gun and controllable source

Abstract

The present disclosure relates to a partition method and device for joint excitation of well gun and controllable source, the method comprising: obtaining elevation data, surface feature data and slope data of a target mountain work area and climbing parameters of a controllable source vehicle; superimposing the elevation data, surface feature data and slope data of the target mountain work area to obtain three-dimensional image data; determining a first connected region on the three-dimensional image data according to the climbing parameters of the controllable source vehicle, and adjusting the first connected region to obtain a second connected region; determining a controllable source construction area and a well gun construction area according to the second connected region on the three-dimensional image data and the unconnected area outside the second connected region, respectively, which can accurately divide the mountain well gun and controllable source construction areas, improve the controllable source mountain construction proportion, ensure the safety of mountain seismic exploration acquisition projects, and achieve the purpose of reducing cost and increasing efficiency.

Description

Technical Field

[0001] This disclosure relates to the field of onshore oil seismic exploration data acquisition technology, and in particular to a zoning method and apparatus for combined excitation of well shot and controlled seismic source. Background Technology

[0002] Seismic exploration is a primary method for oil and gas exploration, and its main tasks include three stages: seismic data acquisition, processing, and interpretation. Seismic data acquisition is the first step in oil and gas seismic exploration engineering, primarily involving the artificial stimulation and reception of seismic waves to obtain raw seismic data containing subsurface geological information. Onshore seismic acquisition work mainly includes three parts: seismic wave stimulation, seismic wave reception, and observation methods. Common onshore seismic sources include controlled-source stimulation and explosive stimulation. Historically, explosive stimulation has been used in the seismic acquisition stage in mountainous areas. When targeting relatively shallow strata, a weaker acquisition scheme and a relatively low borehole density make data acquisition relatively easy, while keeping acquisition costs relatively controllable.

[0003] In recent years, with the increasing depth of exploration targets and growing interference from external human factors, acquiring deep and ultra-deep data has become increasingly difficult. Employing more robust seismic acquisition schemes and significantly increased shot density are currently the most effective methods for obtaining deep and ultra-deep data. However, significantly increased shot density means drilling more excitation wells, leading to higher costs. To ensure data quality and reduce acquisition costs, while also leveraging the advantages of low-frequency controlled-source seismic imaging of deep structures, lower-cost controlled-source excitation methods, particularly those used in mountainous terrain, can be employed to reduce the well-shot ratio.

[0004] With the development and advancement of UAV technology, UAV aerial surveying has been widely applied in the field of seismic exploration. Related technologies utilize UAV aerial survey imagery and DEM (digital terrain model) data to guide indoor reconnaissance, geomorphological zoning, indoor physical point layout, creation of driving route maps, and production management. However, this method lacks precision, and the creation of driving route maps is overly simplistic. Related technologies also propose using slope, contour lines, and bulldozer tracks to create acquisition route maps, providing the first detailed introduction to the process of designing field construction paths using aerial survey data. Further related technologies propose a "road construction first, then layout" route design operation process and a method for acquisition path design based on high-definition geographic imagery and high-precision elevation data obtained from UAV aerial surveying, solving the problem of uneven distribution of excitation points in desert areas. However, high-definition aerial survey imagery data for mountainous areas has not been effectively utilized. Summary of the Invention

[0005] To solve the above-mentioned technical problems, or at least partially solve them, embodiments of this disclosure provide a zoning method and apparatus for joint excitation of well shot and controllable seismic source.

[0006] In a first aspect, embodiments of this disclosure provide a zoning method for joint excitation of well shots and controlled seismic sources, comprising:

[0007] Acquire elevation data, surface feature data, and slope data of the target mountainous work area, as well as the climbing parameters of the controllable seismic source vehicle;

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

[0009] The first connected region on the three-dimensional image data is determined based on the climbing parameters of the controllable seismic source vehicle, and the first connected region is adjusted to obtain the second connected region.

[0010] Based on the second connected region and the non-connected regions outside the second connected region in the 3D image data, the controllable seismic source construction area and the well-blasting construction area are determined respectively.

[0011] In one possible implementation, the climbing parameters of the controllable seismic source vehicle are obtained through the following steps:

[0012] Obtain the tire radius of the controllable seismic source vehicle;

[0013] The maximum driving force of the controllable vibration source vehicle is determined based on the tire radius and engine parameters.

[0014] The climbing parameters of the controllable seismic source vehicle are determined based on the maximum driving force and climbing resistance coefficient of the vehicle.

[0015] In one possible implementation, the tire radius of the controllable vibration source vehicle is obtained using the following expression:

[0016]

[0017] Where R is the tire radius of the controlled seismic source vehicle in meters, D is the tire width of the controlled seismic source vehicle in inches, a is the tire aspect ratio of the controlled seismic source vehicle, and d is the wheel rim diameter of the controlled seismic source vehicle in inches.

[0018] In one possible implementation, the maximum driving force of the controllable vibration source vehicle is determined based on the tire radius and engine parameters of the vehicle using the following expression:

[0019]

[0020] Where C is the maximum driving force of the controllable vibration source vehicle, in Newtons (N), N is the maximum torque of the engine, in Newton-meters (N·m), B is the gear ratio, and W is the engine tail gear.

[0021] In one possible implementation, the climbing parameters of the controllable seismic source vehicle are determined based on the maximum driving force and climbing resistance coefficient of the controllable seismic source vehicle using the following expression:

[0022]

[0023] Where, α max The maximum climbing angle of the controllable seismic source vehicle is given by , k is the climbing resistance coefficient, and Z is the sum of the weights of the controllable seismic source vehicle and the driver, in tons.

[0024] In one possible implementation, determining the first connected region on the three-dimensional image data based on the climbing parameters of the controllable seismic source vehicle includes:

[0025] The region with a surface dip angle less than or equal to the maximum slope angle of the controllable seismic source and a minimum width greater than or equal to a preset width threshold is designated as the first connected region in the 3D image data.

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

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

[0028] Where A is the contact area between the tires of the controllable seismic source vehicle and the ground surface, θ is the ground surface inclination angle, m is the mass of the controllable seismic source vehicle, g is the gravitational acceleration, and h is the height of the center of gravity of the controllable seismic source vehicle.

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

[0030] For the first target unconnected area, if the number of controllable seismic source excitation points deployed after road construction in the first target unconnected area exceeds the first preset number threshold, the first target unconnected area is regarded as the area to be road constructed. The surface dip angle of the first target unconnected area is less than or equal to the maximum climbing angle of the controllable seismic source and the minimum width is less than the preset width threshold.

[0031] The first connected region, which is surrounded by non-connected regions, is considered the third connected region.

[0032] For the second target unconnected area, if the two first connected areas are on the same horizontal plane, the second target unconnected area is regarded as the road to be repaired area. If the two first connected areas are not on the same horizontal plane, the second target unconnected area with a slope less than a preset angle threshold is regarded as the road to be repaired area. The surface dip angle of the second target unconnected area is greater than the maximum climbing angle of the controllable seismic source 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 excitation points of the controllable seismic source both exceed the second preset number threshold.

[0033] The area to be repaired is included in the first connected region, and the third connected region is separated from the first connected region to obtain the second connected region.

[0034] In one possible implementation, for a second target disconnected region where the slope is less than a preset angle threshold and the two first connected regions are not on the same horizontal plane, the road repair length of the second target disconnected region is calculated through the following steps:

[0035] Calculate the elevation difference between the two road construction connection points at the top and bottom of the slope in the disconnected region of the second objective;

[0036] The following expression can be used to calculate the longest horizontal distance between two points:

[0037]

[0038] Where, α max For the maximum climbing angle of the controllable seismic source vehicle, i max This represents the actual maximum gradient of the controllable seismic source after considering the friction coefficients of different terrains. H is the elevation difference between the top and bottom of the slope, and L is the longest horizontal distance perpendicular to the slope direction.

[0039] The length of the road along the slope can be calculated using the following expression, based on the calculated elevation difference and horizontal distance:

[0040]

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

[0042] In one possible implementation, determining the controllable seismic source construction area and the well-blasting construction area based on the second connected region and the non-connected region outside the second connected region in the three-dimensional image data includes:

[0043] The second connected region and the non-connected region outside the second connected region in the three-dimensional image data are respectively regarded as the controllable source construction area and the well shot construction area, and the construction boundary between the well shot and the controllable source is formed according to the controllable source construction area and the well shot construction area.

[0044] The well shot and the controllable seismic source are deployed with specific data points according to the designed observation system;

[0045] At the boundary between the well shot and the controlled seismic source, using the principle of excitation point offset, and under the condition that regular offset is allowed, the points in the well shot construction area are offset and designed as controlled seismic source excitation points, resulting in an optimized construction boundary between the well shot and the controlled seismic source vehicle, which is used to re-divide the controlled seismic source construction area and the well shot construction area on the three-dimensional image data.

[0046] Secondly, embodiments of this disclosure provide a zoned device for combined excitation of well shot and a controllable seismic source, comprising:

[0047] The acquisition module is used to acquire elevation data, surface feature data, and slope data of the target mountainous work area, as well as the climbing parameters of the controllable seismic source vehicle.

[0048] The overlay module is used to overlay elevation data, surface feature data and slope data of the target mountainous work area to obtain three-dimensional image data;

[0049] The adjustment module is used to determine the first connected region on the three-dimensional image data based on the climbing parameters of the controllable seismic source vehicle, and adjust the first connected region to obtain the second connected region;

[0050] The determination module is used to determine the controllable seismic source construction area and the well-shot construction area based on the second connected region and the non-connected region outside the second connected region in the three-dimensional image data.

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

[0052] Memory, used to store computer programs;

[0053] The processor, when executing the program stored in memory, implements the aforementioned partitioning method for joint excitation of well shots and controllable seismic sources.

[0054] Fourthly, embodiments of this disclosure provide a computer-readable storage medium having a computer program stored thereon, characterized in that the computer program, when executed by a processor, implements the aforementioned partitioning method for joint excitation of well shot and controllable seismic source.

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

[0056] The zoning method for combined well-blasting and controlled-source seismic excitation described in this embodiment acquires elevation data, surface feature data, and slope data of the target mountainous work area, as well as the climbing parameters of the controlled-source vehicle. The elevation data, surface feature data, and slope data of the target mountainous work area are overlaid to obtain three-dimensional image data. A first connected region on the three-dimensional image data is determined based on the climbing parameters of the controlled-source vehicle, and the first connected region is adjusted to obtain a second connected region. The controlled-source construction area and the well-blasting construction area are determined based on the second connected region and the non-connected areas outside the second connected region on the three-dimensional image data, respectively. This method can accurately divide the well-blasting and controlled-source construction areas in mountainous areas, increase the proportion of controlled-source construction in mountainous areas, ensure the safety of seismic exploration and acquisition projects in mountainous regions, and achieve the goal of cost reduction and efficiency improvement. Attached Figure Description

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

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

[0059] Figure 1 This schematically illustrates a flow chart of a zoning method for joint excitation of well shot and a controllable seismic source according to an embodiment of the present disclosure;

[0060] Figure 2 A schematic diagram of three-dimensional surface feature data according to an embodiment of the present disclosure is shown.

[0061] Figure 3 A schematic diagram of three-dimensional elevation data according to an embodiment of the present disclosure is shown;

[0062] Figure 4 A schematic diagram illustrating three-dimensional slope data according to an embodiment of the present disclosure is shown.

[0063] Figure 5 A schematic diagram of superimposed three-dimensional image data according to an embodiment of the present disclosure is shown.

[0064] Figure 6 This schematically illustrates a first connected region in superimposed three-dimensional image data according to an embodiment of the present disclosure;

[0065] Figure 7 A schematic diagram illustrating a region where a mountain gully narrows and obstructs connectivity according to an embodiment of the present disclosure;

[0066] Figure 8 A schematic diagram illustrating local obstacles within a large connected region according to an embodiment of the present disclosure is shown.

[0067] Figure 9 This schematic diagram illustrates a road construction trajectory for steep slopes designed with maximum gradient according to an embodiment of the present disclosure.

[0068] Figure 10 The diagram illustrates the partitioning of well shots and controllable seismic sources according to embodiments of the present disclosure.

[0069] Figure 11 A schematic diagram illustrates the structural block diagram of a zoning device for combined excitation of well shot and a controllable seismic source according to an embodiment of the present disclosure;

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

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

[0072] See Figure 1 The embodiments of this disclosure provide a zoning method for joint excitation of well shot and controlled seismic source, including the following steps:

[0073] S1 acquires elevation data, surface feature data, and slope data of the target mountainous work area, as well as the climbing parameters of the controllable seismic source vehicle.

[0074] like Figure 2 As shown, the surface feature data is high-precision imagery obtained through low-altitude drone aerial photography, image processing software, and indoor registration and correction. The image displays surface feature information, particularly information on different landform types in mountainous areas, ground equipment, and vehicle routes along old survey lines in the work area. The surface feature data resolution is greater than or equal to 0.1m x 0.1m. Figure 3 As shown, the elevation data is high-precision elevation imagery obtained through low-altitude aerial photography by a drone and processed by image processing software. The image shows the undulations of the ground surface, and the elevation data resolution is greater than or equal to 1m x 1m. Figure 4 As shown, the slope data is calculated from high-precision elevation image data to reflect the slope magnitude reflecting the undulation of the ground surface. The steepness of the mountains can be seen from the figure.

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

[0076] In this embodiment, the overlay can be performed as follows: the first layer is elevation data, with transparency set to 0%; the second layer is surface feature data, with transparency set to 0%; the third layer is slope data, divided into the following slope segments: below the maximum ramp angle of the controllable seismic source; between the maximum ramp angle of the controllable seismic source and 30 degrees; and above 30 degrees. Slopes below the maximum ramp angle of the controllable seismic source are set to colorless; for the other slope segments, the colors are set from light to dark, from smallest to largest, with the transparency of all slope segments set to 30%. This overlay method takes into account the information from each layer, forming a 3D high-definition image with a three-dimensional and intuitive feel.

[0077] See Figure 5 The vehicle trajectory can be overlaid as the fourth layer with the first, second, and third layers to obtain three-dimensional image data. The bottom layer of the three-dimensional image data is the elevation data layer, the second layer is the surface feature data overlaid on the elevation data layer, the third layer is the slope layer, and the fourth layer is the vehicle trajectory. The vehicle trajectory can be formed by picking the coordinates of gravel roads and cement roads on gentle slopes in mountain gullies and surface feature data.

[0078] S3. Determine the first connected region on the three-dimensional image data based on the climbing parameters of the controllable seismic source vehicle, and adjust the first connected region to obtain the second connected region.

[0079] In this embodiment, step S3, determining the first connected region on the three-dimensional image data based on the climbing parameters of the controllable seismic source vehicle, includes:

[0080] The region with a surface dip angle less than or equal to the maximum climbing angle of the controllable seismic source and a minimum width greater than or equal to a preset width threshold is designated as the first connected region in the 3D image data. Typically, the width of the controllable seismic source vehicle is 3.5m; to ensure safe passage of the controllable seismic source, the preset width threshold can be 4m.

[0081] like Figure 6 As shown, in mountain gullies, areas with continuous sampling points whose slope is equal to or whose width is greater than or equal to the maximum climbing angle of the controllable seismic source are formed into the first connected zone, and controllable seismic source construction is carried out.

[0082] In this embodiment, the region with a surface dip angle less than or equal to the maximum ramp angle of the controllable seismic source and a minimum width greater than or equal to a preset width threshold is defined as the first connected region in the three-dimensional image data, including:

[0083] Establish a 0.5m×0.5m grid, where the coordinate system of the grid data points is the Beijing 54 coordinate system, and the coordinate data includes station number, east coordinate, north coordinate, and ground elevation;

[0084] Create data points on each grid node;

[0085] The slope value of the high-resolution image is appended to each data point;

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

[0087] Data points with slopes less than the maximum climbing angle of the controllable seismic source are extracted in the longitudinal (lateral) direction respectively. Eight consecutive points (i.e., 4m) are set as the area that the controllable seismic source can pass through, and are set as the first connected region.

[0088] The location of discontinuous points or fewer than 8 consecutive excitation points (i.e., 4m) is set as the road repair point. Further judgment is needed on whether road repair is required, thereby further expanding the connectivity area. The road repair point is selected to accurately delineate the boundary between mountain well shots and controllable seismic sources.

[0089] In this embodiment, the maximum climbing angle refers to the maximum slope that the controllable seismic source can climb on a smooth road surface when driven by maximum driving force. It represents the ultimate climbing capability of the controllable seismic source vehicle. Connecting regions with slopes less than the maximum climbing angle of the controllable seismic source constitutes the controllable seismic source construction area; this region is called the connected zone.

[0090] S4. Based on the second connected region and the non-connected region outside the second connected region in the three-dimensional image data, the controllable seismic source construction area and the well-blasting construction area are determined respectively.

[0091] In this embodiment, the controlled source construction area and the well-blasting construction area are used in well-seismic co-mining. Well-seismic co-mining refers to the construction of a seismic acquisition project using a combination of well-blasting and controlled source excitation, which can also be called well-seismic co-excitation.

[0092] In this embodiment, the climbing ability of the controllable seismic source vehicle in mountainous areas is closely related to its maximum driving force, weight, tire width, and radius. The maximum driving force of the seismic source vehicle is closely related to its power system, engine power, drive unit, and transmission system. Generally, the greater the power, the more powerful the drive unit, and the more efficient the transmission system, the greater the maximum driving force of the seismic source vehicle, and its climbing ability in mountainous areas will correspondingly improve. In this embodiment, in step S1, the climbing parameters of the controllable seismic source vehicle are obtained through the following steps:

[0093] Obtain the tire radius of the controllable seismic source vehicle;

[0094] The maximum driving force of the controllable vibration source vehicle is determined based on the tire radius and engine parameters.

[0095] The climbing parameters of the controllable seismic source vehicle are determined based on the maximum driving force and climbing resistance coefficient of the vehicle.

[0096] In this embodiment, the tire radius of the controllable vibration source vehicle is obtained using the following expression:

[0097]

[0098] Where R is the tire radius of the controlled seismic source vehicle in meters, D is the tire width in inches, a is the tire aspect ratio, and d is the wheel diameter in inches. Typically, the tire width of the seismic source vehicle is 43 inches, the aspect ratio is 0.4767, and the wheel diameter is 25 inches. 1 inch = 25.4 millimeters.

[0099] In this embodiment, the maximum driving force of the controllable vibration source vehicle is determined based on the tire radius and engine parameters of the vehicle using the following expression:

[0100]

[0101] Where C is the maximum driving force of the controllable vibration source vehicle, in Newtons (N), N is the maximum engine torque, in Newton-meters (N·m), B is the gear ratio, and W is the engine final drive. Typically, the engine's maximum torque is 2000 N·m, the gear ratio is 4.5, and the engine final drive is 8.42.

[0102] In this embodiment, a slope resistance coefficient k is set, which is related to factors such as surface lithology and surface softness. Slope climbing tests are conducted on different mountain terrains to obtain the slope resistance coefficients of the controllable seismic source on different surfaces. Based on these slope resistance coefficients, combined with the maximum driving force and weight of the controllable seismic source, the slope climbing parameters of the controllable seismic source vehicle are determined using the following expression, based on the maximum driving force and slope resistance coefficient:

[0103]

[0104] Where, α max The maximum climbing angle of the controllable seismic source vehicle is given by denoted by denoted by k, the climbing resistance coefficient is given by denoted by Z, and the weight of the controllable seismic source vehicle and the driver is given by denoted by Z in tons. Typically, Z is 28.1 tons, and the force of 1 kg is 9.8 Newtons.

[0105] In this embodiment, the dynamic testing method is used to calculate the controllable seismic source uphill resistance coefficient k for different ground surfaces:

[0106] The first step is to prepare the testing equipment and surface samples. Prepare a controllable seismic source vehicle as the testing vehicle and install suitable testing equipment, such as force sensors and accelerometers. Simultaneously, prepare surface samples of different surface types, such as loess, gravel, and sandstone, ensuring their area is large enough to guarantee the accuracy and reliability of the test.

[0107] The second step is to conduct field tests. The test vehicle is driven onto different surface samples to conduct motion tests. By controlling parameters such as vehicle speed, acceleration, and tire pressure, the vehicle's motion on different surfaces can be recorded, including acceleration, speed, and tire force. Simultaneously, the climbing angle of the seismic source vehicle on different surfaces is recorded.

[0108] The third step is to process the test data. The obtained test data is processed, such as calculating parameters like friction, speed, and acceleration of the vehicle on different ground surfaces. Statistical analysis methods are used to obtain an estimated value for the controllable seismic source's climbing resistance coefficient.

[0109] The fourth step is to calculate the ramp resistance coefficient of the controllable seismic source. Using the calculation formula, the parameters obtained from the tests are substituted into the calculation to obtain estimated values ​​of the ramp resistance coefficient for different surface types. The ramp resistance coefficient is obtained through the following expression:

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

[0111] Where A is the contact area between the tires of the controllable seismic source vehicle and the ground surface, θ is the ground surface inclination angle, m is the mass of the controllable seismic source vehicle, g is the gravitational acceleration, and h is the height of the center of gravity of the controllable seismic source vehicle.

[0112] Based on multiple rounds of experimental data processing and calculation results, the slope resistance coefficients of controllable seismic sources on mountain loess, grassland, gravel, and rock surfaces were obtained, as shown in Table 1 below.

[0113] Table 1

[0114]

[0115]

[0116] In this embodiment, step S3, adjusting the first connected region to obtain the second connected region, includes:

[0117] For the first target unconnected area, if the number of controllable seismic source excitation points deployed after road construction in the first target unconnected area exceeds the first preset number threshold, the first target unconnected area is regarded as the area to be road constructed. The surface dip angle of the first target unconnected area is less than or equal to the maximum climbing angle of the controllable seismic source and the minimum width is less than the preset width threshold.

[0118] The first connected region, which is surrounded by non-connected regions, is considered the third connected region.

[0119] For the second target unconnected area, if the two first connected areas are on the same horizontal plane, the second target unconnected area is regarded as the road to be repaired area. If the two first connected areas are not on the same horizontal plane, the second target unconnected area with a slope less than a preset angle threshold is regarded as the road to be repaired area. The surface dip angle of the second target unconnected area is greater than the maximum climbing angle of the controllable seismic source 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 excitation points of the controllable seismic source both exceed the second preset number threshold.

[0120] The area to be repaired is included in the first connected region, and the third connected region is separated from the first connected region to obtain the second connected region.

[0121] In practical applications, the primary target non-connected area can be a region where a mountain gully narrows, obstructing connectivity. If a gully or other connecting area narrows to less than 4 meters, a controllable seismic source cannot pass through. It is necessary to determine whether road construction is needed to connect the area, thereby defining the extent of the connecting zone. The specific method is as follows:

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

[0123] Create data points on each grid node;

[0124] The slope value of the high-resolution image is appended to each data point;

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

[0126] Data points with slopes less than the maximum climbing angle of the controllable seismic source are extracted in the longitudinal (lateral) direction respectively. Eight consecutive points (i.e., 4m) are set as the area that the controllable seismic source can pass through, and are set as the first connected region.

[0127] Locations with discontinuities or fewer than 8 consecutive excitation points (i.e., 4m) are designated as road repair points. Further assessment is needed to determine if road repair is necessary, thereby expanding the connectivity zone. Road repair points are optimized to precisely delineate the boundary between mountain well-blasting and controlled seismic sources. The determination of road repair points primarily considers the cost of repair. Firstly, the area blocked by the road repair point, with a slope angle smaller than the maximum climb angle of the controlled seismic source, must be sufficiently large to ensure that at least 10 excitation points can be deployed after repair. If this condition is met, the road repair point is designated as a repairable point, and the blocked area is designated as a connectivity zone, considered a controlled seismic source construction area. If the condition is not met, the road repair point is designated as an unrepairable point, automatically closed, and the blocked area is automatically designated as a well-blasting construction area.

[0128] See Figure 7 The base image is Figure 5 The layer has two large, separated connected areas separated by narrow gullies. The location of the road construction point is determined and connected.

[0129] In practical applications, the first connected region surrounded by unconnected areas can be mountainous terrain with steep slopes, meaning that a partially connected region exists within a large unconnected area, but cannot be directly connected. In large, steeply sloping mountainous areas, there may be localized flat areas with gentle slopes where construction is possible using a controllable seismic source. Such terrain is clearly unsuitable for road construction and cannot be connected to adjacent connected regions. This type of terrain is designated as an obstacle area, with obstacles drawn. No road construction points are set within the obstacle area, thus closing the connected region. Well firing points are permitted to be deployed within the obstacle area.

[0130] In practical applications, obstacles with steep slopes that obstruct connectivity and are surrounded by connected regions are considered non-connected regions. That is, within the area completely enclosed by the connected region, there are local terrains with steep slopes, such as steep riverbanks, dams, ditches, isolated high and steep mountains, tall buildings, etc. (a) obstacles with a longitudinal or transverse width of less than 250m; b) obstacles with both longitudinal and transverse widths of less than 500m. It is impossible to set up excitation points in these areas. Such terrains are designated as obstacle regions, and the obstacles are drawn. The obstacle regions are contained within the connected regions. No road construction points are set up within the obstacles, and controllable seismic source excitation points are not allowed to be set up.

[0131] See Figure 8 Local obstacles within a large connected region, the base map is Figure 5 The layer defines the area within the connected region where controllable seismic sources cannot be deployed as an obstacle area. Road construction points are not set within the obstacles, and controllable seismic source excitation points are not allowed to be deployed.

[0132] In practical applications, the second non-connected area can be used to isolate adjacent large connected areas. Two large connected areas (each with at least 20 firing points) are separated by an obstacle with a large slope (greater than the maximum climbing angle of the controllable source), resulting in a situation where they cannot be connected. ① If the two connected areas are on the same horizontal plane, they can be connected according to the vehicle trajectory. 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), with the top of the slope being relatively flat, they cannot be connected. It is necessary to determine whether they can be connected. First, determine the slope range of the steep slope separating the two connected areas: a. If the slope of the barrier is greater than 30 degrees, it is considered a point where a road cannot be built, and the connected area in this area is automatically isolated and judged as unconnectable. The connected area on the mountaintop is automatically closed and set as a mountain well-blasting construction area, where well-blasting firing points can be deployed. b. If the slope of the barrier is less than 30 degrees, a road construction route can be designed for this steep slope. Among them, the slope-direction connectivity design is used to connect areas that are separated by steep mountain slopes.

[0133] In this embodiment, for terrain where two connecting areas not on the same horizontal plane are separated by a steep mountain slope with a slope less than 30 degrees, requiring road construction for connection, a targeted slope-direction connection design is needed. A lateral route can be designed to be constructed along the slope at a certain angle (β). This reduces the climbing angle of the controllable seismic source vehicle along the road construction direction. The shortest design route is chosen so that the climbing angle of the controllable seismic source vehicle along the road construction direction is lower than its maximum climbing angle, ensuring the controllable seismic source can pass. The maximum climbing angle (i) of the controllable seismic source vehicle is then utilized. max The following steps are used to calculate the road length. First, calculate the elevation difference (H) between the two road connection points at the top and bottom of the slope. Second, calculate the longest horizontal distance (L) between the two points using the following expression. Finally, calculate the road length T along the slope based on the calculated elevation difference and horizontal distance using the following expression:

[0134]

[0135]

[0136] Where, α max For the maximum climbing angle of the controllable seismic source vehicle, i max It is the actual maximum gradient of the controllable seismic source after taking into account 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, and T is the shortest road length along the slope.

[0137] Therefore, as long as it can be guaranteed This allows for the design of the longest possible T-shaped road, ensuring that after road construction, the ramp angle of the controllable seismic source along the construction direction is less than the maximum ramp angle of the controllable seismic source, allowing the controllable seismic source to pass through and thus connecting the two connected regions, further refining the extent of the connected regions. Conversely, even if the designed road route is the longest possible, it still cannot meet the requirements. This condition indicates that the controllable seismic source cannot pass through after the road is repaired. The road repair point can be cancelled, the connecting area on the mountaintop can be closed, and it can be set as a mountain well-blasting construction area, where well-blasting trigger points can be deployed.

[0138] See Figure 9 The road construction trajectory on steep slopes was designed based on the maximum gradient, and the base map is as follows. Figure 5 The layer is used to design the gentlest road construction trajectory on the steep mountain slope between two large connected areas, ensuring the connection between the two connected areas and ensuring that the controllable seismic source can carry out normal construction.

[0139] In this embodiment, step S4, determining the controllable seismic source construction area and the well-blasting construction area based on the second connected region and the non-connected regions outside the second connected region in the three-dimensional image data, includes:

[0140] The second connected region and the non-connected region outside the second connected region in the three-dimensional image data are respectively regarded as the controllable source construction area and the well shot construction area, and the construction boundary between the well shot and the controllable source is formed according to the controllable source construction area and the well shot construction area.

[0141] The well shot and the controllable seismic source are deployed with specific data points according to the designed observation system;

[0142] At the boundary between the well shot and the controlled seismic source, using the principle of excitation point offset, under the condition that regular offset is allowed, points in the well shot construction area are offset and designed as excitation points of the controlled seismic source. This increases the number of controlled seismic source construction areas and excitation points, reduces the boundary of well shot construction in mountainous areas, and obtains an optimized construction boundary between the well shot and the controlled seismic source vehicle, thus achieving a fine division of the well shot and controlled seismic source construction areas.

[0143] See Figure 10 The zoning is divided into well-blasting and controlled-source zones. The black area is the well-blasting construction area in mountainous terrain, where well-blasting is used for excitation. The remaining areas are the controlled-source construction areas, where controlled-source excitation is used.

[0144] This disclosure presents a zoning method for the combined excitation of well drilling and controlled seismic sources. It proposes a mountain well-seismic design method based on 3D stereoscopic overlay high-definition imagery, a calculation method for the "slope climbing ability of controlled seismic sources in mountainous areas," a design method for controlled seismic source connectivity zones, and a concept for indoor design of road construction points. It also presents a novel method and process for precise well-seismic zoning in mountainous areas, replacing the traditional "human-led" field reconnaissance followed by design. This method ensures accurate well-seismic zoning and high precision in excitation point design. Furthermore, addressing the issue of precise well-seismic design in mixed well-seismic mining operations in mountainous areas, a new well-seismic zoning design method and process are proposed. This method utilizes surface feature data (accuracy ≥ 0.1m), high-precision elevation data (accuracy ≥ 1m), and slope data to form a 3D high-definition image. Combined with the controllable seismic source connectivity zones, bulldozer road construction tracks, and reconnaissance tracks, the controllable seismic source construction area in mountainous areas is precisely delineated, and well-seismic boundaries are determined.

[0145] This disclosed zoning method for combined well-shot and controlled-source seismic acquisition addresses the problems of high cost, low efficiency, and poor data quality associated with simple well-shot seismic acquisition in mountainous areas. It utilizes surface feature data, elevation data, slope data, and vehicle travel routes, combined with the mountain climbing ability of the controlled-source seismic acquisition and the aspect connectivity of large-area connected areas in the piedmont zone that are interrupted by steep mountain slopes, to perform zoning for combined well-shot and controlled-source seismic acquisition in mountainous areas. This method effectively supports the economy, safety, and environmental friendliness of complex mountainous seismic acquisition projects, significantly reducing costs and increasing efficiency, thus improving the economic feasibility of seismic acquisition projects in mountainous areas. The method proposed in this invention has been verified in a mountainous seismic acquisition project in western China.

[0146] See Figure 11 The embodiments of this disclosure provide a zoning device for combined excitation of well shot and controlled seismic source, comprising:

[0147] The acquisition module 11 is used to acquire elevation data, surface feature data, and slope data of the target mountainous work area, as well as the climbing parameters of the controllable seismic source vehicle.

[0148] The overlay module 12 is used to overlay the elevation data, surface feature data and slope data of the target mountainous work area to obtain three-dimensional image data;

[0149] The adjustment module 13 is used to determine the first connected region on the three-dimensional image data based on the climbing parameters of the controllable seismic source vehicle, and adjust the first connected region to obtain the second connected region.

[0150] The determination module 14 is used to determine the controllable seismic source construction area and the well-blasting construction area based on the second connected region and the non-connected region outside the second connected region in the three-dimensional image data.

[0151] The zoning device for combined well drilling and controlled seismic source excitation disclosed herein accurately delineates the construction areas of well drilling and controlled seismic source vehicles in well-seismic mixed mining projects, increases the proportion of controlled seismic source construction in mountainous areas, improves acquisition efficiency, and reduces acquisition costs.

[0152] This disclosure discloses a zoning device for the combined excitation of well shot and controlled seismic source to produce 3D stereoscopic high-definition images, experimentally quantifies the slope climbing ability of the controlled seismic source, establishes and delineates connected zones, designs slope-direction connectivity of connected zones separated by steep mountain slopes, and precisely delineates the construction areas of well shot and controlled seismic source. It focuses on improving efficiency and effectiveness in the field, and studies how to implement methods for dividing the location of well shot excitation in mountainous areas and methods to increase the proportion of controlled seismic sources in mountainous areas.

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

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

[0155] In the above embodiments, any and more of the acquisition module 11, superposition module 12, adjustment module 13, and determination module 14 can be combined into one module, or any one of these modules can be split into multiple modules. Alternatively, at least part of the functionality of one or more of these modules can be combined with at least part of the functionality of other modules and implemented in one module. At least one of the acquisition module 11, superposition module 12, adjustment module 13, and determination module 14 can be at least partially implemented as hardware circuitry, 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-package, an application-specific integrated circuit (ASIC), or implemented in hardware or firmware by any other reasonable means of integrating or packaging the circuitry, or implemented in software, hardware, or firmware, or in any appropriate combination of any of these three implementation methods. Alternatively, at least one of the acquisition module 11, superposition module 12, adjustment module 13, and determination module 14 can be at least partially implemented as a computer program module, which can perform corresponding functions when the computer program module is run.

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

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

[0158] When processor 1110 executes the program stored in memory 1130, it implements the following partitioning method for joint excitation of well shots and controlled seismic sources:

[0159] Acquire elevation data, surface feature data, and slope data of the target mountainous work area, as well as the climbing parameters of the controllable seismic source vehicle;

[0160] The elevation data, surface feature data and slope data of the target mountainous work area are overlaid to obtain three-dimensional image data;

[0161] The first connected region on the three-dimensional image data is determined based on the climbing parameters of the controllable seismic source vehicle, and the first connected region is adjusted to obtain the second connected region.

[0162] Based on the second connected region and the non-connected regions outside the second connected region in the 3D image data, the controllable seismic source construction area and the well-blasting construction area are determined respectively.

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

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

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

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

[0167] Embodiments of this disclosure also provide a computer-readable storage medium. The computer-readable storage medium stores a computer program that, when executed by a processor, implements the zoning method for joint excitation of well shots and controlled seismic sources as described above.

[0168] The computer-readable storage medium may be included in the device / apparatus described in the above embodiments; or it may exist independently and not assembled into the device / apparatus. The computer-readable storage medium carries one or more programs that, when executed, implement the zoning method for joint excitation of well shot and controlled seismic source according to the embodiments of this disclosure.

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

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

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

Claims

1. A zoning method for combined excitation of well shots and controlled seismic sources, characterized in that, The method includes: Acquire elevation data, surface feature data, and slope data of the target mountainous work area, as well as the climbing parameters of the controllable seismic source vehicle; The elevation data, surface feature data and slope data of the target mountainous work area are overlaid to obtain three-dimensional image data; The first connected region on the three-dimensional image data is determined based on the climbing parameters of the controllable seismic source vehicle, and the first connected region is adjusted to obtain the second connected region. Based on the second connected region and the non-connected region outside the second connected region in the three-dimensional image data, the controllable seismic source construction area and the well-blasting construction area are determined respectively. The step of determining the first connected region on the three-dimensional image data based on the climbing parameters of the controllable seismic source vehicle includes: The region with a surface dip angle less than or equal to the maximum slope angle of the controllable seismic source and a minimum width greater than or equal to a preset width threshold is designated as the first connected region in the three-dimensional image data. The adjustment of the first connected component to obtain the second connected component includes: For the first target unconnected area, if the number of controllable seismic source excitation points deployed after road construction in the first target unconnected area exceeds the first preset number threshold, the first target unconnected area is regarded as the area to be road constructed. The surface dip angle of the first target unconnected area is less than or equal to the maximum climbing angle of the controllable seismic source and the minimum width is less than the preset width threshold. The first connected region, which is surrounded by non-connected regions, is considered the third connected region. For the second target unconnected area, if the two first connected areas are on the same horizontal plane, the second target unconnected area is regarded as the road to be repaired area. If the two first connected areas are not on the same horizontal plane, the second target unconnected area with a slope less than a preset angle threshold is regarded as the road to be repaired area. The surface dip angle of the second target unconnected area is greater than the maximum climbing angle of the controllable seismic source 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 excitation points of the controllable seismic source both exceed the second preset number threshold. The area to be repaired is included in the first connected region, and the third connected region is separated from the first connected region to obtain the second connected region.

2. The method according to claim 1, characterized in that, The climbing parameters of the controllable seismic source vehicle are obtained through the following steps: Obtain the tire radius of the controllable seismic source vehicle; The maximum driving force of the controllable vibration source vehicle is determined based on the tire radius and engine parameters. The climbing parameters of the controllable seismic source vehicle are determined based on the maximum driving force and climbing resistance coefficient of the vehicle.

3. The method according to claim 2, characterized in that, The tire radius of the controllable vibration source vehicle can be obtained using the following expression: Where R is the tire radius of the controlled seismic source vehicle, in meters, and D is the tire width of the controlled seismic source vehicle, in inches. The aspect ratio of the tires of the controllable vibration source vehicle The diameter of the wheel hub of the controllable seismic source vehicle is in inches.

4. The method according to claim 2, characterized in that, The maximum driving force of the controllable vibration source vehicle is determined using the following expression, based on the tire radius and engine parameters: Where C is the maximum driving force of the controllable vibration source vehicle, in Newtons (N), N is the maximum torque of the engine, in Newton-meters (N·m), B is the gear ratio, and W is the engine tail gear.

5. The method according to claim 2, characterized in that, The climbing parameters of the controllable seismic source vehicle are determined based on its maximum driving force and climbing resistance coefficient using the following expressions: in, The maximum climbing angle of the controllable seismic source vehicle is given by , k is the climbing resistance coefficient, and Z is the sum of the weights of the controllable seismic source vehicle and the driver, in tons.

6. The method according to claim 2, characterized in that, The climbing resistance coefficient can be obtained using the following expression: Where A is the contact area between the tires of the controllable seismic source vehicle and the ground surface, θ is the ground surface inclination angle, m is the mass of the controllable seismic source vehicle, g is the gravitational acceleration, and h is the height of the center of gravity of the controllable seismic source vehicle.

7. The method according to claim 1, characterized in that, For a second target disconnected region where the two first connected regions are not on the same horizontal plane and the slope is less than a preset angle threshold, the road repair length of the second target disconnected region is calculated using the following steps: Calculate the elevation difference between the two road construction connection points at the top and bottom of the slope in the disconnected region of the second objective; The following expression can be used to calculate the longest horizontal distance between two points: in, The maximum climbing angle for the controllable seismic source vehicle. This is the actual maximum gradient of the controllable seismic source after taking into account the friction coefficient of different terrains. This is the height difference between the top and bottom of the slope. It is the longest horizontal distance in the direction perpendicular to the slope. The length of the road along the slope can be calculated using the following expression, based on the calculated elevation difference and horizontal distance: in, This represents the shortest road length along the slope.

8. The method according to claim 1, characterized in that, The determination of the controllable seismic source construction area and the well-blasting construction area based on the second connected region and the non-connected regions outside the second connected region in the three-dimensional image data includes: The second connected region and the non-connected region outside the second connected region in the three-dimensional image data are respectively regarded as the controllable source construction area and the well shot construction area, and the construction boundary between the well shot and the controllable source is formed according to the controllable source construction area and the well shot construction area. The well shot and the controllable seismic source are deployed with specific data points according to the designed observation system; At the boundary between the well shot and the controlled seismic source, using the principle of excitation point offset, and under the condition that regular offset is allowed, the points in the well shot construction area are offset and designed as controlled seismic source excitation points, resulting in an optimized construction boundary between the well shot and the controlled seismic source vehicle, which is used to re-divide the controlled seismic source construction area and the well shot construction area on the three-dimensional image data.

9. A zoned excitation device for combined well shot and controllable seismic source, characterized in that, include: The acquisition module is used to acquire elevation data, surface feature data, and slope data of the target mountainous work area, as well as the climbing parameters of the controllable seismic source vehicle. The overlay module is used to overlay elevation data, surface feature data and slope data of the target mountainous work area to obtain three-dimensional image data; The adjustment module is used to determine the first connected region on the three-dimensional image data based on the climbing parameters of the controllable seismic source vehicle, and adjust the first connected region to obtain the second connected region; The determination module is used to determine the controllable seismic source construction area and the well-blasting construction area based on the second connected region and the non-connected region outside the second connected region in the three-dimensional image data, respectively. The step of determining the first connected region on the three-dimensional image data based on the climbing parameters of the controllable seismic source vehicle includes: The region with a surface dip angle less than or equal to the maximum slope angle of the controllable seismic source and a minimum width greater than or equal to a preset width threshold is designated as the first connected region in the three-dimensional image data. The adjustment of the first connected component to obtain the second connected component includes: For the first target unconnected area, if the number of controllable seismic source excitation points deployed after road construction in the first target unconnected area exceeds the first preset number threshold, the first target unconnected area is regarded as the area to be road constructed. The surface dip angle of the first target unconnected area is less than or equal to the maximum climbing angle of the controllable seismic source and the minimum width is less than the preset width threshold. The first connected region, which is surrounded by non-connected regions, is considered the third connected region. For the second target unconnected area, if the two first connected areas are on the same horizontal plane, the second target unconnected area is regarded as the road to be repaired area. If the two first connected areas are not on the same horizontal plane, the second target unconnected area with a slope less than a preset angle threshold is regarded as the road to be repaired area. The surface dip angle of the second target unconnected area is greater than the maximum climbing angle of the controllable seismic source 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 excitation points of the controllable seismic source both exceed the second preset number threshold. The area to be repaired is included in the first connected region, and the third connected region is separated from the first connected region to obtain the second connected region.

10. An electronic device, characterized in that, It includes a processor, a communication interface, a memory, and a communication bus, wherein the processor, the communication interface, and the memory communicate with each other through the communication bus; Memory, used to store computer programs; The processor, when executing a program stored in memory, implements the partitioning method for joint excitation of well shot and controllable seismic source as described in any one of claims 1-8.

11. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the zoning method for joint excitation of well shot and controllable seismic source as described in any one of claims 1-8.

Images