Method and device for determining explosive column

Abstract

The present application provides a method and device for determining an explosive column. The method includes: constructing a geometric model based on the well depth and well diameter of a target well, assembling the geometric model based on the explosion parameters of different types of explosives to form at least two sets of explosion models, where each set of explosion models configured with an explosive column contains explosives of at least two explosive types, respectively performing explosion simulations on at least two sets of explosion models to obtain the energy distribution corresponding to each set of explosion models, for each set of explosion models, determining the objective function value corresponding to the effective area according to the energy distribution corresponding to the explosion model, and determining the target explosive column applicable to the target well according to the objective function values corresponding to at least two sets of explosion models. Through the present application, an explosive column applicable to single-well seismic exploration can be determined, greatly saving costs, and capable of enhancing the seismic energy excited, thereby enhancing the signal-to-noise ratio of seismic data.

Description

Technical Field

[0001] The present application relates to the technical field of seismic exploration, and in particular to a method and device for determining an explosive column. Background Art

[0002] Seismic exploration technology is a method of finding oil and gas reservoirs or other exploration targets by artificially stimulating seismic waves and receiving the propagation information of seismic waves on the ground or underground to identify geological structures, strata distribution, etc. In complex mountainous areas, since large vehicle-mounted controllable seismic sources are difficult to reach, explosives are a common means of artificially stimulating seismic waves.

[0003] At present, in production, the method of exciting explosives in multiple wells on a plane is usually used to excite seismic waves, that is, the single well explosive source is split into multiple wells on a plane, the same type of explosives is used in each well, and multiple wells are combined for excitation to change the spatial distribution of energy propagating outward when the source is excited, so that it is focused directly downward. In this method of exciting seismic waves, 2-5 wells need to be drilled simultaneously for each shot of seismic data, which has the problem of high cost. Summary of the invention

[0004] The present application provides a method and device for determining an explosive column to solve the problem of high cost of using explosives to excite seismic waves.

[0005] In a first aspect, the present application provides a method for determining an explosive column, which is applicable to single well seismic exploration, comprising:

[0006] Constructing a geometric model according to the depth and diameter of the target well;

[0007] Assembling geometric models based on explosion parameters of different types of explosives to form at least two sets of explosion models, wherein the explosive columns configured in each set of explosion models contain explosives of at least two types of explosives;

[0008] Perform explosion simulations on at least two sets of explosion models respectively to obtain energy distribution corresponding to each set of explosion models;

[0009] For each set of explosion models, the objective function value corresponding to the effective area is determined according to the energy distribution corresponding to the explosion model. The objective function value is used to reflect the relationship between the energy of the effective area and the energy of the ineffective area. The effective area is used to represent the effective range of the seismic wave generated by the explosive column as the source to propagate toward the target layer;

[0010] According to the objective function values ​​corresponding to at least two sets of explosion models, a target explosive column suitable for the target well is determined, and the target explosive column is the explosive column to be selected included in the at least two sets of explosion models.

[0011] Optionally, the energy distribution includes a particle pressure distribution, and determining the objective function value corresponding to the effective area according to the energy distribution corresponding to the explosion model includes:

[0012] Determine the effective area according to the seismic common shot gather data and the maximum offset information in the seismic exploration observation system; determine the particle pressure value of the particle pressure distribution corresponding to the explosion model within the effective area; obtain the objective function value corresponding to the explosion model according to the particle pressure value within the effective area.

[0013] Optionally, obtaining the objective function value corresponding to the explosion model according to the particle pressure value within the effective area includes: obtaining the objective function value corresponding to the explosion model according to the following formula:

[0014]

[0015] Wherein, taking the position of the seismic source as the pole and the ray perpendicular to the ground as the polar axis to establish a polar coordinate system, is the polar angle, r is the polar radius, α1 and α2 are the polar angles corresponding to the left and right boundaries of the effective area, r1 and r2 are the upper and lower limits of the polar radius, is the particle pressure value at the corresponding position, and F is the objective function value corresponding to the explosion model.

[0016] Optionally, determining the target explosive applicable to the target well according to the objective function values corresponding to at least two sets of explosion models includes: determining the maximum function value among the objective function values corresponding to at least two sets of explosion models; determining that the explosive column configured in the explosion model corresponding to the maximum function value is the target explosive column applicable to the target well.

[0017] Optionally, assembling geometric models based on the explosion parameters of different types of explosives to form at least two sets of explosion models includes: based on the explosion parameters of different types of explosives, screwing different types of explosives at different positions of the explosive column to obtain at least two types of explosive columns; using at least two types of explosive columns to assemble geometric models to form at least two sets of explosion models.

[0018] Optionally, constructing a geometric model according to the well depth and well diameter of the target well includes:

[0019] In a preset surrounding rock medium, construct a target well with the well depth and well diameter to obtain a geometric model.

[0020] In a second aspect, the present application provides an explosive column determination device applicable to single-well seismic exploration, including:

[0021] A construction module, configured to construct a geometric model according to the well depth and well diameter of the target well;

[0022] A formation module, configured to assemble geometric models based on the explosion parameters of different types of explosives to form at least two sets of explosion models, and each set of explosion models configured with an explosive column contains explosives of at least two types of explosives;

[0023] An acquisition module, configured to perform explosion simulations on at least two sets of explosion models respectively, and obtain the energy distribution corresponding to each set of explosion models;

[0024] A determination module, configured to, for each set of explosion models, determine the objective function value corresponding to the effective area according to the energy distribution corresponding to the explosion model, where the objective function value is used to reflect the relationship between the energy in the effective area and the energy in the ineffective area, and the effective area is used to represent the effective range in which the seismic wavefront generated with the explosive column as the seismic source propagates towards the target layer;

[0025] A processing module, configured to determine a target explosive column applicable to the target well according to the objective function values corresponding to at least two sets of explosion models, where the target explosive column is a candidate explosive column included in at least two sets of explosion models.

[0026] Optionally, the energy distribution includes the particle pressure distribution. The determination module is specifically configured to: determine the effective area according to the seismic common shot gather data and the maximum offset information in the seismic exploration observation system; determine the particle pressure values in the effective area of the particle pressure distribution corresponding to the explosion model; and obtain the objective function value corresponding to the explosion model according to the particle pressure values in the effective area.

[0027] Optionally, the determination module is specifically configured to: obtain the objective function value corresponding to the explosion model according to the following formula;

