Method and system for predicting propagation direction of seismic wave field excited by high-speed rail seismic source

Abstract

The invention discloses a method and a system for predicting the propagation direction of a seismic wave field excited by a high-speed rail seismic source. The method comprises the following steps: simplifying the high-speed rail seismic source into vertical downward impact force of wheels on the ground; taking the Ricker wavelet as a deformation function of force, summing all wheels of the high-speed rail train to obtain a high-speed rail sub-seismic source function, and performing Fourier transform to obtain an amplitude spectrum; the method comprises the following steps: modeling uniformly arranged sleepers or piers of a high-speed rail line into a uniform linear array, exciting seismic waves at each array element in the uniform linear array, and summing seismic wave signals of each array element at a far-field receiving point to obtain an array factor; a radiation factor formula of the P wave and the S wave of the seismic wave field of the high-speed rail seismic source is obtained by combining the radiation pattern factors of the seismic P wave and the S wave excited by the single-point concentrated force; drawing a theoretical interference directional diagram of seismic waves; and establishing a model for simulation, comparing a theoretically predicted interference pattern with a simulation data wave field snapshot, and keeping the predicted direction consistent with the simulated wave field direction, thereby verifying the rationality and correctness of the proposed method.

Description

Background Art

[0002] The high - speed rail seismic source is a complex mobile combined seismic source, and its repeatability, stability, and low cost are beneficial to high - speed rail seismology research. At present, there are huge challenges in using the seismic wave field of the high - speed rail seismic source to detect the medium around the high - speed rail line. One of the key points lies in how to predict the directional characteristics of the seismic wave field excited by the high - speed rail seismic source. Since it is impossible to deploy a large - scale geophone array along the high - speed rail line to observe the directional characteristics of the seismic wave field excited by the high - speed rail seismic source, it is necessary to give the directional characteristics of the seismic wave field excited by the high - speed rail seismic source according to the high - speed rail train parameters and structural parameters such as the subgrade. The existing predictions of the directional characteristics of the seismic wave field excited by the high - speed rail seismic source mainly include: Prior Art 1: Deploying a geophone array along the high - speed rail line to simplify the wave field propagation direction For example, through the cross - correlation method, extracting the seismic wave signal propagating along the viaduct as a characteristic quantity of the high - speed rail seismic wave field. According to the characteristics of the coherent phase stationary point region, a method for separating the virtual seismic wave field excited by a single pier is used.

[0003] Disadvantages of Prior Art 1: When collecting data, geophones are only arranged along the line, and the directivity of the interference wave field of the high - speed rail seismic source is not considered when analyzing the data.

[0004] Prior Art 2: Using a simple method of representing wave fronts and wave crests with concentric circles to obtain the basic distribution of the interference wave field.

[0005] For example, by drawing concentric circles when the high - speed rail wheel passes through the excitation position, using a simple method of representing wave fronts and wave crests with concentric circles to obtain the basic distribution of the interference wave field.

[0006] Disadvantages of Prior Art 2: Simply regarding the seismic waves excited by the high - speed rail seismic source as concentric circles, without considering the influence of the spatial distribution of the high - speed rail seismic source on the wave field interference phenomenon. Summary of the Invention

[0007] The technical problem to be solved by the present invention is to provide a method and system for predicting the propagation direction of the seismic wave field excited by a high - speed rail seismic source in view of the deficiencies in the above - mentioned prior art, so as to solve the technical problems of lacking consideration of the directivity of the interference wave field of the high - speed rail seismic source and the influence of the spatial distribution of the high - speed rail seismic source on the wave field interference phenomenon in the existing methods.

[0008] The present invention adopts the following technical solutions: A method for predicting the propagation direction of the seismic wave field excited by a high - speed rail seismic source, comprising the following steps: Simplify the interaction between the high - speed rail wheel and the ground into a vertically downward impact force of the high - speed rail wheel on the ground, obtain the theoretical ground force function, and decompose the theoretical force function to obtain the theoretical force function and the amplitude spectrum; Let the Ricker wavelet be the deformation function of the force, convolve it with the obtained theoretical force function and amplitude spectrum to get the high-speed rail source function, and then perform Fourier transform on the high-speed rail source function to obtain the high-speed rail source function and amplitude spectrum; Construct a uniform linear array based on the sleepers or piers of the high-speed rail line, excite seismic wave signals at each element in the uniform linear array, and sum all the seismic wave signals to obtain the array factor; Based on the obtained high-speed rail source function and amplitude spectrum, combined with the seismic P-wave and S-wave radiation pattern factors excited by a single-point concentrated force, multiply them by the obtained array factor to get the radiation factors of the P-wave and S-wave of the high-speed rail source seismic wave field; Draw the theoretical interference direction diagram of the seismic wave according to the array factor in the obtained radiation factor to complete the prediction of the propagation direction of the seismic wave field.

[0009] Preferably, the amplitude spectrum of the high-speed rail source function is:

[0010] where, N t is the number of high-speed rail carriages, v t is the running speed of the high-speed rail, C is a constant related to the mass of the high-speed rail, is the Fourier transform form of the function, is the Fourier transform form of the function, is the Fourier transform form of the function, is the distance between the second set of wheel pairs and the first set of wheel pairs in a single carriage, is the distance between the third set of wheel pairs and the first set of wheel pairs in a carriage, is the length of a single high-speed rail carriage, is the frequency.

[0011] Preferably, perform Fourier transform on the high-speed rail source function to obtain the high-speed rail source function and amplitude spectrum

[0012] where, d(k) is the k distance from the group of wheel pairs to the first group of wheel pairs,

[0013] Preferably, perform Fourier transform on the high-speed rail source function to obtain the high-speed rail source function and amplitude spectrum as follows:

[0014]

[0015] Among them, is the high-speed rail point source function at a single sleeper or pier, is the Ricker wavelet, is a constant related to the mass of the high-speed rail, is the structure function of the second set of wheel pairs related to the first set of wheel pairs in a single carbody, is the structure function of the third set of wheel pairs related to the first set of wheel pairs in a single carbody, is the periodic input function related to the number of high-speed rail carriages, is the distance between the second set of wheel pairs and the first set of wheel pairs in a single carbody, is the running speed of the high-speed rail, is the angular frequency, is the distance between the third set of wheel pairs and the first set of wheel pairs in a single carbody, is the number of high-speed rail carriages, is the length of a single high-speed rail carbody, is the Fourier transform form of the Ricker wavelet function.

