A method and system for obtaining optimal triggers

Abstract

The application provides a method and system for obtaining optimal excitation factors, and belongs to the field of oil and gas geophysical exploration. The method for obtaining optimal excitation factors adopts a frequency scanning mode to obtain a signal-to-noise ratio, then utilizes the signal-to-noise ratio to obtain a plurality of signal-to-noise ratio curves, and finds an intersection area of the plurality of signal-to-noise ratio curves, and finally determines optimal excitation factors according to the intersection area; the optimal excitation factors include optimal well depth and optimal dosage. The optimal excitation factors obtained by the application can simultaneously consider the signal-to-noise ratio and the resolution. Through actual data verification, the method has good effect, can well obtain the optimal excitation factors, and thus improves the quality of seismic data.

Description

Technical Field

[0001] This invention belongs to the field of oil and gas geophysical exploration, specifically relating to a method and system for obtaining optimal excitation factors. Background Technology

[0002] In seismic exploration, triggering factors include well depth and charge quantity. Currently, in the initial stages of seismic exploration, it is necessary to select an optimal triggering factor through experiments. Only the optimal triggering factor can guarantee the best quality of the obtained seismic data. However, current methods are relatively simple, including qualitative and quantitative analysis. Qualitative analysis involves placing different triggering factors in individual shots and visually identifying the best one; this method is heavily influenced by subjective factors. Quantitative analysis analyzes the frequency and signal-to-noise ratio of data from several triggering factors, ultimately relying on human judgment to determine the optimal factor; this method also carries a degree of subjectivity and error. Summary of the Invention

[0003] The purpose of this invention is to solve the problems existing in the prior art and provide a method and system for obtaining the optimal excitation factors, which can simultaneously take into account the signal-to-noise ratio and resolution, obtain the optimal excitation factors for seismic exploration, and improve the quality of seismic data.

[0004] This invention is achieved through the following technical solution:

[0005] In a first aspect, the present invention provides a method for obtaining optimal excitation factors. The method uses a frequency division scanning method to obtain the signal-to-noise ratio (SNR), then uses the SNR to obtain multiple SNR curves, and finds the intersection region of the multiple SNR curves. Finally, the optimal excitation factors are determined based on the intersection region. The optimal excitation factors include: optimal well depth and optimal charge.

[0006] A further improvement of the present invention is that:

[0007] The method includes:

[0008] (1) Conduct activation tests to obtain single-shot data at different well depths under a set charge amount;

[0009] (2) The optimal well depth is obtained by using single shot data at different well depths;

[0010] (3) At the optimal well depth, different amounts of explosives were used to conduct activation tests to obtain single-shot data under different amounts of explosives.

[0011] (4) The optimal amount of gunpowder is obtained by using single-shot data under different amounts of gunpowder.

[0012] A further improvement of the present invention is that:

[0013] The operation of step (1) includes:

[0014] (11) Determine the set amount of pesticide at the test point based on the amount of pesticide used in the surrounding work areas;

[0015] (12) Select a test point in the work area and conduct activation tests at different well depths using a set amount of explosives to obtain single-shot data at different well depths.

[0016] A further improvement of the present invention is that:

[0017] The operation of step (2) includes:

[0018] (21) The single-shot data at different well depths are divided into frequency and scanned to obtain the filtered single-shot data;

[0019] (22) The signal-to-noise ratio is obtained by using the filtered data from a single shot.

[0020] (23) Draw the signal-to-noise ratio diagram for the well depth;

[0021] (24) The optimal well depth is obtained by using the signal-to-noise ratio plot of the well depth.

[0022] A further improvement of the present invention is that:

[0023] The operation of step (23) includes:

[0024] On a coordinate graph with well depth as the horizontal axis and signal-to-noise ratio as the vertical axis, first plot the points of signal-to-noise ratio within different frequency filtering ranges at each well depth. Then, connect the points of signal-to-noise ratio at different well depths within the same frequency filtering range with straight lines to form the signal-to-noise ratio curve of that frequency filtering range, thus obtaining the signal-to-noise ratio map of the well depth.

[0025] A further improvement of the present invention is that:

[0026] The operation of step (24) includes:

[0027] In the signal-to-noise ratio (SNR) plot of well depth, all the SNR curves intersect in an intersection region, and the well depth closest to the intersection region is the optimal well depth.

[0028] A further improvement of the present invention is that:

[0029] The operation of step (3) includes:

[0030] Activation tests were conducted at the optimal well depth using progressively increasing amounts of propellant to obtain single-shot data for different propellant amounts.

