A multi-shot transient surface wave exploration method, device and storage medium

Abstract

The application provides a multi-excitation-point transient surface wave exploration method, equipment and a storage medium. The method comprises the following steps: exciting multiple seismic sources one by one, collecting multiple sets of common-source measurement data by a detector and performing screening, selecting common-source measurement data of the detector within a set offset distance range, taking one of the screened detectors as a central channel, taking the detectors within a set adjacent range as adjacent channels, and obtaining surface wave dispersion spectrum by using the common-source measurement data of the central channel and the adjacent channels to obtain the surface wave dispersion spectrum of the central channel; in this way, the surface wave dispersion spectrum of multiple detectors at different seismic sources is obtained in sequence, and the surface wave dispersion spectrum is superimposed; the underground transverse wave velocity structure of each detector is obtained based on the surface wave dispersion spectrum in sequence, so that the underground transverse wave velocity profile in the measurement range is obtained to complete the exploration. In addition, a multi-excitation-point transient surface wave exploration equipment and a storage medium are also provided. The imaging effect of the surface wave dispersion spectrum obtained by the method is greatly improved, and the exploration result is more accurate.

Description

Technical Field

[0001] This invention belongs to the field of geophysical exploration technology. More specifically, it relates to a method, equipment, and storage medium for multi-excitation point transient surface wave exploration. Background Technology

[0002] Active surface wave exploration (APW) is a seismic exploration method that uses artificial seismic sources such as hammering, falling weights, and explosives to generate surface waves of a certain frequency range on the Earth's surface, mainly Rayleigh-type surface waves. The shallow shear wave velocity is estimated by inverting the dispersion characteristics of the surface waves received by the geophone, thereby obtaining the subsurface shear wave velocity structure. This method has the advantages of high resolution, strong layering ability, convenient construction, and no influence from high-velocity layers. It is widely used in shallow surface exploration, especially in complex engineering geological exploration, and has a very broad application prospect.

[0003] However, conventional transient surface wave methods typically use single-sided or double-sided firing to excite seismic source signals. This often results in stronger energy in seismic traces closer to the source and weaker energy in traces farther away. This affects the overall dispersion spectrum imaging quality when calculating the dispersion spectrum, thus impacting the accuracy of dispersion curve acquisition and inversion. Furthermore, in terms of data processing, conventional transient surface wave exploration methods only calculate one inversion result per arrangement, representing the velocity structure beneath the center of that arrangement. A two-dimensional shear wave velocity profile is then obtained through lateral interpolation of inversion results from multiple arrangements. Since the inversion result of a single arrangement represents the comprehensive reflection of the entire subsurface strata, the final result has a low lateral resolution, making it difficult to meet the requirements of refined exploration in complex environments. Summary of the Invention

[0004] To address at least one of the deficiencies or improvement needs of the prior art, this invention provides a multi-excitation point transient surface wave exploration method, equipment, and storage medium, which can be used for exploration in complex environments, is simple and efficient to implement, and provides high resolution exploration results.

[0005] To achieve the above objectives, according to the first aspect of the present invention, the present invention adopts the following technical solution.

[0006] A multi-excitation-point transient surface wave exploration method, the method comprising:

[0007] Multiple seismic sources are excited one by one, and the excitation of each seismic source generates a set of seismic source signals; the detector collects each of the seismic source signals one by one to generate multiple sets of resonant source measurement data;

[0008] Each set of resonance source measurement data is filtered, wherein the filtering is performed on a portion of the resonance source measurement data whose offset meets the set conditions, and the portion of the resonance source measurement data comes from at least one of the corresponding detectors;

[0009] One of the detectors corresponding to the partial resonance source measurement data is selected as the center channel, and the remaining detectors within a set range adjacent to the center channel are selected as adjacent channels; the surface wave dispersion spectrum of the center channel is calculated using the resonance source measurement data of the center channel and the adjacent channels.

[0010] In this way, multiple detectors generate multi-channel surface wave dispersion spectra, with each source corresponding to one such multi-channel surface wave dispersion spectrum;

[0011] Determine whether a detector obtains at least two center channel surface wave dispersion spectra; if so, calculate by superimposing all the center channel surface wave dispersion spectra obtained by the detector to obtain the surface wave dispersion spectrum of the detector; if not, use the obtained center channel surface wave dispersion spectrum as the surface wave dispersion spectrum of the detector.