[0016] Preferably, the Fourier transform of the Ricker wavelet function is:

[0017] Among them, is the Fourier transform form of the Ricker wavelet function, is the main frequency, is the angular frequency. Preferably, the array factor is:

[0018] Among them, is the array factor affecting the wave field direction, is the frequency, is the adjacent element spacing, is the angle between the wave propagation direction from the source to the receiving point and the running direction of the high-speed rail, is the wave speed, is the running speed of the high-speed rail, is the number of array elements.

[0019] Preferably, the seismic waves excited at each array element received by the far-field receiving point are expressed as:

[0020] Among them, is the time variable of the source function, is the wave path from the source to the receiving point.

[0021] Preferably, the radiation factors of the P-wave and S-wave in the seismic wave field of the high-speed rail source are:

[0022] Among them, is the radiation pattern factor of the seismic P-wave excited by the high-speed rail sub-source, is the radiation pattern factor of the P-wave, and are the array factors affecting the wave field direction, is the angle between the wave propagation direction and the direction of the concentrated force application, is the frequency, is the interval between adjacent array elements, is the wave velocity of the P-wave, is the running speed of the high-speed rail, is the angle between the wave propagation direction from the source to the receiving point and the running direction of the high-speed rail, is the number of array elements, is the wave velocity of the S-wave, and are the radiation pattern factors of the P-wave and S-wave respectively, is the radiation pattern factor of the seismic P-wave excited by the high-speed rail sub-source.

[0023] Preferably, the radiation pattern factors of the seismic P-wave and S-wave excited by a single-point concentrated force are: .

[0025] Secondly, the embodiment of the present invention provides a prediction system for the propagation direction of the seismic wave field excited by a high-speed rail source, including: The decomposition module simplifies the interaction between the high-speed rail wheel and the ground into a vertically downward impact force of the high-speed rail wheel on the ground, obtains the theoretical force function of the ground, and decomposes the theoretical force function to obtain the theoretical force function and the amplitude spectrum; The convolution module makes the Ricker wavelet the deformation function of the force, convolves it with the obtained theoretical force function and amplitude spectrum to obtain the high-speed rail sub-source function, and then performs Fourier transform on the high-speed rail sub-source function to obtain the high-speed rail sub-source function and the amplitude spectrum; The summation module constructs a uniform linear array based on the sleepers or piers of the high-speed rail line, excites seismic wave signals at each array element in the uniform linear array, and sums all the seismic wave signals to obtain the array factor; The factor module multiplies the high-speed rail sub-source function and the amplitude spectrum, combines the radiation pattern factors of the seismic P-wave and S-wave excited by a single-point concentrated force, and multiplies them with the array factor to obtain the radiation factors of the P-wave and S-wave in the seismic wave field of the high-speed rail source; A prediction module draws the theoretical interference direction diagram of seismic waves based on the array factor in the obtained radiation factor, and completes the prediction of the propagation direction of the seismic wave field.

[0026] In a third aspect, a computer device includes a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the steps of the method for predicting the propagation direction of the seismic wave field excited by the high-speed rail seismic source are implemented.

[0027] In a fourth aspect, an embodiment of the present invention provides a computer-readable storage medium, including a computer program, and when the computer program is executed by a processor, the steps of the method for predicting the propagation direction of the seismic wave field excited by the high-speed rail seismic source are implemented.

[0028] In a fifth aspect, a chip includes a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the steps of the method for predicting the propagation direction of the seismic wave field excited by the high-speed rail seismic source are implemented.

[0029] In a sixth aspect, an embodiment of the present invention provides an electronic device, including a computer program, and when the computer program is executed by the electronic device, the steps of the method for predicting the propagation direction of the seismic wave field excited by the high-speed rail seismic source are implemented.

[0030] Compared with the prior art, the present invention has at least the following beneficial effects: A method for predicting the propagation direction of the seismic wave field excited by a high-speed rail seismic source includes a verification model designed according to the disclosed prediction analysis method. The high-speed rail seismic source is simplified into a vertical downward impact force of the wheel on the ground, and the Ricker wavelet is used as the deformation function of the force. The high-speed rail sub-seismic source function is obtained by summing all the wheels of the high-speed rail train and Fourier transform is performed to obtain the amplitude spectrum. Then, the sleepers or piers evenly arranged on the high-speed rail line are modeled into a uniform linear array, the array factor is calculated, and combined with the radiation pattern factors of the seismic P-wave and S-wave excited by a single-point concentrated force, and multiplied by the obtained array factor to obtain the radiation factor formulas of the P-wave and S-wave of the high-speed rail seismic source seismic wave field. The theoretical interference direction diagram of the seismic wave is drawn according to the formula. Compared with the existing method for analyzing the characteristics of the high-speed rail seismic source wave field, the influence of the directivity of the high-speed rail seismic source interference wave field and the influence of the spatial distribution of the high-speed rail seismic source on the wave field interference phenomenon are considered.

[0031] Further, the high-speed rail seismic source is simplified into a vertical downward impact force of the wheel on the ground; the Ricker wavelet is used as the deformation function of the force, and the high-speed rail sub-seismic source function is obtained by summing all the wheels of the high-speed rail train and Fourier transform is performed to obtain the amplitude spectrum.

[0032] Furthermore, the sleepers or piers arranged uniformly along the high-speed railway line are modeled as a uniform linear array. Due to the interaction between the wheels and the rails, seismic waves are excited at each element in the uniform linear array. A far-field receiving point is set, and the seismic wave signals of each element are summed at the far-field receiving point to obtain the array factor.

[0033] Furthermore, by combining the seismic P-wave and S-wave radiation pattern factors excited by a single-point concentrated force, and multiplying them with the obtained array factor, the radiation factor formulas of the P-wave and S-wave of the high-speed railway seismic source wave field are derived.

[0034] Furthermore, according to the array factor in the formula, the theoretical interference direction diagram of the seismic wave is drawn to facilitate subsequent comparative analysis.

[0035] It can be understood that the beneficial effects of the above-mentioned second to sixth aspects can be referred to the relevant descriptions in the above-mentioned first aspect, and will not be elaborated here.

[0036] In summary, the present invention makes innovations in the existing analysis methods of the characteristics of high-speed railway seismic source wave fields, gives a theoretical formula for a method for predicting the propagation direction of seismic waves excited by high-speed railway seismic sources. On this basis, a model is established for simulation, and the theoretical direction diagram is compared with the wave field snapshots of the simulation data. The predicted direction is consistent with the simulated wave field, verifying the correctness and rationality of the present method.

[0037] Next, through the drawings and embodiments, the technical solutions of the present invention will be further described in detail. BRIEF DESCRIPTION OF THE DRAWINGS

[0038] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings required to be used in the embodiments of the present application will be briefly introduced below. Obviously, the following described drawings are only some embodiments of the present application. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts.

[0039] Figure 1 is the flow chart of the present invention; FIG. 2 is the time-domain waveform diagram and frequency spectrum diagram of the high-speed railway sub-seismic source, where (a) is the time-domain waveform diagram; (b) is the frequency spectrum diagram; FIG. 3 is a schematic diagram of the amplitude size of the array factor, where (a) is the element spacing d taking 1 m; (b) is the element spacing d taking 2 m: (c) is the element spacing d taking 5 m; FIG. 4 is the theoretical prediction direction diagram, where (a) is the Z-axis section; (b) is the X-axis section; Figure 5 is the radiation pattern diagram of seismic waves excited by a single-point concentrated force source; Figure 6 is a snapshot of the simulated wave field, where (a) is the Z-axis cross-section; (b) is the X-axis cross-section; Figure 7 is a combined diagram of the theoretical radiation pattern and the snapshot of the simulated wave field, where (a) is the Z-axis cross-section; (b) is the X-axis cross-section; Figure 8 Schematic diagram of a computer device provided by an embodiment of the present invention; Figure 9 Block diagram of an electronic device provided by an embodiment of the present invention.

[0040] Among them, 60. Computer device; 61. Processor; 62. Memory; 63. Computer program; 600. Electronic device; 610. Processing unit; 620. Storage unit; 6201. Random access storage unit; 6202. Cache storage unit; 6203. Read-only storage unit; 6204. Program / utilities; 6205. Program module; 630. Bus; 640. Display unit; 650. Input / output interface; 660. Network adapter; 700. External device. Detailed implementation manners

[0041] Next, the technical solutions in the embodiments of the present invention will be clearly and completely described in conjunction with the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are some, but not all, of the embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0042] In the description of the present invention, it should be understood that the terms "include" and "comprise" indicate the presence of the described features, wholes, steps, operations, elements, and / or components, but do not exclude the presence or addition of one or more other features, wholes, steps, operations, elements, components, and / or their combinations.

[0043] It should also be understood that the terms used in the specification of the present invention are only for the purpose of describing specific embodiments and are not intended to limit the present invention. As used in the specification of the present invention and the appended claims, unless the context clearly indicates otherwise, the singular forms "a", "an", and "the" are intended to include the plural forms.

[0044] It should be further understood that the term " / and / or" used in the specification of the present invention and the appended claims refers to any combination and all possible combinations of one or more of the associated listed items, and includes these combinations. For example, A and / or B can represent: A exists alone, A and B exist simultaneously, and B exists alone. In addition, the character " / " in the present invention generally represents an "or" relationship between the preceding and following related objects.

[0045] It should be understood that although terms such as first, second, and third may be used in the embodiments of the present invention to describe preset ranges, etc., these preset ranges should not be limited to these terms. These terms are only used to distinguish the preset ranges from each other. For example, without departing from the scope of the embodiments of the present invention, the first preset range may also be referred to as the second preset range, and similarly, the second preset range may also be referred to as the first preset range.

[0046] Depending on the context, the word "if" as used herein may be interpreted as "when" or "while" or "in response to determining" or "in response to detecting". Similarly, depending on the context, the phrase "if determined" or "if detected (stated condition or event)" may be interpreted as "when determined" or "in response to determining" or "when detected (stated condition or event)" or "in response to detecting (stated condition or event)".

[0047] Various structural schematic diagrams according to the disclosed embodiments of the present invention are shown in the drawings. These figures are not drawn to scale, where for the purpose of clear expression, some details are enlarged and some details may be omitted. The shapes of various regions and layers shown in the figures and their relative sizes and positional relationships are only exemplary, and in practice, there may be deviations due to manufacturing tolerances or technical limitations, and those skilled in the art can design regions / layers with different shapes, sizes, and relative positions according to actual needs.

[0048] The present invention provides a method for predicting the propagation direction of seismic wave fields excited by high-speed rail sources. The sleepers or piers arranged uniformly along the high-speed rail line are modeled as a uniform linear array. According to the total length of the train, the running speed of the train, the sleeper spacing or pier spacing, and the seismic wave speed in the medium, the array factor can be calculated. Combining the radiation pattern factors of seismic P-waves and S-waves excited by a single-point concentrated force and multiplying them with the obtained array factor, the radiation factor formulas of P-waves and S-waves of the high-speed rail source seismic wave field are obtained, and the direction of the seismic wave field excited by the high-speed rail source is predicted.

[0049] Embodiment 1 A method for predicting the propagation direction of seismic wave fields excited by high-speed rail sources according to the present invention includes the following steps: S1. Simplify the interaction between the high-speed rail wheels and the ground into a vertically downward impact force of the high-speed rail wheels on the ground to obtain a theoretical ground force function, and decompose the theoretical force function to obtain the theoretical force function and the amplitude spectrum; The theoretical force function is:

[0050] Wherein, d(k) is the distance from the k th wheel pair to the first wheel pair, vt is the running speed of the high - speed train, N t is the number of carriages of the high - speed train.

[0051] Among them, p ( t ) can also be expressed as:

[0052] Among them,

[0053] Among them, d 1 , d 2 are respectively the distance between the second set of wheel sets and the first set of wheel sets and the distance between the third set of wheel sets and the first set of wheel sets in one carriage. Performing Fourier transform on this formula gives:

[0054] Therefore, the theoretical force - bearing function and amplitude spectrum are:

[0055] S2. Let the Ricker wavelet be the deformation function of the force, convolve it with the theoretical force - bearing function to obtain the high - speed train sub - source function, and then perform Fourier transform on the high - speed train sub - source to obtain the high - speed train sub - source function and amplitude spectrum; The Ricker wavelet is expressed as:

[0056] Among them, f 0 is the main frequency.

[0057] Performing Fourier transform on the Ricker wavelet function gives:

[0058] Convolving the Ricker wavelet function with the theoretical force - bearing function to obtain the high - speed train sub - source function and amplitude spectrum as:

[0059]

[0060] S3. Model the sleepers or piers arranged uniformly on the high - speed rail line as a uniform linear array. The interaction between the wheels and the rails generates vibrations, and the vibrations are transmitted to the ground through the sleepers or piers, thereby exciting seismic waves. For the seismic waves excited at each array element in the uniform linear array modeled in the present invention, set a far - field receiving point, and sum the seismic wave signals of each array element at the far - field receiving point to obtain the array factor; The seismic waves received at the far-field receiving points are expressed as:

[0061] where v m is the wave velocity, v t is the traveling speed of the high-speed rail, D is the wave path from the seismic source to the receiving point, d is the adjacent array element interval, N is the number of array elements, is the angle between the wave propagation direction from the seismic source to the receiving point and the running direction of the high-speed rail. Performing Fourier transform on this formula gives:

[0062] Using the superposition principle, summing and taking the modulus of the seismic wave signals of each array element gives:

[0063] Therefore, the array factor is obtained:

[0064] S4. Combining the seismic P-wave and S-wave radiation pattern factors excited by a single-point concentrated force in the seismic source propagation theory, multiplying them with the obtained array factor to obtain the radiation factor formulas for the P-wave and S-wave of the high-speed rail seismic source seismic wave field; The seismic P-wave and S-wave radiation pattern factors excited by a single-point concentrated force are:

[0065] where is the wave propagation direction and the direction of the concentrated force application, is the P-wave radiation pattern factor, is the S-wave radiation pattern factor.

[0066] Multiplying the P-wave and S-wave radiation pattern factors with the array factor gives the radiation factor formulas for the P-wave and S-wave of the high-speed rail seismic source excited seismic wave field:

[0067] S5. Draw the theoretical interference direction diagram of the seismic waves according to the array factor in the formula; S6. Establish a finite element model and perform simulation; compare the theoretical interference direction diagram with the simulation data wave field snapshots to verify the rationality and correctness of the proposed method.

[0068] Those skilled in the art can understand that various aspects of the present invention can be implemented as a system, method, or program product. Therefore, various aspects of the present invention can be specifically implemented in the following forms, namely: a complete hardware implementation, a complete software implementation (including firmware, microcode, etc.), or an implementation combining hardware and software aspects, which can be collectively referred to as "circuit", "module", or "platform".

[0069] Embodiment 2 The present invention provides a prediction system for the propagation direction of seismic wave fields excited by high-speed rail seismic sources. This system can be used to implement the above-mentioned method for predicting the propagation direction of seismic wave fields excited by high-speed rail seismic sources. Specifically, the prediction system for the propagation direction of seismic wave fields excited by high-speed rail seismic sources includes a decomposition module, a convolution module, a summation module, a factor module, and a prediction module.

[0070] Among them, the decomposition module simplifies the interaction between the high-speed rail wheel and the ground into a vertically downward impact force of the high-speed rail wheel on the ground, obtains the theoretical ground force function, and decomposes the theoretical ground force function to obtain the theoretical ground force function and the amplitude spectrum. The convolution module takes the Ricker wavelet as the deformation function of force, convolves it with the obtained theoretical ground force function and amplitude spectrum to obtain the high-speed rail sub-seismic source function, and then performs a Fourier transform on the high-speed rail sub-seismic source function to obtain the high-speed rail sub-seismic source function and the amplitude spectrum. The summation module constructs a uniform linear array based on the sleepers or piers of the high-speed rail line, excites seismic wave signals at each element of the uniform linear array, and sums all the seismic wave signals to obtain the array factor. The factor module, based on the high-speed rail sub-seismic source function and the amplitude spectrum, combines the radiation pattern factors of seismic P-waves and S-waves excited by a single-point concentrated force, and multiplies them by the array factor to obtain the radiation factors of seismic P-waves and S-waves in the high-speed rail seismic source wave field. The prediction module draws the theoretical interference direction diagram of seismic waves according to the array factor in the obtained radiation factor, and completes the prediction of the propagation direction of the seismic wave field.