[0031] A further improvement of the present invention is that:

[0032] The operation of step (4) includes:

[0033] (41) The single-shot data under different amounts of gunpowder are divided into frequency scanning filters to obtain the filtered single-shot data;

[0034] (42) The signal-to-noise ratio is obtained by using the filtered data from a single shot.

[0035] (43) Draw a signal-to-noise ratio diagram of drug dosage;

[0036] (44) The optimal drug dosage is obtained by using the signal-to-noise ratio plot of the drug dosage.

[0037] A further improvement of the present invention is that:

[0038] The operation of step (43) includes:

[0039] On a coordinate graph with drug dosage as the horizontal axis and signal-to-noise ratio as the vertical axis, first plot the points of signal-to-noise ratio in different frequency filtering ranges for each drug dosage. Then, connect the points of each signal-to-noise ratio for different drug dosages in the same frequency filtering range with straight lines to form the signal-to-noise ratio curve of that frequency filtering range, thus obtaining the signal-to-noise ratio graph for the drug dosage.

[0040] A further improvement of the present invention is that:

[0041] The operation of step (44) includes:

[0042] In the signal-to-noise ratio (SNR) graph of drug dosage, all SNR curves intersect in an intersection region, and the drug dosage closest to the intersection region is the optimal drug dosage.

[0043] A further improvement of the present invention is that:

[0044] The frequency filter used in the frequency division scanning filter in steps (21) and (41) is:

[0045] 10-20Hz, 20-40Hz, 30-60Hz, 40-80Hz, 50-100Hz.

[0046] A second aspect of the present invention provides a system for obtaining optimal excitation factors, the system comprising:

[0047] The first single-shot data acquisition unit is used to acquire single-shot data at different well depths with a set charge amount obtained after the activation test.

[0048] Optimal well depth acquisition unit: connected to the first single shot data acquisition unit, used to obtain the optimal well depth using single shot data at different well depths;

[0049] The second single-shot data acquisition unit is connected to the optimal well depth acquisition unit and is used to acquire single-shot data at different charge levels obtained after the activation test at the optimal well depth.

[0050] Optimal charge acquisition unit: connected to the second single-shot data acquisition unit, used to obtain the optimal charge using single-shot data under different charge rates.

[0051] A third aspect of the present invention provides a computer-readable storage medium storing at least one computer-executable program, which, when executed by the computer, causes the computer to perform the steps in the method for obtaining the optimal stimulating factors described above.

[0052] Compared with the prior art, the beneficial effects of the present invention are:

[0053] The optimal excitation factor obtained using this invention can simultaneously achieve a good balance between signal-to-noise ratio and resolution. Verification with actual data shows that the method of this invention is effective in obtaining the optimal excitation factor, thereby improving the quality of seismic data. Attached Figure Description

[0054] Figure 1 A flowchart illustrating the steps of the method of this invention;

[0055] Figure 2 Signal-to-noise ratio diagram of well depth in embodiments of the present invention;

[0056] Figure 3 Signal-to-noise ratio diagram of drug dosage in embodiments of the present invention. Detailed Implementation

[0057] The present invention will now be described in further detail with reference to the accompanying drawings:

[0058] This invention is applicable to seismic acquisition work in petroleum exploration and provides a method for obtaining optimal excitation factors. The method uses frequency-division scanning to obtain the signal-to-noise ratio (SNR), then uses the SNR to generate multiple SNR curves, and identifies the intersection region of these curves. This intersection region can simultaneously consider both SNR and resolution. Using the excitation factors determined based on this intersection region as the optimal excitation factors for this seismic exploration is scientifically sound and can effectively improve data quality.

[0059] The invention first selects a test point in the work area to conduct an activation test. The activation factors are designed in a progressive manner. For example, the activation well depth is designed to be 1, 2, 3, 5, 7, and 9 m below the high-speed layer, and the activation dosage is designed to be progressive, such as 4, 5, 6, 7, 8, 9, and 10 kg.

[0060] In practice, when conducting well depth tests, the same amount of explosive is always used. However, testing each amount of explosive would be too costly and unnecessary. Therefore, this invention first uses a pre-set amount of explosive for testing. Generally, the test point uses the amount of explosive used in the surrounding work area as the pre-set amount. For example, if the surrounding work area uses 8 kg of explosive, then the pre-set amount for this test point is also 8 kg, i.e., 8 kg of explosive is used for the excitation test. Through the excitation test, single-shot data at different well depths are obtained, and the optimal well depth under the pre-set amount of explosive is obtained using the single-shot data at different well depths. Then, at the optimal well depth, different amounts of explosive are used for testing to obtain the optimal amount of explosive. The optimal well depth and the optimal amount of explosive constitute the optimal excitation factor.