[0012] In this manner, the surface wave dispersion spectrum of all the detectors is obtained sequentially;

[0013] Based on the surface wave dispersion spectrum, the subsurface shear wave velocity structure of each of the detectors is obtained sequentially, thereby obtaining the subsurface shear wave velocity profile within the measurement range to complete the exploration.

[0014] Furthermore, the calculation method for the dispersion spectrum of the surface wave of the center channel is as follows: the phase tilt superposition calculation is performed on the resonant source measurement data generated by one source of the center channel and the adjacent channel.

[0015] The resonance source measurement data includes vibration measurement data and the coordinate position of the detector.

[0016] Furthermore, the formula for calculating the surface wave dispersion spectrum of the center channel detector is as follows:

[0017]

[0018] Among them, A i (f) represents the normalized amplitude spectrum of the vibration measurement data within the set range; d c d represents the offset distance of the central channel; i The offset distance between adjacent lanes is represented by L; the radius of the set range is represented by D1 and D2, which both represent the range values ​​of the offset distance of the center lane that meet the set conditions.

[0019] Furthermore, the resonance source measurement data is obtained by acquiring the vibration source signal using the detector and preprocessing the vibration source signal;

[0020] The preprocessing includes noise suppression, bad sector processing, filtering, and gain compensation.

[0021] Furthermore, the detector is a distributed cableless seismic detector; the detectors are arranged in two ways: with equal spacing and with unequal spacing.

[0022] Furthermore, the underground shear wave velocity structure is obtained by inverting the dispersion curve of each of the central channel detectors based on the surface wave dispersion spectrum.

[0023] The inversion includes the least squares method, simulated annealing method, and genetic algorithm.

[0024] Furthermore, the underground shear wave velocity profile is obtained by performing lateral two-dimensional interpolation on the underground shear wave velocity structure of each of the central channel geophones.

[0025] Furthermore, when the central channel is located at the beginning or end of the arrangement, the surface wave dispersion spectrum of the central channel is calculated after extending the arrangement by adding a detector based on the set range.

[0026] According to a second aspect of the present invention, a multi-excitation point transient surface wave exploration device is also provided, comprising at least one processing unit and at least one storage unit, wherein the storage unit stores a computer program that, when executed by the processing unit, causes the processing unit to perform the steps of any of the methods described above.

[0027] According to a third aspect of the invention, a storage medium is also provided that stores a computer program executable by an access authentication device, which, when run on the access authentication device, causes the access authentication device to perform the steps of any of the methods described above.

[0028] In summary, compared with the prior art, the above-described technical solutions conceived by this invention can achieve the following beneficial effects:

[0029] (1) The present invention provides a multi-excitation point transient surface wave exploration method, equipment, and storage medium. When exciting the seismic source signal, multiple seismic sources are used. According to the actual site conditions, multiple seismic source points are set up along the survey line and excited sequentially to collect multiple sets of co-source measurement data. By screening each set of co-source measurement data, a geophone that meets the set conditions is obtained. Then, using this geophone as the center channel, the surface wave dispersion spectrum of the center channel is calculated with the co-source measurement data of adjacent channels. The present invention selects a geophone with moderate source energy and does not need to calculate all geophones sequentially. This is more conducive to high-quality imaging of the surface wave dispersion spectrum and also improves the exploration efficiency to a certain extent.

[0030] (2) The present invention provides a multi-excitation point transient surface wave exploration method, equipment, and storage medium. It calculates the surface wave dispersion spectrum of the central channel by superimposing phase tilt data obtained from the same source excitation of the central channel and adjacent channels. Furthermore, if a detector obtains at least two central channel surface wave dispersion spectra based on different sources, the superimposed central channel surface wave dispersion spectra are used to obtain a more accurate surface wave dispersion spectrum for that detector. Based on this, the subsurface shear wave velocity structure is obtained by inverting channel by channel. Because channel data with moderate energy is selected for segmented surface wave dispersion spectrum calculation, and the central channel surface wave dispersion spectra are superimposed, the imaging effect of the final surface wave dispersion spectrum is greatly improved, thereby increasing the picking accuracy of the dispersion curve and the accuracy of the inversion. This results in a more accurate subsurface shear wave velocity profile within the final measurement range, meeting the requirements for refined exploration in complex environments. Attached Figure Description

