Three-component artificial vibration synthesis method and device

Abstract

The invention provides a three-component artificial vibration synthesis method and device. The three-component artificial vibration synthesis method comprises the following steps: determining model parameter estimation values in three directions according to seismic parameters; determining a main energy group wavelet packet spectrum and a secondary energy group wavelet packet spectrum corresponding to the estimated values of the model parameters in the three directions; synthesizing the main energy group wavelet packet spectrum and the secondary energy group wavelet packet spectrum to obtain a total wavelet packet energy spectrum in three directions; wavelet packet inverse transformation is carried out on the total wavelet packet energy spectrum, a three-direction seismic oscillation time history is obtained to generate three-component artificial seismic waves, the distribution law of actually measured earthquakes is met, and reliable artificial seismic wave input can be provided for three-dimensional anti-seismic analysis of engineering structures.

Description

Technical Field

[0001] The present invention relates to the technical field of seismic wave synthesis, and in particular to a three-component artificial seismic motion synthesis method and device. Background Art

[0002] Currently, many regions are rich in hydropower resources but experience frequent seismic activity and complex geological structures. This poses numerous challenges to the construction of high-dam hydropower projects and other water conservancy facilities. Failure of these structures in an earthquake poses a safety hazard to downstream facilities. Therefore, scientific seismic safety analysis and assessment of these structures is crucial. Performance-based seismic engineering design often requires a large number of seismic wave inputs for time-history analysis to understand the structural response patterns under specific earthquake conditions. However, when conducting seismic design, measured ground motion records that meet design standards are often insufficient. Seismic network construction is still underway, and problems such as short observation histories and incomplete station site data exist. Due to the lack of usable measured ground motion records, artificial seismic waves are widely used in time-history analysis for seismic design of structures.

[0003] In addition, observations of several strong earthquakes in recent years have shown that the vertical component of strong earthquake motion near the earthquake source has a large amplitude and great potential destructive power. The vertical earthquake load has a significant impact on the dynamic response of dam structures, seismic slopes, and other structures. For seismic analysis of complex engineering structures, the role of multidimensional seismic waves must be considered. However, most current artificial seismic wave technologies use a method that modifies white noise through a filter to obtain a random process in the time domain, which cannot effectively quantify the joint time-frequency characteristics of seismic waves. Current simulation methods are mostly aimed at fitting response spectra or measured seismic motions, and cannot be applied to engineering design scenarios where a large number of seismic motions that conform to the measured distribution patterns under earthquake scenarios must be determined. Most artificial seismic wave technologies can only simulate the horizontal seismic wave time history of a single site and cannot simultaneously obtain three interrelated artificial seismic waves. Summary of the Invention

[0004] The main purpose of the embodiments of the present invention is to provide a three-component artificial seismic motion synthesis method and device to conform to the distribution law of measured earthquakes and provide reliable artificial seismic wave input for three-dimensional seismic analysis of engineering structures.

[0005] To achieve the above objectives, an embodiment of the present invention provides a three-component artificial ground motion synthesis method, comprising:

[0006] Determine the estimated values ​​of the model parameters in three directions according to the earthquake parameters;

[0007] Determine the primary energy group wavelet packet spectrum and the secondary energy group wavelet packet spectrum corresponding to the estimated values ​​of the model parameters in three directions;

[0008] synthesizing the wavelet packet spectrum of the primary energy group and the wavelet packet spectrum of the secondary energy group to obtain a total wavelet packet energy spectrum in three directions;

[0009] The total wavelet packet energy spectrum is subjected to an inverse wavelet packet transform to obtain a three-directional earthquake motion time history to generate a three-component artificial seismic wave.

[0010] In one embodiment, the earthquake parameters include target scene ground motion parameters, model coefficients, and a model parameter covariance matrix;

[0011] Determining the estimated values ​​of the model parameters in three directions according to the earthquake parameters includes:

[0012] Determining the model parameter median according to the target scene seismic parameter and the model coefficient;

[0013] The estimated values ​​of the model parameters in three directions are determined according to the model parameter median and the model parameter covariance matrix.

[0014] In one embodiment, determining the estimated values ​​of the model parameters in three directions according to the model parameter median and the model parameter covariance matrix includes:

[0015] Decomposing the model parameter covariance matrix to obtain a model parameter triangular matrix;

[0016] Obtaining model parameter residual random values ​​according to the model parameter triangular matrix and standard normal distribution random values;

[0017] The estimated values ​​of the model parameters in three directions are determined according to the model parameter median and the model parameter residual random value.

[0018] In one embodiment, the model parameter estimates include energy parameters, primary parameters, secondary parameters, and distribution random parameters;

[0019] The primary energy group wavelet packet spectrum and the secondary energy group wavelet packet spectrum corresponding to the three-direction model parameter estimation values ​​are determined as follows:

[0020] Determine the corresponding main energy group wavelet packet spectrum according to the energy parameter and the main parameter;

[0021] The corresponding secondary energy group wavelet packet spectrum is determined according to the secondary parameter and the distribution random parameter.

[0022] An embodiment of the present invention further provides a three-component artificial ground motion synthesis device, comprising:

[0023] A model parameter estimation module is used to determine the estimated values ​​of the model parameters in three directions according to the earthquake parameters;

[0024] A wavelet packet spectrum module is used to determine the primary energy group wavelet packet spectrum and the secondary energy group wavelet packet spectrum corresponding to the estimated values ​​of the model parameters in three directions;