[0061] like Figure 1 As shown, embodiments of the method of the present invention are as follows:

[0062] Example 1

[0063] (1) Conduct activation tests to obtain single-shot data at different well depths under a set charge amount;

[0064] (11) Determine the set amount of pesticide at the test point based on the amount of pesticide used in the surrounding work areas;

[0065] (12) Select a test point in the work area and conduct tests at different well depths using a set amount of explosives to obtain single-shot data at different well depths. The single-shot data refers to the seismic wave information obtained, including waveform, amplitude, frequency and other information.

[0066] The experimental factors were designed in a progressive manner: the well depth was designed in a progressive manner based on the depth of the near-surface high-velocity layer top interface. For example, if the depth of the high-velocity layer top interface is h, then the designed well depths would be h+1, h+2, h+3, h+4, h+5, h+6, h+7, h+8, h+9... In this embodiment, the excitation well depth was designed as 1, 2, 3, 5, 7, and 9 meters below the high-velocity layer (the well depth was designed below the high-velocity layer to avoid the shielding effect of the high-velocity layer interface on seismic waves).

[0067] (2) Obtain the optimal well depth using single-shot data at different well depths;

[0068] (21) The single-shot data at different well depths were subjected to frequency division scanning filtering to obtain the filtered single-shot data:

[0069] Frequency-division filtering is performed on single-shot data for each well depth. Preferably, frequency-division filtering is performed at 10-20Hz, 20-40Hz, 30-60Hz, 40-80Hz, and 50-100Hz to obtain filtered single-shot data. These frequency filtering settings are commonly used in the geophysical exploration industry. More frequency-division filtering can be performed as needed, but these are the most prevalent in the industry.

[0070] (22) The signal-to-noise ratio is obtained using the filtered data from a single shot:

[0071] The filtered data from each individual shot were extracted into 40-50ms segments based on the target layer location and analyzed to determine the signal-to-noise ratio.

[0072] The method of opening a time window is an existing technology, which is briefly described below:

[0073] The detonation point of the explosive is 1-9m below the high-velocity layer. The target layer is located at a depth of 3000-4000m, which corresponds to 2-3s in a single-shot record. Since the vertical axis of a single-shot record is time, recording the travel time of seismic waves as they propagate underground and reflect back to the surface, time represents the depth. The method for opening a time window based on the target layer is as follows: the reflected wave from the target layer is a wave group. A rectangular window is drawn according to the location of the wave group, containing a portion of the reflected wave from the target layer.

[0074] The method for determining the signal-to-noise ratio is an existing technology, which is briefly described below:

[0075] The signal-to-noise ratio (SNR) is calculated by first opening a time window (time window) at the target layer location in the seismic data. Specifically, a 40-50ms segment of data from 20-30 channels (the general range is around 20-30 channels, but can be adjusted according to actual needs) is extracted for analysis. (20-30 channels represent the horizontal range, and 40-50ms represents the vertical range; the chosen range should adequately encompass the reflected waves from the target layer.) This data is then subjected to a Fourier transform to the frequency domain. The energy at a frequency of 15-50Hz (because the frequency of reflected seismic waves from the target layer is generally 15-50Hz) is defined as the effective signal energy, and the energy from other regions is defined as noise energy. The ratio of the effective signal energy to the noise energy is calculated, yielding the SNR. The 15-50Hz frequency range is used for SNR calculation, while the preceding frequency division filtering is used to convert the single-shot data to different frequency ranges. For example, first convert the single-shot data to a frequency range of 30-60Hz, and then calculate the target layer signal-to-noise ratio within that frequency range.

[0076] Specifically, the data from each individual shot is divided into frequency filters of 10-20Hz, 20-40Hz, 30-60Hz, 40-80Hz, and 50-100Hz. Then, a signal-to-noise ratio is calculated for each frequency filter of each individual shot.

[0077] In this way, the signal-to-noise ratio of the target layer within each frequency filter setting at each well depth can be calculated.

[0078] (23) Draw the signal-to-noise ratio diagram for the well depth:

[0079] The signal-to-noise ratio (SNR) of the well depth is plotted using the SNR of the same frequency filter settings, as shown below:

[0080] On a coordinate graph with well depth on the x-axis and signal-to-noise ratio (SNR) on the y-axis, first, plot the SNR points for different frequency filter levels at each well depth. Then, connect the SNR points for different well depths within the same frequency filter level with straight lines to form the SNR curve for that frequency filter level, thus obtaining the SNR map for the well depth. The number of curves corresponds to the number of frequency filter levels, with each curve representing one frequency filter level.

[0081] For example, the signal-to-noise ratio (SNR) of well depth data in the 10-20Hz range can be plotted as an SNR curve according to the order of 1, 2, 3, 5, 7, and 9m below the high-speed layer. Each curve represents the data of a frequency filter range, and the curve shows the SNR of all well depths at that frequency filter range corresponding to the same dosage. Figure 2 The five curves in the figure correspond to the signal-to-noise ratios at five frequency filtering settings (10-20Hz, 20-40Hz, 30-60Hz, 40-80Hz, and 50-100Hz) at well depths of 1, 2, 3, 5, 7, and 9, respectively, with an 8kg dosage of reagent. For example, Figure 2 In the 10-20Hz curve, the dots correspond to the signal-to-noise ratios at depths of wells 1, 2, 3, 5, 7, and 9 with an 8kg dosage of reagent. Connecting these dots with straight lines in sequence yields the signal-to-noise ratio curve for that frequency filtering setting.

[0082] (24) The optimal well depth is obtained using the signal-to-noise ratio plot of the well depth:

[0083] In the signal-to-noise ratio (SNR) plot of well depth, all the SNR curves will intersect in an intersection region. The well depth closest to the intersection region is the most reasonable well depth, i.e., the optimal well depth.

[0084] Specifically, Figure 2 In the illustrated embodiment, the five signal-to-noise ratio (SNR) curves roughly intersect at the same location, but not completely at the same point; that is, they intersect within an intersection region. Within this region, the SNR curves are closest to each other, representing the point where the curves intersect most densely. Figure 2In the illustrated embodiment, a well depth of 3m is the closest to the confluence area, therefore 3m is the most reasonable well depth, i.e., the optimal well depth. This yields the optimal well depth for the set dosage (8 kg in this embodiment). The well depth must be an integer. If the optimal well depth is chosen to be 3.2m or 2.9m, it would be impractical in actual field drilling. Therefore, a general range is sufficient; the integer well depth closest to the confluence area is the optimal well depth.

[0085] (3) At the optimal well depth, different charge concentrations were used to conduct activation tests to obtain single-shot data under different charge concentrations:

[0086] In this embodiment, the activation charge dosage was designed to be 4, 5, 6, 7, 8, 9, and 10 kg. A test point was selected in the work area, and activation tests were conducted at the optimal well depth using charge dosages of 4, 5, 6, 7, 8, 9, and 10 kg respectively to obtain single-shot data under different charge dosages.

[0087] (4) Obtain the optimal charge rate using single-shot data under different charge rates:

[0088] The processing method in step (4) is similar to that in step (2), as follows:

[0089] (41) The single-shot data under different charges were subjected to frequency division scanning filtering to obtain the filtered single-shot data:

[0090] Frequency division filtering is performed on the single-shot data for each charge. Preferably, frequency division filtering is performed according to the following filter levels: 10-20Hz, 20-40Hz, 30-60Hz, 40-80Hz, and 50-100Hz, to obtain the filtered single-shot data.

[0091] (42) The signal-to-noise ratio is obtained using the filtered data from a single shot:

[0092] The filtered data from each individual shot were extracted into 40-50ms segments based on the target layer location and analyzed to determine the signal-to-noise ratio.

[0093] The method of opening a time window is an existing technology, which is briefly described below:

[0094] The detonation point of the explosive is 3 meters below the high-velocity layer. The target layer is located at a depth of 3000-4000 meters underground, which corresponds to the 2-3 second mark in a single-shot record. Since the vertical axis of a single-shot record is time, recording the travel time of seismic waves as they propagate underground and reflect back to the surface, time represents the depth. The method for opening a time window based on the target layer is as follows: the reflected wave from the target layer is a wave group. A rectangular window is drawn according to the location of the wave group, containing a portion of the reflected wave from the target layer.

[0095] The method for determining the signal-to-noise ratio is an existing technology, which is briefly described below:

[0096] The signal-to-noise ratio (SNR) is calculated by first opening a time window (time window) at the target layer location in the seismic data. Specifically, a 40-50ms segment of data from 20-30 channels (the general range is around 20-30 channels, but can be adjusted according to actual needs) is extracted for analysis. (20-30 channels represent the horizontal range, and 40-50ms represents the vertical range; the chosen range should adequately encompass the reflected waves from the target layer.) This data is then subjected to a Fourier transform to the frequency domain. The energy at a frequency of 15-50Hz (because the frequency of reflected seismic waves from the target layer is generally 15-50Hz) is defined as the effective signal energy, and the energy from other regions is defined as noise energy. The ratio of the effective signal energy to the noise energy is calculated, yielding the SNR. The 15-50Hz frequency range is used for SNR calculation, while the preceding frequency division filtering is used to convert the single-shot data to different frequency ranges. For example, first convert the single-shot data to a frequency range of 30-60Hz, and then calculate the target layer signal-to-noise ratio within that frequency range.

[0097] Specifically, the data from each individual shot is divided into frequencies of 10-20Hz, 20-40Hz, 30-60Hz, 40-80Hz, and 50-100Hz. Then, a signal-to-noise ratio is calculated for each frequency division of each individual shot.

[0098] In this way, the signal-to-noise ratio of the target layer within each frequency filter setting under each drug dosage can be calculated.

[0099] (43) Draw the signal-to-noise ratio diagram of drug dosage:

[0100] A signal-to-noise ratio (SNR) plot of the drug dosage is drawn using the SNR of the same filter setting, as shown below:

[0101] On a coordinate graph with drug dosage as the horizontal axis and signal-to-noise ratio as the vertical axis, first plot the points of signal-to-noise ratio under different frequency filtering levels for each drug dosage. Then, connect the points of each signal-to-noise ratio under different drug dosages within the same frequency filtering level with straight lines to form the signal-to-noise ratio curve of that frequency filtering level, thus obtaining the signal-to-noise ratio graph for the drug dosage. The number of curves corresponds to the number of frequency filtering levels, and each curve corresponds to one frequency filtering level.

[0102] For example, the signal-to-noise ratio (SNR) of well depth data in the 10-20Hz range can be plotted according to the order of reagent dosage: 4, 5, 6, 7, 8, 9, 10 kg. Each SNR curve represents the data of a frequency filter range, and the curve shows the SNR of all reagent dosages at the same well depth within that frequency filter range. Figure 3The five curves in the image correspond to the signal-to-noise ratio curves at five frequency filtering levels (10-20Hz, 20-40Hz, 30-60Hz, 40-80Hz, and 50-100Hz) at a well depth of 3 meters, with dosages of 4, 5, 6, 7, 8, 9, and 10 kg of reagent. For example, Figure 3 In the 10-20Hz curve, the dots correspond to the signal-to-noise ratios of 4, 5, 6, 7, 8, 9, and 10 kg of reagent at a well depth of 3 meters. Connecting these dots with straight lines in sequence forms the signal-to-noise ratio curve for that frequency filtering setting.

[0103] (44) The optimal drug dosage is obtained using the signal-to-noise ratio plot of the drug dosage:

[0104] In the signal-to-noise ratio (SNR) graph of drug dosage, the various SNR curves will converge in an intersection region. The drug dosage closest to the intersection region is the most reasonable drug dosage, i.e., the optimal drug dosage.

[0105] Specifically, Figure 3 In the illustrated embodiment, the five signal-to-noise ratio (SNR) curves roughly intersect at the same location, but not completely at the same point; that is, they intersect within an intersection region, where the SNR curves are closest to each other. Specifically... Figure 3 In the illustrated embodiment, the 6 kg position represents the dosage closest to the junction area; therefore, 6 kg is the most reasonable dosage, i.e., the optimal dosage. This yields the optimal dosage for the optimal well depth.

[0106] Thus, the optimal well depth for this embodiment is 3 meters, the optimal dosage is 6 kg, and the optimal activating factor is obtained.

[0107] The present invention also provides a system for obtaining optimal excitation factors, an embodiment of which is as follows:

[0108]

Example 2

[0109] The system includes:

[0110] The first single-shot data acquisition unit is used to acquire single-shot data at different well depths with a set charge amount obtained after the activation test.

[0111] Optimal well depth acquisition unit: connected to the first single shot data acquisition unit, used to obtain the optimal well depth using single shot data at different well depths;

[0112] The second single-shot data acquisition unit is connected to the optimal well depth acquisition unit and is used to acquire single-shot data at different charge levels obtained after the activation test at the optimal well depth.

[0113] Optimal charge acquisition unit: connected to the second single-shot data acquisition unit, used to obtain the optimal charge using single-shot data under different charge rates.