[0031] To more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings used in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0032] Figure 1 This is a schematic diagram of the process flow of a multi-excitation point transient surface wave exploration method, equipment, and storage medium according to the present invention.

[0033] Figure 2 This is a schematic diagram of the detector arrangement and data processing of a multi-excitation point transient surface wave exploration method, equipment, and storage medium according to the present invention. Detailed Implementation

[0034] To make the objectives, technical solutions, and advantages of this invention clearer, the invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the invention. Furthermore, the technical features involved in the various embodiments of this invention described below can be combined with each other as long as they do not conflict with each other.

[0035] Figure 1 This is a schematic diagram of the process flow for a multi-excitation-point transient surface wave exploration method, equipment, and storage medium provided by the present invention. This invention significantly improves the imaging effect of the surface wave dispersion spectrum, making the subsurface shear wave velocity profile within the measurement range more accurate and the results more precise, meeting the requirements for refined exploration in complex environments. The method includes:

[0036] Multiple seismic sources are excited one by one, and the excitation of each seismic source generates a set of seismic source signals; the detector collects each seismic source signal one by one to generate multiple sets of resonant source measurement data;

[0037] Each set of resonance source measurement data is filtered out as follows: partial resonance source measurement data whose offset meets the set conditions, and partial resonance source measurement data comes from at least one corresponding detector.

[0038] One of the detectors corresponding to a portion of the resonance source measurement data is selected as the center channel, and the remaining detectors within a set range adjacent to the center channel are selected as adjacent channels; the surface wave dispersion spectrum of the center channel is calculated using the resonance source measurement data of the center channel and the adjacent channels.

[0039] In this way, multiple detectors generate multi-channel surface wave dispersion spectra, with each source corresponding to a multi-channel surface wave dispersion spectrum;

[0040] Determine whether a detector obtains at least two center channel surface wave dispersion spectra; if so, superimpose and calculate all the center channel surface wave dispersion spectra obtained by the above detector to obtain the surface wave dispersion spectrum of the detector; if not, use the obtained center channel surface wave dispersion spectrum as the surface wave dispersion spectrum of the detector.

[0041] In this manner, the surface wave dispersion spectrum of all center channel detectors is obtained sequentially;

[0042] The subsurface shear wave velocity structure of each detector is obtained sequentially based on the surface wave dispersion spectrum, thereby obtaining the subsurface shear wave velocity profile within the measurement range to complete the exploration.

[0043] First, multiple seismic sources are excited one by one, and the excitation of each source generates a set of seismic source signals. A geophone then collects each seismic source signal to generate multiple sets of resonance source measurement data. Each resonance source measurement data set includes vibration measurement data and the coordinate position of the geophone.

[0044] Specifically, based on the actual scenario conditions, multiple seismic sources are deployed at multiple locations along the survey line and excited sequentially. Multiple sets of resonant source measurement data are obtained using the vibration source signals generated by the detector seismic sources and the coordinate positions of the detectors. The vibration measurement data are then obtained after preprocessing the vibration source signals. Preprocessing includes, but is not limited to, noise suppression, bad channel processing, filtering, and gain compensation.

[0045] Existing technologies typically employ multiple geophones connected in series via cables to acquire vibration source signals. These geophones need to be arranged in a straight line with a fixed spacing, usually 1 meter, 2 meters, or 5 meters. However, in complex environments, obstacles such as buildings, ditches, and vehicles are unavoidable, making it impossible to arrange all geophones in a straight line with equal spacing. Therefore, the geophones in this invention preferably utilize distributed cableless seismic geophones. These geophones can be arranged arbitrarily along the survey line, with the arrangement including both equal and unequal spacing, depending on the specific conditions of the actual scenario.

[0046] This invention utilizes distributed cableless seismic detectors to acquire the vibration source signals and coordinate positions of the seismic sources. The vibration source signals are then preprocessed to obtain vibration measurement data. It is worth noting that the vibration source signal is a waveform vibration signal transmitted from the seismic source, which is artificially generated by methods such as hammering or dropping weights to induce seismic waves. Multiple seismic sources are used, and they are activated sequentially. Each time a seismic source is activated, the detector acquires data. By activating all seismic sources sequentially, each detector generates multiple sets of co-source measurement data.