[0025] A total wavelet packet energy spectrum module is used to synthesize the wavelet packet spectrum of the primary energy group and the wavelet packet spectrum of the secondary energy group to obtain a total wavelet packet energy spectrum in three directions;

[0026] The three-component artificial seismic wave module is used to perform wavelet packet inverse transformation on the total wavelet packet energy spectrum to obtain three-directional seismic motion time history to generate three-component artificial seismic waves.

[0027] In one embodiment, the earthquake parameters include target scene ground motion parameters, model coefficients, and a model parameter covariance matrix;

[0028] The model parameter estimation module includes:

[0029] a model parameter median unit, configured to determine the model parameter median according to the target scene seismic parameter and the model coefficient;

[0030] The model parameter estimation unit is used to determine the model parameter estimation values ​​in three directions according to the model parameter median and the model parameter covariance matrix.

[0031] In one embodiment, the model parameter estimation unit includes:

[0032] A decomposition subunit, configured to decompose the model parameter covariance matrix to obtain a model parameter triangular matrix;

[0033] A model parameter residual random value subunit is used to obtain a model parameter residual random value according to the model parameter triangular matrix and a standard normal distribution random value;

[0034] The model parameter estimation value subunit is used to determine the model parameter estimation values ​​in three directions according to the model parameter median and the model parameter residual random value.

[0035] In one embodiment, the model parameter estimates include energy parameters, primary parameters, secondary parameters, and distribution random parameters;

[0036] The wavelet packet spectrum module includes:

[0037] A main energy group wavelet packet spectrum unit is used to determine the corresponding main energy group wavelet packet spectrum according to the energy parameter and the main parameter;

[0038] The secondary energy group wavelet packet spectrum unit is used to determine the corresponding secondary energy group wavelet packet spectrum according to the secondary parameter and the distribution random parameter.

[0039] An embodiment of the present invention further provides an electronic device comprising a memory, a processor, and a computer program stored in the memory and running on the processor. When the processor executes the computer program, the steps of the three-component artificial ground motion synthesis method are implemented.

[0040] An embodiment of the present invention further provides a computer-readable storage medium storing a computer program. When the computer program is executed by a processor, the steps of the three-component artificial ground motion synthesis method are implemented.

[0041] An embodiment of the present invention further provides a computer program product, including a computer program / instruction, which implements the steps of the three-component artificial ground motion synthesis method when executed by a processor.

[0042] The three-component artificial seismic motion synthesis method and device of the embodiment of the present invention first determine the estimated values ​​of the model parameters in three directions according to the earthquake parameters to obtain the corresponding main energy group wavelet packet spectrum and the secondary energy group wavelet packet spectrum, and then synthesize the wavelet packet spectrum and perform the wavelet packet inverse transform to obtain the three-directional seismic motion time history to generate a three-component artificial seismic wave, which conforms to the distribution law of measured earthquakes and can provide reliable artificial seismic wave input for the three-dimensional seismic resistance analysis of engineering structures. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0044] Figure 1 is a flow chart of a three-component artificial ground motion synthesis method according to an embodiment of the present invention;

[0045] Figure 2 is a flow chart of a three-component artificial earthquake motion synthesis method according to another embodiment of the present invention;

[0046] Figure 3 is a flow chart of S101 in an embodiment of the present invention;

[0047] Figure 4 This is a flow chart of S202 in an embodiment of the present invention;

[0048] Figure 5 This is a flow chart of S102 in an embodiment of the present invention;

[0049] Figure 6 Schematic diagram of wavelet packet spectrum of main energy groups in an embodiment of the present invention;

[0050] Figure 7is a schematic diagram of the wavelet packet spectrum of the secondary energy group in an embodiment of the present invention;

[0051] Figure 8 Schematic diagram of the total wavelet packet spectrum in an embodiment of the present invention;

[0052] Figure 9 is a schematic diagram of the time history of artificial earthquake motion in the first horizontal direction in an embodiment of the present invention;

[0053] Figure 10 is a schematic diagram of the time history of artificial earthquake motion in the second horizontal direction in an embodiment of the present invention;

[0054] Figure 11 1 is a schematic diagram of the time history of artificial earthquake motion in the vertical direction according to an embodiment of the present invention;

[0055] Figure 12 is a structural block diagram of a three-component artificial ground motion synthesis device according to an embodiment of the present invention;

[0056] Figure 13 This is a schematic block diagram of the system structure of the electronic device 9600 according to an embodiment of the present application. DETAILED DESCRIPTION

[0057] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of the present invention.

[0058] Those skilled in the art will appreciate that the embodiments of the present invention may be implemented as a system, apparatus, device, method, or computer program product. Therefore, the present disclosure may be implemented in the following forms: entirely in hardware, entirely in software (including firmware, resident software, microcode, etc.), or in a combination of hardware and software.

[0059] To address the problem that existing artificial seismic wave technologies can only simulate the horizontal seismic wave time history of a single site and cannot simultaneously obtain three-component, interrelated artificial seismic waves, the present invention utilizes wavelet packet analysis to simulate three-component artificial seismic waves under given seismic conditions, providing reliable artificial seismic wave input for three-dimensional seismic analysis of engineering structures. The present invention is described in detail below with reference to the accompanying drawings.

[0060] Figure 1 4 is a flow chart of a three-component artificial ground motion synthesis method according to an embodiment of the present invention. Figure 2 FIG. 1 is a flow chart of a three-component artificial ground motion synthesis method according to another embodiment of the present invention. Figure 1 and Figure 2As shown in FIG, the three-component artificial ground motion synthesis method includes:

[0061] S101: Determine estimated values ​​of model parameters in three directions according to earthquake parameters.

[0062] The earthquake parameters include target scene seismic parameters, model coefficients and model parameter covariance matrix.

[0063] Figure 3 This is a flow chart of S101 in an embodiment of the present invention. Figure 3 As shown, S101 includes:

[0064] S201: Determine the model parameter median according to the target scene seismic parameter and the model coefficient.

[0065] In one embodiment, the median value of the model parameter is determined by the following formula:

[0066]

[0067] Among them, Y is the median of the model parameters, representing 13 model parameters, which are defined by wavelet packet decomposition of seismic waves; α, β1-β6 are model coefficients, and the 13 model parameters correspond to different model coefficients; M is the moment magnitude, R hyp is the focal distance, that is, the distance from the field point to the focal location of the fault; R rup is the fault distance, i.e. the shortest distance from the field point to the fault plane where the rupture occurs; h is the distance parameter considering distance saturation, V S30 is the site condition of the site, that is, the equivalent shear wave velocity in the soil layer 30m below the surface.

[0068] S202: Determine estimated values ​​of the model parameters in three directions according to the model parameter median and the model parameter covariance matrix.

[0069] Figure 4 This is a flow chart of S202 in an embodiment of the present invention. Figure 4 As shown, S202 includes:

[0070] S301: Decomposing the model parameter covariance matrix to obtain a model parameter triangular matrix.

[0071] S302: Obtaining model parameter residual random values ​​according to the model parameter triangular matrix and standard normal distribution random values.

[0072] In one embodiment, the model parameter residual random values ​​are as follows:

[0073] ε=L T y;

[0074] Σ=L T L;

[0075] Among them, ε is the residual random value of the model parameter, the vector y is 36 independent standard normal distribution random values, the matrix Σ is the given model parameter covariance matrix, and the matrix L T is the triangular matrix obtained by performing Cholesky decomposition on the covariance matrix Σ.

[0076] S303: Determine estimated values ​​of the model parameters in three directions according to the model parameter median and the model parameter residual random value.

[0077] In one embodiment, the estimated values ​​of the model parameters are as follows:

[0078] Z=μ+ε;

[0079] Among them, μ is the vector composed of the median value Y of the model parameters in sequence, and Z is the estimated value vector of the model parameters used to generate artificial seismic waves, including two horizontal directions and one vertical direction, for a total of three groups of model parameters

[0080] S102: Determine the primary energy group wavelet packet spectrum and the secondary energy group wavelet packet spectrum corresponding to the estimated values ​​of the model parameters in three directions.

[0081] The model parameter estimates include energy parameters, main parameters, secondary parameters and distribution random parameters. For example, two energy parameters E acc 、E(a) major ; Four time domain related parameters E(t) major and E(t) minor , E(f) major and E(f) minor ; Four frequency domain related parameters S(t) major and S(t) minor , characterizes the standard deviation of the wavelet packet coefficients in the time domain; S(f) major and S(f) minor ; Two time-frequency domain correlation parameters ρ(t,f) major and ρ(t,f) minor ; a distribution random parameter S(ξ). The main parameters include E(t) major 、E(f) major 、S(t) major 、S(f) major 、ρ(t,f) major ; Secondary parameters include E(t) minor 、E(f) minor 、S(t) minor 、S(f) minor 、ρ(t,f) minor .

[0082] For example, take a set of model parameter estimates vector Z H1 Perform dimensional reduction, for the correlation parameter ρ(t,f) major and ρ(t,f) minor , there are the following transformations:

[0083]

[0084] For all except the correlation parameter ρ(t,f) major ,ρ(t,f) minor The remaining parameters except the random parameter S(ξ) of the sum distribution are exponentially expressed with the natural logarithm e as the base.

[0085] The total wavelet packet spectrum is divided into the main energy group wavelet packet spectrum and the secondary energy group wavelet packet spectrum. The total amplitude of the main energy group wavelet packet spectrum is 70% of the total wavelet packet spectrum amplitude of the earthquake motion, and the total amplitude of the secondary energy group wavelet packet spectrum is 30% of the total wavelet packet spectrum amplitude of the earthquake motion. The total wavelet packet spectrum amplitude of the earthquake motion corresponds to the model parameter E acc .

[0086] Figure 5 This is a flow chart of S102 in an embodiment of the present invention. Figure 5 As shown, S102 includes:

[0087] S401: Determine a corresponding main energy group wavelet packet spectrum according to the energy parameter and the main parameter.

[0088] In one embodiment, the wavelet packet spectrum of the main energy group is as follows:

[0089]

[0090] in, is the wavelet packet coefficient distributed in the time-frequency domain (wavelet packet spectrum of the main energy group), the superscript i represents the position of the wavelet packet on the frequency axis in the wavelet packet spectrum, the subscript k represents the position on the time axis, and the subscript j is the number of wavelet packet decomposition layers. maj is the square value of the wavelet packet spectrum coefficient of the main energy group t k,maj and f i,maj are the position of the wavelet packet on the time axis and the frequency axis in the wavelet packet spectrum respectively; is the mean of the square values ​​of the wavelet packet spectral coefficients of the main energy group, which is the model parameter E(a) major ;M maj and Σ maj as follows:

[0091] M maj =[E(lnt k,maj )E(lnf i,maj )]

[0092]