[0028]

[0029] Wherein, taking the position of the seismic source as the pole and the ray perpendicular to the ground as the polar axis to establish a polar coordinate system, is the polar angle, r is the polar radius, α1 and α2 are the polar angles corresponding to the left and right boundaries of the effective area, r1 and r2 are the upper and lower limits of the polar radius, is the particle pressure value at the corresponding position, and F is the objective function value corresponding to the explosion model.

[0030] Optionally, the processing module is specifically configured to: determine the maximum function value among the objective function values corresponding to at least two sets of explosion models; and determine the explosive column configured in the explosion model corresponding to the maximum function value as the target explosive column applicable to the target well.

[0031] Optionally, the formation module is specifically configured to: based on the explosion parameters of different types of explosives, screw different types of explosives at different positions of the explosive column to obtain at least two explosive columns; and use at least two explosive columns to assemble a geometric model to form at least two sets of explosion models.

[0032] Optionally, the construction module is specifically configured to: construct a target well with a well depth and a well diameter in a preset surrounding rock medium to obtain a geometric model.

[0033] In a third aspect, the present application provides an explosive column determination device, including: a memory and a processor;

[0034] The memory is used to store program instructions.

[0035] The processor is used to call the program instructions in the memory to execute the explosive column determination method as described in the first aspect of this application.

[0036] In a fourth aspect, this application provides a computer-readable storage medium. Computer program instructions are stored in the computer-readable storage medium. When the computer program instructions are executed, the explosive column determination method as described in the first aspect of this application is implemented.

[0037] In a fifth aspect, this application provides a computer program product, including a computer program. When the computer program is executed by a processor, the explosive column determination method as described in the first aspect of this application is implemented.

[0038] The explosive column determination method and device provided by this application construct a geometric model according to the well depth and well diameter of the target well, assemble the geometric model based on the explosion parameters of different types of explosives to form at least two sets of explosion models. Each set of explosion models configured with an explosive column contains explosives of at least two explosive types. The explosion simulation is respectively carried out on at least two sets of explosion models to obtain the energy distribution corresponding to each set of explosion models. For each set of explosion models, the objective function value corresponding to the effective area is determined according to the energy distribution corresponding to the explosion model. According to the objective function values corresponding to at least two sets of explosion models, the target explosive column applicable to the target well is determined. Compared with the existing planar multi-well combined excitation technology, the explosive column determined by this application is applicable to single-well seismic exploration. For each shot excitation, only one excitation well needs to be drilled. Therefore, the cost can be greatly saved. Moreover, the explosive column contains explosives of at least two explosive types, making the excitation energy more focused in the direction of propagation towards the target layer, capable of enhancing the excitation seismic energy, and further enhancing the signal-to-noise ratio of seismic data. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] In order to more clearly illustrate the technical solutions in the embodiments of this application or the prior art, the following will briefly introduce the drawings required for use in the description of the embodiments or the prior art. Obviously, the drawings in the following description are some embodiments of this application. For those of ordinary skill in the art, other drawings can also be obtained based on these drawings without creative efforts.

[0040] Figure 1 It is a schematic diagram of an application scenario provided by an embodiment of this application;

[0041] Figure 2 It is a flowchart of the explosive column determination method provided by an embodiment of this application;

[0042] Figure 3Schematic diagram of the measured stress-strain curve of sandstone provided by an embodiment of the present application;

[0043] Figure 4 Schematic diagram of the geometric model provided by an embodiment of the present application;

[0044] Figure 5 Schematic diagram of the explosive column provided by an embodiment of the present application;

[0045] Figure 6 Schematic diagram of the time-varying curve of the particle pressure at different positions within the range of 0 - 3 m directly below the bottom of the well provided by an embodiment of the present application;

[0046] Figure 7a Schematic diagram of the effective area provided by an embodiment of the present application;

[0047] Figure 7b Schematic diagram of the effective area provided by another embodiment of the present application;

[0048] Figure 8 Flowchart of the method for determining the explosive column provided by another embodiment of the present application;

[0049] Figure 9a Schematic diagram of the spatial distribution of the particle pressure when different types of explosive combinations at a total explosive excitation mass of 18 kg propagate out of the plastic zone at the tail of the vibration provided by an embodiment of the present application;

[0050] Figure 9b Schematic diagram of the spatial distribution of the particle pressure when different types of explosive combinations at a total explosive excitation mass of 18 kg propagate out of the plastic zone at the tail of the vibration provided by another embodiment of the present application;

[0051] Figure 9c Schematic diagram of the spatial distribution of the particle pressure when different types of explosive combinations at a total explosive excitation mass of 18 kg propagate out of the plastic zone at the tail of the vibration provided by another embodiment of the present application;

[0052] Figure 9d Schematic diagram of the spatial distribution of the particle pressure when different types of explosive combinations at a total explosive excitation mass of 18 kg propagate out of the plastic zone at the tail of the vibration provided by another embodiment of the present application;

[0053] Figure 9e Schematic diagram of the spatial distribution of the particle pressure when different types of explosive combinations at a total explosive excitation mass of 18 kg propagate out of the plastic zone at the tail of the vibration provided by another embodiment of the present application;

[0054] Figure 9f Schematic diagram of the spatial distribution of the particle pressure when different types of explosive combinations at a total explosive excitation mass of 18 kg propagate out of the plastic zone at the tail of the vibration provided by another embodiment of the present application;

[0055] Figure 10 Schematic diagram of the particle pressure distribution curve at a certain moment after the explosive column provided in an embodiment of the present application is excited;

[0056] Figure 11 Schematic diagram of the objective function values corresponding to six sets of explosion models provided in an embodiment of the present application;

[0057] Figure 12 Schematic diagram of the seismic geological model provided in an embodiment of the present application;

[0058] Figure 13a Forward modeling of the acoustic wave equation provided in an embodiment of the present application Figure 12 Schematic diagram of the single-shot record of the model;

[0059] Figure 13b Forward modeling of the acoustic wave equation provided in another embodiment of the present application Figure 12 Schematic diagram of the single-shot record of the model;

[0060] Figure 14a Schematic diagram of the field test seismic record provided in an embodiment of the present application;

[0061] Figure 14b Schematic diagram of the field test seismic record provided in another embodiment of the present application;

[0062] Figure 15 Schematic diagram of the structure of the explosive column determination device provided in an embodiment of the present application;