[0047] Then, the measurement data of each set of resonance sources are filtered to include: partial resonance source measurement data whose offset meets the set conditions, and partial resonance source measurement data comes from at least one corresponding detector.

[0048] Since the locations of the seismic source and the geophone significantly affect the calculation results of the surface wave dispersion spectrum, if the geophone is too close to the source and the seismic trace energy is too strong, the acquired source signal will be too strong; if the geophone is too far from the source and the seismic trace energy is too weak, the acquired source signal will be too weak or even nonexistent. Both excessively strong and weak source signals will affect the overall imaging quality of the surface wave dispersion spectrum, thus impacting the accuracy of dispersion curve acquisition and inversion. Therefore, by screening the resonance source measurement data for each group, geophones with moderately strong source signals are selected for calculation, resulting in higher imaging quality of the obtained surface wave dispersion spectrum.

[0049] Each set of resonance source measurement data is screened. The main screening method is to select detectors with appropriate vibration source signals by filtering the resonance source measurement data whose offset meets the set conditions. It should be noted that the resonance source measurement data comes from at least one corresponding detector. The set conditions are also measured by the strength of the vibration source signal collected by the detector in continuous testing.

[0050] Secondly, one of the detectors corresponding to a portion of the resonance source measurement data is selected as the center channel, and the remaining detectors within a set range adjacent to the center channel are selected as adjacent channels; the surface wave dispersion spectrum of the center channel is calculated using the resonance source measurement data of the center channel and the adjacent channels. Preferably, the surface wave dispersion spectrum of the center channel is calculated by performing phase tilt superposition of the resonance source measurement data generated by the same set of sources in the center channel and the adjacent channels.

[0051] It should be noted that the center channel is selected from the geophones corresponding to the portion of the resonant source measurement data whose offset meets the set conditions, while adjacent channels must be the remaining geophones within the set range adjacent to the center channel. The set conditions refer to the range of offset values ​​for the center channel, which must be within a suitable range of seismic channel energy; the set range is a circle with radius L centered on the center channel geophone. Radius L is calculated based on multiple experiments, and the evaluation criterion is obtaining a high-quality surface wave dispersion spectrum for the center channel.

[0052] In this way, multiple detectors generate multi-channel surface wave dispersion spectra, with each seismic source corresponding to one multi-channel surface wave dispersion spectrum.

[0053] Specifically, to generate a seismic source, the phase tilt of the co-source measurement data of one detector meeting the set conditions and the detectors in its adjacent channels is superimposed to obtain the surface wave dispersion spectrum of the center channel. This process is repeated to obtain the multi-channel surface wave dispersion spectra of all detectors meeting the set conditions. Then, generating another seismic source yields another set of co-source measurement data, resulting in another set of multi-channel surface wave dispersion spectra for detectors meeting the set conditions. This process is repeated for all the set seismic sources, resulting in multiple multi-channel surface wave dispersion spectra. Multiple seismic sources correspond to multiple multi-channel surface wave dispersion spectra.

[0054] Next, taking one detector as a unit, it is determined whether a detector corresponds to obtaining at least two center channel surface wave dispersion spectra; if so, the above at least two center channel surface wave dispersion spectra are superimposed and calculated to obtain the surface wave dispersion spectrum of this detector; if not, the obtained center channel surface wave dispersion spectrum is used as the surface wave dispersion spectrum of this detector.

[0055] In other words, when a detector acquires and calculates at least two center channel surface wave dispersion spectra, the surface wave dispersion spectrum of the detector is obtained by superimposing all the center channel surface wave dispersion spectra acquired by the detector.

[0056] In this way, the surface wave dispersion spectrum of all detectors is obtained sequentially.

[0057] Preferably, the specific formula for calculating the surface wave dispersion spectrum of the detector is as follows:

[0058]

[0059] Among them, A i (f) represents the normalized amplitude spectrum of vibration measurement data within the set range; d c Indicates the offset distance of the center track; d i D1 represents the offset distance between adjacent lanes; L represents the radius of the set range; D1 and D2 both represent the range values ​​of the offset distance of the center lane that meet the set conditions, that is, the offset distance of the center lane must satisfy D1≤|d c |≤D2, D1 and D2 are also calculated based on multiple experiments, and the evaluation criterion is to obtain a high-quality multi-segment surface wave dispersion spectrum.