[0093] Among them, E(lnt k,maj ) and E(lnf i,maj ) respectively correspond to the model parameters E(t) major and E(f) major ; S(lnt k,maj ) and S(lnf i,maj ) respectively correspond to the model parameters S(t) major and S(f) major ;cov(lnt k,maj ,lnf i,maj ) through ρ(t,f) major Calculated.

[0094] S402: Determine a corresponding secondary energy group wavelet packet spectrum according to the secondary parameter and the distribution random parameter.

[0095] In one embodiment, the wavelet packet spectrum of the secondary energy group is as follows:

[0096] X k =ln(t k ),Y i =ln(f i );

[0097]

[0098] Among them, t k and f i Respectively represent the position of the wavelet packet in the time axis and frequency axis of the wavelet packet spectrum; E(X), S(X), E(Y), S(Y), ρ(X,Y) correspond to E(t) in the model parameters respectively. minor ,S(t) minor , E(f) minor , S(f) minor and ρ(t,f) minor ,ξ k,i It is a standard lognormally distributed random number with the model parameter S(ξ) as the standard deviation.

[0099] S103: synthesizing the wavelet packet spectrum of the primary energy group and the wavelet packet spectrum of the secondary energy group to obtain a total wavelet packet energy spectrum in three directions.

[0100] In one embodiment, the total wavelet packet energy spectrum is as follows:

[0101]

[0102] in, is the total wavelet packet energy spectrum.

[0103] S104: performing an inverse wavelet packet transform on the total wavelet packet energy spectrum to obtain three-directional seismic motion time histories to generate three-component artificial seismic waves.

[0104] In specific implementation, the wavelet packet coefficients (total wavelet packet energy spectrum) Randomly assign symbols, and then perform inverse wavelet packet transform on the wavelet packet spectrum using the following formula:

[0105]

[0106] in, is the wavelet packet basis function, 2 N is the number of data points in the time series, and x(t) is the simulated earthquake motion time history. For the other two groups of earthquake motion components, the model parameters are and Z V , repeat the above steps, and finally obtain the artificial seismic wave time histories (three-directional seismic motion time histories) corresponding to the three sets of model parameters.

[0107] Figure 1 The execution subject of the three-component artificial earthquake motion synthesis method shown can be a computer. Figure 1 As can be seen from the process shown, the three-component artificial seismic motion synthesis method of the embodiment of the present invention first determines the estimated values ​​of the model parameters in three directions according to the seismic parameters to obtain the corresponding main energy group wavelet packet spectrum and the secondary energy group wavelet packet spectrum, and then synthesizes the wavelet packet spectrum and performs inverse wavelet packet transform to obtain the three-directional seismic motion time history to generate a three-component artificial seismic wave.

[0108] The specific embodiments of the present invention are as follows:

[0109] 1. Set the earthquake parameters of the near-field large earthquake scenario M=7, R rup =10km,R hyp =10km,V S30 =760m / s, and the number of three-component artificial seismic wave simulations is 100.

[0110] 2. According to the model parameter median prediction formula, calculate the median values ​​of the thirteen model parameters for the horizontal and vertical components of seismic motion.

[0111] 3. Calculate the corresponding number of model parameter residual random values ​​according to the correlation matrix, and add them to the corresponding model parameter median to obtain three sets of model parameter estimates, including two horizontal directions and one vertical direction, a total of three sets of model parameters.

[0112] 4. First, take the estimated values ​​of the model parameters in the first horizontal direction to calculate the wavelet packet spectrum of the main energy group, where the wavelet packet spectrum amplitude is independent of the wavelet packet position. The wavelet packet spectrum amplitude is randomly generated by exponential distribution, and its energy sum is 70% E acc . Figure 6Schematic diagram of the wavelet packet spectrum of the main energy groups in the embodiment of the present invention. Figure 6 As shown, the wavelet packet spectrum position is randomly generated by log-normal distribution. Figure 7 Schematic diagram of the wavelet packet spectrum of the secondary energy group in the embodiment of the present invention. Figure 7 As shown, the wavelet packet spectrum of the secondary energy group is calculated using the estimated values ​​of the model parameters. First, the wavelet packet spectrum value is obtained according to the binary log-normal distribution, and then fine-tuned according to the randomness parameter. Finally, the total energy of the wavelet packet spectrum is adjusted to 30%.

[0113] 5. Figure 8 Schematic diagram of the total wavelet packet spectrum in an embodiment of the present invention. Figure 9 : is a schematic diagram of the time history of artificial earthquake motion in the first horizontal direction in an embodiment of the present invention. Figure 8 and Figure 9 As shown in Figure 1, the wavelet packet spectrum of the primary energy group and the wavelet packet spectrum of the secondary energy group are combined to obtain the total wavelet packet spectrum. The total wavelet packet spectrum is subjected to inverse wavelet packet transform to generate the artificial ground motion time history.

[0114] 6. Figure 10 1 is a schematic diagram of the time history of artificial earthquake motion in the second horizontal direction in an embodiment of the present invention. Figure 11 FIG is a schematic diagram of the vertical artificial earthquake time history in an embodiment of the present invention. Figure 10-11 As shown, by selecting the model parameters of the second horizontal direction and the vertical direction and repeating steps 4-5, the corresponding artificial seismic waves can be obtained.

[0115] 7. Repeat steps 3-6 to generate 100 sets of three-component artificial seismic waves.