[0071] Embodiment 3 The present invention provides a terminal device, which includes a processor and a memory. The memory is used to store a computer program, and the computer program includes program instructions. The processor is used to execute the program instructions stored in the computer storage medium. The processor may be a Central Processing Unit (CPU), or may also be other general-purpose processors, Graphics Processing Unit (GPU), Tensor Processing Unit (TPU), Digital Signal Processor (DSP), Application Specific Integrated Circuit (ASIC), Field-Programmable Gate Array (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. It is the computing core and control core of the terminal, and is suitable for implementing one or more instructions. Specifically, it is suitable for loading and executing one or more instructions to implement the corresponding method flow or corresponding function. The processor described in the embodiments of the present invention can be used for the operation of the method for predicting the propagation direction of seismic wave fields excited by high-speed rail seismic sources, including: Simplify the interaction between the high-speed rail wheel and the ground into a vertically downward impact force of the high-speed rail wheel on the ground to obtain a theoretical ground force function, decompose the theoretical force function to obtain the theoretical force function and the amplitude spectrum; let the Ricker wavelet be the deformation function of the force, convolve it with the theoretical force function and the amplitude spectrum to obtain the high-speed rail sub-seismic source function, and then perform Fourier transform on the high-speed rail sub-seismic source function to obtain the high-speed rail sub-seismic source function and the amplitude spectrum; construct a uniform linear array based on the sleepers or piers of the high-speed rail line, excite seismic wave signals at each element in the uniform linear array, and sum all the seismic wave signals to obtain the array factor; based on the high-speed rail sub-seismic source function and the amplitude spectrum, combine the radiation pattern factors of seismic P-waves and S-waves excited by a single-point concentrated force, multiply them by the array factor to obtain the radiation factors of seismic P-waves and S-waves in the high-speed rail seismic source wave field; draw the theoretical interference direction diagram of the seismic waves according to the array factor in the radiation factor to complete the prediction of the propagation direction of the seismic wave field.

[0072] Please refer to Figure 8, the terminal device is a computer device. The computer device 60 in this embodiment includes: a processor 61, a memory 62, and a computer program 63 stored in the memory 62 and executable on the processor 61. When the computer program 63 is executed by the processor 61, it implements the method for predicting the propagation direction of the seismic wave field excited by the high-speed rail seismic source in the embodiment. To avoid repetition, it will not be elaborated here one by one. Alternatively, when the computer program 63 is executed by the processor 61, it implements the functions of each model / unit in the system for predicting the propagation direction of the seismic wave field excited by the high-speed rail seismic source in the embodiment. To avoid repetition, it will not be elaborated here one by one.

[0073] The computer device 60 can be a computing device such as a desktop computer, a notebook, a palm computer, and a cloud server. The computer device 60 may include, but is not limited to, a processor 61 and a memory 62. Those skilled in the art can understand that Figure 8 This is only an example of the computer device 60 and does not constitute a limitation on the computer device 60. It may include more or fewer components than shown in the figure, or combine some components, or different components. For example, the computer device may also include input / output devices, network access devices, buses, etc.

[0074] The so-called processor 61 may be a central processing unit (CPU), or may also be other general-purpose processors, a graphics processing unit (GPU), a tensor processing unit (TPU), a digital signal processor (DSP), an application specific integrated circuit (ASIC), a field-programmable gate array (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor may be a microprocessor or the processor may also be any conventional processor, etc.

[0075] The memory 62 may be an internal storage unit of the computer device 60, such as the hard disk or memory of the computer device 60. The memory 62 may also be an external storage device of the computer device 60, such as a plug-in hard disk equipped on the computer device 60, a smart media card (SMC), a secure digital (SD) card, a flash card, etc.

[0076] Further, the memory 62 may also include both internal storage units of the computer device 60 and external storage devices. The memory 62 is used to store computer programs and other programs and data required by the computer device. The memory 62 may also be used to temporarily store data that has been output or is to be output.

[0077] Please refer to Figure 9 , the terminal device is the electronic device 600, and the electronic device 600 is presented in the form of a general-purpose computing device. The components of the electronic device may include but are not limited to: at least one processing unit 610, at least one storage unit 620, a bus 630 connecting different platform components (including the storage unit 620 and the processing unit 610), a display unit 640, etc.

[0078] Among them, the storage unit stores program codes, and the program codes can be executed by the processing unit 610, so that the processing unit 610 executes the steps according to various exemplary embodiments of the present invention described in the above method part of this specification. For example, the processing unit 610 can execute steps as shown in Figure 1 .

[0079] The storage unit 620 may include a readable medium in the form of a volatile storage unit, such as a random access storage unit (RAM) 6201 and / or a cache storage unit 6202, and may further include a read-only storage unit (ROM) 6203.

[0080] The storage unit 620 may also include a program / utility 6204 having a set (at least one) of program modules 6205. Such program modules 6205 include but are not limited to: an operating system, one or more application programs, other program modules, and program data. Each or some combination of these examples may include the implementation of a network environment.

[0081] The bus 630 may represent one or more of several types of bus structures, including a memory bus or a memory controller, a peripheral bus, a graphics acceleration port, a processing unit, or a local bus using any one of the multiple bus structures.

[0082] The electronic device 600 can also communicate with one or more external devices 700 (such as a keyboard, a pointing device, a Bluetooth device, etc.), and can also communicate with one or more devices that enable a user to interact with the electronic device 600, and / or communicate with any device that enables the electronic device 600 to communicate with one or more other computing devices (such as a router, a modem). Such communication can be carried out through the input / output interface 650. Moreover, the electronic device 600 can also communicate with one or more networks (such as a local area network, a wide area network, and / or a public network, such as the Internet) through the network adapter 660. The network adapter 660 can communicate with other modules of the electronic device 600 through the bus 630. It should be understood that, although not shown in the figure, other hardware and / or software modules can be used in combination with the electronic device 600, including but not limited to: microcode, device drivers, redundant processing units, external disk drive arrays, RAID systems, tape drives, and data backup storage platforms, etc.

[0083] Embodiment 4 The present invention also provides a storage medium, specifically a computer-readable storage medium, which is a memory device in a terminal device and is used to store programs and data. It can be understood that the computer-readable storage medium here can include both the built-in storage medium in the terminal device and, of course, the extended storage medium supported by the terminal device. It can be any tangible medium that contains or stores a program, and this program can be used by or in combination with an instruction execution system, device, or component. The computer-readable storage medium provides a storage space, and this storage space stores the operating system of the terminal. And, in this storage space, one or more instructions suitable for being loaded and executed by the processor are also stored. These instructions can be one or more computer programs (including program codes). It should be noted that more specific examples of the computer-readable storage medium here include: an electrical connection having one or more wires, a portable disk, a hard disk, a random access memory, a read-only memory, an erasable programmable read-only memory, an optical fiber, a portable compact disk read-only memory, an optical storage device, a magnetic storage device, or any suitable combination of the above.

[0084] The computer-readable storage medium also includes a data signal propagated in a baseband or as part of a carrier wave, which carries the readable program code. Such a propagated data signal can take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination of the above. The readable storage medium can also be any readable medium other than the readable storage medium, which can send, propagate, or transmit a program for use by or in conjunction with an instruction execution system, apparatus, or device. The program code contained on the readable storage medium can be transmitted using any appropriate medium, including but not limited to wireless, wired, optical fiber, radio frequency, etc., or any suitable combination of the above.

[0085] The program code for performing the operations of the present invention can be written in any combination of one or more programming languages. The programming languages include object-oriented programming languages such as Java, C++, etc., and also include conventional procedural programming languages such as the "C" language or similar programming languages. The program code can be executed entirely on the user's computing device, partially on the user's device, executed as an independent software package, partially on the user's computing device and partially on a remote computing device, or entirely on a remote computing device or server. In the case of a remote computing device, the remote computing device can be connected to the user's computing device through any type of network, including a local area network or a wide area network, or can be connected to an external computing device (e.g., by using an Internet service provider to connect through the Internet).

[0086] One or more instructions stored in the computer-readable storage medium can be loaded and executed by a processor to implement the corresponding steps of the method for predicting the propagation direction of seismic wave fields excited by high-speed rail sources in the above embodiments; one or more instructions in the computer-readable storage medium are loaded and executed by the processor as follows: Simplify the interaction between the high-speed rail wheel and the ground into a vertically downward impact force of the high-speed rail wheel on the ground, obtain the theoretical ground force function, decompose the theoretical force function to obtain the theoretical force function and the amplitude spectrum; let the Ricker wavelet be the deformation function of the force, convolve it with the theoretical force function and the amplitude spectrum to obtain the high-speed rail sub-source function, and then perform a Fourier transform on the high-speed rail sub-source function to obtain the high-speed rail sub-source function and the amplitude spectrum; construct a uniform linear array based on the sleepers or piers of the high-speed rail line, excite seismic wave signals at each element of the uniform linear array, and sum all the seismic wave signals to obtain the array factor; based on the high-speed rail sub-source function and the amplitude spectrum, combine the radiation pattern factors of seismic P-waves and S-waves excited by a single-point concentrated force, and multiply them by the array factor to obtain the radiation factors of seismic P-waves and S-waves of the high-speed rail source seismic wave field; draw the theoretical interference direction diagram of the seismic waves according to the array factor in the radiation factor to complete the prediction of the propagation direction of the seismic wave field.

[0087] In each of the embodiments provided in this application, the databases involved may include at least one of a relational database and a non-relational database. The non-relational database may include a distributed database based on blockchain, etc., and is not limited thereto. In each of the embodiments provided in this application, the processors involved may be general-purpose processors, central processors, graphics processors, digital signal processors, programmable logic devices, data processing logics based on quantum computing, etc., and are not limited thereto.

[0088] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are some, but not all, of the embodiments of the present invention. Usually, the components described and shown in the accompanying drawings here can be arranged and designed in various different configurations. Therefore, the following detailed description of the embodiments of the present invention provided in the drawings is not intended to limit the scope of the claimed present invention, but merely represents selected embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts fall within the protection scope of the present invention.

[0089] Simulation Experiment Design 1. Experimental Purpose Verify the accuracy of the prediction method for the propagation direction of seismic wave fields of high-speed rail vibration sources, analyze the influence of array parameters (such as sleeper spacing, number of array elements) and geological conditions on direction prediction, and evaluate the applicability of the model in complex scenarios.

[0090] 2. Experimental Setup 2.1 Parameter Configuration

[0091] 3. Simulation Process Step 1: Generate the theoretical source function and array factor Synthesize the Ricker wavelet time-domain signal according to the high-speed rail impact force model, and verify the consistency between its frequency-domain amplitude spectrum and the theoretical force function.

[0092] Construct uniform linear arrays with different element spacings, and calculate the variation of the array factor with the azimuth angle.

[0093] Step 2: Full waveform simulation of elastic wave fields In a uniform half-space and a layered medium model, simulate the seismic wave fields (P waves and S waves) excited by a single array element (sleeper) and a uniform linear array respectively, and record the data of surface and subsurface geophones.

[0094] Step 3: Direction Prediction and Verification Extract the radiation factor based on the simulation data and draw the theoretical interference direction diagram.

[0095] Compare the actual wavefield propagation direction (through wavefront arrival time / polarization analysis) with the predicted direction, and calculate the angular deviation.

[0096] 4. Key experimental results 4.1 Verification of direction prediction accuracy

[0097] The error between the predicted direction and the measured direction is ≤2°, verifying the effectiveness of the method; the error slightly increases in layered media, and the geological model coupling algorithm needs to be further optimized.

[0098] Please refer to Figure 2, which is the time-domain waveform diagram and spectrum diagram of the high-speed rail sub-source. Referring to the Chinese CRH380 high-speed rail, each carriage has two bogies, and each bogie has two sets of wheel pairs, so a single carriage has four sets of wheel pairs. Take 2.5 m, Take 14.5 m Take 80 m / s, Take 8, Figure 2 can be obtained when taking 25 m. The number of carriages can be seen from the time-domain waveform diagram. The spectral lines in the spectrum diagram have the characteristics of equally spaced narrow-band discreteness, which is consistent with the measured data. This shows that the selection of the high-speed rail sub-source in the present invention is reasonable.

[0099] Please refer to Figures 3(a)-(c), which are respectively the d array factor magnitude diagrams when the element spacing takes 1 m, 2 m, and 5 m, where the vehicle speed is 80 m / s and the wave speed is 1206 m / s. N is determined by the ratio of the total length of the high-speed rail vehicle to the element spacing d Here, the total length of the high-speed rail vehicle is 200 m. In Figure 3(a), there is only one bright curve. This indicates that when the element spacing d = 1 m and the frequency is lower than 75 Hz, the magnitude of the array factor at any angle is small. Further, when the element spacing d = 1 m and the frequency is lower than 75 Hz, the wavefield excited by the high-speed train source has no obvious interference phenomenon in the entire plane space. In Figures 3(b) and 3(c), as the element spacing increases, the number of bright curves in the figure also increases. This shows that when the element spacing is large, the wavefield excited by the high-speed train source is more likely to show a significant interference phenomenon in the plane space. The present invention analyzes the factors affecting the interference direction of the high-speed rail source wavefield through the specific array factor formula, and the parametric analysis makes the analysis of the high-speed rail source wavefield direction more perfect and more reasonable for direction prediction.

[0100] Please refer to Fig. 4(a) for the theoretical predicted radiation pattern in the Z-axis cross-section, and Fig. 4(b) for the theoretical predicted radiation pattern in the X-axis cross-section. and The representations in polar coordinates. The parameters are: d = 2m, v p = 2010 m / s, v s = 1206 m / s, v t = 80 m / s, N is determined by the ratio of the total length of the high-speed train to the element spacing d Here, the total length of the high-speed train is 200 m. In the Z-axis cross-section At this time so there is no P-wave theoretically in the Z-axis cross-section. In the X-axis cross-section so there are both P-waves and S-waves theoretically in the X-axis cross-section. In Fig. 4(a), the direction of the main lobe of the S-wave interference is at 127.5°, and in Fig. 4(b), the direction of the main lobe of the S-wave interference is at 128.5°, and the direction of the main lobe of the P-wave interference is at 156°. The present invention predicts the direction of the high-speed train seismic source wave field through parameterization of theoretical formulas, and uses different cross-sections to represent the theoretical predicted radiation patterns of the plane wave fields in each space, and the results are intuitive and clear.

[0101] Please refer to Figure 5 for the radiation pattern of the seismic wave excited by a single-point concentrated force source, which is obtained when the direction of the concentrated force is the negative Z-axis direction. The radiation patterns of the P-waves and S-waves excited by the single-point concentrated force source are drawn in the figure. The present invention is based on the source theory and takes into account the radiation pattern excited by the single-point force source, making the prediction method of the present invention more in line with the actual situation.

[0102] Please refer to Fig. 6(a) for the snapshot of the simulated wave field in the Z-axis cross-section, and Fig. 6(b) for the snapshot of the simulated wave field in the X-axis cross-section. The size of the simulation model is length × width × height = 800 m × 600 m × 400 m, v p = 2010 m / s, v s = 1206 m / s, v t = 80 m / s, the load interval (element spacing) d = 2m, and the load application direction is the negative Z-axis direction. There are obvious and stable main lobes in both Fig. 6(a) and Fig. 6(b). The wave field on the positive half-axis of the Z-axis in Fig. 6(b) is stronger than that on the negative half-axis, which is the manifestation of the load direction (the load application direction is the negative Z-axis direction). The present invention establishes a simulation model to verify the theoretical prediction direction, further demonstrating the persuasiveness of the prediction method of the present invention.

[0103] Please refer to Fig. 7(a), which is a combined diagram of the theoretical direction pattern in the Z-axis cross-section and the simulated wave field snapshot, and Fig. 7(b), which is a combined diagram of the theoretical direction pattern in the X-axis cross-section and the simulated wave field snapshot. It can be seen that whether it is the Z-axis cross-section or the X-axis cross-section, the predicted direction is consistent with the wave field direction of the simulation, thus verifying the rationality and correctness of the proposed method.

[0104] In summary, for the method and system for predicting the propagation direction of seismic wave fields excited by high-speed rail seismic sources of the present invention, the entire planar space and the distribution of the high-speed rail seismic sources themselves are considered, high-speed rail sub-seismic sources are reasonably selected, and the wave field direction of the high-speed rail seismic sources is predicted through parameterization of theoretical formulas. The theoretical prediction direction diagrams of the wave fields in each planar space are represented by different cross-sections. Considering the radiation pattern excited by a single-point force source, a simulation model is established to verify the theoretical prediction direction. The predicted direction is consistent with the wave field direction of the simulation, thus verifying the rationality and correctness of the proposed method.

[0105] Those skilled in the art can clearly understand that, for the convenience and brevity of description, only the above division of each functional unit and module is used as an example. In actual applications, the above functions can be allocated to different functional units and modules according to needs, that is, the internal structure of the device is divided into different functional units or modules to complete all or part of the functions described above. Each functional unit and module in the embodiment can be integrated into a processing unit, or each unit can exist physically alone, or two or more units can be integrated into one unit. The above integrated unit can be implemented in the form of hardware or in the form of a software functional unit. In addition, the specific names of each functional unit and module are only for the convenience of mutual distinction and do not limit the protection scope of the present application. The specific working process of the units and modules in the above system can refer to the corresponding process in the foregoing method embodiment and will not be elaborated here.

[0106] In the above embodiments, the descriptions of each embodiment have their own emphases. For the parts not detailed or recorded in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0107] Those of ordinary skill in the art can realize that the units and algorithm steps of each example described in combination with the embodiments disclosed in the present invention can be implemented by electronic hardware, or by a combination of computer software and electronic hardware. Whether these functions are executed in a hardware or software manner depends on the specific application and design constraints of the technical solution. Professional technicians can use different methods to implement the described functions for each specific application, but such implementation should not be considered to exceed the scope of the present invention.

[0108] In the embodiments provided by the present invention, it should be understood that the disclosed device / terminal and method can be implemented in other ways. For example, the device / terminal embodiments described above are merely illustrative. For example, the division of the modules or units is only a logical function division. In actual implementation, there may be other division methods. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the displayed or discussed coupling, direct coupling, or communication connection between each other can be through some interfaces. The indirect coupling or communication connection of the device or unit can be in electrical, mechanical, or other forms.

[0109] The units described as separate components may or may not be physically separated. The components displayed as units may or may not be physical units, that is, they can be located in one place or distributed to multiple network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0110] In addition, each functional unit in various embodiments of the present invention can be integrated into a processing unit, or each unit can exist physically alone, or two or more units can be integrated into one unit. The above integrated units can be implemented in the form of hardware or in the form of software functional units.