[0060] When the center channel is located at the beginning or end of the array, the surface wave dispersion spectrum of the center channel is calculated after extending the array by adding a detector based on a set range.

[0061] Because the detectors at both ends of the array are unique, when calculating the dispersion spectrum of the center channel wave at both ends, there may be insufficient data channels on one side. That is, when the center channel is located at the beginning or end of the array, the detector at the beginning or end only has an adjacent channel on one side within the set range, which cannot accurately represent the dispersion spectrum of the center channel wave at the beginning or end, resulting in an inaccurate underground shear wave velocity structure. Therefore, it is necessary to appropriately extend or move the array and repeat a certain number of channels. This involves setting up several additional detectors in front of the beginning detector along the survey line to obtain sufficient channel data. These newly added detectors only collect data for use by the beginning or end detectors; their center channel wave dispersion spectrum is not used.

[0062] Finally, the subsurface shear wave velocity structure of each detector is obtained sequentially based on the surface wave dispersion spectrum, thereby obtaining the subsurface shear wave velocity profile within the measurement range to complete the exploration.

[0063] Specifically, the dispersion curve of the detector is first obtained based on its surface wave dispersion spectrum. Then, the dispersion curve is inverted to obtain the subsurface shear wave velocity structure of the detector. The inversion process includes methods such as least squares, simulated annealing, and genetic algorithms. Next, the subsurface shear wave velocity structure obtained through channel-by-channel inversion is interpolated laterally in two dimensions to obtain a two-dimensional shear wave velocity profile. This provides the subsurface shear wave velocity profile within the measurement range, enabling exploration in complex environments.

[0064] To gain a more detailed and clear understanding of the present invention, in conjunction with Figure 2 Further explanation. For example... Figure 2 The diagram shown is a schematic representation of the detector arrangement and data processing of a multi-excitation point transient surface wave exploration method, equipment, and storage medium according to the present invention.

[0065] The geophones are arranged linearly along the survey line. Due to some obstructions, the spacing between the geophones is unequal. The seismic source is located to the left of the first geophone and at a certain distance from it. It is worth noting that when there are no obstructions, the geophones are preferably arranged at equal intervals; when there are obstructions, they are arranged at unequal intervals, excluding the obstructions.

[0066] First, the first seismic source is excited, and the offset distance satisfies D1≤|d c The detector with an offset of |d1≤|d2 is selected based on the condition that the first detector meets the condition. c If |≤D2, then the first detector can be the center channel. With the first detector as the center channel, the remaining detectors within a radius L centered on the first detector are the adjacent channels. Here, L is a radius within a predetermined range determined experimentally; preferably, such as... Figure 2 As shown, the second and third detectors are adjacent channels. By performing phase tilt superposition calculation on the resonance source measurement data of the center channel and the two adjacent channels, the surface wave dispersion spectrum of the center channel of the first detector is obtained.

[0067] Since the first geophone has an adjacent channel on only one side, a certain number of channels are repeated when the arrangement is appropriately extended or moved. That is, several more geophones are set up to the left of the first geophone along the survey line to obtain sufficient channel data. Based on the L value, two geophones, A and B, are preferably set up. The resonance source measurement data of the first geophone, the second geophone, the third geophone, and the newly added geophones A and B are used to perform phase tilt superposition calculations to obtain the center channel wave dispersion spectrum of the first geophone.

[0068] Next, taking the second detector as the center channel and the detectors within a radius of L as adjacent channels, that is, the first detector, the third detector, the fourth detector and the newly added detector B are all adjacent channels. The resonant source measurement data of these five detectors, including the center channel and the adjacent channels, are superimposed by phase tilt calculation to obtain the surface wave dispersion spectrum of the center channel of the second detector.

[0069] This process continues until the offset satisfies D1≤|d cAfter all center channels of |≤D2 have been detected, the multi-channel surface wave dispersion spectrum corresponding to the first source is obtained. The multi-channel surface wave dispersion spectrum includes the center channel surface wave dispersion spectrum generated by each detector.