[0116] In summary, the three-component artificial seismic motion synthesis method of the embodiment of the present invention can simulate the three-component artificial seismic waves under the earthquake conditions of a given scenario, generate randomly distributed model parameter estimates based on the prediction model and correlation structure of the model parameters, use the main energy group wavelet packet spectrum and the secondary energy group wavelet packet spectrum to synthesize the total wavelet packet spectrum, and generate the seismic motion time history through the inverse wavelet packet transform. The generated artificial seismic motion conforms to the distribution law of measured earthquakes, providing reliable artificial seismic wave input for the three-dimensional seismic analysis of engineering structures.

[0117] Based on the same inventive concept, an embodiment of the present invention also provides a three-component artificial seismic motion synthesis device. Since the principle of solving the problem by this device is similar to that of the three-component artificial seismic motion synthesis method, the implementation of the three-component artificial seismic motion synthesis device can refer to the implementation of the three-component artificial seismic motion synthesis method, and the repeated parts will not be repeated.

[0118] Figure 12 FIG. 1 is a structural block diagram of a three-component artificial ground motion synthesis device according to an embodiment of the present invention. Figure 12 As shown, the three-component artificial ground motion synthesis device includes:

[0119] A model parameter estimation module is used to determine the estimated values ​​of the model parameters in three directions according to the earthquake parameters;

[0120] A wavelet packet spectrum module is used to determine the primary energy group wavelet packet spectrum and the secondary energy group wavelet packet spectrum corresponding to the estimated values ​​of the model parameters in three directions;

[0121] A total wavelet packet energy spectrum module is used to synthesize the wavelet packet spectrum of the primary energy group and the wavelet packet spectrum of the secondary energy group to obtain a total wavelet packet energy spectrum in three directions;

[0122] The three-component artificial seismic wave module is used to perform wavelet packet inverse transformation on the total wavelet packet energy spectrum to obtain three-directional seismic motion time history to generate three-component artificial seismic waves.

[0123] In one embodiment, the earthquake parameters include target scene ground motion parameters, model coefficients, and a model parameter covariance matrix;

[0124] The model parameter estimation module includes:

[0125] a model parameter median unit, configured to determine the model parameter median according to the target scene seismic parameter and the model coefficient;

[0126] The model parameter estimation unit is used to determine the model parameter estimation values ​​in three directions according to the model parameter median and the model parameter covariance matrix.

[0127] In one embodiment, the model parameter estimation unit includes:

[0128] A decomposition subunit, configured to decompose the model parameter covariance matrix to obtain a model parameter triangular matrix;

[0129] A model parameter residual random value subunit is used to obtain a model parameter residual random value according to the model parameter triangular matrix and a standard normal distribution random value;

[0130] The model parameter estimation value subunit is used to determine the model parameter estimation values ​​in three directions according to the model parameter median and the model parameter residual random value.

[0131] In one embodiment, the model parameter estimates include energy parameters, primary parameters, secondary parameters, and distribution random parameters;

[0132] The wavelet packet spectrum module includes:

[0133] A main energy group wavelet packet spectrum unit is used to determine the corresponding main energy group wavelet packet spectrum according to the energy parameter and the main parameter;

[0134] The secondary energy group wavelet packet spectrum unit is used to determine the corresponding secondary energy group wavelet packet spectrum according to the secondary parameter and the distribution random parameter.

[0135] In summary, the three-component artificial seismic motion synthesis method and device of the embodiment of the present invention first determine the estimated values ​​of the model parameters in three directions according to the seismic parameters to obtain the corresponding main energy group wavelet packet spectrum and the secondary energy group wavelet packet spectrum, and then synthesize the wavelet packet spectrum and perform the wavelet packet inverse transform to obtain the three-directional seismic motion time history to generate a three-component artificial seismic wave, which conforms to the distribution law of measured earthquakes and can provide reliable artificial seismic wave input for the three-dimensional seismic resistance analysis of engineering structures.

[0136] Figure 13 Schematic block diagram of the system structure of the electronic device 9600 according to an embodiment of the present application. Figure 13 As shown, the electronic device 9600 may include a central processing unit 9100 and a memory 9140; the memory 9140 is coupled to the central processing unit 9100. It is worth noting that the Figure 13 is exemplary; other types of structures may also be used to supplement or replace this structure to implement telecommunication functions or other functions.

[0137] In one embodiment, the three-component artificial ground motion synthesis method function may be integrated into the central processing unit 9100. The central processing unit 9100 may be configured to perform the following control:

[0138] Determine the estimated values ​​of the model parameters in three directions according to the earthquake parameters;

[0139] Determine the primary energy group wavelet packet spectrum and the secondary energy group wavelet packet spectrum corresponding to the estimated values ​​of the model parameters in three directions;

[0140] synthesizing the wavelet packet spectrum of the primary energy group and the wavelet packet spectrum of the secondary energy group to obtain a total wavelet packet energy spectrum in three directions;

[0141] The total wavelet packet energy spectrum is subjected to an inverse wavelet packet transform to obtain a three-directional earthquake motion time history to generate a three-component artificial seismic wave.

[0142] From the above description, it can be seen that the three-component artificial seismic motion synthesis method provided in this application first determines the estimated values ​​of the model parameters in three directions according to the earthquake parameters to obtain the corresponding main energy group wavelet packet spectrum and the secondary energy group wavelet packet spectrum, and then synthesizes the wavelet packet spectrum and performs an inverse wavelet packet transform to obtain the three-directional seismic motion time history to generate a three-component artificial seismic wave, which conforms to the distribution law of measured earthquakes and can provide reliable artificial seismic wave input for three-dimensional seismic analysis of engineering structures.