[0063] Figure 16 Schematic diagram of the structure of the explosive column determination device provided in another embodiment of the present application;

[0064] Figure 17 Schematic diagram of the structure of the explosive column determination device provided in another embodiment of the present application. Detailed implementation manners

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

[0066] Seismic exploration technology is a method of artificially generating seismic waves and receiving the propagation information of seismic waves on the ground or downhole to identify geological structures, stratigraphic distributions, etc., and to search for oil and gas reservoirs or other exploration targets. In complex mountainous terrain areas, since large vehicle-mounted vibrators cannot reach, explosives are a common means of artificially generating seismic waves. Usually, one or several boreholes are drilled, explosive columns with the same drug type are placed at the bottom of the well, and the construction process of detonating with electronic detonators is used. In such areas, due to the complex surface and underground geological structures, the reflected signal energy of the data obtained by single-well excitation is weak, and the signal-to-noise ratio is extremely low. To obtain greater reflected signal energy, usually the method of increasing the amount of explosives is chosen. Since the diameter of the explosive column is fixed and the density is fixed, the increase in the amount of explosives can only be achieved by increasing the length of the explosive column. Taking 30 kg of explosives as an example, the length of the explosive column can reach 7.5 m, while the diameter is only 8 cm. The slender shape causes the energy to preferentially propagate downward to the side after the detonator detonates, rather than downward as expected in seismic exploration. Changing the spatial distribution pattern of the energy propagating outward during the source excitation to make it focus downward can only be achieved through a combination method.

[0067] Currently, in production, usually the method of combining and exciting explosives with multiple wells in a plane is adopted to generate seismic waves, that is, the single-well explosive source is split into multiple wells in the plane, and the same type of explosives is used in each well, and multiple wells are combined and excited to change the spatial distribution pattern of the energy propagating outward during the source excitation to make it focus downward. In this method of generating seismic waves, to obtain one shot of seismic data, 2 - 5 wells need to be drilled simultaneously, which has the problem of high cost.

[0068] Based on the above, the present application provides a method and device for determining an explosive column. By installing different types of explosives in segments in a single well and exciting them simultaneously, the excitation energy is more focused in the direction of propagation towards the target layer, and the seismic excitation energy can be improved without increasing the field construction cost, thereby improving the signal-to-noise ratio of seismic data. The specific implementation process of determining the explosive column with different types of explosives in a single well can refer to the solutions of the following embodiments.

[0069] Figure 1 The following is a schematic diagram of an application scenario provided by an embodiment of the present application. As Figure 1 shown, in this application scenario, the explosive column 120 is placed at the bottom of the well 110, detonated at the top with an electronic detonator, the generated seismic waves are reflected back to the surface by the underground layer interface 150, and the geophone 130 receives the seismic wave signal reflected back to the surface and transmits the received seismic wave signal to the computer device 140 to obtain seismic data.

[0070] It should be noted that Figure 1 is only a schematic diagram of an application scenario provided by an embodiment of the present application. The embodiments of the present application do not limit Figure 1 the equipment included in Figure 1Limit the positional relationship between devices. For example, in the Figure 1 shown application scenario, it may further include a data storage device, which may be an external memory relative to the computer device 120 or an internal memory integrated in the computer device 120.

[0071] Figure 2 The flowchart of the explosive column determination method provided by an embodiment of the present application. The method of this embodiment can be applied to an electronic device, which can be a terminal device, a server, or a server cluster, etc. The terminal device can be, for example, a mobile phone, a tablet computer, a laptop computer, a desktop computer, etc. As Figure 2 shown, the method of this embodiment includes:

[0072] S201. Construct a geometric model according to the well depth and well diameter of the target well.

[0073] In this embodiment, the well depth and well diameter of the target well can be input by the user to the electronic device executing the method embodiment of the present application, or sent by other devices to the electronic device executing the method embodiment of the present application. The present application embodiment does not limit its acquisition method. Exemplarily, the well depth of the target well is 30m and the well diameter is 10cm. After obtaining the well depth and well diameter of the target well, a geometric model can be constructed.

[0074] Optionally, this step may include: constructing a target well with a well depth and a well diameter in a preset surrounding rock medium to obtain a geometric model. It can be understood that the settings of the well depth and / or well diameter are related to the parameters of the preset surrounding rock medium, and different preset surrounding rock media correspond to different well depths and / or well diameters. Further, the parameters of the preset surrounding rock medium may include, but are not limited to: the density, elastic modulus, tangent modulus, Poisson's ratio, yield stress information, etc. of the surrounding rock. Exemplarily, the preset surrounding rock medium is, for example, sandstone. Table 1 shows the measured sandstone parameter data provided by an embodiment of the present application.

[0075] Table 1 Measured sandstone parameter data

[0076]

[0077]

[0078] Figure 3 The schematic diagram of the measured stress-strain curve of sandstone provided by an embodiment of the present application. Figure 3 Quantitatively describes the process in which the tested sandstone undergoes elastic deformation -> plastic deformation -> fracture under axial compression. Among them, the ordinate is usually represented by differential stress, which is the pressure difference between the axial and radial directions. The mechanical properties of the surrounding rock are described by a piecewise linear plastic model, that is, the measured stress-strain curve is piecewise fitted (as Figure 3 shown).

[0079] After obtaining the well depth and well diameter of the target well, a target well with the well depth and well diameter is constructed in a preset surrounding rock medium to obtain a geometric model. Exemplarily, Figure 4 As shown in the schematic diagram of the geometric model provided by an embodiment of the present application, Figure 4 as shown, a cylindrical well with a well depth of 30 m and a well diameter of 10 cm is constructed at the center of a cylindrical surrounding rock medium with a diameter of 60 m and a height of 60 m, that is, the geometric model. Among them, an explosive column can be placed at the bottom of the well, and the top of the explosive column is filled with surrounding rock medium.

[0080] S202. Assemble the geometric model based on the explosion parameters of different types of explosives to form at least two sets of explosion models, and each set of explosion models configured with an explosive column contains explosives of at least two types of explosives.

[0081] Exemplarily, the explosion parameters may include the charge density, detonation velocity, detonation pressure, and initial internal energy parameters of different types of explosives. Table 2 shows the explosion parameters corresponding to five types of explosives provided by an embodiment of the present application. As shown in Table 2, different types of explosives follow the following equation of state to convert chemical energy into mechanical energy:

[0082]