[0070] A second seismic source is then excited, and the detectors that meet the conditions are selected to obtain the multi-channel surface wave dispersion spectrum corresponding to the second seismic source. This process is repeated to excite all seismic sources in sequence, resulting in multiple multi-channel surface wave dispersion spectra.

[0071] Next, taking the first detector as a unit, based on multiple multi-channel surface wave dispersion spectra, it is determined whether the first detector has obtained at least two center channel surface wave dispersion spectra: if so, then all the center channel surface wave dispersion spectra obtained by the first detector are superimposed and calculated to finally obtain a more accurate surface wave dispersion spectrum of the first detector; if not, it means that the first detector has only obtained one center channel surface wave dispersion spectrum, then the only center channel surface wave dispersion spectrum obtained by this detector is used as the surface wave dispersion spectrum of this detector.

[0072] Then, based on the surface wave dispersion spectrum of each detector, the dispersion curve of the corresponding detector is obtained, and the dispersion curve is inverted channel by channel to obtain the underground shear wave velocity structure of each detector. Figure 2 In the data graph, the horizontal axis represents speed and the vertical axis represents depth; the greater the depth, the greater the speed.

[0073] Finally, the subsurface shear wave velocity structure of each detector obtained by inversion is subjected to lateral two-dimensional interpolation to obtain a two-dimensional shear wave velocity profile, thereby obtaining the subsurface shear wave velocity profile within the measurement range, thus completing the exploration in complex environments.

[0074] According to a second aspect of the present invention, a multi-source transient surface wave exploration device is also provided, comprising at least one processing unit and at least one storage unit, wherein the storage unit stores a computer program, and when the computer program is executed by the processing unit, the processing unit performs the steps of the above-described method.

[0075] According to a third aspect of the invention, the invention also provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the above-described method. The computer-readable storage medium may include, but is not limited to, any type of disk, including floppy disks, optical disks, DVDs, CD-ROMs, microdrives, as well as magneto-optical disks, ROMs, RAMs, EPROMs, EEPROMs, DRAMs, VRAMs, flash memory devices, magnetic cards or optical cards, nanosystems (including molecular memory ICs), or any type of medium or device suitable for storing instructions and / or data.

[0076] It should be noted that, for the sake of simplicity, the foregoing method embodiments are all described as a series of actions. However, those skilled in the art should understand that the present invention is not limited to the described order of actions, as some steps can be performed in other orders or simultaneously according to the present invention. Furthermore, those skilled in the art should also understand that the embodiments described in the specification are preferred embodiments, and the actions and modules involved are not necessarily essential to the present invention.

[0077] In the above embodiments, the descriptions of each embodiment have different focuses. For parts not described in detail in a certain embodiment, please refer to the relevant descriptions in other embodiments.

[0078] In the several embodiments provided by this invention, it should be understood that the disclosed apparatus can be implemented in other ways. For example, the apparatus embodiments described above are merely illustrative; for instance, the division of units is only a logical functional division, and in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Furthermore, the coupling or direct coupling or communication connection shown or discussed may be through some service interface; the indirect coupling or communication connection between devices or units may be electrical or other forms.

[0079] The units described as separate components may or may not be physically separate. 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 can be selected to achieve the purpose of this embodiment according to actual needs.

[0080] Furthermore, the functional units in the various embodiments of the present invention can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit. The integrated unit can be implemented in hardware or as a software functional unit.

[0081] If the integrated unit is implemented as a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage device (CMD). Based on this understanding, the technical solution of this invention, in essence, or the part that contributes to the prior art, or all or part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a memory and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods of the various embodiments of this invention. The aforementioned memory includes various media capable of storing program code, such as USB flash drives, read-only memory (ROM), random access memory (RAM), portable hard drives, magnetic disks, or optical disks.

[0082] Those skilled in the art will understand that all or part of the steps in the various methods of the above embodiments can be implemented by a program instructing related hardware. The program can be stored in a computer-readable storage medium, which may include: a flash drive, a read-only memory (ROM), a random access memory (RAM), a magnetic disk, or an optical disk, etc.