[0114] The present invention also provides a computer-readable storage medium, embodiments of which are as follows:

[0115]

Example 3

[0116] The computer-readable storage medium stores at least one computer-executable program, which, when executed by the computer, causes the computer to perform the steps in the method for obtaining the optimal stimulating factors described above.

[0117] Finally, it should be noted that the above technical solution is only one embodiment of the present invention. For those skilled in the art, based on the application methods and principles disclosed in the present invention, it is easy to make various types of improvements or modifications, and not limited to the methods described in the above specific embodiments of the present invention. Therefore, the methods described above are only preferred and have no limiting significance.

Claims

1. A method of obtaining optimal priming factors, characterized by: The method comprises the following steps: Step 1, performing excitation test to obtain single-shot data at different well depths under a set explosive amount, comprising: Step 11, determining the set explosive amount of the test point according to the explosive amount used in the surrounding work area; Step 12, selecting a test point in the work area, and performing excitation test at different well depths under the set explosive amount to obtain single-shot data at different well depths; Step 2, obtaining the optimal well depth by using the single-shot data at different well depths, comprising: Step 21, performing frequency scanning filtering on the single-shot data at different well depths to obtain filtered single-shot data; Step 22, obtaining signal-to-noise ratio by using the filtered single-shot data; Step 23, drawing a well depth signal-to-noise ratio graph: on a coordinate graph with well depth as the horizontal coordinate and signal-to-noise ratio as the vertical coordinate, first draw the signal-to-noise ratio points in different frequency filter ranges at each well depth, then connect the signal-to-noise ratio points at different well depths in the same frequency filter range with straight lines to form the signal-to-noise ratio curve of the frequency filter range, and obtain the well depth signal-to-noise ratio graph; Step 24, obtaining the optimal well depth by using the well depth signal-to-noise ratio graph: in the well depth signal-to-noise ratio graph, the signal-to-noise ratio curves intersect in an intersection area, and the well depth closest to the intersection area is the optimal well depth; Step 3, performing excitation test at the optimal well depth under different explosive amounts to obtain single-shot data under different explosive amounts; Step 4, obtaining the optimal explosive amount by using the single-shot data under different explosive amounts, comprising: Step 41, performing frequency scanning filtering on the single-shot data under different explosive amounts to obtain filtered single-shot data; Step 42, obtaining signal-to-noise ratio by using the filtered single-shot data; Step 43, drawing an explosive amount signal-to-noise ratio graph: on a coordinate graph with explosive amount as the horizontal coordinate and signal-to-noise ratio as the vertical coordinate, first draw the signal-to-noise ratio points in different frequency filter ranges under each explosive amount, then connect the signal-to-noise ratio points under different explosive amounts in the same frequency filter range with straight lines to form the signal-to-noise ratio curve of the frequency filter range, and obtain the explosive amount signal-to-noise ratio graph; Step 44, obtaining the optimal explosive amount by using the explosive amount signal-to-noise ratio graph: in the explosive amount signal-to-noise ratio graph, the signal-to-noise ratio curves intersect in an intersection area, and the explosive amount closest to the intersection area is the optimal explosive amount.

2. The method of obtaining optimal trigger factors according to claim 1, wherein: The operation of step 3 comprises: Performing excitation test at the optimal well depth under progressively increasing explosive amounts to obtain single-shot data under different explosive amounts.

3. The method of obtaining optimal trigger factors as claimed in claim 1 wherein: The frequency filter ranges used in the frequency scanning filtering in steps 21 and 41 are: 10-20 Hz, 20-40 Hz, 30-60 Hz, 40-80 Hz, and 50-100 Hz.

4. A system for obtaining optimal stimuli based on the method of any of claims 1-3, characterized by: The system comprises: A first single-shot data collection unit for collecting single-shot data at different well depths under a set explosive amount obtained after excitation test; An optimal well depth acquisition unit connected with the first single-shot data collection unit and configured to obtain the optimal well depth by using the single-shot data at different well depths; A second single-shot data collection unit connected with the optimal well depth acquisition unit and configured to collect single-shot data under different explosive amounts obtained after excitation test at the optimal well depth; An optimal charge amount acquisition unit, connected with the second single-shot data acquisition unit, is configured to obtain an optimal charge amount by using single-shot data under different charge amounts.

5. A computer-readable storage medium, characterized in that: The computer readable storage medium stores at least one program executable by the computer, and the at least one program, when executed by the computer, causes the computer to perform the steps in the method for obtaining an optimal excitation factor according to any one of claims 1-3.

Images