[0143] In another embodiment, the three-component artificial seismic motion synthesis device can be configured separately from the central processing unit 9100. For example, the three-component artificial seismic motion synthesis device can be configured as a chip connected to the central processing unit 9100, and the function of the three-component artificial seismic motion synthesis method can be realized through the control of the central processing unit.

[0144] like Figure 13 As shown, the electronic device 9600 may further include: a communication module 9110, an input unit 9120, an audio processor 9130, a display 9160, and a power supply 9170. It is worth noting that the electronic device 9600 does not necessarily have to include Figure 13 In addition, the electronic device 9600 may also include all components shown in Figure 13 For components not shown, reference may be made to the prior art.

[0145] like Figure 13 As shown, the central processing unit 9100 is sometimes also referred to as a controller or operation control, and may include a microprocessor or other processor device and / or logic device. The central processing unit 9100 receives input and controls the operation of various components of the electronic device 9600.

[0146] Memory 9140 can be, for example, one or more of a cache, flash memory, hard drive, removable media, volatile memory, non-volatile memory, or other suitable devices. It can store the aforementioned failure-related information and also store programs that execute the relevant information. The CPU 9100 can execute the programs stored in memory 9140 to implement information storage or processing.

[0147] The input unit 9120 provides input to the central processing unit 9100. The input unit 9120 may be, for example, a keypad or touch input device. The power supply 9170 is used to provide power to the electronic device 9600. The display 9160 is used to display objects such as images and text. The display may be, for example, an LCD display, but is not limited thereto.

[0148] The memory 9140 may be a solid-state memory, such as a read-only memory (ROM), a random access memory (RAM), or a SIM card. Alternatively, it may be a memory that retains information even when power is off, can be selectively erased, and is provided with more data. Examples of such memory are sometimes referred to as EPROMs. The memory 9140 may also be some other type of device. The memory 9140 includes a buffer 9141 (sometimes referred to as a buffer memory). The memory 9140 may include an application / function storage unit 9142 for storing application programs and function programs or processes for executing the operation of the electronic device 9600 by the central processing unit 9100.

[0149] The memory 9140 may also include a data storage unit 9143 for storing data, such as contacts, digital data, pictures, sounds, and / or any other data used by the electronic device. The driver storage unit 9144 of the memory 9140 may include various driver programs for communication functions of the electronic device and / or for executing other functions of the electronic device (such as messaging applications, address book applications, etc.).

[0150] The communication module 9110 is a transmitter / receiver 9110 that transmits and receives signals via an antenna 9111. The communication module (transmitter / receiver) 9110 is coupled to the central processor 9100 to provide input signals and receive output signals, which may be the same as in a conventional mobile communication terminal.

[0151] Based on different communication technologies, multiple communication modules 9110 can be provided in the same electronic device, such as a cellular network module, a Bluetooth module, and / or a wireless local area network module. The communication module (transmitter / receiver) 9110 is also coupled to a speaker 9131 and a microphone 9132 via an audio processor 9130 to provide audio output via the speaker 9131 and receive audio input from the microphone 9132, thereby implementing common telecommunication functions. The audio processor 9130 may include any suitable buffer, decoder, amplifier, etc. Furthermore, the audio processor 9130 is also coupled to the central processing unit 9100, enabling local recording via the microphone 9132 and playback of stored audio via the speaker 9131.

[0152] An embodiment of the present invention also provides a computer-readable storage medium capable of implementing all steps of the three-component artificial ground motion synthesis method in the above embodiment, where the execution subject is a server or a client. The computer-readable storage medium stores a computer program. When the computer program is executed by a processor, all steps of the three-component artificial ground motion synthesis method in the above embodiment are implemented. For example, when the processor executes the computer program, the following steps are implemented:

[0153] Determine the estimated values ​​of the model parameters in three directions according to the earthquake parameters;

[0154] Determine the primary energy group wavelet packet spectrum and the secondary energy group wavelet packet spectrum corresponding to the estimated values ​​of the model parameters in three directions;

[0155] synthesizing the wavelet packet spectrum of the primary energy group and the wavelet packet spectrum of the secondary energy group to obtain a total wavelet packet energy spectrum in three directions;

[0156] The total wavelet packet energy spectrum is subjected to an inverse wavelet packet transform to obtain a three-directional earthquake motion time history to generate a three-component artificial seismic wave.

[0157] In summary, the computer-readable storage medium of an embodiment of the present invention first determines the estimated values ​​of the model parameters in three directions according to the earthquake parameters to obtain the corresponding wavelet packet spectrum of the main energy group and the wavelet packet spectrum of the secondary energy group, and then synthesizes the wavelet packet spectrum and performs an inverse wavelet packet transform to obtain the three-directional seismic motion time history to generate a three-component artificial seismic wave, which conforms to the distribution law of measured earthquakes and can provide reliable artificial seismic wave input for three-dimensional seismic analysis of engineering structures.

[0158] An embodiment of the present invention also provides a computer program product capable of implementing all steps of the three-component artificial ground motion synthesis method in the above embodiment, where the execution subject is a server or a client. The computer program product includes a computer program / instructions. When the computer program / instructions are executed by a processor, all steps of the three-component artificial ground motion synthesis method in the above embodiment are implemented. For example, when the processor executes the computer program, the following steps are implemented:

[0159] Determine the estimated values ​​of the model parameters in three directions according to the earthquake parameters;

[0160] Determine the primary energy group wavelet packet spectrum and the secondary energy group wavelet packet spectrum corresponding to the estimated values ​​of the model parameters in three directions;