[0083] The foregoing description is merely an exemplary embodiment of this disclosure and should not be construed as limiting the scope of this disclosure. Any equivalent changes and modifications made in accordance with the teachings of this disclosure shall still fall within the scope of this disclosure. Those skilled in the art will readily conceive of other embodiments of this disclosure upon considering the specification and practicing the disclosure herein. This invention is intended to cover any variations, uses, or adaptations of this disclosure that follow the general principles of this disclosure and include common knowledge or customary techniques in the art not described herein. The specification and embodiments are to be considered exemplary only, and the scope and spirit of this disclosure are defined by the claims.

[0084] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0085] Those skilled in the art will readily understand that the above description is merely a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, and improvements made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

Claims

1. A multi-excitation-point transient surface wave exploration method, characterized in that, The method includes: Multiple seismic sources are excited one by one, and the excitation of each seismic source generates a set of seismic source signals; detectors arranged along the measurement line collect the seismic source signals one by one to generate multiple sets of co-source measurement data; Each set of resonance source measurement data is filtered, wherein the filtering is performed on a portion of the resonance source measurement data whose offset meets the set conditions, and the portion of the resonance source measurement data comes from at least one of the corresponding detectors; One of the detectors corresponding to the partial resonance source measurement data is selected as the center channel, and the remaining detectors within the adjacent set range of the center channel are selected as adjacent channels. The resonance source measurement data generated by one source of the center channel and the adjacent channels are superimposed by phase tilt to obtain the surface wave dispersion spectrum of the center channel. The set range is a circle with a radius of L centered on the center channel detector. When the center channel is located at the beginning or end of the arrangement, the surface wave dispersion spectrum of the center channel is calculated after extending the arrangement by adding detectors based on the set range. In this way, multiple detectors generate multi-channel surface wave dispersion spectra, with each source corresponding to one such multi-channel surface wave dispersion spectrum; Determine whether a detector obtains at least two center channel surface wave dispersion spectra; if so, calculate by superimposing all the center channel surface wave dispersion spectra obtained by the detector to obtain the surface wave dispersion spectrum of the detector; if not, use the obtained center channel surface wave dispersion spectrum as the surface wave dispersion spectrum of the detector. In this manner, the surface wave dispersion spectrum of all the detectors is obtained sequentially; Based on the surface wave dispersion spectrum, the underground shear wave velocity structure of each of the detectors is obtained sequentially, and the underground shear wave velocity structure of each of the central channel detectors is interpolated laterally in two dimensions to obtain the underground shear wave velocity profile within the measurement range, thus completing the exploration.

2. The multi-excitation-point transient surface wave exploration method as described in claim 1, characterized in that, The resonance source measurement data includes vibration measurement data and the coordinate position of the detector.

3. The multi-excitation point transient surface wave exploration method as described in claim 2, characterized in that, The formula for calculating the surface wave dispersion spectrum of the center channel detector is as follows: in, This represents the normalized amplitude spectrum of the vibration measurement data within the defined range; This indicates the offset distance of the central channel; This indicates the offset distance between the adjacent channels; This represents the radius of the defined range; , All of these represent the range of offset distances of the center track that meet the set conditions.

4. The multi-excitation point transient surface wave exploration method as described in claim 2, characterized in that, The resonance source measurement data is obtained by acquiring the vibration source signal using the detector and then preprocessing the vibration source signal. The preprocessing includes noise suppression, bad sector processing, filtering, and gain compensation.

5. The multi-excitation point transient surface wave exploration method as described in claim 1, characterized in that, The detector is a distributed cableless seismic detector; the detectors are arranged in two ways: with equal spacing and with unequal spacing.

6. The multi-excitation point transient surface wave exploration method as described in claim 1, characterized in that, The underground shear wave velocity structure is obtained by inverting the dispersion curve of each of the central channel detectors based on the surface wave dispersion spectrum. The inversion includes least squares method, simulated annealing method and genetic algorithm.

7. A multi-excitation-point transient surface wave exploration device, characterized in that, It includes at least one processing unit and at least one storage unit, wherein the storage unit stores a computer program that, when executed by the processing unit, causes the processing unit to perform the steps of the method according to any one of claims 1-6.

8. A multi-excitation-point transient surface wave exploration storage medium, characterized in that, It stores a computer program executable by an access authentication device, which, when run on the access authentication device, causes the access authentication device to perform the steps of the method described in any one of claims 1-6.

Images