CN114199489A - Method and device for adjusting natural seismic wave fitting response spectrum based on digital filtering - Google Patents
Method and device for adjusting natural seismic wave fitting response spectrum based on digital filtering Download PDFInfo
- Publication number
- CN114199489A CN114199489A CN202111508287.5A CN202111508287A CN114199489A CN 114199489 A CN114199489 A CN 114199489A CN 202111508287 A CN202111508287 A CN 202111508287A CN 114199489 A CN114199489 A CN 114199489A
- Authority
- CN
- China
- Prior art keywords
- spectrum
- frequency
- natural
- response spectrum
- time course
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M7/00—Vibration-testing of structures; Shock-testing of structures
- G01M7/02—Vibration-testing by means of a shake table
- G01M7/022—Vibration control arrangements, e.g. for generating random vibrations
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a method and a device for adjusting a natural seismic wave fitting reaction spectrum based on digital filtering, wherein the method comprises the following steps: determining the number of frequency bands and the frequency range contained in each frequency band by taking the frequency point of the target response spectrum as a node, decomposing the natural seismic motion record into wavelets with different frequency bands by using a digital filtering technology, and reconstructing each wavelet to obtain an initial time course required by iteration; and for the initial time course, iteratively and gradually adjusting the amplitude coefficient of each wavelet by using an influence matrix method until the fitting precision of the corresponding time course response spectrum and the target response spectrum meets the requirement. The method overcomes the problem of insufficient preservation of seismic oscillation natural characteristics in the current reaction spectrum fitting technology, has high iterative calculation efficiency and high matching precision, and is beneficial to improving the power time-course response analysis efficiency of the engineering structure and obtaining an analysis result with higher reliability.
Description
Technical Field
The invention relates to a structural seismic design and analysis method, in particular to a method for adjusting a natural seismic wave fitting reaction spectrum.
Background
In anti-seismic specifications and standards, power time-course response analysis is often required to be carried out on important engineering structures, and the power time-course response analysis needs proper design seismic motion input. There are three methods of adjusting and generating design seismic motions that are currently in common use: the first is to amplitude-modulate the historical earthquake motion as the input earthquake motion, the second is to use the pure artificial synthetic earthquake motion as the earthquake motion input, and the third is to use the numerical method to adjust the natural earthquake motion to match the design response spectrum and then use it as the earthquake motion input. Because the number of strong earthquakes and the regional distribution of historical records are limited, fewer original earthquake records can be used for constructing a site in actual engineering; the pure artificial synthesis earthquake motion completely utilizes the artificial wave spectrum to represent earthquake waves, so that the action effect of the earthquake motion is difficult to simulate really; therefore, the method for adjusting natural earthquake motion by using a numerical method to fully match or envelop the design response spectrum and meet the seismic motion engineering characteristics required by the specification is a more reasonable method, and the method is also the key focus content of theoretical research in the field of current seismic engineering and the problem to be solved urgently. However, the problems of long calculation time consumption, low fitting precision and the like often exist in the prior art, and the requirements of engineering structure seismic design and analysis cannot be met.
Disclosure of Invention
The purpose of the invention is as follows: aiming at the defects of the prior art, the invention aims to provide a method for adjusting the fitting reaction spectrum of the natural seismic waves based on digital filtering, which can fully retain the characteristics of the natural seismic motion record and has higher calculation efficiency and precision.
Another object of the present invention is to provide a corresponding apparatus for adjusting the fitted response spectrum of natural seismic waves based on digital filtering.
The technical scheme is as follows: in a first aspect, a method for adjusting a fitted response spectrum of natural seismic waves based on digital filtering includes the following steps:
(1) selecting a natural earthquake motion time course A (T) based on the attribute of a target reaction spectrum, wherein the duration length of the time course is T, namely T belongs to [0, T ], and the time interval is dt;
(2) based on target response spectra
Determining M frequency calculation points as f1；f2；...；fMDetermining the number of frequency bands as M-1 and the frequency intervals as (- ∞, f)2],[f2,f3]...,[fM-2,fM-1],[fM-1,+∞)；
(3) Performing band-pass filtering on the natural earthquake motion time course by using a digital filtering technology and taking the upper limit and the lower limit of a 2 nd to an M-2 nd frequency interval as cut-off frequencies, wherein the 1 st frequency band is f2Low pass filtering for cut-off frequency, the M-1 th band being defined by fM-1High-pass filtering is carried out for the cut-off frequency, and M-1 basis functions g are obtained by decomposing the natural earthquake motion time coursen(t), n is 1, 2.., M-1, wherein the frequency range contained in the basis function covers the frequency range possessed by the original natural earthquake;
(4) determining amplitude coefficients constituting a basis function of the time course of the natural earthquake
Reconstructing to obtain an initial time course
And an initial time course A(0)(t) has a reaction spectrum of S(0)；
(5) And carrying out iterative operation on the initial time interval by using an influence matrix method until the fitting accuracy of the time interval response spectrum and the target response spectrum meets the requirement, thereby obtaining the seismic motion time interval fitted with the target response spectrum.
Further, in the step (1), a natural earthquake motion time interval corresponding to the target response spectrum is selected according to the field type, the seismic source characteristics and the response spectrum shape.
Further, a finite impulse response FIR digital filter based on a Kaiser window is used for digital filtering in the step (3).
Further, the step (4) of performing an iterative operation on the initial time interval by using an influence matrix method includes: based on an initial time course A(0)(t) calculating the nth basis function gn(t) for the mth frequency point fmThe contribution of the reaction spectrum value is marked as an influence factor ImnAccording to ImnConstructing an influence matrix
For the polarity of the time course response spectrum at the mth frequency point, the amplitude coefficient of each basis function is adjusted step by step through iteration until the time course response spectrum S(0)Spectrum S of reaction with targetTUntil the matching accuracy of the data processing system meets the requirement.
Further, the influence factor ImnCalculated as follows:
wherein, ω ismCircle frequency omega corresponding to mth frequency pointm＝2πfm；hm(t) is a unit pulse function; tau ismThe moment of the acceleration peak response of the mth frequency single degree of freedom.
In a second aspect, an apparatus for adjusting a fitted response spectrum of natural seismic waves based on digital filtering is provided, which includes:
the natural earthquake motion time course selection module is used for selecting a natural earthquake motion time course A (T) based on the attribute of the target reaction spectrum, the duration length of the time course is T, namely T belongs to [0, T ], and the time interval is dt;
a target response spectrum band determination module for determining a target response spectrum based on the target response spectrum
Determining M frequency calculation points as f1；f2；...；fMDetermining the number of frequency bands as M-1 and the frequency intervals as (- ∞, f)2],[f2,f3]...,[fM-2,fM-1],[fM-1,+∞)；
A digital filtering decomposition module for performing band-pass filtering on the natural earthquake motion time course by using the digital filtering technology and taking the upper and lower limits of the 2 nd to M-2 nd frequency intervals as cut-off frequencies, wherein the 1 st frequency band is f2Low pass filtering for cut-off frequency, the M-1 th band being defined by fM-1High-pass filtering is carried out for the cut-off frequency, and M-1 basis functions g are obtained by decomposing the natural earthquake motion time coursen(t), n is 1, 2.., M-1, wherein the frequency range contained in the basis function covers the frequency range possessed by the original natural earthquake;
an initial time-course reconstruction module for determining the amplitude coefficient of the basis function of the natural seismic motion time-course
Reconstructing to obtain an initial time course
And an initial time course A(0)(t) has a reaction spectrum of S(0)；
And the iterative fitting module is used for carrying out iterative operation by utilizing an influence matrix method according to the initial time interval until the fitting precision of the time interval response spectrum and the target response spectrum meets the requirement, so that the seismic motion time interval fitted with the target response spectrum is obtained.
In a third aspect, a computer device is provided, comprising:
one or more processors;
a memory; and
one or more programs stored in the memory and configured to be executed by the one or more processors, the programs when executed by the processors implementing the steps of the method for adjusting a natural seismic wave fit response spectrum based on digital filtering according to the first aspect of the invention.
Has the advantages that: compared with the prior art, the invention has technical changes. The invention introduces a digital filtering technology to carry out filtering decomposition on the natural recording earthquake motion time course to obtain a basis function, and reconstructs the original earthquake motion time course; and then the earthquake time course response spectrum and the target response spectrum are gradually fitted according to the iterative calculation of an influence matrix method so as to meet the precision requirement required by the standard or specific engineering. Compared with the prior art, the method has high calculation efficiency, and the matching precision is higher, and simultaneously, the natural earthquake motion characteristic of the earthquake time course is kept as much as possible.
Drawings
FIG. 1 is a general flow chart of a method for adjusting a fitted response spectrum of natural seismic waves based on digital filtering according to an embodiment of the invention;
FIG. 2 is an example of an initial iterative time course response spectrum and a target response spectrum of an embodiment of the present invention;
fig. 3 is an example of intermediate results and final results obtained by iteration in the embodiment of the present invention.
Detailed Description
The technical scheme of the invention is further explained by combining the attached drawings.
In one embodiment, the R.G.1.60 design spectrum is selected as a target response spectrum, and the seismic motion input time course matched with the target response spectrum is generated by using the method. As shown in fig. 1, a method for generating designed seismic motion based on filtering technology includes the following steps:
(1) and selecting a natural earthquake motion time course which has the same field type and earthquake source characteristics as the target field and has the frequency component distribution of the response spectrum similar to the target response spectrum.
The response spectrums have similar shapes, for example, the target spectrum has a higher middle part and lower ends, and the seismic motion time-course response spectrum also has a higher middle part and lower ends. In this embodiment, the acceleration time interval of the Mammoth Lake earthquake recorded by the Long Valley station is selected as the original natural earthquake motion time interval, and the comparison between the initial time interval and the target response spectrum is shown in fig. 2. The duration of the time interval is T-30 s and the time interval dt-0.005 s. Maximum and minimum frequencies involved in r.g.1.60 design spectraAt a frequency of fmin＝0.1Hz，fmax100Hz, i.e. the target response spectrum has a calculated frequency range of [0.1,100]Hz。
(2) And determining the number of frequency bands and the frequency range contained in each frequency band according to the target response spectrum frequency point.
The number of frequency points of the target reaction spectrum R.G.1.60 design spectrum is M301, and the frequency range is [0.1,100]]Frequency points in Hz are uniformly distributed in logarithmic coordinates, i.e. f1＝0.1000,f2＝0.1023,...,f300＝97.7237,f301100.0000. Therefore, 300 frequency band intervals are determined to be (- ∞, f)2],...,[f299,f300],[f300,+∞)。
(3) And decomposing the natural earthquake into wavelets containing various frequency bands by using a digital filtering technology, and reconstructing to obtain an initial time course required by iteration.
In this embodiment, a finite impulse response filter with a Kaiser window function of 900 orders is used to perform filter decomposition on the historical seismic acceleration in 300 divided frequency bands, wherein the 2 nd to 299 th frequency bands are respectively subjected to band-pass filtering with the upper and lower limits of the interval as cut-off frequencies to respectively obtain basis functions g2(t),g2(t),...,g299(t) 1 st frequency band with f2Low pass filtering the cut-off frequency to obtain the basis function g1(t) 300 th band with f300High-pass filtering the cut-off frequency to obtain a basis function g300(t) finally obtaining filtering time intervals g respectively covering 300 frequency band intervalsn(t),n＝1,2,...,300。
The amplitude coefficient of the original natural seismic oscillation reconstructed by 300 basis functions obtained by filtering is obtained by solving by utilizing a generalized minimum residual error method
The reconstructed time course is recorded as
The reaction spectrum of the reconstructed time course is S(0)。
(5) And performing iterative fitting on the target reaction spectrum by using an influence matrix method.
Assuming that the relative error between the time-course reaction spectrum and the target spectrum at each frequency point is not more than eta which is 15%, adjusting the amplitude coefficients of all basis functions within the frequency range of 0.1,100 Hz in each iteration by using an influence matrix method, gradually fitting the time-course reaction spectrum obtained in each iteration to the target reaction spectrum, and stopping the iteration until the maximum relative error between the time-course reaction spectrum and the target spectrum is not more than a threshold eta.
(6) And outputting a time interval when the relative error meets the threshold requirement.
The intermediate result and the final result of the iterative process are shown in fig. 3, the obtained acceleration response spectrum of the time course and the target response spectrum meet the required matching precision, and the method enables the time course response spectrum to be uniformly and consistently close to the target response spectrum.
Based on the same technical concept of the above method embodiments, in another embodiment, there is provided an apparatus for adjusting a natural seismic wave fitting response spectrum based on digital filtering, including:
the natural earthquake motion time course selection module is used for selecting a natural earthquake motion time course A (T) based on the attribute of the target reaction spectrum, the duration length of the time course is T, namely T belongs to [0, T ], and the time interval is dt;
a target response spectrum band determination module for determining a target response spectrum based on the target response spectrum
Determining M frequency calculation points as f1；f2；...；fMDetermining the number of frequency bands as M-1 and the frequency intervals as (- ∞, f)2],[f2,f3],...,[fM-2,fM-1],[fM-1,+∞)；
A digital filtering decomposition module for performing band-pass filtering on the natural earthquake motion time course by using the digital filtering technology and taking the upper and lower limits of the 2 nd to M-2 nd frequency intervals as cut-off frequencies, wherein the 1 st frequency band is f2Low pass filtering for cut-off frequency, the M-1 th band being defined by fM-1High-pass filtering is carried out for the cut-off frequency, and M-1 basis functions g are obtained by decomposing the natural earthquake motion time coursen(t), n ═ 1, 2.., M-1, where the basis functions compriseThe frequency range covers the frequency range of the original natural earthquake motion;
an initial time-course reconstruction module for determining the amplitude coefficient of the basis function of the natural seismic motion time-course
Reconstructing to obtain an initial time course
And an initial time course A(0)(t) has a reaction spectrum of S(0)；
And the iterative fitting module is used for carrying out iterative operation by utilizing an influence matrix method according to the initial time interval until the fitting precision of the time interval response spectrum and the target response spectrum meets the requirement, so that the seismic motion time interval fitted with the target response spectrum is obtained.
The natural earthquake motion time interval selection module selects a natural earthquake motion time interval corresponding to the target response spectrum according to the field type, the earthquake source characteristics and the response spectrum shape.
In this embodiment, the digital filter decomposition module uses a finite impulse response FIR digital filter based on a Kaiser window.
As a preferred embodiment, the iterative fitting module includes:
an influence matrix construction unit for constructing an influence matrix based on the initial time interval A(0)(t) calculating the nth basis function gn(t) for the mth frequency point fmThe contribution of the reaction spectrum value is marked as an influence factor ImnAccording to ImnConstructing an influence matrix
The polarity of the time course response spectrum at the mth frequency point;
an iterative calculation unit for gradually adjusting the amplitude coefficient of each basis function through iteration until the time course response spectrum S(0)Spectrum S of reaction with targetTUntil the matching precision of the target meets the requirement; and
and the time interval output unit is used for calculating to obtain a corresponding speed time interval and a corresponding displacement time interval through the amplitude coefficient of the acceleration time interval meeting the matching precision, and outputting acceleration, speed and displacement time interval data meeting the requirements.
Wherein the influence factor ImnCalculated as follows:
wherein, ω ismCircle frequency omega corresponding to mth frequency pointm＝2πfm；hm(t) is a unit pulse function; tau ismThe moment of the acceleration peak response of the mth frequency single degree of freedom.
It should be understood that the apparatus for adjusting a natural seismic wave fitting response spectrum based on digital filtering provided in this embodiment may implement all technical solutions in the above method embodiments, functions of each functional module may be specifically implemented according to the method in the above method embodiments, and specific implementation processes not described in detail in this embodiment may refer to relevant descriptions in the above embodiments.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, apparatus, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus, devices (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solutions of the present invention and not for limiting the same, and although the present invention is described in detail with reference to the above embodiments, those of ordinary skill in the art should understand that: modifications and equivalents may be made to the embodiments of the invention without departing from the spirit and scope of the invention, which is to be covered by the claims.
Claims (10)
1. A method for adjusting a natural seismic wave fitting reaction spectrum based on digital filtering is characterized by comprising the following steps:
(1) selecting a natural earthquake motion time course A (T) based on the attribute of a target reaction spectrum, wherein the duration length of the time course is T, namely T belongs to [0, T ], and the time interval is dt;
(2) based on target response spectra
Determining M frequency calculation points as f1；f2；...；fMDetermining the number of frequency bands as M-1 and the frequency intervals as (- ∞, f)2],[f2,f3]...,[fM-2,fM-1],[fM-1,+∞)；
(3) Performing band-pass filtering on the natural earthquake motion time course by using a digital filtering technology and taking the upper limit and the lower limit of a 2 nd to an M-2 nd frequency interval as cut-off frequencies, wherein the 1 st frequency band is f2Low pass filtering for cut-off frequency, the M-1 th band being defined by fM-1High-pass filtering is carried out for the cut-off frequency, and M-1 basis functions g are obtained by decomposing the natural earthquake motion time coursen(t), n is 1, 2.., M-1, wherein the frequency range contained in the basis function covers the frequency range possessed by the original natural earthquake;
(4) determining amplitude coefficients constituting a basis function of the time course of the natural earthquake
Reconstructing to obtain an initial time course
And an initial time course A(0)(t) has a reaction spectrum of S(0)；
(5) And carrying out iterative operation on the initial time interval by using an influence matrix method until the fitting accuracy of the time interval response spectrum and the target response spectrum meets the requirement, thereby obtaining the seismic motion time interval fitted with the target response spectrum.
2. The method for adjusting the fitted response spectrum of natural seismic waves based on digital filtering according to claim 1, wherein: and (2) selecting a natural earthquake motion time interval corresponding to the target reaction spectrum according to the field type, the earthquake source characteristics and the reaction spectrum shape in the step (1).
3. The method for adjusting the fitted response spectrum of natural seismic waves based on digital filtering according to claim 1, wherein: and (3) when the digital filtering is carried out in the step (3), a finite impulse response FIR digital filter based on a Kaiser window is used.
4. The method for adjusting the fitted response spectrum of natural seismic waves based on digital filtering according to claim 1, wherein: the step (4) of performing iterative operation on the initial time interval by using an influence matrix method comprises the following steps: based on an initial time course A(0)(t) calculating the nth basis function gn(t) for the mth frequency point fmThe contribution of the reaction spectrum value is marked as an influence factor ImnAccording to ImnConstructing an influence matrix
For the polarity of the time course response spectrum at the mth frequency point, the amplitude coefficient of each basis function is adjusted step by step through iteration until the time course response spectrum S(0)Spectrum S of reaction with targetTUntil the matching accuracy of the data processing system meets the requirement.
5. The method for adjusting the fitted response spectrum of natural seismic waves based on digital filtering according to claim 4, wherein: the influencing factor ImnCalculated as follows:
wherein, ω ismCircle frequency omega corresponding to mth frequency pointm＝2πfm；hm(t) is a unit pulse function; tau ismThe moment of the acceleration peak response of the mth frequency single degree of freedom.
6. An apparatus for adjusting a fitted response spectrum of natural seismic waves based on digital filtering, comprising:
the natural earthquake motion time course selection module is used for selecting a natural earthquake motion time course A (T) based on the attribute of the target reaction spectrum, the duration length of the time course is T, namely T belongs to [0, T ], and the time interval is dt;
a target response spectrum band determination module for determining a target response spectrum based on the target response spectrum
Determining M frequency calculation points as f1；f2；...；fMDetermining the number of frequency bands as M-1 and the frequency intervals as (- ∞, f)2],[f2,f3]...,[fM-2,fM-1],[fM-1,+∞)；
A digital filtering decomposition module for performing band-pass filtering on the natural earthquake motion time course by using the digital filtering technology and taking the upper and lower limits of the 2 nd to M-2 nd frequency intervals as cut-off frequencies, wherein the 1 st frequency band is f2Low pass filtering for cut-off frequency, the M-1 th band being defined by fM-1High-pass filtering is carried out for the cut-off frequency, and M-1 basis functions g are obtained by decomposing the natural earthquake motion time coursen(t), n is 1,2, M-1, wherein the frequency range contained in the basis function covers the frequency range possessed by the original natural earthquake motion;
an initial time-course reconstruction module for determining the amplitude coefficient of the basis function of the natural seismic motion time-course
Reconstructing to obtain an initial time course
And an initial time course A(0)(t) has a reaction spectrum of S(0)；
And the iterative fitting module is used for carrying out iterative operation by utilizing an influence matrix method according to the initial time interval until the fitting precision of the time interval response spectrum and the target response spectrum meets the requirement, so that the seismic motion time interval fitted with the target response spectrum is obtained.
7. The apparatus for adjusting the fitted response spectrum of natural seismic waves based on digital filtering according to claim 6, wherein: and the natural earthquake motion time interval selection module selects a natural earthquake motion time interval corresponding to the target response spectrum according to the field type, the earthquake source characteristics and the response spectrum shape.
8. The apparatus for adjusting the fitted response spectrum of natural seismic waves based on digital filtering according to claim 6, wherein: the digital filter decomposition module uses a finite impulse response FIR digital filter based on a Kaiser window.
9. The apparatus for adjusting the fitted response spectrum of natural seismic waves based on digital filtering according to claim 6, wherein: the iterative fitting module comprises:
an influence matrix construction unit for constructing an influence matrix based on the initial time interval A(0)(t) calculating the nth basis function gn(t) for the mth frequency point fmThe contribution of the reaction spectrum value is marked as an influence factor ImnAccording to ImnConstructing an influence matrix
The polarity of the time course response spectrum at the mth frequency point;
an iterative calculation unit for gradually adjusting the amplitude coefficient of each basis function through iteration until the time course response spectrum S(0)Spectrum S of reaction with targetTUntil the matching precision of the target meets the requirement;
and the time interval output unit is used for calculating to obtain a corresponding speed time interval and a corresponding displacement time interval through the amplitude coefficient of the acceleration time interval meeting the matching precision, and outputting acceleration, speed and displacement time interval data meeting the requirements.
10. A computer device, comprising:
one or more processors;
a memory; and
one or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the one or more processors, which when executed by the processors, implement the steps of the method of adjusting natural seismic wave fitting response spectra based on digital filtering according to any of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111508287.5A CN114199489B (en) | 2021-12-10 | 2021-12-10 | Method and device for adjusting natural seismic wave fitting reaction spectrum based on digital filtering |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111508287.5A CN114199489B (en) | 2021-12-10 | 2021-12-10 | Method and device for adjusting natural seismic wave fitting reaction spectrum based on digital filtering |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114199489A true CN114199489A (en) | 2022-03-18 |
| CN114199489B CN114199489B (en) | 2023-05-09 |
Family
ID=80652168
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111508287.5A Active CN114199489B (en) | 2021-12-10 | 2021-12-10 | Method and device for adjusting natural seismic wave fitting reaction spectrum based on digital filtering |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114199489B (en) |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH11118939A (en) * | 1997-10-17 | 1999-04-30 | Toshiba Corp | Earthquake observation system |
| EP1089507A2 (en) * | 1999-09-28 | 2001-04-04 | Siemens Information and Communication Networks S.p.A. | Method for the calculation of the coefficients of a polyphase FIR filter |
| CN102193107A (en) * | 2010-03-05 | 2011-09-21 | 西安石油大学 | Method for separating and denoising seismic wave field |
| US20120163122A1 (en) * | 2010-12-23 | 2012-06-28 | Ahmet Kemal Ozdemir | Removing Noise from a Seismic Measurement |
| US20140278382A1 (en) * | 2005-02-23 | 2014-09-18 | Digital Intelligence, L.L.C. | Signal decomposition, analysis and reconstruction using high-resolution filter banks and component tracking |
| US20140288841A1 (en) * | 2013-03-19 | 2014-09-25 | Westerngeco L.L.C. | Removing noise from a seismic measurement |
| CN106093701A (en) * | 2016-06-06 | 2016-11-09 | 国家电网公司 | A kind of cable fault signal detecting method based on empirical mode decomposition filtering |
| CN106096579A (en) * | 2016-06-22 | 2016-11-09 | 天津理工大学 | A kind of method of ECG signal processing |
| CN107144880A (en) * | 2017-05-12 | 2017-09-08 | 招商局重庆交通科研设计院有限公司 | A kind of seismic wave wave field separation method |
| US20180024263A1 (en) * | 2016-07-25 | 2018-01-25 | Schlumberger Technology Corporation | Seismic Spectral Balancing |
| CN110020400A (en) * | 2019-01-07 | 2019-07-16 | 河海大学 | A kind of adjustment seismic wave and goal response compose accurate matched influence matrix method |
| CN110069836A (en) * | 2019-04-03 | 2019-07-30 | 河海大学 | A kind of height frequency range is alternately and target composes matched improvement influence matrix method |
-
2021
- 2021-12-10 CN CN202111508287.5A patent/CN114199489B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH11118939A (en) * | 1997-10-17 | 1999-04-30 | Toshiba Corp | Earthquake observation system |
| EP1089507A2 (en) * | 1999-09-28 | 2001-04-04 | Siemens Information and Communication Networks S.p.A. | Method for the calculation of the coefficients of a polyphase FIR filter |
| US20140278382A1 (en) * | 2005-02-23 | 2014-09-18 | Digital Intelligence, L.L.C. | Signal decomposition, analysis and reconstruction using high-resolution filter banks and component tracking |
| CN102193107A (en) * | 2010-03-05 | 2011-09-21 | 西安石油大学 | Method for separating and denoising seismic wave field |
| US20120163122A1 (en) * | 2010-12-23 | 2012-06-28 | Ahmet Kemal Ozdemir | Removing Noise from a Seismic Measurement |
| US20140288841A1 (en) * | 2013-03-19 | 2014-09-25 | Westerngeco L.L.C. | Removing noise from a seismic measurement |
| CN106093701A (en) * | 2016-06-06 | 2016-11-09 | 国家电网公司 | A kind of cable fault signal detecting method based on empirical mode decomposition filtering |
| CN106096579A (en) * | 2016-06-22 | 2016-11-09 | 天津理工大学 | A kind of method of ECG signal processing |
| US20180024263A1 (en) * | 2016-07-25 | 2018-01-25 | Schlumberger Technology Corporation | Seismic Spectral Balancing |
| CN107144880A (en) * | 2017-05-12 | 2017-09-08 | 招商局重庆交通科研设计院有限公司 | A kind of seismic wave wave field separation method |
| CN110020400A (en) * | 2019-01-07 | 2019-07-16 | 河海大学 | A kind of adjustment seismic wave and goal response compose accurate matched influence matrix method |
| CN110069836A (en) * | 2019-04-03 | 2019-07-30 | 河海大学 | A kind of height frequency range is alternately and target composes matched improvement influence matrix method |
Non-Patent Citations (11)
| Title |
|---|
| 冯雪玲等: "2008年新疆乌恰M_S6.8地震和M_S6.3余震的重叠地震波形分析", 《内陆地震》 * |
| 吴琛: "小波分解法在多自由度线性体系地震反应分析中的应用", 《福建工程学院学报》 * |
| 徐龙军;赵国臣;谢礼立;: "基于分量分离方法的地震动反应谱" * |
| 文建波,周进雄,张陵: "设计地震波的小波包多尺度调整" * |
| 李安: "体表电位标测中的关键技术――心电信号的数字滤波", 《南昌航空工业学院学报》 * |
| 樊剑;吕超;张辉;: "地震波时频特征及与结构地震响应的关系" * |
| 武安绪;林向东;穆会泳;赵桂儒;柴金翼;吴培稚;: "EMD新技术在数字波形预处理中的初步应用" * |
| 盛涛等: "拟合D-V-A联合谱及各项峰值的人造地震动", 《振动.测试与诊断》 * |
| 许立英;刘宁;吴应雄;: "基于抗震分析的长周期地震动的界定与选取" * |
| 赵志航;车伟;冯群;: "基于移动平均分解的近断层地震动弹性反应谱研究" * |
| 陈可洋;范兴才;吴清岭;陈树民;李来林;刘振宽;王建民;关昕;: "基于扩散滤波的多尺度分解和重构方法及应用初探" * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114199489B (en) | 2023-05-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107272066B (en) | A kind of noisy seismic signal first-arrival traveltime pick-up method and device | |
| CN110598166B (en) | Wavelet denoising method for adaptively determining wavelet layering progression | |
| CN106997767A (en) | Method of speech processing and device based on artificial intelligence | |
| CN104807534B (en) | Equipment eigentone self study recognition methods based on on-line vibration data | |
| KR20140079369A (en) | System and method of processing a sound signal including transforming the sound signal into a frequency-chirp domain | |
| CN108334987B (en) | Wavelet decomposition-neural network-based sea wave height prediction method | |
| CN107480391B (en) | Near-fault non-stationary seismic oscillation simulation method based on data driving | |
| KR20050074451A (en) | Method and system for processing a sound field representation | |
| CN114578417B (en) | Pulse type earthquake motion fitting method based on wavelet adjustment | |
| CN106443787A (en) | Prestack seismic gather noise suppression method and device | |
| CN109658935A (en) | The generation method and system of multichannel noisy speech | |
| CA2462325A1 (en) | Reverberation sound generating apparatus | |
| Lee | Improvement of decomposing results of empirical mode decomposition and its variations for sea-level records analysis | |
| CN114199489A (en) | Method and device for adjusting natural seismic wave fitting response spectrum based on digital filtering | |
| WO2020057286A1 (en) | Wave field forward modeling method and device | |
| CN106950597A (en) | The mixing source data separation method filtered based on three sides | |
| Webster | A random number generator for ocean noise statistics | |
| CN113866827B (en) | Interpretation velocity modeling seismic imaging method, system, medium and equipment | |
| CN110376642B (en) | Three-dimensional seismic velocity inversion method based on conical surface waves | |
| JP2000296118A (en) | Method and device for analyzing living body signal | |
| Miguel et al. | Reconstruction of the heart rate signal using the dominant frequencies | |
| JP3984131B2 (en) | Waveform restoration apparatus and waveform restoration method | |
| CN114063167B (en) | Seismic data processing method for suppressing exogenous interference | |
| CN113945974B (en) | Seismic data reconstruction method, device, equipment and medium | |
| CN117497194B (en) | Biological information processing method and device, electronic equipment and storage medium |
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 |