[0111] If the integrated module / 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, to implement all or part of the processes in the above method embodiments of the present invention, it can also be completed by a computer program instructing relevant hardware. The computer program can be stored in a computer-readable storage medium. When the computer program is executed by a processor, the steps of the above method embodiments can be implemented. Among them, the computer program includes computer program code, and the computer program code can be in the form of source code, object code, executable file, or some intermediate form, etc. The computer-readable medium can include: any entity or device capable of carrying the computer program code, recording medium, USB flash drive, mobile hard disk, magnetic disk, optical disc, computer memory, read-only memory (ROM), random access memory (RAM), electrical carrier signal, telecommunication signal, and software distribution medium, etc. It should be noted that the content included in the computer-readable medium can be appropriately increased or decreased according to the requirements of legislation and patent practice in the jurisdiction. For example, in some jurisdictions, according to legislation and patent practice, the computer-readable medium does not include electrical carrier signals and telecommunication signals.

[0112] This application is described with reference to the flowcharts and / or block diagrams of methods, apparatuses, and computer program products according to embodiments of the present application. It should be understood that each flow and / or block in the flowchart and / or block diagram, and the combination of flows and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing devices to generate a machine, such that the instructions executed by the processor of the computer or other programmable data processing devices generate means for implementing the functions specified in one flow Figure 1 one flow or multiple flows and / or blocks Figure 1 or means for implementing the functions specified in one block or multiple blocks.

[0113] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, such that the instructions stored in the computer-readable memory generate a manufactured article including instruction means that implement the functions specified in one flow Figure 1 one flow or multiple flows and / or blocks Figure 1 or means for implementing the functions specified in one block or multiple blocks.

[0114] These computer program instructions can also be loaded onto a computer or other programmable data processing device, such that a series of operation steps are executed on the computer or other programmable device to generate a computer-implemented process, and thus the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in one flow Figure 1 one flow or multiple flows and / or blocks Figure 1 or means for implementing the functions specified in one block or multiple blocks.

[0115] The above content is only to illustrate the technical idea of the present invention and cannot be used to limit the protection scope of the present invention. Any modification made on the basis of the technical solution according to the technical idea proposed by the present invention falls within the protection scope of the claims of the present invention.

Claims

1. A method for predicting the propagation direction of seismic wave fields excited by high-speed railway sources, characterized in that: The following steps are involved: The interaction between the high-speed rail wheel and the ground is simplified into the vertical downward impact force of the high-speed rail wheel on the ground, and the theoretical force function of the ground is obtained. The theoretical force function is decomposed to obtain the theoretical force function and amplitude spectrum; The Ricker wavelet is set as the deformation function of the force, and is convolved with the obtained theoretical force function and amplitude spectrum to obtain the high-speed rail sub-source function, and then the high-speed rail sub-source function is Fourier transformed to obtain the high-speed rail sub-source function and amplitude spectrum; A uniform linear array is constructed based on the sleepers or bridge piers of the high-speed railway line. A seismic wave signal is excited at each array element in the uniform linear array, and the array factor is obtained by summing up all the seismic wave signals. Based on the obtained high-speed railway sub-source function and amplitude spectrum, combined with the single-point concentrated force excitation seismic P-wave and S-wave radiation pattern factors, multiplied by the obtained array factor, the radiation factors of the P-wave and S-wave seismic wave field of the high-speed railway source are obtained; The theoretical interference direction diagram of seismic waves is drawn according to the array factor in the obtained radiation factor to complete the prediction of the propagation direction of the seismic wave field.