[0161] synthesizing the wavelet packet spectrum of the primary energy group and the wavelet packet spectrum of the secondary energy group to obtain a total wavelet packet energy spectrum in three directions;

[0162] The total wavelet packet energy spectrum is subjected to an inverse wavelet packet transform to obtain a three-directional earthquake motion time history to generate a three-component artificial seismic wave.

[0163] In summary, the computer program product of the embodiment of the present invention first determines the estimated values ​​of the model parameters in three directions according to the earthquake parameters to obtain the corresponding wavelet packet spectrum of the main energy group and the wavelet packet spectrum of the secondary energy group, and then synthesizes the wavelet packet spectrum and performs an inverse wavelet packet transform to obtain the three-directional seismic motion time history to generate a three-component artificial seismic wave, which conforms to the distribution law of measured earthquakes and can provide reliable artificial seismic wave input for the three-dimensional seismic analysis of engineering structures.

[0164] The various embodiments in this specification are described in a progressive manner. Similar portions between the various embodiments can be referenced to each other. Each embodiment focuses on the differences between the other embodiments. In particular, the hardware + program embodiments are generally similar to the method embodiments, so their description is relatively simple. For relevant portions, refer to the description of the method embodiments.

[0165] The foregoing description of this specification describes specific embodiments. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recited in the claims can be performed in an order different from that described in the embodiments and still achieve the desired results. Furthermore, the processes depicted in the accompanying drawings do not necessarily require the specific order shown or the sequential order to achieve the desired results. In certain embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0166] Although the present application provides method operation steps as described in the embodiments or flowcharts, more or fewer operation steps may be included based on conventional or non-creative work. The order of steps listed in the embodiments is only one way of executing the steps among many steps and does not represent the only execution order. When the actual device or client product is executed, it can be executed in sequence or in parallel according to the method shown in the embodiments or the drawings (for example, in a parallel processor or multi-threaded processing environment).

[0167] Although the present specification embodiment provides the method operation steps as described in the embodiment or flow chart, more or less operation steps may be included based on conventional or non-creative means. The order of steps listed in the embodiment is only one way in the order of execution of many steps and does not represent a unique execution order. When the device or terminal product in practice is executed, it can be performed in sequence or in parallel according to the method shown in the embodiment or the accompanying drawings (such as a parallel processor or a multi-threaded processing environment, or even a distributed data processing environment). The term "comprise", "comprising" or any other variant thereof is intended to cover non-exclusive inclusion, so that the process, method, product or equipment including a series of elements not only include those elements, but also include other elements not clearly listed, or also include elements inherent to such process, method, product or equipment. In the absence of more restrictions, it is not excluded that there are other identical or equivalent elements in the process, method, product or equipment including the elements.

[0168] For the convenience of description, the above devices are described in terms of functions divided into various modules. Of course, when implementing the embodiments of this specification, the functions of each module can be implemented in the same or multiple software and / or hardware, or the modules that implement the same function can be implemented by a combination of multiple sub-modules or sub-units, etc. The device embodiments described above are merely schematic. For example, the division of the units is only a logical function division. In actual implementation, there may be other division methods, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be an indirect coupling or communication connection through some interfaces, devices or units, which can be electrical, mechanical or other forms.

[0169] Those skilled in the art will also appreciate that, in addition to implementing the controller in pure computer-readable program code, it is entirely possible to implement the same functionality by logically programming the method steps in the form of logic gates, switches, application-specific integrated circuits, programmable logic controllers, embedded microcontrollers, and the like. Therefore, such a controller can be considered a hardware component, and the devices included therein for implementing various functions can also be considered structures within the hardware component. Alternatively, the devices for implementing various functions can be considered both software modules implementing the method and structures within the hardware component.

[0170] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowcharts and / or block diagrams, as well as combinations of processes and / or blocks in the flowcharts and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowcharts and / or block diagrams. Figure 1 a process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0171] These computer program instructions may 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, so that the instructions stored in the computer readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 a process or multiple processes and / or boxes Figure 1 The function specified in one or more boxes.

[0172] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing the instructions executed on the computer or other programmable device for implementing the process. Figure 1 a process or multiple processes and / or boxes Figure 1 The steps for the function specified in one or more boxes.

[0173] In a typical configuration, a computing device includes one or more processors (CPUs), input / output interfaces, network interfaces, and memory.

[0174] Memory may include non-permanent storage in a computer-readable medium, random access memory (RAM) and / or non-volatile memory in the form of read-only memory (ROM) or flash RAM. Memory is an example of a computer-readable medium.

[0175] Computer-readable media includes permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. The information can be computer-readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disc read-only memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassettes, magnetic disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined herein, computer-readable media does not include transitory computer-readable media (transitory media), such as modulated data signals and carrier waves.

[0176] Those skilled in the art will appreciate that the embodiments of this specification may be provided as methods, systems, or computer program products. Therefore, the embodiments of this specification may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Furthermore, the embodiments of this specification may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to magnetic disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0177] The various embodiments of this specification may be described in the general context of computer-executable instructions executed by a computer, such as program modules. Generally, program modules include routines, programs, objects, components, data structures, and the like that perform specific tasks or implement specific abstract data types. The embodiments of this specification may also be practiced in distributed computing environments where tasks are performed by remote processing devices connected via a communications network. In a distributed computing environment, program modules may be located in local and remote computer storage media, including storage devices.

