CN114578417A - Pulse type earthquake motion fitting method based on wavelet adjustment - Google Patents
Pulse type earthquake motion fitting method based on wavelet adjustment Download PDFInfo
- Publication number
- CN114578417A CN114578417A CN202210229458.9A CN202210229458A CN114578417A CN 114578417 A CN114578417 A CN 114578417A CN 202210229458 A CN202210229458 A CN 202210229458A CN 114578417 A CN114578417 A CN 114578417A
- Authority
- CN
- China
- Prior art keywords
- wavelet
- formula
- pulse
- time course
- seismic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000001228 spectrum Methods 0.000 claims abstract description 83
- 230000001133 acceleration Effects 0.000 claims abstract description 82
- 230000004044 response Effects 0.000 claims abstract description 40
- 230000010355 oscillation Effects 0.000 claims abstract description 22
- 238000000354 decomposition reaction Methods 0.000 claims description 22
- 238000004088 simulation Methods 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000000694 effects Effects 0.000 claims description 9
- 238000006073 displacement reaction Methods 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 5
- 238000009795 derivation Methods 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- 238000000605 extraction Methods 0.000 claims 1
- 230000009466 transformation Effects 0.000 abstract description 6
- 238000006243 chemical reaction Methods 0.000 description 9
- 230000007246 mechanism Effects 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000013016 damping Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000002194 synthesizing effect Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000005253 gamme decay Effects 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000001308 synthesis method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a pulse type earthquake motion fitting method based on a wavelet adjustment method, which comprises the following steps: fitting a pulse type earthquake motion wavelet function; substituting the field condition information to obtain a target acceleration response spectrum; obtaining the mother wave to calculate the acceleration response spectrum; adjusting the mother wave by adopting wavelet transformation based on the target acceleration response spectrum to obtain wavelets; acquiring a father wave by adopting an equivalent pulse velocity model; respectively solving wavelet coefficients of the wavelet and the parent wave; adjusting the wavelet coefficient of the wavelet to obtain an adjusted wavelet coefficient; and reconstructing the adjusted wavelet coefficient to obtain the near fault pulse seismic oscillation which is more in line with the actual situation.
Description
Technical Field
The invention provides an impulse type earthquake motion fitting method based on a wavelet adjustment method, and belongs to the field of impulse type earthquake motion fitting.
Background
The pulse type earthquake motion generated at the present stage is earthquake motion synthesized based on Fourier transform, and the earthquake motion can realize non-stationarity in a time domain through uniform modulation on amplitude in the time domain, but cannot realize non-stationarity in a frequency domain. The uniform modulation is realized by multiplying an envelope function by an acceleration time course, and the method cannot simulate the seismic motion containing velocity pulses. The time course of seismic oscillation is a time-frequency non-stationary random process with abrupt change characteristics. The seismic motion synthesis method based on Fourier transform cannot ensure that the synthesized seismic motion is time-frequency non-stationary. In order to be able to simulate the non-stationary nature of the time-frequency behavior of seismic oscillations, new signal processing tools must be sought. Therefore, the invention provides a pulse type seismic oscillation fitting method based on a wavelet adjustment method.
With the continuous development of national economic and technological levels, more and more large-sized buildings are built as lifelines in the near fault area and even across the fault in recent years. The natural vibration period of such a building (structure) is often larger than that of a general building (structure), and is easily affected by low-frequency excitation. The impulse seismic motion is a seismic motion containing a large amount of low-frequency components, and therefore, it is necessary to consider the influence of the impulse seismic motion on the structure. So far, the seismic record of the near fault pulse type seismic motion is insufficient, and for the area lacking the seismic record, no measured seismic record can be used as seismic motion input when the seismic design of a large building is carried out. But the artificial earth vibration can generate a series of earthquake motion meeting the condition of the proposed site in a short time, which can effectively make up for the deficiency. Therefore, it is necessary to study a method of fitting a near-fault pulse type seismic oscillation. The artificial synthetic seismic motion can not only obtain seismic waves meeting various site conditions, but also promote the development of structural seismic design to a certain extent. Therefore, it is very important to research a method capable of synthesizing earthquake motion which is as consistent as possible with the actual situation.
The method comprises the following steps of Wangyuan navigation, near fault area division and near fault velocity pulse type earthquake motion simulation [ D ]. southwest transportation university 2015:
(1) determining a target reaction spectrum S (T);
the target acceleration response spectrum was determined according to 5.2.1 in road engineering seismic Specification (JTG B02-2013), and is of the form shown below:
in the formula: t isg-a characteristic period of the reaction spectrum;
t-structure natural vibration period;
s (T) -target response spectrum;
Smax-reaction spectrum maximum.
(2) Calculating a target power spectrum S (omega);
according to the relation between the reaction spectrum and the power spectrum conversion which is proposed currently, the form is shown as the following formula.
Acceleration response spectrum from target
In the formula: sa T(xi, ω) -target acceleration response spectrum;
s (ω) -power spectrum;
ξ -damping ratio;
omega-structural natural frequency of vibration;
t-structure natural vibration period;
gamma-probability of surging.
(3) Calculating a Fourier amplitude spectrum A (omega);
the mathematical relationship between the power spectrum and the Fourier amplitude spectrum is shown in the formula.
In the formula: a (omega) -Fourier amplitude spectrum
(4) Calculating a phase spectrum;
generating a phase spectrum uniformly distributed over-2 pi to 2 pi
(5) Obtaining a steady acceleration time course by IFFT;
the earthquake motion X is simulated by adopting a trigonometric series model shown as the following formulaa。
Xa=g(t)·xa(t)
ωk-the circular frequency;
A(ωk) -an amplitude spectrum;
n-number of frequency separation points in the frequency domain of the target response spectrum;
g (t) -an intensity envelope function;
xa(t) -smooth acceleration time course;
Xaacceleration time course
Calculating to obtain the time course of steady acceleration by using fast Fourier inverse transformation
(6) Multiplying by an envelope function to generate a non-stationary acceleration time interval;
to convert a stationary acceleration time course into a non-stationary acceleration time course, the generated stationary acceleration time course is non-uniformly modulated by multiplying by an intensity envelope function. The mathematical expression of the envelope function is shown as
Wherein: t is1-the start time of the plateau;
T2-end time of the plateau;
c-rate of change of decay phase.
(7) Fitting and iterating with a target acceleration response spectrum to obtain an acceleration time course;
in order to match the response spectrum of the non-stationary acceleration time course with the target acceleration response spectrum, iterative calculation needs to be performed on the non-stationary acceleration time course, so that the generated response spectrum of the non-stationary acceleration time course approaches the target spectrum.
The specific iterative calculation steps are as follows:
a. calculating the ratio of the target acceleration response spectrum to the non-stationary acceleration time course;
b. performing fast Fourier transform calculation on the non-stationary acceleration time course to obtain a phase spectrum and a Fourier amplitude spectrum, and multiplying the ratio calculated in the step a by the Fourier amplitude spectrum;
c. combining the adjusted Fourier amplitude spectrum with the phase spectrum, and calculating to obtain a stable acceleration time interval through fast Fourier inverse transformation;
d. multiplying by an intensity envelope function to calculate a non-stationary acceleration time course;
e. calculating the average relative error Em between the response spectrum of the non-stationary acceleration time course and the target acceleration response spectrum, judging whether the average relative error Em is less than 5%, if the average relative error Em is less than 5%, stopping iterative calculation, and otherwise returning to the step a.
After n times of iterative computation, the average relative error between the reaction spectrum of the non-stationary acceleration time course and the target acceleration reaction spectrum is less than 5 percent, and the iterative computation is stopped. Finally obtaining a calculated acceleration reaction spectrum curve and a target acceleration reaction spectrum curve
(8) Fourier transform is carried out on the acceleration time course finally generated in the step (7) to obtain a Fourier spectrum, and the frequency is smaller than frIs set to 0
(9) Combining the Fourier amplitude spectrum finally generated in the step (8) with the phase spectrum generated in the step (4), and calculating and generating a high-frequency acceleration time interval through inverse fast Fourier transform
(10) Simulating a low-frequency speed pulse time course by adopting an equivalent speed pulse model:
where ω (t) is an envelope function and is calculated by the following equation (3-9).
Wherein:
Vp-peak value of pulse
fpFrequency, 1/Tp,TpIs a pulse period
t1-the moment of occurrence of the peak of the cosine function
Self-oscillation period of T-structure
Gamma-decay rate
t0-time of occurrence of peak of envelope function
The low-frequency speed pulse time interval is derived to obtain a low-frequency acceleration time interval
(11) High and low frequency components are superposed to generate near fault pulse type earthquake motion
And (4) superposing the high-frequency acceleration time interval obtained in the step (9) and the low-frequency pulse acceleration time interval obtained in the step (10) to generate the near-fault pulse type seismic oscillation.
The literature [ wangynavigation, 2015] method for synthesizing seismic oscillation has the following disadvantages:
(1) the impulse type earthquake motion generated by the harmonic method based on Fourier transform cannot reflect the frequency domain non-stationarity of the earthquake motion. In the simulation of the high-frequency components, fourier transform methods are used in the generation of the acceleration time intervals, which are only suitable for processing stationary signals. However, the seismic waves are generally non-stationary signals, and the frequency of the seismic waves changes along with the time change, so that the frequency domain non-stationarity of the seismic motion cannot be reflected by the Fourier transform, and therefore the pulse type seismic motion generated by the harmonic wave method based on the Fourier transform cannot completely meet the seismic motion of the actual situation.
(2) The single equivalent velocity pulse model is used for simulating the low-frequency component of the artificial near fault pulse type earthquake motion, and the influence of different fault fracture mechanisms is not considered. A great deal of research shows that different fault fracture mechanisms can generate different types of impulse type earthquake motion, so that the low-frequency components of the artificial near fault impulse type earthquake motion simulated by using a single equivalent velocity impulse model are not comprehensive enough, and the earthquake motion generated by using the single equivalent velocity impulse model cannot completely meet the earthquake motion of the actual situation.
Disclosure of Invention
The invention provides a pulse type seismic motion fitting method based on a wavelet adjustment method, and the pulse type seismic motion time course obtained by fitting by the method can be more reasonably matched with the actual near fault seismic motion characteristic.
Step S1.1: determining a scale vector S from equation (1)j
In the formula:
Sj-scale vector
n0Controlling the scale vector SjParameters of the range
m-scale vector SjNumber of points
In the formula:
Zeta-wavelet basis function attenuation coefficient
Omega-wavelet basis function frequency coefficient
pi-time vector
tkWavelet basis function time parameter
Step S2: according to the impulse type earthquake motion wavelet basis function obtained in the step S1Substituting the field condition information to obtain a target acceleration response spectrum Sa(Tj)]Target spectra
Step S3: calculating acceleration response spectrum (S) of original seismic oscillation, namely mother wave MWa(Tj)]Calculating the spectrum
Step S4: obtaining an adjustment coefficient gamma according to the formula (3)
In the formula: gamma-adjustment factor
[Sa(Tj)]Target spectraTarget acceleration response spectrum
[Sa(Tj)]Calculating the spectrumCalculating an acceleration response spectrum
Step S5: obtaining a target acceleration response spectrum [ S ] according to the formula (4)a(Tj)]Target spectraAnd calculating acceleration response spectrum Sa(Tj)]Calculating the spectrumError E ofrror。
In the formula:
Errorerror of target acceleration response spectrum from calculated acceleration response spectrum
m-scale vector SjNumber of points
[Sa(Tj)]Target spectrumTarget acceleration reactionStress spectrum
[Sa(Tj)]Calculating the spectrumCalculating an acceleration response spectrum
Step S6: obtaining wavelet coefficients C (s, p) according to equation (5)
In the formula:
C(sj,pi) -wavelet coefficients
Delta t-time interval of seismic oscillation acceleration time course
tkWavelet basis function time parameter
Sj-scale vector
piTime vector f (t)k) Primary seismic oscillation
m-scale vector SjNumber of points
M-number of data points of original seismic oscillation
Step S7: obtaining detail function D (s, t) in wavelet transform according to equation (6)
In the formula:
D(sj,tk) -detail function
Δp-wavelet function time interval
Step S8: obtaining a wavelet SW according to equation (7)
In the formula:
SW-wavelet
Step S9: the average relative error E between the calculated acceleration response spectrum and the target acceleration response spectrum is obtained in step S5rrorAnd judges whether it is less than 5%. If the average relative error is greater than 5%, returning to step S4; and if the average relative error is less than or equal to 5%, outputting the wavelets SW of the pulse seismic motion time course.
Step S10: the pulse period T is obtained according to the formula (8)pThe peak value V of the pulse is obtained from the formula (9)pThe peak time t of the pulse is obtained from the equation (10)1,V
ln(Tp)=-6.45+1.11MwFormula (8)
ln(Vp)=3.680+0.065Mw+0.025ln (R) formula (9)
ln(tl,v)=1.35Mw-6.88 formula (10)
Wherein: t isp-pulse period
Vp-peak value of pulse
t1,V-pulse peak time
Distance of R-fault
MwMoment-magnitude
Step S11: speed time course v for simulating sliding impulse type near fault impulse type earthquake to provide near fault impulse type earthquakegAAs shown in equation (11).
Method for simulating forward directivity effect to provide time course v of near fault pulse type seismic velocitygBAs shown in equation (12).
νgB(t)=Vpsin(ωpt) 0≤t≤TpFormula (12)
Wherein:
νgAsimulation sliding impulse type near fault impulse type seismic velocity time course
νgBSimulation of forward directional effect near fault pulse type seismic velocity time course
Vp-peak value of pulse
TpPulse period
ωpPulse frequency of fromp=2π/TpDetermining
Step S12: velocity time course v obtained in step S11gADerivation is carried out to obtain acceleration time course alphagAAs shown in formula (13). Velocity time course v obtained in step S11gAIntegral calculation to obtain displacement time interval dgA. As shown in equation (14).
In the formula: t ispDetermined by dividing the maximum of equations (11) and (14), i.e.:
velocity time interval v obtained in step S11gBDerivation is carried out to obtain acceleration time course alphagBAs shown in formula (15). Velocity time course v obtained in step S11gBIntegral calculation to obtain displacement time interval dgBAs shown in formula (16).
αgB(t)=ωpVpcos(ωpt),0≤t≤TpFormula (15)
In the formula:
αgAsimulation sliding impulse type near fault impulse type earthquake acceleration time course
dgASimulation sliding impulse type near fault impulse type earthquake displacement time course
Alpha gB-simulated forward directional effect near fault pulse type seismic acceleration time course
dgBSimulation forward directivity effect near fault pulse type seismic displacement time course
Vp-peak value of pulse
Tp-pulse period
ωpPulse frequency of fromp=2π/TpDetermining
TpDetermined by dividing the maximum of equations (12) and (16), i.e.:
step S13: father wave FW for obtaining pulse type earthquake motion time range
Step S14: calculating the pulse period Tp
ln(Tp)=-6.45+1.11MwFormula (17)
In the formula: t isp-pulse period
Step S15: determining the number of layers n of a wavelet decomposition
In the formula: f. ofpPulse frequency
number of layers for n-wavelet transform decomposition
foriginalFrequency band range of the signal
Delta t-time interval of seismic dynamic acceleration time interval
Step S16: decomposing the parent wave into n layers by multi-scale discrete wavelet transform as shown in formula (18)
In the formula: FW-father wave
FWcD1,…,FWcDn-each frequency band corresponding to the detail component of the parent wave after wavelet decomposition
FWcAnFrequency band corresponding to approximate component of parent wave after wavelet decomposition
n-number of layers of wavelet decomposition
Step S17: the parent wave wavelet coefficient LW is determined by step S16cAn
LWcAn=FWcAn
In the formula:
LWcAn-wavelet coefficient of parent wave
Step S18: decomposing the wavelet into n layers by multi-scale discrete wavelet transform as shown in formula (18)
In the formula: SW-wavelet
SWcD1,…,SWcDn-each frequency band corresponding to a detail component of a wavelet decomposition
SWcAn-bands corresponding to approximate components of wavelets after wavelet decomposition
n-number of layers of wavelet decomposition
Step S19: wavelet coefficients HW of wavelets are determined by step S18cAn
HWcAn=SWcAn
In the formula:
SWcAnwavelet coefficients of wavelets
Step S20: solving wavelet adjustment coefficient beta according to equation (20)coef
βcoef=FWcAn/SWcAnFormula (20)
In the formula: beta is acoefWavelet adjustment coefficient
Step S21: obtaining the adjusted wavelet coefficient SW 'according to the formula (21)'cAn
SW′cAn=βcoef·SWcAnFormula (21)
In the formula: SW'cAn-adjusted wavelet coefficients
Step S22: the adjusted wavelet coefficient SW 'is obtained according to the formula (21)'cAnReplacing wavelet coefficients SW after wavelet SW decompositioncAnPerforming wavelet change reconstruction to obtain the near fault pulse type seismic oscillation AW
AW-near fault pulse type seismic motion.
(1) The impulse type earthquake motion generated based on the wavelet transformation can reflect the frequency domain non-stationarity of the earthquake motion. The wavelet transformation can acquire frequency by constructing wavelet basis functions with different characteristics to realize non-stationary signal processing, so that the frequency domain non-stationarity of earthquake motion is reflected, and the pulse earthquake motion generated based on the wavelet method is more suitable for the earthquake motion of actual conditions.
(2) The invention provides an impulse type earthquake motion fitting method based on a wavelet adjustment method for the first time. The study on the seismic dynamic acceleration time-history curve generated based on the wavelet adjustment method finds that the original time domain characteristics of the mother wave in the time domain are reserved after the mother wave is adjusted through discrete wavelet transformation. In the frequency domain, the spectral characteristics of impulsive seismic motion combine the characteristics of parent and parent waves. Therefore, the characteristics of the pulse type earthquake motion based on the wavelet adjustment method and the actual near fault earthquake motion are more consistent.
(3) The invention provides an optimal number of decomposition layers in the generation of a near fault impulse type earthquake. There is a correlation between the number of decomposition layers and the pulse period. The invention passes through the sum pulse period TpHarmonic magnitude MwThe statistical analysis shows that the pulse period and the magnitude have strong correlation, and the invention is just about the pulse period TpHarmonic vibration level MwThe statistical rule obtained by least square fitting is more reasonable and closer to the real pulse period value. Therefore, it isAccording to the invention, more accurate wavelet coefficients are obtained by determining the optimal decomposition layer number, and seismic oscillation more in line with actual conditions can be generated.
(4) The invention uses different equivalent velocity pulse models to generate the father wave of the near fault pulse type earthquake motion, and considers the influence of different fault fracture mechanisms. It has been found through extensive research that different fault-breaking mechanisms produce different types of impulse-type seismic motion. Therefore, the earthquake motion generated by using different equivalent velocity pulse models can better accord with the earthquake motion of the earthquake of the actual situation.
Drawings
FIG. 1 is a flow chart of the present invention.
FIG. 2 is a time course curve of the acceleration of the parent wave.
FIG. 3 is a graph of parent wave acceleration time course.
Fig. 4 is a graph of a near fault pulse type acceleration time course.
Detailed Description
Seismic motion fitting parameters:
(1) example parameters are as follows: the engineering field category is II type field; the earthquakes are grouped into a third group (Tg ═ 0.45 s); the seismic fortification intensity is 7 degrees; the moment-vibration level is 6.5; the structural damping ratio is 0.05; the 50-year overrun probability is 10% (recurrence period 475 years), and the corresponding design basic earthquake motion peak acceleration is 0.15 g; the class A highway bridge E1 is used for seismic action and is in the form of near fault pulse seismic motion.
The mother wave is an El Mayor-Cucapah _ Mexico wave. The specific information is shown in table 1, and the time course curve of the acceleration of the parent wave is shown in fig. 1.
TABLE 1
The father wave is determined by an equivalent velocity pulse model, and the values of model parameters are as follows: t20 s, Vp=60cm/s、t0=5s、t15.6s, γ is 3; the acceleration time course curve of the parent wave is shown in figure 2.
(2) And (3) fitting: taking the approximate fault pulse type seismic oscillation of an equivalent velocity pulse model of a synthetic sliding impact effect as an example;
the results obtained by the calculations performed in the above steps S1 to S18 are shown in FIG. 3.
Wherein, fig. 2 is a time course curve of the acceleration of the mother wave.
Fig. 3 is a parent wave acceleration time course curve.
Fig. 4 is a graph of a near fault pulse type acceleration time course.
(1) The invention provides a near fault seismic motion fitting method based on a wavelet adjustment method, and the method can better inherit the time domain and frequency domain characteristics of parent waves and parent waves of the fitted seismic motion and is more consistent with the characteristics of the actual near fault seismic motion (steps S1-S22).
(2) The invention provides the method for generating the optimal decomposition layer number in the near fault impulse type earthquake, which can calculate more accurate wavelet coefficient and generate earthquake motion more conforming to the actual situation (step S15).
(3) The invention adjusts the seismic oscillation based on the target response spectrum and applies the continuous wavelet transform method, can realize the fitting of the target response spectrum, and can better reserve the time-frequency characteristic of the in-situ vibration (step S9).
Claims (1)
1. A pulse type seismic oscillation fitting method based on a wavelet adjustment method is characterized by comprising the following steps:
Step S1.1: determining a scale vector S from equation (1)j
In the formula:
Sj-scale vector
n0Controlling the scale vector SjParameters of the range
m-scale vector SjNumber of points
In the formula:
Zeta-wavelet basis function attenuation coefficient
Omega-wavelet basis function frequency coefficient
pi-time vector
tkWavelet basis function time parameter
Step S2: according to the impulse type seismic motion wavelet basis function obtained in the step S1Substituting the field condition information to obtain a target acceleration response spectrum Sa(Tj)]Target spectra
Step S3: calculating acceleration response spectrum (S) of original seismic oscillation, namely mother wave MWa(Tj)]Calculating the spectrum
Step S4: an adjustment coefficient gamma is obtained from the equation (3)
In the formula: gamma-adjustment factor
[Sa(Tj)]Target spectraTarget acceleration response spectrum
[Sa(Tj)]Calculating the spectrum-calculating a sumVelocity response spectrum
Step S5: obtaining a target acceleration response spectrum [ S ] according to the formula (4)a(Tj)]Target spectraAnd calculating acceleration response spectrum Sa(Tj)]Calculating the spectrumError E ofrror;
In the formula:
Errorerror of target acceleration response spectrum from calculated acceleration response spectrum
m-scale vector SjNumber of points
[Sa(Tj)]Target spectraTarget acceleration response spectrum
[Sa(Tj)]Calculating the spectrumCalculating an acceleration response spectrum
Step S6: obtaining wavelet coefficients C (s, p) according to equation (5)
In the formula:
C(sj,pi) -wavelet coefficients
Delta t-time interval of seismic oscillation acceleration time course
tkWavelet basis function time parameter
Sj-scale vector
pi-time vector
f(tk) Primary seismic oscillation
m-scale vector SjNumber of points
M-number of data points of original seismic oscillation
Step S7: obtaining detail function D (s, t) in wavelet transform according to equation (6)
In the formula:
D(sj,tk) -detail function
Δp-wavelet function time interval
Step S8: obtaining a wavelet SW according to equation (7)
In the formula:
SW-wavelet
Step S9: the average relative error E between the calculated acceleration response spectrum and the target acceleration response spectrum is obtained in step S5rrorAnd judging whether the content is less than 5%; if the average relative error is greater than 5%, returning to step S4; if the average relative error is less than or equal to 5%, outputting wavelets SW of the pulse type earthquake motion time course;
step S10: the pulse period T is obtained according to the formula (8)pThe peak value V of the pulse is obtained from the formula (9)pThe peak time t of the pulse is obtained from the equation (10)1,V
ln(Tp)=-6.45+1.11MwFormula (8)
ln(Vp)=3.680+0.065Mw+0.025ln (R) formula (9)
ln(tl,v)=1.35Mw-6.88 formula (10)
Wherein: t isp-pulse period
Vp-peak value of pulse
t1,V-pulse peak time
Distance of R-fault
MwMoment-magnitude
Step S11: analog slide-impact type near fault pulse type earthquake extraction near faultVelocity time course v of pulse type earthquakegAAs shown in formula (11);
method for simulating forward directivity effect to provide time course v of near fault pulse type seismic velocitygBAs shown in formula (12);
νgB(t)=Vpsin(ωpt)0≤t≤Tpformula (12)
Wherein:
νgAsimulation sliding impulse type near fault impulse type seismic velocity time course
νgBSimulation of forward directional effect near fault pulse type seismic velocity time course
Vp-peak value of pulse
Tp-pulse period
ωpPulse frequency of fromp=2π/TpDetermining
Step S12: velocity time course v obtained in step S11gADerivation is carried out to obtain acceleration time course alphagAAs shown in formula (13); velocity time course v obtained in step S11gAIntegral calculation to obtain displacement time interval dgA(ii) a As shown in formula (14);
in the formula: t ispDetermined by dividing the maximum of equations (11) and (14), i.e.:
velocity time course v obtained in step S11gBDerivation is carried out to obtain acceleration time course alphagBAs shown in formula (15); velocity time course v obtained in step S11gBIntegral calculation to obtain displacement time interval dgBAs shown in formula (16);
αgB(t)=ωpVpcos(ωpt),0≤t≤Tpformula (15)
In the formula:
αgAsimulation sliding impulse type near fault impulse type earthquake acceleration time course
dgASimulation sliding impulse type near fault impulse type earthquake displacement time course
αgBSimulation of forward directional effect near fault impulse type seismic acceleration time course
dgBSimulation of forward directional effect near fault pulse type seismic displacement time course
Vp-peak value of pulse
Tp-pulse period
ωpPulse frequency, from ωp=2π/TpDetermining
TpDetermined by dividing the maximum of equations (12) and (16), i.e.:
step S13: obtaining parent wave FW of pulse type earthquake motion time
Step S14: calculating the pulse period Tp
ln(Tp)=-6.45+1.11MwFormula (17)
In the formula: t isp-pulse period
Step S15: determining the number of layers n of a wavelet decomposition
In the formula: f. ofpPulse frequency
number of layers for n-wavelet transform decomposition
foriginalFrequency band range of the signal
Delta t-time interval of seismic oscillation acceleration time course
Step S16: decomposing the parent wave into n layers by multi-scale discrete wavelet transform as shown in formula (18)
In the formula: FW-father wave
FWcD1,…,FWcDn-each frequency band corresponding to the detail component of the parent wave after wavelet decomposition
FWcAnFrequency band corresponding to approximate component of parent wave after wavelet decomposition
n-number of layers of wavelet decomposition
Step S17: the parent wave wavelet coefficient LW is determined by step S16cAn
LWcAn=FWcAn
In the formula:
LWcAn-wavelet coefficient of parent wave
Step S18: decomposing the wavelet into n layers by multi-scale discrete wavelet transform as shown in formula (18)
In the formula: SW-wavelet
SWcD1,…,SWcDn-each frequency band corresponding to a detail component of a wavelet decomposition
SWcAn-bands corresponding to approximate components of wavelets after wavelet decomposition
n-number of layers of wavelet decomposition
Step S19: wavelet coefficients HW of wavelets are determined by step S18cAn
HWcAn=SWcAn
In the formula:
SWcAnwavelet coefficients of wavelets
Step S20: solving wavelet adjustment coefficient beta according to equation (20)coef
βcoef=FWcAn/SWcAnFormula (20)
In the formula: beta is acoefWavelet adjustment coefficient
Step S21: obtaining the adjusted wavelet coefficient SW 'according to the formula (21)'cAn
SW′cAn=βcoef·SWcAnFormula (21)
In the formula: SW'cAn-adjusted wavelet coefficients
Step S22: the adjusted wavelet coefficient SW 'is obtained according to the formula (21)'cAnReplacing wavelet coefficients SW after wavelet SW decompositioncAnPerforming wavelet change reconstruction to obtain the near fault pulse type seismic oscillation AW
AW-near fault pulse type seismic motion.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210229458.9A CN114578417B (en) | 2022-03-09 | 2022-03-09 | Pulse type earthquake motion fitting method based on wavelet adjustment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210229458.9A CN114578417B (en) | 2022-03-09 | 2022-03-09 | Pulse type earthquake motion fitting method based on wavelet adjustment |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114578417A true CN114578417A (en) | 2022-06-03 |
CN114578417B CN114578417B (en) | 2022-09-20 |
Family
ID=81779146
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210229458.9A Active CN114578417B (en) | 2022-03-09 | 2022-03-09 | Pulse type earthquake motion fitting method based on wavelet adjustment |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114578417B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115236732A (en) * | 2022-07-19 | 2022-10-25 | 福州大学 | Near-fault seismic oscillation fitting method based on phase spectrum reconstruction |
CN115903017A (en) * | 2022-10-17 | 2023-04-04 | 西南交通大学 | Near fault pulse type seismic oscillation parameterization method |
CN115951396A (en) * | 2022-12-30 | 2023-04-11 | 北京建筑大学 | Method and device for acquiring seismic motion recorded data, electronic equipment and storage medium |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002168964A (en) * | 2000-11-30 | 2002-06-14 | Tokyo Electric Power Co Inc:The | Earthquake motion amplitude estimation method, device thereof and storage medium storing program thereof |
CN110389379A (en) * | 2019-07-12 | 2019-10-29 | 中国地震局地球物理研究所 | The near-fault ground motion Acceleration time course approximating method of ground permanent displacement can be characterized |
CN110794459A (en) * | 2019-11-14 | 2020-02-14 | 福州大学 | Fitting method for seabed near fault seismic oscillation |
CN111487676A (en) * | 2020-03-31 | 2020-08-04 | 中国地震局工程力学研究所 | Method for synthesizing earthquake motion field by applying support vector machine, principal component analysis and particle swarm optimization in machine learning |
-
2022
- 2022-03-09 CN CN202210229458.9A patent/CN114578417B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002168964A (en) * | 2000-11-30 | 2002-06-14 | Tokyo Electric Power Co Inc:The | Earthquake motion amplitude estimation method, device thereof and storage medium storing program thereof |
CN110389379A (en) * | 2019-07-12 | 2019-10-29 | 中国地震局地球物理研究所 | The near-fault ground motion Acceleration time course approximating method of ground permanent displacement can be characterized |
CN110794459A (en) * | 2019-11-14 | 2020-02-14 | 福州大学 | Fitting method for seabed near fault seismic oscillation |
CN111487676A (en) * | 2020-03-31 | 2020-08-04 | 中国地震局工程力学研究所 | Method for synthesizing earthquake motion field by applying support vector machine, principal component analysis and particle swarm optimization in machine learning |
Non-Patent Citations (3)
Title |
---|
JIAN FAN 等: "Simulation of artificial near-fault ground motions based on the S-transform", 《2008 CONGRESS ON IMAGE AND SIGNAL PROCESSING》 * |
王珊 等: "基于小波函数的地震动反应谱与峰值位移拟合方法", 《地震学报》 * |
陈笑宇 等: "近断层地震动脉冲特性研究综述", 《工程力学》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115236732A (en) * | 2022-07-19 | 2022-10-25 | 福州大学 | Near-fault seismic oscillation fitting method based on phase spectrum reconstruction |
CN115903017A (en) * | 2022-10-17 | 2023-04-04 | 西南交通大学 | Near fault pulse type seismic oscillation parameterization method |
CN115903017B (en) * | 2022-10-17 | 2024-03-15 | 西南交通大学 | Near fault pulse type earthquake motion parameterization method |
CN115951396A (en) * | 2022-12-30 | 2023-04-11 | 北京建筑大学 | Method and device for acquiring seismic motion recorded data, electronic equipment and storage medium |
CN115951396B (en) * | 2022-12-30 | 2024-02-20 | 北京建筑大学 | Method and device for acquiring earthquake motion record data, electronic equipment and storage medium |
Also Published As
Publication number | Publication date |
---|---|
CN114578417B (en) | 2022-09-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114578417B (en) | Pulse type earthquake motion fitting method based on wavelet adjustment | |
Pestiaux et al. | Paleoclimatic variability at frequencies ranging from 1 cycle per 10 000 years to 1 cycle per 1000 years: Evidence for nonlinear behaviour of the climate system | |
Huang et al. | A confidence limit for the empirical mode decomposition and Hilbert spectral analysis | |
CN113806845B (en) | POD interpolation-based non-stationary wind field efficient simulation method | |
CN114442153B (en) | Near fault earthquake motion fitting method | |
Broccardo et al. | Simulation of near-fault ground motions using frequency-domain discretization | |
CN109002678B (en) | Artificial earthquake simulation method based on hydroelectric engineering hydraulic structure earthquake-resistant design | |
Li et al. | Simulation and generation of spectrum-compatible ground motions based on wavelet packet method | |
Shahjouei et al. | Synthetic seismograms using a hybrid broadband ground‐motion simulation approach: Application to central and eastern United States | |
Wang et al. | Tidal modulation of surface gravity waves in the Gulf of Maine | |
CN108169794A (en) | Nearly tomography Artificial Seismic Wave and the matched method of adjustment of Code response spectrum | |
Shama | Spectrum compatible earthquake ground motions by morlet wavelet | |
Temperton | Some experiments in dynamic initialization for a simple primitive equation model | |
Chalikov | Numerical Modeling of Sea Waves | |
Doron et al. | A comparison of grid‐point and spectral methods in a meteorological problem | |
Samelson et al. | Currents forced by stochastic winds with meridionally varying amplitude | |
CN117252013B (en) | Method for constructing broadband limited seismic source model | |
CN115236732A (en) | Near-fault seismic oscillation fitting method based on phase spectrum reconstruction | |
CN114966835A (en) | Near-fault seismic motion fitting method based on improved time domain superposition method | |
Mashayekhi et al. | Non-Stationary Spectral Matching by Unconstrained Optimization and Discrete Wavelet Transformation 2 | |
Bender et al. | Wave modeling with SWAN+ ADCIRC for the South Carolina coastal storm surge study | |
CN116558753A (en) | Long-period earthquake motion simulation generation method | |
Litvenko et al. | Numerical stochastic models of sea surface undulation and extreme ocean waves | |
Nguyen et al. | Establishment of Artificial Accelerogram for Structural Analysis according to the Earthquake Characteristics of Vietnam | |
Guo et al. | Estimation of wave velocity and analysis of dispersion characteristics based on wavelet decomposition |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |