Detection method, device and equipment with multiple dimensions and multiple visual angles

Abstract

The embodiment of the invention relates to the technical field of detection imaging, and provides a multi-dimensional and multi-view detection method, device and equipment, and the method comprises the steps: obtaining all single-transmitting-end response signals of a target detection region, and obtaining a total-transmitting-end total response data set; decomposing the total response data set of all transmitting ends into a plurality of sub-data sets, and setting a corresponding detection view angle; calculating the time translation size of the response signal of each single transmitting end according to a preset time domain translation rule; translating each single transmitting end response signal in a time domain according to the time translation size; superposing the translated single transmitting end response signals to obtain a synthetic response signal of each sub-data set; and obtaining a multi-dimensional multi-view total response data set of the target detection area according to the synthesized response signal of each sub-data set. According to the embodiment of the invention, the signal-to-noise ratio of the data of the target area can be improved, so that the resolution of detection imaging is improved.

Description

Technical Field

[0001] The embodiments of this specification relate to the field of detection imaging technology, and in particular to a detection method, device, and equipment with multi-dimensional and multi-viewing angles. Background Art

[0002] In active detection imaging technology, whether it is seismic wave methods in geological exploration or ultrasound, electromagnetic waves (such as MRI) and other methods in medical imaging, the core process relies on the basic framework of "signal transmission-reception-processing-imaging". As the exploration level continues to deepen, the detailed evaluation of the target area medium becomes increasingly urgent. For example, thin layer identification in oil and gas exploration, detection of hidden disaster-causing geological bodies in coal mine safety monitoring, and accurate imaging of Moho surface boundaries in deep seismic exploration. However, existing detection imaging technology has the following bottlenecks in data processing and imaging accuracy:

[0003] In land seismic exploration, multiple-coverage seismic acquisition technology is often used as a conventional means of improving the signal-to-noise ratio of seismic data. However, traditional multiple-coverage acquisition technology has limited adaptability to specific geological conditions commonly encountered in field operations, such as areas with dense rivers and lakes, target areas with dramatic mountainous terrain, and residential areas in urban seismic exploration. These areas often have irregularly arranged seismic observation systems, resulting in limited coverage times for larger target areas, making it difficult to achieve ideal exploration results. For small geological bodies, such as small faults and small-diameter collapse columns, the reflected energy of reflected waves passing through these bodies is very weak, and data from a single perspective or limited angles makes it difficult to accurately focus on the target reflected waves, resulting in blurred imaging and insufficient noise suppression.

[0004] For medical ultrasound imaging, ultrasonic phased array technology transmits sound waves through a combination of array elements. However, existing methods often treat all array elements as a single sound source, failing to fully utilize the multi-angle echo information generated by a single element. For deep tissue or strongly scattering media, the echo superposition of a single synthetic perspective is susceptible to sidelobe interference, which affects the resolution of detailed lesions.

[0005] Therefore, a new detection method is urgently needed to improve the signal-to-noise ratio and resolution of complex medium data to meet the development needs of high-precision detection technology. Summary of the Invention

[0006] In response to the above-mentioned problems in the prior art, the purpose of the embodiments of this specification is to provide a detection method, device and equipment with multi-dimensional and multi-perspective capabilities to solve the problems of low signal-to-noise ratio and low resolution when using existing detection methods to image complex media.

[0007] In order to solve the above technical problems, the specific technical solutions of the embodiments of this specification are as follows:

[0008] On the one hand, embodiments of this specification provide a multi-dimensional and multi-perspective detection method, the method comprising:

[0009] Obtain all single transmitter response signals in the target detection area to obtain the total response data set of all transmitters;

[0010] Decomposing the total transmitter response data set into a plurality of sub-data sets, and setting a corresponding detection viewing angle for each sub-data set, wherein each sub-data set includes one or more single transmitter response signals;

[0011] Calculate the time shift of each single transmitter response signal in each sub-data set at the corresponding detection viewing angle according to a preset time domain shift rule;

[0012] Shifting each single transmitter response signal in each sub-data set in the time domain according to the time shift size;

[0013] Superimposing the response signals of each single transmitter in each sub-data set after translation to obtain a synthetic response signal of each sub-data set;

[0014] Obtaining a multi-dimensional and multi-view total response data set of the target detection area according to the synthetic response signal of each sub-data set;

[0015] Imaging processing is performed using the multi-dimensional and multi-view total response data set.

[0016] Furthermore, the time domain translation rule includes a first translation rule and a second translation rule;

[0017] The calculating, according to a preset time domain shift rule, the time shift size of each single transmitting end response signal in each sub-data set at the corresponding detection viewing angle includes:

[0018] determining a size of the target detection area;

[0019] Determining whether the size of the target detection area is greater than a preset threshold;

[0020] If yes, then calculating the time shift of each single transmitting end response signal in each sub-data set at the corresponding detection viewing angle according to the first shift rule;

[0021] If not, the time shift of each single transmitting end response signal in each sub-data set at the corresponding detection viewing angle is calculated according to the second shift rule.

[0022] Furthermore, the calculating, according to the first translation rule, the time translation size of each single transmitter response signal in each sub-data set at the corresponding detection viewing angle includes:

[0023] Calculating the head-to-tail transmitter distance of each sub-data set;

[0024] Determining a detection viewing angle direction according to the detection viewing angle of each sub-data set;

[0025] The time shift of the response signal of each single transmitter in each sub-data set is calculated based on the head-to-tail transmitter distance, detection viewing angle, detection viewing direction and seismic wave velocity of each sub-data set.

[0026] Furthermore, the time shift of the response signal of each single transmitter in each sub-data set is calculated based on the head-to-tail transmitter spacing, the detection viewing angle, the detection viewing direction, and the seismic wave velocity of each sub-data set, including:

[0027] The time shift of the response signal of each single transmitter in each sub-data set is calculated using the following formula:

[0028] ;

[0029] in, Indicates the time shift of the response signal of the mth single transmitter in the jth sub-data set, represents the distance between the first and last transmitters in the jth sub-dataset, Represents the detection viewing angle of the jth sub-dataset, when When , it means that the detection angle is in the direction of increasing offset; when When , it means that the detection viewing angle is in the direction of decreasing offset, and c represents the seismic wave velocity.

[0030] Furthermore, the calculating, according to the second shift rule, the time shift of each single-transmitter response signal in each sub-data set at the corresponding detection viewing angle includes:

[0031] Determine the wave focusing point based on the spatial position and depth of the target detection area;

[0032] Calculating the distance from each transmitting end in each sub-data set to the wave focusing point;

[0033] Determining a detection viewing angle direction according to the detection viewing angle of each sub-data set;

[0034] The time shift of the response signal of each single transmitter in each sub-data set is calculated based on the distance from each transmitter in each sub-data set to the wave focusing point, the detection viewing angle direction and the seismic wave velocity.

[0035] Furthermore, the time shift of the response signal of each single transmitter in each sub-data set is calculated based on the distance from each transmitter in each sub-data set to the wave focusing point, the detection viewing angle direction, and the seismic wave velocity, including:

[0036] The time shift of the response signal of each single transmitter in each sub-data set is calculated using the following formula:

[0037] ;

[0038] in, Indicates the time shift of the response signal of the mth single transmitter in the jth sub-data set, represents the distance from the first transmitter to the wave focusing point in the jth sub-data set, represents the distance from the tail emission end to the wave focusing point in the jth sub-data set, represents the distance from the i-th transmitter to the wave focusing point in the j-th sub-data set, Represents the detection viewing angle of the jth sub-dataset, when When , it means that the detection angle is in the direction of increasing offset; when When , it means that the detection viewing angle is in the direction of decreasing offset, and c represents the seismic wave velocity.

[0039] Furthermore, the calculating of the distance from each transmitting end in each sub-data set to the wave focusing point includes:

[0040] Calculate the first distance from the central transmitting end to the wave focusing point in each sub-data set;

[0041] Calculating a second distance from the first transmitter to the central transmitter and a third distance from the first transmitter to each transmitter in each sub-data set;

[0042] The distance from each transmitting end in each sub-data set to the wave focusing point is calculated according to the first distance, the second distance, and the third distance.

[0043] Furthermore, calculating the distance from each transmitting end in each sub-data set to the wave focusing point based on the first distance, the second distance, and the third distance includes:

[0044] The distance from each transmitter to the wave focusing point in each sub-data set is calculated using the following formula:

[0045] ;

[0046] ;

[0047] ;

[0048] in, represents the distance from the first transmitter to the wave focusing point in the jth sub-data set, represents the distance from the tail emission end to the wave focusing point in the jth sub-data set, represents the distance from the i-th transmitter to the wave focusing point in the j-th sub-data set, represents the detection viewing angle of the j-th sub-dataset, represents the first distance from the central transmitter to the wave focusing point in the jth sub-dataset, represents the second distance from the first transmitter to the center transmitter in the jth sub-data set, Represents the third distance from the first transmitter to the i-th transmitter in the j-th sub-data set.

[0049] On the other hand, an embodiment of this specification provides a detection device with multi-dimensional and multi-viewing angles, the device comprising:

[0050] An acquisition module is used to acquire the response signals of all single transmitters in the target detection area and obtain a total response data set of all transmitters;

[0051] a decomposition module, configured to decompose the total transmitter response data set into a plurality of sub-data sets and set a corresponding detection viewing angle for each sub-data set, wherein each sub-data set includes one or more single transmitter response signals;

[0052] A calculation module, configured to calculate the time shift of each single transmitter response signal in each sub-data set at a corresponding detection viewing angle according to a preset time domain shift rule;

[0053] A translation module, configured to translate each single transmitter response signal in each sub-data set in the time domain according to the time translation size;

[0054] A superposition module is used to superimpose the response signals of each single transmitter in each sub-data set after translation to obtain a synthetic response signal of each sub-data set;

[0055] a total data set determination module, configured to obtain a multi-dimensional and multi-view total response data set of the target detection area according to the synthetic response signal of each sub-data set;

[0056] A processing module is used to perform imaging processing using the multi-dimensional and multi-view total response data set.

[0057] On the other hand, an embodiment of this specification further provides a computer device, including a memory, a processor, and a computer program stored in the memory, wherein when the computer program is run by the processor, the computer program executes instructions of any one of the above methods.

[0058] On the other hand, an embodiment of the present specification further provides a computer-readable storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor of a computer device, executes instructions of any one of the above methods.

[0059] On the other hand, the embodiments of this specification further provide a computer program product, which, when executed by a processor of a computer device, executes instructions of any one of the above methods.

[0060] Using the above-mentioned technical solution, the embodiments of this specification provide a multi-dimensional, multi-perspective detection method. When decomposing the total response dataset of the entire transmitter, a different detection perspective is set for each sub-dataset, so that each sub-dataset corresponds to an independent perspective dimension value. Time domain translation rules are used to simulate the time difference in signal propagation under different perspectives (such as wave velocity differences in geological exploration and acoustic path differences in medical ultrasound). This allows the translated signals to focus on the event axis of the target area when superimposed, enhancing the energy of the effective signal (such as reflected waves and echoes) and suppressing noise in the signal, thereby improving the signal-to-noise ratio of the target area data. A multi-dimensional, multi-perspective total dataset is generated based on the synthesized response signals of each sub-dataset. Because the multi-dimensional, multi-perspective total dataset contains information from different angle domains, it provides a more complete input for subsequent offset imaging, enabling the imaging process to utilize multi-angle constraints, thereby improving the imaging accuracy of complex media (such as geological layered structures and heterogeneous regions of human tissue).

[0061] The above description is only an overview of the technical solutions of some embodiments of this specification. In order to more clearly understand the technical means of some embodiments of this specification, they can be implemented in accordance with the contents of the specification. In order to make the above and other purposes, features and advantages of the embodiments of this specification more obvious and easy to understand, the following specifically cites preferred embodiments and provides detailed descriptions in conjunction with the accompanying drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0062] In order to more clearly illustrate the embodiments of this specification or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are only some embodiments of this specification. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0063] Figure 1 A schematic diagram showing the steps of a multi-dimensional and multi-viewing detection method in some embodiments of this specification is shown;

[0064] FIG2( a ) shows a schematic diagram of a source end synthesizing a multi-dimensional and multi-view wave-focused data set in some embodiments of this specification;

[0065] FIG2( b ) shows a schematic diagram of a receiving end synthesizing a multi-dimensional and multi-view wave-focused data set in some embodiments of this specification;

[0066] FIG2( c ) shows a schematic diagram of geometric seismological implementation of synthesizing multi-dimensional and multi-view wave-focused datasets in some embodiments of this specification;

[0067] FIG2( d ) shows a schematic diagram of a high-coverage observation shot gather after wave focusing in a target area in some embodiments of this specification;

[0068] Figure 3 A schematic diagram showing the steps of calculating the time shift of the response signal of each single transmitting end in some embodiments of this specification is shown;

[0069] FIG4( a ) shows a schematic diagram of beam directionality for a large target area in some embodiments of this specification;

[0070] FIG4( b ) shows a schematic diagram of beam directional focusing in a small target area in some embodiments of this specification;

[0071] Figure 5 A schematic diagram showing the steps of calculating the time shift size of each single transmitting end response signal according to the second shift rule in some embodiments of this specification is shown;

[0072] Figure 6 Shows a schematic diagram of a single-layer geological model and seismic data acquisition in some embodiments of this specification;

[0073] Figure 7 A schematic diagram showing low signal-to-noise ratio seismic data of a single-layer geological model in some embodiments of this specification;

[0074] FIG8( a ) is a schematic diagram showing a delay time scheme for recording each single shot of a focused wave at 700 meters in the vertical direction in some embodiments of this specification;

[0075] FIG8( b ) shows a schematic diagram of a focused wave detection area at 700 meters in the vertical direction in some embodiments of this specification;

[0076] FIG8( c ) shows a concentrated wave recording diagram at 700 meters in the vertical direction in some embodiments of this specification;

[0077] FIG8( d ) shows a reverse time migration imaging result of the focused wave detection area at 700 meters in the vertical direction in some embodiments of this specification;

[0078] FIG9( a ) shows a time delay scheme for recording each single shot of a wave focused at 808.3 meters after deflecting 30° in the direction of increasing offset in some embodiments of this specification;

[0079] FIG9( b ) shows a schematic diagram of a focused wave detection area at 808.3 meters after deflecting 30° toward the direction of increasing offset in some embodiments of this specification;

[0080] FIG9( c ) shows a concentrated wave recording at 808.3 meters after deflection by 30° in the direction of increasing offset in some embodiments of this specification;

[0081] FIG9( d ) shows the reverse time migration imaging result of the focused wave detection area at 808.3 meters after deflection by 30° in the direction of increasing offset distance in some embodiments of this specification;

[0082] FIG10( a ) shows a time delay scheme for recording each single shot of a wave focused at 808.3 meters after deflecting 30° toward the direction of decreasing offset in some embodiments of this specification;

[0083] FIG10( b ) shows a schematic diagram of a focused wave detection area at 808.3 meters after deflecting 30° toward the direction of decreasing offset in some embodiments of this specification;

[0084] FIG10( c ) shows a wave gathering record at 808.3 meters after deflecting 30° toward the direction of decreasing offset in some embodiments of this specification;

[0085] FIG10( d ) shows the reverse time migration imaging result of the focused wave detection area at 808.3 meters after deflection by 30° in the direction of decreasing offset in some embodiments of this specification;

[0086] Figure 11 A schematic structural diagram of a multi-dimensional and multi-viewing detection device in some embodiments of this specification is shown;

[0087] Figure 12 A schematic structural diagram of a computer device in this specification is shown.

[0088] Description of the accompanying symbols:

[0089] 1101. Get module;

[0090] 1102. Decomposition module;

[0091] 1103. Calculation module;

[0092] 1104, translation module;

[0093] 1105, superposition module;

[0094] 1106. Total data set determination module;

[0095] 1107, processing module;

[0096] 1202. Computer equipment;

[0097] 1204, processor;

[0098] 1206. Memory;

[0099] 1208, driving mechanism;

[0100] 1210, input / output module;

[0101] 1212. Input device;

[0102] 1214. Output device;

[0103] 1216. Presentation equipment;

[0104] 1218. Graphical User Interface;

[0105] 1220, network interface;

[0106] 1222, communication link;

[0107] 1224. Communication bus. DETAILED DESCRIPTION

[0108] The following will be combined with the drawings in the embodiments of this specification to clearly and completely describe the technical solutions in the embodiments of this specification. Obviously, the embodiments described are only part of the embodiments of this specification, not all of the embodiments. Based on the embodiments in this specification, all other embodiments obtained by ordinary technicians in this field without making any creative efforts are within the scope of protection of this specification.

[0109] It should be noted that the terms "first," "second," and the like in this specification, the claims, and the accompanying drawings are used to distinguish similar objects and are not necessarily used to describe a particular order or precedence. It should be understood that the terms used in this manner are interchangeable where appropriate, so that the embodiments of this specification described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having," and any variations thereof, are intended to cover non-exclusive inclusions. For example, a process, method, apparatus, product, or device comprising a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units that are not explicitly listed or that are inherent to such processes, methods, products, or devices.

[0110] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, storage, and display, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties. Furthermore, the acquisition, storage, use, and processing of data in the technical solutions described in the embodiments of this application comply with relevant regulations.

[0111] In order to solve the above problems, the embodiments of this specification provide a detection method with multi-dimensional and multi-viewing angles. Figure 1 It is a schematic diagram of the steps of a multi-dimensional and multi-perspective detection method provided in the embodiment of this specification. This specification provides the method operation steps described in the embodiment or flowchart, but it may include more or fewer operation steps based on conventional or non-creative labor. The order of steps listed in the embodiment is only one way of executing the steps among many steps, and does not represent the only execution order. When the actual system or device product is executed, it can be executed in the order or in parallel according to the method shown in the embodiment or the accompanying drawings. Specifically, Figure 1 As shown, the method may include:

[0112] S101: Acquire all single transmitter response signals in the target detection area to obtain a total response data set of all transmitters;

[0113] S102: Decomposing the total transmitter response data set into a plurality of sub-data sets, and setting a corresponding detection viewing angle for each sub-data set, wherein each sub-data set includes one or more single transmitter response signals;

[0114] S103: Calculating the time shift of each single transmitter response signal in each sub-data set at a corresponding detection viewing angle according to a preset time domain shift rule;

[0115] S104: Shifting each single transmitter response signal in each sub-data set in the time domain according to the time shift size;

[0116] S105: superimposing the response signals of each single transmitting end in each sub-data set after translation to obtain a composite response signal of each sub-data set;

[0117] S106: Obtaining a multi-dimensional and multi-view total response data set of the target detection area according to the synthetic response signal of each sub-data set;

[0118] S107: Perform imaging processing using the multi-dimensional and multi-view total response data set.

[0119] Using the above-mentioned technical solution, the embodiments of this specification provide a multi-dimensional, multi-perspective detection method. When decomposing the total response dataset of the entire transmitter, a different detection perspective is set for each sub-dataset, so that each sub-dataset corresponds to an independent perspective dimension value. Time domain translation rules are used to simulate the time difference in signal propagation under different perspectives (such as wave velocity differences in geological exploration and acoustic path differences in medical ultrasound). This allows the translated signals to focus on the event axis of the target area when superimposed, enhancing the energy of the effective signal (such as reflected waves and echoes) and suppressing noise in the signal, thereby improving the signal-to-noise ratio of the target area data. A multi-dimensional, multi-perspective total dataset is generated based on the synthesized response signals of each sub-dataset. Because the multi-dimensional, multi-perspective total dataset contains information from different angle domains, it provides a more complete input for subsequent offset imaging, enabling the imaging process to utilize multi-angle constraints, thereby improving the imaging accuracy of complex media (such as geological layered structures and heterogeneous regions of human tissue).

[0120] The embodiments of this specification have no special dependence on signal acquisition equipment and detection systems, and can be applied to the fields of inter-well seismic exploration, VSP seismic exploration, tunnel seismic exploration, surface or underground passive source seismic exploration, three-dimensional land and offshore seismic exploration, tunnel seismic advance detection, urban seismic exploration, etc. The embodiments of this specification are also applicable to other types of detection methods based on wave propagation theory, including but not limited to geological radar, medical imaging (ultrasonic method or electromagnetic wave method) and various non-destructive testing methods in material structure health monitoring. The purpose of the embodiments of this specification is to generate a multi-dimensional and multi-perspective data set based on existing data without changing the traditional data acquisition method. This method can increase the coverage of data in any target area, thereby improving the signal-to-noise ratio of the data in that area. The technical principles in this specification are explained below using the geological exploration application scenario as an example.

[0121] In geological exploration, when a wave propagates downward from a source (shot point) and reaches an interface, according to Huygens' principle, each point on this interface can be considered a secondary point source. The wavefield excited by this secondary point source then propagates to each receiving point (detector point). Simply put, ignoring noise, the wavefield recorded by the receiving point is the sum of the wavefield from the source point and the wavefields of all secondary point sources on a specific interface in the subsurface. Because the wavepath from the source to the receiving point is reversible, multidimensional, multi-view datasets can be generated at both the source end (common-detector gathers) and the receiving end (common-shot gathers). The basic principles for obtaining multidimensional, multi-view datasets at the source and receiver ends are explained below with reference to Figures 2(a) and 2(b), respectively.

[0122] Figure 2(a) illustrates the basic concept of synthesizing a multi-dimensional, multi-perspective, wave-focused data set at the source. Assume that the detection system consists of three shots and one receiver. After each shot excites a seismic wave, the wavefield originates from the shot, reflects from a subsurface reflection point, and is received by the receiver, forming a common-receiver gather. A common-receiver gather is a collection of seismic wave records received at a particular receiver when all different shots are fired. Each data track in the gather corresponds to an independent shot excitation event and contains information such as the time, amplitude, and phase of the wavefield's arrival at the receiver, reflecting the wave propagation path from the shot to the receiver (including direct waves, reflected waves, and refracted waves). When multiple shots excite the same receiver, each shot corresponds to an independent "shot → reflection point → receiver" path, resulting in multiple coverage of the subsurface reflection points. By performing time-domain shifting (extending the shot to a virtual shot) and in-phase superposition (wave focusing) on the common-receiver gather, the signal at the target reflection point can be enhanced and random noise suppressed. Due to the reversibility of the wave path, the path of the wave from the shot point through the reflection point to the receiver point is exactly the same as the path of the wave from the receiver point through the reflection point to the shot point, only the propagation direction is opposite. Therefore, the wavefield recorded at the receiver point can be regarded as a "virtual source", and the forward propagation from "shot point excitation → reflection point → receiver point reception" is converted into the reverse propagation from "receiver point as virtual source excitation → reflection point → shot point reception". Because the underground interface has multiple secondary point sources, the original common receiver point gather contains the superposition information of multiple reflection points. By calculating the Green's function from the target reflection point back to the shot point, as shown by the dashed arrow in Figure 2(a), the record of each shot point is converted into reverse propagation data from "reflection point to shot point". At this point, the trace data of each shot point no longer contains the contributions of other reflection points and corresponds only to a single target reflection point, thus narrowing the detection target to a single reflection point and forming a focused gather, namely the new common receiver point gather. In this case, each trace in the gather corresponds to the reverse propagation from the same reflection point to a different shot point. The new common-receiver gathers are then time-shifted to extend the surface shot to the location of the underground virtual shot, as shown by the shaded triangle in Figure 2(a). Finally, all traces in the new common-receiver gathers are stacked in phase (i.e., wave focusing), achieving a three-fold coverage observation of the target reflection point. Through these steps, the common-receiver gathers are transformed from a mixed observation of multiple reflection points to a precise three-fold coverage observation of a single reflection point by utilizing the reversibility of the wave path to separate the target, extending the virtual shot to construct a multi-dimensional perspective, and stacking and accumulating enhanced signals, thereby improving the signal-to-noise ratio of the seismic data.

[0123] Similarly, Figure 2(b) illustrates the concept of synthesizing a multidimensional, multi-view dataset at the receiving end. Assuming the detection system consists of one shot and three receivers, the original common-shot gather contains three data traces, recording the mixed wavefield from shot excitation to the subsurface interface (including multiple reflection points) and then received at each receiver. From each trace in the common-shot gather, a Green's function containing only the upgoing wave at the target reflection point is extracted, as shown by the dashed arrows in Figure 2(b). This yields a new common-shot gather. In this case, each trace corresponds only to the upgoing propagation from the target reflection point to different receivers, narrowing the observation target to a single reflection point observed by the minimum shot moment receiver, as shown in Figure 2(b). The new common-shot gather is then extended in the time domain to move the ground receiver to the location of a virtual underground receiver, as shown in the shaded upper triangle in Figure 2(b). Each virtual receiver is located at a different receiving direction from the target reflection point. This is equivalent to creating three virtual receivers at different locations below the target reflection point, receiving the reflected wave from the target reflection point from multiple angles, thus forming a multidimensional observation of the reflection point. Finally, the amplitude of the target point reflected wave signal can be enhanced by superimposing the records of each receiving point in the new common shot point gather in phase, and three coverage observations can be achieved. From the above analysis, it can be seen that n shots and 1 detection point or 1 shot and n detection points can achieve n coverage observations. If it is n shots and n detection points, it can achieve n coverage observations. 2 In traditional exploration, coverage is linearly determined by the number of shot-detection pairs. However, the embodiments of this specification utilize multidimensional synthesis at the source and receiver ends (e.g., combined processing of common receiver gathers and common shot gathers) to generate multiple virtual observations (virtual shot points / virtual receiver points) for each physical shot-detection pair, further increasing the effective coverage. Without changing traditional data acquisition methods, the embodiments of this specification generate multidimensional, multi-view datasets based on existing data. This method increases the coverage of data in any target layer area, thereby improving the signal-to-noise ratio of the data in that area, achieving a new high-coverage, high-resolution exploration method.

[0124] The aforementioned approach to synthesizing multidimensional, multi-perspective datasets from either the source or receiver is based on the perspective of wave seismology. Based on classical principles of geometric seismology, multidimensional, multi-perspective wave-focused datasets can also be obtained. As shown in Figure 2(c), the approach to achieving multidimensional, multi-perspective wave-focused datasets based on geometric seismic principles is relatively simple and direct. All shot gathers are shifted in the time domain according to a regular time interval. For example, from the 1st to the Nth shot, each is shifted by 0-n time sampling intervals, as shown in Figure 2(d). These shot gathers are then superimposed to obtain a high-coverage observation shot gather after wave focusing in the target area.

[0125] In some embodiments of this specification, in step S102, the signal-to-noise ratio (SNR) of each single-shot record can be first evaluated, and then the original total shot gather can be divided into several sub-shot gathers based on the principle of uniform SNR distribution. The purpose of uniform SNR distribution is to avoid SNR polarization within the sub-shot gathers. If a sub-shot gather contains too many low-SNR signals, the noise suppression effect after superposition is poor, which becomes an imaging shortcoming. If a sub-shot gather contains only high-SNR signals, resources are wasted and the SNR cannot be further improved. Therefore, uniform SNR distribution ensures that the SNR improvement efficiency of each sub-dataset after superposition is consistent, and the resulting multi-dimensional, multi-view data has balanced reliability. First arrival picking or reflected wave phase axis tracking is performed on the single-shot records to obtain the signal peak amplitude at the target reflection point. A time window without valid signal is selected, and the noise standard deviation of the data within this time window is calculated. The SNR of each single-shot record can be calculated based on the signal amplitude and noise standard deviation. All single-shot records are sorted in descending order of signal-to-noise ratio (SNR), and the number of sub-shot sets is determined. The single-shot records are then allocated to the sub-shot sets according to a round-robin allocation strategy. For example, starting with channel 1 (highest SNR), channel 1 is allocated to sub-shot set 1, channel 2 to sub-shot set 2, and so on. The Kth channel is then allocated to sub-shot set K. The K+1th channel is then returned to sub-shot set 1, the K+2th channel to sub-shot set 2, and so on. This ensures that each sub-shot set gradually includes high, medium, and low SNR signals. In other embodiments, single-shot records can be allocated to sub-shot sets according to an allocation strategy such as stratified sampling.

[0126] When sub-shot collections are time-domain shifted, different time interval rules will produce high-coverage observation data for target areas of different sizes and azimuth depths. The time interval rules based on which the original shot collections are time-domain shifted are described in detail below. In some embodiments of this specification, the time domain shift rules include a first shift rule and a second shift rule. Figure 3 As shown, the time shift of each single transmitter response signal in each sub-data set at the corresponding detection viewing angle is calculated according to a preset time domain shift rule, including:

[0127] S301: Determine the size of the target detection area;

[0128] S302: Determine whether the size of the target detection area is greater than a preset threshold;

[0129] S303: If yes, then calculating the time shift of each single transmitter response signal in each sub-data set at the corresponding detection viewing angle according to the first shift rule;

[0130] S304: If not, then calculating the time shift of each single-transmitter response signal in each sub-data set at the corresponding detection viewing angle according to the second shift rule.

[0131] It can be understood that the size of the target area has different effects on wave propagation, resulting in different levels of difficulty in distinguishing it in the imaging profile. When the target area is large (significantly larger than the wavelength), it can generally be accurately identified. However, due to the influence of external environments (such as topography and urban agglomerations), the number of coverage times for larger target areas is limited, and the effectiveness of multiple coverage techniques in improving the signal-to-noise ratio is affected to varying degrees, which in turn affects the accuracy of subsequent migration imaging. When the target area is small (less than or equal to the wavelength), or even smaller than the vertical and lateral resolution limits, the reflected energy of the wave passing through such geological bodies is very weak, making it difficult to accurately identify such geological bodies directly from the imaging results. Therefore, the embodiments of this specification design different detection methods for different sizes of exploration target areas. When the exploration target area is large and multiple coverage observations are difficult to achieve, a method of forming multi-dimensional, multi-viewpoint shot gathers in the target area can be used to improve the signal-to-noise ratio of the seismic data in the target area. When the exploration target geological body is small, a multi-dimensional, multi-viewpoint focused wave gather can be formed in the target area to increase the reflected energy.

[0132] As an embodiment of this specification, a target detection area is determined in the exploration area, and all single shot records in the detection area are obtained. Assuming that the total number of single shot records is n, the i-th shot record can be recorded as , where t represents the travel time, and the original total gun set is Next, the original total shot set is decomposed into N sub-shot sets, each of which corresponds to a detection angle, and each sub-shot set contains M single-shot records. Then the original total shot set can be recorded as , where the j-th sub-shot set can be expressed as . General and Cannon Collection Each single shot record is translated in the time dimension according to certain rules. Assuming that the j-th sub-shot set The time shift of the mth shot is If the observed target area is larger than a preset threshold, the time shift of each single transmitter response signal in each sub-dataset at the corresponding detection viewing angle is calculated according to the first shift rule, thereby obtaining a multi-view recording. In some embodiments of this specification, the preset threshold is determined based on the wavelength of the seismic wave excited by the shot point. In other embodiments, it can also be determined based on, for example, the offset coverage range.

[0133] As an embodiment of the present specification, calculating the time shift of each single transmitter response signal in each sub-data set at a corresponding detection viewing angle according to the first shift rule includes:

[0134] Calculating the head-to-tail transmitter distance of each sub-data set;

[0135] Determining a detection viewing angle direction according to the detection viewing angle of each sub-data set;

[0136] The time shift of the response signal of each single transmitter in each sub-data set is calculated based on the head-to-tail transmitter distance, detection viewing angle, detection viewing direction and seismic wave velocity of each sub-data set.

[0137] As an embodiment of this specification, the first translation rule may be expressed using the following formula:

[0138] (1)

[0139] in, Indicates the time shift of the response signal of the mth single transmitter in the jth sub-data set, represents the distance between the first and last transmitters in the jth sub-dataset, Represents the detection viewing angle of the jth sub-dataset, when When , it means that the detection angle is in the direction of increasing offset; when When , it means that the detection viewing angle is in the direction of decreasing offset, and c represents the seismic wave velocity.

[0140] As shown in Figure 4(a), s represents the sub-shot set The distance between the 1st shot and the Mth shot, the z-axis represents the depth direction, the x-axis represents the survey line direction, the right is positive, that is, the direction of increasing offset, the angle between the detection direction and the z-axis is , whose spatial orientation can be decomposed into vertical components: , that is, the distance along the z-axis, and the horizontal component: , that is, the offset direction. hour, , x is positive, corresponding to the increase of right offset, when hour, , x is negative, corresponding to a decrease in the right offset. It reflects the difference in wavefront arrival time caused by the position differences of different shot points. By applying different time delays or advances to each single shot record, the seismic waves from the target detection perspective can be superimposed in phase during synthesis, thus forming a focused beam.

[0141] When the observed target area is smaller than a preset threshold, it indicates a small area. Since the reflected energy when the wave passes through such an area is very weak, it is difficult to accurately identify such a target directly from the imaging results. Therefore, the reflected energy can be increased by generating regional data of multi-dimensional, multi-perspective focused waves in the target area. In some embodiments of this specification, when the observed target area is smaller than a preset threshold, the time shift of each single transmitter response signal in each sub-data set at the corresponding detection angle is calculated according to a second shift rule, thereby generating a multi-perspective record.

[0142] As an example of this specification, refer to Figure 5 As shown, the calculating of the time shift of each single transmitter response signal in each sub-data set at the corresponding detection viewing angle according to the second shift rule includes:

[0143] S501: Determine a wave focusing point according to the spatial position and depth of the target detection area.

[0144] It's understandable that in seismic exploration or array signal processing, by time-delaying or time-leading signals from different transmitters, the signals are superimposed in phase within the target area, thereby enhancing the signal at a specific detection angle and suppressing noise from other angles. The convergence point can be understood as a focal point within the target detection area, where signals from all transmitters must be superimposed in phase to enhance the signal. As shown in Figure 4(b), the convergence point location P can be determined based on the spatial location of the detection area, such as the horizontal coordinates of the geological structure, and the depth, such as the burial depth of the target layer.

[0145] S502: Calculate the distance from each transmitting end in each sub-data set to the wave focusing point.

[0146] In some embodiments of this specification, the following formula is used to calculate the distance from each transmitting end in each sub-data set to the wave focusing point:

[0147] (2)

[0148] (3)

[0149] (4)

[0150] in, represents the distance from the first transmitter to the wave focusing point in the jth sub-data set, represents the distance from the tail emission end to the wave focusing point in the jth sub-data set, represents the distance from the i-th transmitter to the wave focusing point in the j-th sub-data set, represents the detection viewing angle of the j-th sub-dataset, represents the first distance from the central transmitter to the wave focusing point in the jth sub-dataset, represents the second distance from the first transmitter to the center transmitter in the jth sub-data set, Represents the third distance from the first transmitter to the i-th transmitter in the j-th sub-data set.

[0151] S503: Determine a detection viewing angle direction according to the detection viewing angle of each sub-data set.

[0152] Detection viewing angle Determines the horizontal direction of the beam, and by adjusting the time shift, the signal energy is focused in this direction. When , it means that the detection angle is in the direction of increasing offset. , it indicates that the detection angle of view is in the direction of decreasing offset.

[0153] S504: Calculate the time shift of the response signal of each single transmitter in each sub-data set according to the distance from each transmitter in each sub-data set to the wave focusing point, the detection viewing angle direction and the seismic wave velocity.

[0154] In some embodiments of this specification, the second translation rule may be expressed using the following formula:

[0155] (5)

[0156] in, Indicates the time shift of the response signal of the mth single transmitter in the jth sub-data set, represents the distance from the first transmitter to the wave focusing point in the jth sub-data set, represents the distance from the tail emission end to the wave focusing point in the jth sub-data set, represents the distance from the i-th transmitter to the wave focusing point in the j-th sub-data set, Represents the detection viewing angle of the jth sub-dataset, when When , it means that the detection angle is in the direction of increasing offset; when When , the detection angle of view is toward decreasing offset, and c represents the seismic wave velocity. Through the above steps, calculating the time shift essentially compensates for the path differences of the signals from each transmitter, allowing the signals to superimpose in phase at the focal point, forming a beam directed toward a specific angle of view, thereby enhancing signal energy in the target area.

[0157] After calculating the time shift of each single transmitter response signal in each sub-data set at the corresponding detection viewing angle, the single transmitter response signal in each sub-data set is shifted in the time domain according to the time shift to obtain the shifted single shot record. , then the single shot records after translation in each sub-dataset are superimposed to obtain the synthetic response signal of each sub-dataset, that is:

[0158] (6)

[0159] Then the total multi-dimensional and multi-view shot collection in the target detection area can be recorded as Finally, the multi-dimensional, multi-view total shot gather can be used for migration or inversion imaging. Thus, by controlling the discrete time intervals of each shot recorded in the sub-shot gather, the embodiments of this specification can radiate seismic wave energy within any detection viewing angle. When necessary, the waves can be directly focused to a specific depth, thereby enhancing the reflected energy within a specific area or small range, thereby achieving improved coverage and resolution.

[0160] The following takes the single-layer horizontal interface geological model as an example to illustrate the specific implementation steps and application effects. Figure 6 As shown, the layer interface is located at , set a total of 7 shots, the shots are arranged at equal intervals, the shot spacing is 60m, and the center shot is located at , the detection points are arranged at , the channel spacing is 30m, and a total of 81 channels are received. The longitudinal wave velocity of the upper medium is , the density is ; The longitudinal wave velocity of the lower medium is , the density is . After the numerical simulation, noise is introduced. Figure 7 The center shot record shows a low signal-to-noise ratio. Figure 7 To simplify the description of the problem, the total shot gather and sub-shot gather in this embodiment are equivalent, that is, This embodiment only presents the multi-dimensional and multi-view focused wave shot gathering formation process of the horizontal layer interface medium. This embodiment is also applicable to more complex underground media.

[0161] Figures 8(a)-8(d) show the realization of seismic beam focusing at a vertical position of 700 m from the central shot point. At this point, the time shift of each shot can be calculated according to formulas (2)-(5), as shown in Figure 8(a). Then, each single shot record is shifted according to the corresponding time shift and superimposed to obtain a high-coverage, high-resolution seismic record of the detection area, as shown in Figure 8(b). Figure 8(c) shows that the reflected wave energy in the target area is significantly enhanced. Finally, the high-coverage, high-resolution seismic record of the detection area is subjected to reverse time migration imaging, resulting in the reverse time migration imaging result shown in Figure 8(d).

[0162] Figures 9(a) to 9(d) show the realization of wave focusing at 808.3 m after the seismic beam is deflected 30° in the direction of increasing offset. =30°, , According to formulas (2)-(5), the time shift of each shot can be calculated, as shown in Figure 9(a). Then, the single shot records are shifted according to the corresponding time shift and superimposed to obtain the high-coverage and high-resolution seismic records of the detection area as shown in Figure 9(b). As can be seen from Figure 9(c), the reflected wave energy in the target area is significantly enhanced. Finally, the high-coverage and high-resolution seismic records of the detection area are subjected to reverse time migration imaging, and the reverse time migration imaging results are obtained as shown in Figure 9(d).

[0163] Figures 10(a) to 10(d) show the realization of wave focusing at 808.3 m after the seismic beam is deflected 30° in the direction of decreasing offset. =-30°, , According to formulas (2)-(5), the time shift size of each shot can be calculated, as shown in Figure 10 (a). Then, each single shot record is shifted according to the corresponding time shift size and superimposed to obtain a high-coverage, high-resolution seismic record of the detection area as shown in Figure 10 (b). From Figure 10 (c), it can be seen that the reflected wave energy of the target area is significantly enhanced. Finally, the high-coverage, high-resolution seismic record of the detection area is subjected to reverse time migration imaging processing to obtain the reverse time migration imaging result diagram shown in Figure 10 (d). Therefore, the high-coverage, high-resolution seismic record obtained in the embodiment of this specification is used to perform migration imaging processing. Since multiple high-coverage, high-resolution seismic records contain information in different angle domains, the imaging process can utilize multi-angle constraints, thereby improving the imaging accuracy of complex media.

[0164] Based on the above-mentioned detection method with multiple dimensions and multiple perspectives, the embodiments of this specification also provide a corresponding detection device with multiple dimensions and multiple perspectives. The device may include a system (including a distributed system), software (application), module, component, server, client, etc. that uses the method described in the embodiments of this specification and is combined with the necessary implementation hardware. Based on the same innovative concept, the devices in one or more embodiments provided in the embodiments of this specification are as described in the following embodiments. Since the implementation scheme of the device to solve the problem is similar to the method, the implementation of the specific device in the embodiments of this specification can refer to the implementation of the aforementioned method, and the repeated parts will not be repeated. As used below, the term "unit" or "module" can be a combination of software and / or hardware that implements the predetermined function. Although the device described in the following embodiments is preferably implemented in software, the implementation of hardware, or a combination of software and hardware, is also possible and conceived.

[0165] Specifically, Figure 11 This is a schematic diagram of the module structure of an embodiment of a multi-dimensional and multi-viewing detection device provided in the embodiment of this specification, with reference to Figure 11As shown, the embodiment of this specification provides a multi-dimensional and multi-viewing detection device including:

[0166] An acquisition module 1101 is configured to acquire response signals from all single transmitters in a target detection area to obtain a total response data set of all transmitters;

[0167] a decomposition module 1102 configured to decompose the total transmitter response data set into a plurality of sub-data sets and set a corresponding detection viewing angle for each sub-data set, wherein each sub-data set includes one or more single transmitter response signals;

[0168] A calculation module 1103 is configured to calculate the time shift of each single transmitter response signal in each sub-data set at a corresponding detection viewing angle according to a preset time domain shift rule;

[0169] A shift module 1104 is configured to shift each single transmitter response signal in each sub-data set in the time domain according to the time shift size;

[0170] A superposition module 1105 is configured to superimpose the response signals of each single transmitting end in each sub-data set after translation to obtain a composite response signal of each sub-data set;

[0171] a total data set determining module 1106, configured to obtain a multi-dimensional and multi-view total response data set of the target detection area according to the synthetic response signal of each sub-data set;

[0172] The processing module 1107 is configured to perform imaging processing using the multi-dimensional and multi-view total response data set.

[0173] The beneficial effects achieved by the device provided in the embodiments of this specification are consistent with the beneficial effects achieved by the above-mentioned method and will not be repeated here.

[0174] Reference Figure 12As shown, based on the above-described multi-dimensional and multi-perspective detection method, one embodiment of this specification further provides a computer device 1202, wherein the above-described method is executed on the computer device 1202. The computer device 1202 may include one or more processors 1204, such as one or more central processing units (CPUs), each of which may implement one or more hardware threads. The computer device 1202 may also include any memory 1206 for storing any type of information, such as code, settings, data, etc. For example, without limitation, the memory 1206 may include any one or more combinations of the following: any type of RAM, any type of ROM, a flash memory device, a hard disk, an optical disk, etc. More generally, any memory may use any technology to store information. Furthermore, any memory may provide volatile or non-volatile retention of information. Furthermore, any memory may represent a fixed or removable component of the computer device 1202. In one embodiment, when the processor 1204 executes associated instructions stored in any memory or combination of memories, the computer device 1202 may perform any operation of the associated instructions. The computer device 1202 also includes one or more drive mechanisms 1208 for interacting with any storage, such as a hard disk drive mechanism, an optical disk drive mechanism, and the like.

[0175] Computer device 1202 may also include an input / output module 1210 (I / O) for receiving various inputs (via input devices 1212) and providing various outputs (via output devices 1214). A specific output mechanism may include a presentation device 1216 and an associated graphical user interface (GUI) 1218. In other embodiments, input / output module 1210 (I / O), input devices 1212, and output devices 1214 may not be included, and the computer device 1202 may simply function as a single computer device on a network. Computer device 1202 may also include one or more network interfaces 1220 for exchanging data with other devices via one or more communication links 1222. One or more communication buses 1224 couple the components described above.

[0176] The communication link 1222 may be implemented in any manner, for example, via a local area network, a wide area network (e.g., the Internet), a point-to-point connection, etc., or any combination thereof. The communication link 1222 may include any combination of hardwired links, wireless links, routers, gateway functions, name servers, etc., governed by any protocol or combination of protocols.

[0177] Corresponding to Figure 1 、 Figure 3 and Figure 5In addition to the method shown, an embodiment of this specification also provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the steps of the above method are executed.

[0178] The embodiment of this specification also provides a computer-readable instruction, wherein when the processor executes the instruction, the program therein causes the processor to execute the following Figure 1 、 Figure 3 and Figure 5 method.

[0179] The embodiment of this specification also provides a computer program product, including at least one instruction or at least one program, which is loaded and executed by a processor to implement the following Figure 1 、 Figure 3 and Figure 5 method.

[0180] It should be understood that in the various embodiments of this specification, the size of the serial numbers of the above-mentioned processes does not mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of this specification.

[0181] It should also be understood that in the embodiments of this specification, the term "and / or" is merely a description of the relationship between associated objects, indicating that three possible relationships exist. For example, "A and / or B" can represent three situations: A exists alone, A and B exist simultaneously, and B exists alone. Furthermore, the character " / " in this specification generally indicates that the associated objects are in an "or" relationship.

[0182] Those skilled in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed in this specification can be implemented in electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the above description has generally described the composition and steps of each example according to function. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of this specification.

[0183] Those skilled in the art will clearly understand that, for the convenience and brevity of description, the specific working processes of the systems, devices and units described above can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0184] In the several embodiments provided in this specification, it should be understood that the disclosed systems, devices and methods can be implemented in other ways. For example, the device embodiments described above are merely illustrative. For example, the division of the units is merely a logical function division. In actual implementation, there may be other division methods, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the mutual coupling or direct coupling or communication connection shown or discussed can be an indirect coupling or communication connection through some interfaces, devices or units, or can be an electrical, mechanical or other form of connection.

[0185] The units described as separate components may or may not be physically separate, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed across multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the embodiments of this specification.

[0186] In addition, the functional units in the various embodiments of this specification may be integrated into a single processing unit, or each unit may exist physically separately, or two or more units may be integrated into a single unit. The aforementioned integrated units may be implemented in the form of hardware or software functional units.

[0187] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this specification, or the part that contributes to the existing technology, or all or part of the technical solution can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes a number of instructions for enabling a computer device (which can be a personal computer, server, or network device, etc.) to execute all or part of the steps of the method described in each embodiment of this specification. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk, and other media that can store program code.

[0188] Specific embodiments are used in this specification to illustrate the principles and implementation methods of this specification. The description of the above embodiments is only used to help understand the methods and core ideas of this specification. At the same time, for those skilled in the art, based on the ideas of this specification, there will be changes in the specific implementation methods and application scope. In summary, the contents of this specification should not be understood as limiting this specification.

Claims

1. A multi-dimensional and multi-view detection method, characterized in that: The method comprises: Obtain all single transmitter response signals in the target detection area to obtain the total response data set of all transmitters; Decomposing the total transmitter response data set into a plurality of sub-data sets, and setting a corresponding detection viewing angle for each sub-data set, wherein each sub-data set includes one or more single transmitter response signals; Calculate the time shift of each single transmitter response signal in each sub-data set at the corresponding detection viewing angle according to a preset time domain shift rule; Shifting each single transmitter response signal in each sub-data set in the time domain according to the time shift size; Superimposing the response signals of each single transmitter in each sub-data set after translation to obtain a synthetic response signal of each sub-data set; Obtaining a multi-dimensional and multi-view total response data set of the target detection area according to the synthetic response signal of each sub-data set; Imaging processing is performed using the multi-dimensional and multi-view total response data set.

2. The method according to claim 1, characterized in that The time domain translation rule includes a first translation rule and a second translation rule; The calculating, according to a preset time domain shift rule, the time shift size of each single transmitting end response signal in each sub-data set at the corresponding detection viewing angle includes: determining a size of the target detection area; Determining whether the size of the target detection area is greater than a preset threshold; If yes, then calculating the time shift of each single transmitting end response signal in each sub-data set at the corresponding detection viewing angle according to the first shift rule; If not, the time shift of each single transmitting end response signal in each sub-data set at the corresponding detection viewing angle is calculated according to the second shift rule.

3. The method according to claim 2, characterized in that The calculating, according to the first translation rule, the time translation size of each single transmitting end response signal in each sub-data set at the corresponding detection viewing angle includes: Calculating the head-to-tail transmitter distance of each sub-data set; Determining a detection viewing angle direction according to the detection viewing angle of each sub-data set; The time shift of the response signal of each single transmitter in each sub-data set is calculated based on the head-to-tail transmitter distance, detection viewing angle, detection viewing direction and seismic wave velocity of each sub-data set.

4. The method according to claim 3, characterized in that The calculating of the time shift of the response signal of each single transmitter in each sub-data set according to the head-to-tail transmitter spacing, the detection viewing angle, the detection viewing direction, and the seismic wave velocity of each sub-data set includes: The time shift of the response signal of each single transmitter in each sub-data set is calculated using the following formula: ； in, Indicates the time shift of the response signal of the mth single transmitter in the jth sub-data set, represents the distance between the first and last transmitters in the jth sub-dataset, Represents the detection viewing angle of the jth sub-dataset, when When , it means that the detection angle is in the direction of increasing offset; when When , it means that the detection viewing angle is in the direction of decreasing offset, and c represents the seismic wave velocity.

5. The method according to claim 2, characterized in that The calculating, according to the second translation rule, the time translation size of each single transmitting end response signal in each sub-data set at the corresponding detection viewing angle includes: Determine the wave focusing point based on the spatial position and depth of the target detection area; Calculating the distance from each transmitting end in each sub-data set to the wave focusing point; Determining a detection viewing angle direction according to the detection viewing angle of each sub-data set; The time shift of the response signal of each single transmitter in each sub-data set is calculated based on the distance from each transmitter in each sub-data set to the wave focusing point, the detection viewing angle direction and the seismic wave velocity.

6. The method according to claim 5, characterized in that The calculating of the time shift of the response signal of each single transmitter in each sub-data set according to the distance from each transmitter in each sub-data set to the wave focusing point, the detection viewing angle direction, and the seismic wave velocity includes: The time shift of the response signal of each single transmitter in each sub-data set is calculated using the following formula: ； in, Indicates the time shift of the response signal of the mth single transmitter in the jth sub-data set, represents the distance from the first transmitter to the wave focusing point in the jth sub-data set, represents the distance from the tail emission end to the wave focusing point in the jth sub-data set, represents the distance from the i-th transmitter to the wave focusing point in the j-th sub-data set, Represents the detection viewing angle of the jth sub-dataset, when When , it means that the detection angle is in the direction of increasing offset; when When , it means that the detection viewing angle is in the direction of decreasing offset, and c represents the seismic wave velocity.

7. The method according to claim 5, characterized in that The calculating the distance from each transmitting end in each sub-data set to the wave focusing point includes: Calculate the first distance from the central transmitting end to the wave focusing point in each sub-data set; Calculating a second distance from the first transmitter to the central transmitter and a third distance from the first transmitter to each transmitter in each sub-data set; The distance from each transmitting end in each sub-data set to the wave focusing point is calculated according to the first distance, the second distance, and the third distance.

8. The method according to claim 7, characterized in that The calculating, based on the first distance, the second distance, and the third distance, the distance from each transmitting end in each sub-data set to the wave focusing point includes: The distance from each transmitter to the wave focusing point in each sub-data set is calculated using the following formula: ； ； ； in, represents the distance from the first transmitter to the wave focusing point in the jth sub-data set, represents the distance from the tail emission end to the wave focusing point in the jth sub-data set, represents the distance from the i-th transmitter to the wave focusing point in the j-th sub-data set, represents the detection viewing angle of the j-th sub-dataset, represents the first distance from the central transmitter to the wave focusing point in the jth sub-dataset, represents the second distance from the first transmitter to the center transmitter in the jth sub-data set, Represents the third distance from the first transmitter to the i-th transmitter in the j-th sub-data set.

9. A multi-dimensional and multi-viewing detection device, characterized in that: The device comprises: An acquisition module is used to acquire the response signals of all single transmitters in the target detection area and obtain a total response data set of all transmitters; a decomposition module, configured to decompose the total transmitter response data set into a plurality of sub-data sets and set a corresponding detection viewing angle for each sub-data set, wherein each sub-data set includes one or more single transmitter response signals; A calculation module, configured to calculate the time shift of each single transmitter response signal in each sub-data set at a corresponding detection viewing angle according to a preset time domain shift rule; A translation module, configured to translate each single transmitter response signal in each sub-data set in the time domain according to the time translation size; A superposition module is used to superimpose the response signals of each single transmitter in each sub-data set after translation to obtain a synthetic response signal of each sub-data set; a total data set determination module, configured to obtain a multi-dimensional and multi-view total response data set of the target detection area according to the synthetic response signal of each sub-data set; A processing module is used to perform imaging processing using the multi-dimensional and multi-view total response data set.

10. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the computer program, the method according to any one of claims 1 to 8 is implemented.

11. A computer-readable storage medium, characterized in that The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 8 is implemented.

12. A computer program product, characterized in that The method comprises at least one instruction or at least one program, wherein the at least one instruction or the at least one program is loaded and executed by a processor to implement the method according to any one of claims 1 to 8.

Images