[0178] The various embodiments in this specification are described in a progressive manner. Similar or identical parts between the various embodiments can be referenced to each other. Each embodiment focuses on the differences from the other embodiments. In particular, since the system embodiments are generally similar to the method embodiments, the description is relatively simple. For relevant parts, reference can be made to the description of the method embodiments. Throughout this specification, reference to the terms "one embodiment," "some embodiments," "example," "specific example," or "some examples" means that the specific features, structures, materials, or characteristics described in conjunction with that embodiment or example are included in at least one embodiment or example of the embodiments in this specification. In this specification, the schematic representations of these terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, those skilled in the art may combine and integrate the different embodiments or examples, and features of different embodiments or examples, described in this specification, without conflict.

[0179] The above description is merely an example of the embodiments of this specification and is not intended to limit the embodiments of this specification. For those skilled in the art, various modifications and variations of the embodiments of this specification are possible. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the embodiments of this specification shall be included within the scope of the claims of the embodiments of this specification.

Claims

1. A three-component artificial ground motion synthesis method, characterized in that: include: Determine the estimated values ​​of the model parameters in three directions according to the earthquake parameters; Determine the primary energy group wavelet packet spectrum and the secondary energy group wavelet packet spectrum corresponding to the estimated values ​​of the model parameters in three directions; synthesizing the wavelet packet spectrum of the primary energy group and the wavelet packet spectrum of the secondary energy group to obtain a total wavelet packet energy spectrum in three directions; The total wavelet packet energy spectrum is subjected to an inverse wavelet packet transform to obtain a three-directional earthquake motion time history to generate a three-component artificial seismic wave.

2. The three-component artificial ground motion synthesis method according to claim 1, characterized in that: The earthquake parameters include target scene earthquake parameters, model coefficients and model parameter covariance matrix; Determining the estimated values ​​of the model parameters in three directions according to the earthquake parameters includes: Determining the model parameter median according to the target scene seismic parameter and the model coefficient; The estimated values ​​of the model parameters in three directions are determined according to the model parameter median and the model parameter covariance matrix.

3. The three-component artificial ground motion synthesis method according to claim 2, characterized in that: Determining the estimated values ​​of the model parameters in three directions according to the model parameter median and the model parameter covariance matrix includes: Decomposing the model parameter covariance matrix to obtain a model parameter triangular matrix; Obtaining model parameter residual random values ​​according to the model parameter triangular matrix and standard normal distribution random values; The estimated values ​​of the model parameters in three directions are determined according to the model parameter median and the model parameter residual random value.

4. The three-component artificial ground motion synthesis method according to claim 1, characterized in that: The model parameter estimates include energy parameters, main parameters, secondary parameters and distribution random parameters; The primary energy group wavelet packet spectrum and the secondary energy group wavelet packet spectrum corresponding to the three-direction model parameter estimation values ​​are determined as follows: Determine the corresponding main energy group wavelet packet spectrum according to the energy parameter and the main parameter; The corresponding secondary energy group wavelet packet spectrum is determined according to the secondary parameter and the distribution random parameter.

5. A three-component artificial earthquake motion synthesis device, characterized in that: include: A model parameter estimation module is used to determine the estimated values ​​of the model parameters in three directions according to the earthquake parameters; A wavelet packet spectrum module is used to determine the primary energy group wavelet packet spectrum and the secondary energy group wavelet packet spectrum corresponding to the estimated values ​​of the model parameters in three directions; A total wavelet packet energy spectrum module is used to synthesize the wavelet packet spectrum of the primary energy group and the wavelet packet spectrum of the secondary energy group to obtain a total wavelet packet energy spectrum in three directions; The three-component artificial seismic wave module is used to perform wavelet packet inverse transformation on the total wavelet packet energy spectrum to obtain three-directional seismic motion time history to generate three-component artificial seismic waves.

6. The three-component artificial earthquake vibration synthesis device according to claim 5, characterized in that: The earthquake parameters include target scene earthquake parameters, model coefficients and model parameter covariance matrix; The model parameter estimation module includes: a model parameter median unit, configured to determine the model parameter median according to the target scene seismic parameter and the model coefficient; The model parameter estimation unit is used to determine the model parameter estimation values ​​in three directions according to the model parameter median and the model parameter covariance matrix.

7. The three-component artificial earthquake vibration synthesis device according to claim 6, characterized in that: The model parameter estimation unit includes: a decomposition subunit, configured to decompose the model parameter covariance matrix to obtain a model parameter triangular matrix; A model parameter residual random value subunit is used to obtain a model parameter residual random value according to the model parameter triangular matrix and a standard normal distribution random value; The model parameter estimation value subunit is used to determine the model parameter estimation values ​​in three directions according to the model parameter median and the model parameter residual random value.

8. The three-component artificial earthquake vibration synthesis device according to claim 5, characterized in that: The model parameter estimates include energy parameters, main parameters, secondary parameters and distribution random parameters; The wavelet packet spectrum module includes: A main energy group wavelet packet spectrum unit is used to determine the corresponding main energy group wavelet packet spectrum according to the energy parameter and the main parameter; The secondary energy group wavelet packet spectrum unit is used to determine the corresponding secondary energy group wavelet packet spectrum according to the secondary parameter and the distribution random parameter.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and running on the processor, characterized in that: When the processor executes the computer program, the steps of the three-component artificial ground motion synthesis method according to any one of claims 1 to 4 are implemented.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the three-component artificial ground motion synthesis method according to any one of claims 1 to 4 are implemented.

11. A computer program product comprising a computer program / instructions, characterized in that When the computer program / instructions are executed by a processor, the steps of the three-component artificial ground motion synthesis method according to any one of claims 1 to 4 are implemented.

Images