2. The method for predicting the propagation direction of the seismic wave field excited by the high-speed railway source according to claim 1 is characterized in that: The amplitude spectrum of the high-speed rail source function is: in, N t is the number of high-speed rail carriages, v t is the high-speed rail speed, C is a constant related to the mass of the high-speed rail, yes The Fourier transform form of the function is, yes The Fourier transform form of the function is, yes The Fourier transform form of the function is, is the distance between the second set of wheels and the first set of wheels in a single carriage. is the distance between the third wheel set and the first wheel set in a carriage. is the length of a single high-speed rail carriage. is the angular frequency.

3. The method for predicting the propagation direction of the seismic wave field excited by the high-speed railway source according to claim 2 is characterized in that: Theoretical ground force function for: in, d(k) It is k The distance from the first wheelset to the first wheelset, is the impulse function.

4. The method for predicting the propagation direction of the seismic wave field excited by the high-speed railway source according to claim 1 is characterized in that: Perform Fourier transform on the high-speed rail source function to obtain the high-speed rail source function and amplitude spectrum , as follows: in, is the high-speed rail source function at a single sleeper or bridge pier, is the Ricker wavelet, is a constant related to the mass of the high-speed rail, is the structural function of the second wheelset in a single carriage related to the first wheelset, is the structural function of the third wheelset in a single carriage related to the first wheelset, is a periodic input function related to the number of high-speed rail carriages, is the distance between the second set of wheels and the first set of wheels in a single carriage. is the high-speed rail speed, is the frequency, is the distance between the third wheel set and the first wheel set in a carriage. is the number of high-speed rail carriages, is the length of a single high-speed rail carriage. is the Fourier transform of the Ricker wavelet function.

5. The method for predicting the propagation direction of the seismic wave field excited by the high-speed railway source according to claim 4 is characterized in that: Perform Fourier transform on Ricker wavelet function: in, is the Fourier transform form of the Ricker wavelet function, The main frequency.

6. The method for predicting the propagation direction of the seismic wave field excited by the high-speed railway source according to claim 1, characterized in that: The array factor is: in, is the array factor that affects the wave field direction, is the frequency, is the spacing between adjacent array elements, It is the angle between the wave propagation direction from the source to the receiving point and the running direction of the high-speed rail. is the wave speed, is the high-speed rail speed, is the number of array elements.

7. The method for predicting the propagation direction of the seismic wave field excited by the high-speed railway source according to claim 6 is characterized in that: The seismic waves excited at each array element received at the far-field receiving point It is expressed as: in, is the time variable of the source function, is the wave path from the source to the receiving point.

8. The method for predicting the propagation direction of the seismic wave field excited by the high-speed railway source according to claim 1, characterized in that: The radiation factors of P and S waves in the high-speed railway source seismic wave field are: in, is the radiation pattern factor of the P-wave excited by the high-speed iron sub-source, is the P-wave radiation pattern factor, and is the array factor that affects the wave field direction, is the angle between the wave propagation direction and the concentrated force direction, is the frequency, is the spacing between adjacent array elements, is the velocity of the P wave, is the high-speed rail speed, It is the angle between the wave propagation direction from the source to the receiving point and the running direction of the high-speed rail. is the number of array elements, is the speed of the S wave, and are the P-wave and S-wave radiation pattern factors, respectively, is the radiation pattern factor of the seismic P wave excited by the high-speed iron sub-source.

9. The method for predicting the propagation direction of the seismic wave field excited by the high-speed railway source according to claim 8, characterized in that: The radiation pattern factors of seismic P and S waves excited by a single point concentrated force are: 。 10. A system for predicting the propagation direction of seismic wave fields excited by high-speed railway sources, characterized in that: include: The decomposition module simplifies the interaction between the high-speed rail wheel and the ground into the vertical downward impact force of the high-speed rail wheel on the ground, obtains the theoretical force function of the ground, and decomposes the theoretical force function to obtain the theoretical force function and amplitude spectrum; Convolution module, let Ricker wavelet be the deformation function of force, convolve with the obtained theoretical force function and amplitude spectrum to obtain the high-speed rail sub-source function, and then perform Fourier transform on the high-speed rail sub-source function to obtain the high-speed rail sub-source function and amplitude spectrum; The summation module constructs a uniform linear array based on the sleepers or bridge piers of the high-speed railway line, excites seismic wave signals at each array element in the uniform linear array, and sums all seismic wave signals to obtain the array factor; The factor module, based on the high-speed railway sub-source function and amplitude spectrum, combines the radiation pattern factors of the seismic P and S waves excited by the single-point concentrated force, and multiplies them with the array factor to obtain the radiation factors of the P and S waves of the high-speed railway source seismic wave field; The prediction module draws the theoretical interference direction diagram of the seismic wave according to the array factor in the obtained radiation factor, and completes the prediction of the propagation direction of the seismic wave field.

Images