[0083] where p is the particle pressure, E is the initial internal energy of the detonation products per unit volume, V is the ratio of the detonation product volume to the initial volume, V usually takes the value of 1, A, B, R1, R2, and ω are undetermined constants, usually obtained from real explosion tests.

[0084] Table 2 Explosion parameters corresponding to five types of explosives

[0085]

[0086]

[0087] After constructing the geometric model, at least two sets of explosion models are formed by assembling the geometric model based on the explosion parameters of different types of explosives, and each set of explosion models configured with an explosive column contains explosives of at least two types of explosives. For example, an explosive column configured in a set of explosion models contains explosives of two types of explosives, namely: TNT explosive and high-density ammonium nitrate explosive.

[0088] It should be added that the explosion parameters of the above different types of explosives can be input by the user to the electronic device executing the embodiment of the method, or sent by other devices to the electronic device executing the embodiment of the method, which is not limited in the embodiments of the present application.

[0089] In some embodiments, this step may include: screwing different types of explosives at different positions of the explosive column based on the explosion parameters of different types of explosives to obtain at least two explosive columns; using at least two explosive columns to assemble a geometric model to form at least two sets of explosion models.

[0090] Exemplarily, Figure 5 is a schematic diagram of an explosive column provided in an embodiment of the present application. As Figure 5 shown, the explosive column is screwed together by explosives of different explosive types, and the explosives of different explosive types respectively correspond to Figure 5 a small section of the explosive column in. The length of each section can be set to 0.5 m, and the mass can be set to 2 kg. The specific number of sections used depends on the total mass of explosives in the production design. Optionally, the same or different types of explosives can be set for each section of the explosive column. If each section is the same, it degenerates into the same type of explosive used in production.

[0091] Among them, using at least two explosive columns to assemble the geometric model as Figure 4 shown, that is, placing each type of explosive column at the position for placing the explosive column in the geometric model as Figure 4 shown to form at least two sets of explosion models.

[0092] S203. Perform explosion simulations on at least two sets of explosion models respectively to obtain the energy distribution corresponding to each set of explosion models.

[0093] In this embodiment, after obtaining at least two sets of explosion models, among them, the surrounding rock parameters of each set of explosion models are, for example, the sandstone parameters shown in Table 1, and the explosion parameters of different types of explosives in each set of explosion models are, for example, the explosion parameters corresponding to the five types of explosives shown in Table 2. Therefore, perform explosion simulations on at least two sets of explosion models respectively to obtain the energy distribution corresponding to each set of explosion models.

[0094] When performing the explosion simulation, from the explosion center outwards, in a short period of time, the surrounding rock will experience the processes of fragmentation, plastic deformation, and elastic deformation. Using whether the particle can return to the equilibrium position after jumping as the judgment criterion, if it can return to the equilibrium position, it indicates that it has entered the elastic region, otherwise it is the plastic region. Exemplarily, Figure 6 is a schematic diagram of the time-varying curve of the particle pressure at different positions within 0 - 3 m directly below the bottom of the well provided in an embodiment of the present application. Figure 6 Shows the time-varying curves of the particle pressure at positions 0 m, 0.5 m, 1 m, 1.5 m, 2 m, 2.5 m, and 3 m directly below the bottom of the well. As Figure 6 shown, before 3 m, the particle still maintains a certain amplitude after forced vibration and cannot return to the equilibrium state. Based on this, the range of the plastic region can be judged: the plastic region is approximately more than 3 m away from the explosion center. When performing energy statistics, it is necessary to ensure that the wave tail of the particle vibration spreads more than 3 m away.

[0095] S204. For each set of explosion models, determine the objective function value corresponding to the effective area according to the energy distribution corresponding to the explosion model.

[0096] Among them, the objective function value is used to reflect the relationship between the energy in the effective area and the energy in the ineffective area. The effective area is used to represent the effective range in which the seismic wavefront generated with the explosive column as the seismic source propagates towards the target layer. As an alternative, the objective function value can be the ratio of the energy in the effective area to the energy in the ineffective area.

[0097] For the effective area, Figure 7a is a schematic diagram of the effective area provided by an embodiment of the present application, Figure 7b is a schematic diagram of the effective area provided by another embodiment of the present application. As Figure 7a shown, the seismic source is excited at point O, corresponding to the signal receiving length of 0 - 1200m on the surface. Only the reflection signals in the BC section on the target layer can be received, and the areas to the right of point B and to the left of point C cannot be received. The reflections at points B and C correspond to Figure 7b the two farthest positions at both ends of the reflection wave of the target layer in Figure 7a . Only the energy propagating within an angle less than α is useful, and this interval is the effective area (in

[0098] S205. Determine the target explosive column applicable to the target well according to the objective function values corresponding to at least two sets of explosion models.

[0099] Among them, the target explosive column is the candidate explosive column included in at least two sets of explosion models.

[0100] It can be understood that the determination of the target explosive column is related to the specific meaning of the objective function value. Exemplarily, when a larger objective function value indicates that the energy of the seismic wave excited in the effective area is greater, the explosive column in the explosion model corresponding to the largest objective function value is used as the target explosive column; or, when a smaller objective function value indicates that the energy of the seismic wave excited in the effective area is greater, the explosive column in the explosion model corresponding to the largest objective function value is used as the target explosive column, and so on. After determining the target explosive column applicable to the target well, place the target explosive column at the bottom of the target well and excite the target explosive column to obtain seismic data.

[0101] The method for determining the explosive column provided in this embodiment constructs a geometric model according to the well depth and well diameter of the target well, assembles the geometric model based on the explosion parameters of different types of explosives to form at least two sets of explosion models. Each explosive column configured in each set of explosion models contains explosives of at least two types of explosives. Respectively perform explosion simulations on at least two sets of explosion models to obtain the energy distribution corresponding to each set of explosion models. For each set of explosion models, determine the objective function value corresponding to the effective area according to the energy distribution corresponding to the explosion model. According to the objective function values corresponding to at least two sets of explosion models, determine the target explosive column applicable to the target well. Compared with the existing planar multi-well combined excitation technology, the explosive column determined through the embodiments of the present application is applicable to single-well seismic exploration. For each shot excitation, only one excitation well needs to be drilled. Therefore, the cost can be greatly saved. Moreover, the explosive column contains explosives of at least two types of explosives, making the excitation energy more focused in the direction of propagation towards the target layer, capable of enhancing the excitation seismic energy, and further enhancing the signal-to-noise ratio of seismic data.

[0102] In the above embodiment, the energy distribution is mentioned. Among them, the energy distribution has various forms. Next, taking the particle pressure distribution as an example to characterize the energy distribution for relevant description. Figure 8 It is a flowchart of the method for determining the explosive column provided in another embodiment of the present application. As Figure 8 shown, the method of the embodiment of the present application may include:

[0103] S801. Construct a geometric model according to the well depth and well diameter of the target well.

[0104] S802. Assemble the geometric model based on the explosion parameters of different types of explosives to form at least two sets of explosion models. Each explosive column configured in each set of explosion models contains explosives of at least two types of explosives.

[0105] In this embodiment, the specific implementation processes of S801 and S802 can refer to the relevant descriptions of the Figure 2 shown embodiment, and will not be elaborated here.

[0106] S803. Respectively perform explosion simulations on at least two sets of explosion models to obtain the particle pressure distribution corresponding to each set of explosion models.

[0107] When the wave tail of the vibration in each set of explosion models propagates out of the plastic zone, record the distribution pattern of the particle pressure at this moment. For example, use a total explosive excitation mass of 18 kg to perform explosion simulations on different explosion models, and record the spatial distribution of the particle pressure when the wave tail of the vibration in each set of explosion models propagates out of the plastic zone; the total explosive amount of 18 kg can be disassembled into 9 segments of 2 kg. Figures 9a to 9f Respectively show schematic diagrams of the spatial distribution of the particle pressure of different types of explosive combinations when the wave tail of the vibration propagates out of the plastic zone under a total explosive excitation mass of 18 kg. Among them, Figure 9aThe different types of explosive combinations are 3 small segments of TNT explosive in the upper part - 6 small segments of high-density ammonium nitrate explosive in the lower part. Figure 9b The different types of explosive combinations are 3 small segments of high-density ammonium nitrate explosive in the upper part - 6 small segments of TNT explosive in the lower part. Figure 9c The different types of explosive combinations are 3 small segments of ammonium nitrate and TNT explosive in the upper part - 6 small segments of high-density ammonium nitrate explosive in the lower part. Figure 9d The different types of explosive combinations are high-density ammonium nitrate explosive in the upper segment - 6 small segments of ammonium nitrate and TNT explosive in the lower part. Figure 9e The different types of explosive combinations are 3 small segments of TNT explosive in the upper part - 6 small segments of ammonium nitrate and TNT explosive in the lower part. Figure 9f The different types of explosive combinations are 3 small segments of ammonium nitrate and TNT explosive in the upper part - 6 small segments of TNT explosive in the lower part. At the same time, referring to Figures 9a to 9f , it can be seen that the energy spatial distributions corresponding to different types of explosive ratios vary greatly, but generally they all have the characteristic of focusing downward.

[0108] After that, the following processing is carried out for each set of explosion models:

[0109] S804. Determine the effective area according to the seismic common shot gather data and the maximum offset information in the seismic exploration observation system.

[0110] Still referring to Figure 7a and Figure 7b , the angle α can be obtained through the following formula 1:

[0111]

[0112] where x is the maximum offset length, and h A is the depth of the target formation. Among them, for the depth h of the target formation A , the stacking velocity v of the target formation can be obtained first through the following formula 2, and then the depth h of the target formation can be obtained according to the stacking velocity v through the following formula 3 A . That is, extract the self-excitation and self-reception point time t A , the reflection time t of the target formation at the maximum offset B from the seismic common shot gather data, obtain the maximum offset information x from the seismic exploration observation system parameters, calculate the stacking velocity v of the target formation, and multiply the stacking velocity by half of the self-excitation and self-reception time to obtain the estimated depth h of the target formation A .

[0113]

[0114] where t A is the self-excitation and self-reception point time extracted from the seismic common shot gather data, t B is the reflection time of the target formation at the maximum offset, x is the maximum offset length obtained from the seismic exploration observation system parameters, and v is the stacking velocity v of the target formation.

[0115] h A = v * t A / 2 Formula Three

[0116] Based on Figure 7a and Figure 7b above, the stacking velocity v can be obtained as 1268 m / s through the above Formula Two, and then the target layer depth h A is 409 m. Comparing Figure 7a and Figure 7b the estimated target layer depth of 409 m and the actual target layer depth of 400 m, the layer depths of the two are close.

[0117] Therefore, after pre-collecting the seismic common shot gather data and the maximum offset information in the seismic exploration observation system information, the effective area can be determined to serve the subsequent calculation of the objective function. Exemplarily, for the parameters corresponding to this calculation example, the incident angle α corresponding to the effective area is about 36.9°, and the integration upper and lower limits α1 = -36.9° and α2 = 36.9° corresponding to the objective function calculation formula Four.

[0118] S805. For each set of explosion models, determine the particle pressure values of the particle pressure distribution corresponding to the explosion model within the effective area.

[0119] After obtaining the particle pressure distribution corresponding to each set of explosion models, the particle pressure values within the effective area can be extracted according to the particle pressure distribution.

[0120] Optionally, use red to mark the particle pressure distribution within the effective area and extract the particle pressure values of all particles in the red area.

[0121] S806. Obtain the objective function value corresponding to the explosion model according to the particle pressure values within the effective area.

[0122] Optionally, obtain the objective function value corresponding to the explosion model according to the following Formula Four.

[0123]

[0124] where, taking the position of the seismic source as the pole and the ray perpendicular to the ground as the polar axis to establish a polar coordinate system, is the polar angle, r is the polar radius, α1 and α2 are the polar angles corresponding to the left and right boundaries of the effective area, r1 and r2 are the upper and lower limits of the polar radius, is the particle pressure value at the corresponding position, and F is the objective function value corresponding to the explosion model.

[0125] The numerator in Equation 4 represents the cumulative energy value of the effective signal interval, and the denominator represents the energy value after subtracting the effective signal interval from the total energy in the full polar angle space within the polar radius range from r1 to r2. Therefore, this objective function reflects the relationship between the energy in the effective area and the energy in the ineffective area, that is, the ratio of the energy in the effective area to the energy in the ineffective area. Figure 10 This is a schematic diagram of the particle pressure distribution curve at a certain moment after the explosive column is excited provided by an embodiment of the present application, as Figure 10 shown. The vibration wavefront corresponds to the r2 polar radius length, the vibration wave tail corresponds to the r1 polar radius length, and the energy statistics of the objective function are carried out within the polar radius interval enclosed by r1 and r2 in Figure 10 . The energy in the effective area corresponds to the Figure 10 shaded part in, and it is necessary to ensure that r1 is greater than the plastic zone range.

[0126] Among them, S804 to S806 are further explanations of the Figure 2 step S804 in.

[0127] The larger the objective function value obtained through Equation 2, the higher the matching degree of the corresponding explosive column and the target well, that is, the better the effect of seismic exploration of the target well using this explosive column. Therefore, select this explosive column as the target explosive column.

[0128] S807. Determine the maximum function value among the objective function values corresponding to at least two sets of explosion models.

[0129] Exemplarily, the objective function values corresponding to three sets of explosion models are obtained. Among them, the objective function value corresponding to the first set of explosion models is 0.4, the objective function value corresponding to the second set of explosion models is 0.6, and the objective function value corresponding to the third set of explosion models is 0.5. Then, determine that the maximum function value among the objective function values corresponding to the three sets of explosion models is 0.6.

[0130] S808. Determine that the explosive column configured in the explosion model corresponding to the maximum function value is the target explosive column suitable for the target well.

[0131] Exemplarily, Figure 11 shows the objective function values corresponding to six sets of explosion models in Figures 9a to 9f respectively. As Figure 11 shown, when the upper 3 small sections of unit explosives are ammonium nitrate explosives, the objective function values are generally larger. Combining Figure 9c and Figure 9f , when the upper 3 small half-sections of explosives are ammonium nitrate explosives, the combined excitation has a better suppression effect on the upward energy, and more energy is distributed in the downward direction. Therefore, the objective function value is high. Therefore, determine that the combination of the upper 3 small sections of ammonium nitrate explosives and the lower 6 small sections of high-density ammonium nitrate explosives with the largest objective function value is the target explosive column suitable for the target well.

[0132] The explosion simulation in the above embodiments calculates the vibration state of particles within a limited space range (generally 20m to 50m) near the explosion point. To verify the energy change in the actual seismic wave propagation scale range, a Figure 12 seismic geological model is designed. This seismic geological model consists of three layers of horizontal layered media, and the layer velocity values are marked in the figure. A layer of high-speed anomaly points is evenly arranged near the surface to simulate the surface scattering source. Figure 13a It is the schematic diagram of the single-shot record of the Figure 12 model obtained by the forward modeling of the acoustic wave equation provided by an embodiment of the present application. Figure 13b It is the schematic diagram of the single-shot record of the Figure 12 model obtained by the forward modeling of the acoustic wave equation provided by another embodiment of the present application. In Figure 13a and Figure 13b , the initial wavelet adopts the stable particle vibration form in the elastic region of the explosion numerical simulation. Figure 13a is the simulation result excited by a single-shaped high-density ammonium nitrate explosive for the whole section. Figure 13b is the simulation result excited by the combination of 3 small sections of ammonium-tnt explosive in the upper part and 6 small sections of high-density ammonium nitrate explosive in the lower part. It can be seen that compared with the single-shaped explosive, the diffracted noise generated by the shallow scattering points excited by the combined explosive type is greatly reduced.

[0133] Figure 14a It is the schematic diagram of the field test seismic record provided by an embodiment of the present application. Figure 14b They are respectively the schematic diagrams of the field test seismic records provided by another embodiment of the present application. Among them, Figure 14a is the seismic record excited by 18Kg of a single long charge of high-density ammonium nitrate explosive. Figure 14b is the single-shot record obtained by splitting the surface adjacent position into two sections of 6kg:12kg (3 small sections in the upper part and 6 small sections in the lower part), with ammonium-tnt in the upper section and high-density ammonium nitrate in the lower section. Referring to Figure 14a and Figure 14b , it can be seen that: compared with the energy of the effective reflection of the single-shaped explosive column, the energy of the effective reflection of the explosive column determined by the embodiments of the present application is stronger, and the overall signal-to-noise ratio of the record is improved.

[0134] In summary, the method and device for determining the explosive column provided in this application install different types of explosives in sections in a single well and simultaneously initiate them. When different types of explosives act on the same surrounding rock, the generated wavelet waveforms and energies are different. This application utilizes the reasonable combination of these differences in space to make the excitation energy more focused in the direction of propagation towards the target layer, which is more conducive to improving the signal-to-noise ratio of the data. Compared with the existing planar multi-well combined excitation technology, only one excitation well needs to be drilled for each shot excitation in this application, and the drilling cost is low; compared with the conventional single-well single explosive type, during the construction process, only different explosive types are screwed at different positions of the explosive column, without adding additional field construction processes; compared with the existing single-well multi-stage multi-charge delayed excitation technology, there is no need for delayed excitation, reducing the control error of the delay time. The method for determining the explosive column provided in this application considers the mechanical action process of multiple different types of explosives and the excitation medium. Only one geometric modeling needs to be completed, and then the explosive type parameters of each small section of the explosive column are continuously changed, and the explosion simulation calculation is repeated. Then, the required statistical information is extracted for calculation and comparison, and the best ratio relationship of different explosive types can be selected according to the size of the objective function value. The manual interaction is less, and the calculation is convenient and fast.

[0135] The following is an embodiment of the device of this application, which can be used to execute the method embodiment of this application. For the details not disclosed in the device embodiment of this application, please refer to the method embodiment of this application.

[0136] Figure 15 The following is a schematic structural diagram of the device for determining the explosive column provided in an embodiment of this application, as Figure 15 shown, the device 1500 for determining the explosive column in this embodiment includes: a construction module 1501, a formation module 1502, an acquisition module 1503, a determination module 1504, and a processing module 1505. Among them:

[0137] The construction module 1501 is used to construct a geometric model according to the well depth and well diameter of the target well.

[0138] The formation module 1502 is used to assemble the geometric model based on the explosion parameters of different types of explosives to form at least two sets of explosion models, and each set of explosion models configured with an explosive column contains at least two types of explosives.

[0139] The acquisition module 1503 is used to perform explosion simulations on at least two sets of explosion models respectively to obtain the energy distribution corresponding to each set of explosion models.

[0140] The determination module 1504 is used to, for each set of explosion models, determine the objective function value corresponding to the effective area according to the energy distribution corresponding to the explosion model. The objective function value is used to reflect the relationship between the energy in the effective area and the energy in the ineffective area, and the effective area is used to represent the effective range of the seismic wave generated by the explosive column as the seismic source propagating towards the target layer.

[0141] A processing module 1505, configured to determine a target explosive column applicable to a target well according to objective function values corresponding to at least two sets of explosion models, where the target explosive column is a candidate explosive column included in at least two sets of explosion models.

[0142] Based on any of the above-described embodiments, the energy distribution may specifically be a particle pressure distribution. In this case, the determination module 1504 may specifically be configured to: determine an effective area according to seismic common shot gather data and maximum offset information in a seismic exploration observation system; determine particle pressure values of the particle pressure distribution corresponding to the explosion model within the effective area; and obtain an objective function value corresponding to the explosion model according to the particle pressure values within the effective area.

[0143] Based on any of the above-described embodiments, the determination module 1504 may specifically be configured to: obtain an objective function value corresponding to the explosion model according to the following formula:

[0144]

[0145] where a polar coordinate system is established with the position of the seismic source as the pole and the ray perpendicular to the ground as the polar axis. is the polar angle, r is the polar radius, α1 and α2 are polar angles corresponding to the left and right boundaries of the effective area, r1 and r2 are the upper and lower limits of the polar radius. is the particle pressure value at the corresponding position, and F is the objective function value corresponding to the explosion model.

[0146] Based on any of the above-described embodiments, the processing module 1505 may specifically be configured to: determine the maximum function value among the objective function values corresponding to at least two sets of explosion models; and determine that the explosive column configured in the explosion model corresponding to the maximum function value is the target explosive column applicable to the target well.

[0147] Based on any of the above-described embodiments, the forming module 1502 may specifically be configured to: screw different types of explosives at different positions of the explosive column based on explosion parameters of different types of explosives to obtain at least two explosive columns; and use the at least two explosive columns to assemble a geometric model to form at least two sets of explosion models.

[0148] Based on any of the above-described embodiments, the construction module 1501 may specifically be configured to: construct a target well with a well depth and a well diameter in a preset surrounding rock medium to obtain a geometric model.

[0149] The device in this embodiment may be used to execute the technical solutions of any of the above-described method embodiments. The implementation principles and technical effects are similar and will not be elaborated herein.

[0150] Figure 16 This is a schematic structural diagram of an explosive column determination device provided in another embodiment of the present application. As Figure 16As shown, the explosive column determination device 1600 of this embodiment includes: a memory 1601 and a processor 1602. Among them, the memory 1601 and the processor 1602 are connected through a bus.

[0151] The memory 1601 is used to store program instructions.

[0152] The processor 1602 is used to call the program instructions in the memory to execute:

[0153] Construct a geometric model according to the well depth and well diameter of the target well; assemble the geometric model based on the explosion parameters of different types of explosives to form at least two sets of explosion models, and each set of explosion models configured with an explosive column contains at least two types of explosives; perform explosion simulations on at least two sets of explosion models respectively to obtain the energy distribution corresponding to each set of explosion models; for each set of explosion models, determine the objective function value corresponding to the effective area according to the energy distribution corresponding to the explosion model, and the objective function value is used to reflect the relationship between the energy in the effective area and the energy in the ineffective area, and the effective area is used to represent the effective range where the seismic wave generated by the explosive column propagates towards the target layer; determine the target explosive column applicable to the target well according to the objective function values corresponding to at least two sets of explosion models, and the target explosive column is the candidate explosive column included in at least two sets of explosion models.

[0154] Based on any of the above - shown embodiments, the energy distribution includes the particle pressure distribution, and the processor 1602 can specifically be used to: determine the effective area according to the seismic common - shot gather data and the maximum offset information in the seismic exploration observation system; determine the particle pressure value of the particle pressure distribution corresponding to the explosion model within the effective area; obtain the objective function value corresponding to the explosion model according to the particle pressure value within the effective area.

[0155] Based on any of the above - shown embodiments, the processor 1602 can specifically be used to: obtain the objective function value corresponding to the explosion model according to the following formula;

[0156]

[0157] Among them, taking the position of the seismic source as the pole and the ray perpendicular to the ground as the polar axis to establish a polar coordinate system, is the polar angle, r is the polar radius, α1 and α2 are the polar angles corresponding to the left and right boundaries of the effective area, r1 and r2 are the upper and lower limits of the polar radius, is the particle pressure value at the corresponding position, and F is the objective function value corresponding to the explosion model.

[0158] Based on any of the above - shown embodiments, the processor 1602 can specifically be used to: determine the maximum function value among the objective function values corresponding to at least two sets of explosion models; determine the explosive column configured in the explosion model corresponding to the maximum function value as the target explosive column applicable to the target well.

[0159] Based on any of the above - shown embodiments, the processor 1602 may specifically be configured to: based on the explosion parameters of different types of explosives, screw - connect different types of explosives at different positions of the explosive column to obtain at least two explosive columns; use at least two explosive columns to assemble a geometric model to form at least two sets of explosion models.

[0160] Based on any of the above - shown embodiments, the processor 1602 may specifically be configured to: in a preset surrounding rock medium, construct a target well with a well depth and a well diameter to obtain a geometric model.

[0161] The device of this embodiment can be used to execute the technical solutions of any of the above - shown method embodiments. The implementation principles and technical effects are similar and will not be elaborated here.

[0162] Figure 17 The following is a schematic structural diagram of an explosive - column determination device provided in another embodiment of the present application. Exemplarily, the explosive - column determination device may be provided as a server or a computer. Refer to Figure 17 , the explosive - column determination device 1700 includes a processing component 1701, which further includes one or more processors, and memory resources represented by a memory 1702 for storing instructions executable by the processing component 1701, such as application programs. The application programs stored in the memory 1702 may include one or more modules each corresponding to a set of instructions. In addition, the processing component 1701 is configured to execute instructions to perform any of the above - mentioned method embodiments.

[0163] The explosive - column determination device 1700 may further include a power supply component 1703 configured to perform power management of the explosive - column determination device 1700, a wired or wireless network interface 1704 configured to connect the explosive - column determination device 1700 to a network, and an input / output (I / O) interface 1705. The explosive - column determination device 1700 may operate based on an operating system stored in the memory 1702, such as Windows ServerTM, Mac OS XTM, UnixTM, LinuxTM, FreeBSDTM or the like.

[0164] The present application also provides a computer program product, including a computer program, which when executed by a processor, implements the solution of the above - mentioned explosive - column determination method.

[0165] The present application also provides a computer - readable storage medium, in which computer - executable instructions are stored. When a processor executes the computer - executable instructions, the above - mentioned explosive - column determination method is implemented.

[0166] The above-mentioned computer-readable storage medium may be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic memory, flash memory, a magnetic disk, or an optical disk. The readable storage medium may be any available medium accessible by a general-purpose or special-purpose computer.

[0167] An exemplary readable storage medium is coupled to a processor, enabling the processor to read information from and write information to the readable storage medium. Of course, the readable storage medium may also be a component of the processor. The processor and the readable storage medium may be located in an application specific integrated circuit (ASIC). Of course, the processor and the readable storage medium may also exist as discrete components in the explosive column determination device.

[0168] Those of ordinary skill in the art can understand that all or part of the steps for implementing the above method embodiments may be completed by hardware related to program instructions. The foregoing program may be stored in a computer-readable storage medium. When the program is executed, it performs the steps including those of the above method embodiments; and the foregoing storage medium includes various media that can store program codes, such as ROM, RAM, magnetic disks, or optical disks.

[0169] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present application and are not intended to limit them; although the present application has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions described in the foregoing embodiments, or perform equivalent replacements on some or all of the technical features; and these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present application.

Claims

1. A method for determining an explosive column, characterized in that, Applicable to single-well seismic exploration, including: Construct a geometric model according to the well depth and well diameter of the target well; Assemble the geometric model based on the explosion parameters of different types of explosives to form at least two sets of explosion models, and each set of explosives columns configured in the explosion models contains explosives of at least two types of explosives; Conduct explosion simulations on the at least two sets of explosion models respectively to obtain the energy distribution corresponding to each set of explosion models; For each set of explosion models, determine the objective function value corresponding to the effective area according to the energy distribution corresponding to the explosion model. The objective function value is used to reflect the relationship between the energy in the effective area and the energy in the ineffective area. The effective area is used to represent the effective range where the seismic wave generated by the explosives column propagates towards the target layer; Determine the target explosives column applicable to the target well according to the objective function values corresponding to the at least two sets of explosion models. The target explosives column is the candidate explosives column included in the at least two sets of explosion models; The energy distribution includes the particle pressure distribution. The determining the objective function value corresponding to the effective area according to the energy distribution corresponding to the explosion model includes: Determine the effective area according to the seismic common shot gather data and the maximum offset information in the seismic exploration observation system; Determine the particle pressure values of the particle pressure distribution corresponding to the explosion model within the effective area; Obtain the objective function value corresponding to the explosion model according to the particle pressure values within the effective area; The obtaining the objective function value corresponding to the explosion model according to the particle pressure values within the effective area includes: Obtain the objective function value corresponding to the explosion model according to the following formula: Wherein, taking the position of the seismic source as the pole and the ray perpendicular to the ground as the polar axis, a polar coordinate system is established. is the polar angle. is the polar radius. and are the polar angles corresponding to the left and right boundaries of the effective area. and are the upper and lower limits of the polar radius. is the particle pressure value at the corresponding position. is the objective function value corresponding to the explosion model.

2. The method according to claim 1, wherein The determining the target explosives column applicable to the target well according to the objective function values corresponding to the at least two sets of explosion models includes: Determine the maximum function value among the objective function values corresponding to the at least two sets of explosion models; Determine that the explosives column configured in the explosion model corresponding to the maximum function value is the target explosives column applicable to the target well.

3. The method according to claim 1 or 2, characterized in that, The assembling the geometric model based on the explosion parameters of different types of explosives to form at least two sets of explosion models includes: Based on the explosion parameters of different types of explosives, screw different types of explosives at different positions of the explosives column to obtain at least two types of explosives columns; Use the at least two types of explosives columns to assemble the geometric model to form at least two sets of explosion models.

4. The method according to claim 1 or 2, characterized in that, The constructing a geometric model according to the well depth and well diameter of the target well includes: In a preset surrounding rock medium, construct a target well with the well depth and the well diameter to obtain the geometric model.

5. An explosive column determination device, characterized in that, Applicable to single-well seismic exploration, including: A construction module for constructing a geometric model according to the well depth and well diameter of the target well; A forming module for assembling the geometric model based on the explosion parameters of different types of explosives to form at least two sets of explosion models, and each set of explosives columns configured in the explosion models contains explosives of at least two types of explosives; An acquisition module for conducting explosion simulations on the at least two sets of explosion models respectively to obtain the energy distribution corresponding to each set of explosion models; A determination module, configured to, for each set of explosion models, determine a target function value corresponding to the effective area according to the energy distribution corresponding to the explosion model, where the target function value is used to reflect the relationship between the energy in the effective area and the energy in the ineffective area, and the effective area is used to represent the effective range in which the seismic wave generated with the explosive column as the seismic source propagates towards the target layer; A processing module, configured to determine a target explosive column applicable to the target well according to the target function values corresponding to the at least two sets of explosion models, where the target explosive column is a candidate explosive column included in the at least two sets of explosion models; The energy distribution includes a particle pressure distribution, and the determination module is specifically configured to: Determine the effective area according to the seismic common shot gather data and the maximum offset information in the seismic exploration observation system; Determine the particle pressure values of the particle pressure distribution corresponding to the explosion model in the effective area; Obtain the target function value corresponding to the explosion model according to the particle pressure values in the effective area; The determination module is specifically configured to: Obtain the target function value corresponding to the explosion model according to the following formula: Wherein, taking the position of the seismic source as the pole and the ray perpendicular to the ground as the polar axis, a polar coordinate system is established. is the polar angle. is the polar radius. and are the polar angles corresponding to the left and right boundaries of the effective area. and are the upper and lower limits of the polar radius. is the particle pressure value at the corresponding position. is the objective function value corresponding to the explosion model.

6. An explosive column determination device, characterized in that, Including: A memory and a processor; The memory is used to store program instructions; The processor is used to call the program instructions in the memory to execute the explosive column determination method according to any one of claims 1-4.

7. A computer-readable storage medium, characterized in that, Computer program instructions are stored in the computer-readable storage medium, and when the computer program instructions are executed, the explosive column determination method according to any one of claims 1 to 4 is implemented.

8. A computer program product comprising a computer program, characterized in that, When the computer program is executed by the processor, the explosive column determination method according to any one of claims 1-4 is implemented.

Images