CN110609323B - Simplified calculation method for site excellent period based on stratum information - Google Patents
Simplified calculation method for site excellent period based on stratum information Download PDFInfo
- Publication number
- CN110609323B CN110609323B CN201910903742.8A CN201910903742A CN110609323B CN 110609323 B CN110609323 B CN 110609323B CN 201910903742 A CN201910903742 A CN 201910903742A CN 110609323 B CN110609323 B CN 110609323B
- Authority
- CN
- China
- Prior art keywords
- seismic wave
- time
- stratum
- site
- excellent period
- 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.)
- Active
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 58
- 238000001228 spectrum Methods 0.000 claims abstract description 31
- 238000005070 sampling Methods 0.000 claims description 27
- 230000015572 biosynthetic process Effects 0.000 claims description 18
- 238000011835 investigation Methods 0.000 claims description 5
- 230000000737 periodic effect Effects 0.000 claims description 4
- 101100328518 Caenorhabditis elegans cnt-1 gene Proteins 0.000 claims description 3
- 238000000034 method Methods 0.000 abstract description 18
- 239000002689 soil Substances 0.000 abstract description 9
- 230000001276 controlling effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000000691 measurement 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
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 discloses a simplified calculation method of a site excellent period based on stratum information, which is characterized in that an estimated value of the excellent period is obtained by applying an excellent period simple algorithm, and all calculation parameters except the site stratum information are set according to the value; taking the field soil as a seismic wave signal filter, calculating a ground seismic time course under the action of pulse waves, and performing Fourier transform on the ground seismic time course to obtain an amplitude spectrum; by judging the relation between the maximum value of the amplitude spectrum and the time (or frequency), a high-precision excellent period calculation value is obtained. The simplified calculation method for the site excellent period based on the stratum information is simple in calculation process and high in calculation precision.
Description
Technical Field
The invention belongs to the technical field of geotechnical engineering investigation, and particularly relates to a simplified calculation method for a site excellent period based on stratum information.
Background
The excellent period refers to a harmonic component of seismic waves which have resonance action with the foundation soil layer. At present, methods for determining the excellent period of a field are roughly classified into two major types, namely a direct measurement method and a wave velocity method. In any method, calculation parameters except formation information need to be input, and even the calculation parameters need to be debugged, so that a proper calculation result can be obtained; and the calculation precision is not high, or the calculation process is complicated.
Disclosure of Invention
The invention aims to provide a simplified calculation method of a site excellent period based on stratum information, which is simple in calculation process and high in calculation precision.
The technical scheme adopted by the invention is that a simplified calculation method of the site excellent period based on stratum information is implemented according to the following steps:
step 9, judging whether the amplitude spectrum | Fu (k) | has a maximum value, if not, expanding nT by 20%, and returning to the step 8; if yes, determining the site excellent period T0。
The invention is also characterized in that:
in step 2, site excellent period estimated value t0The formula is as follows:
t0=4(h1/v1+h2/v2+……+hn/vn) (1)。
in step 3, the seismic wave duration tlThe formula is as follows:
tl≥10 t0 (2)
the time sampling interval Δ t is calculated by the following equation (3):
Δt<0.02min(h1/v1,h2/v2,……,hn/vn) (3)。
in step 4, sampling point number nt of ground seismic wave holding time and sampling point number nt of seismic wave propagation time in stratum iiThe sampling point number nT of the Fourier cycle spectrum of the ground seismic wave is calculated by the following formula:
nt=Int(tl/Δt) (4)
nti=Int(hi/vi/Δt) (5)
nT=Int(t0/Δt) (6)
wherein Int is an integer function, i is 1, 2.
In step 5, the reflection coefficient RiThe formula is as follows:
Ri=(ρi+1vi+1-ρivi)/(ρi+1vi+1+ρivi) Wherein, i is 1, 2, a. (7)
In step 6, seismic wave time course u0 dThe formula is as follows:
u0 d(j)=-u1 p(j-nt1)
u1 d(j)=(1+R1)u0 d(j-nt1)-R1u2 p(j-nt2)
u1 p(j)=R1u0 d(j-nt1)+(1-R1)u2 p(j-nt2)
u2 d(j)=(1+R2)u1 d(j-nt2)-R2u3 p(j-nt3)
u2 p(j)=R2u1 d(j-nt2)+(1-R2)u3 p(j-nt3)
......
ui d(j)=(1+Ri)ui-1 d(j-nti)-Riui+1 p(j-nti+1)
ui p(j)=Riui-1 d(j-nti)+(1-Ri)ui+1 p(j-nti+1)
......
un-1 d(j)=(1+Rn-1)un-2 d(j-ntn-1)-Rn-1un p(j-ntn)
un-1 p(j)=Rn-1un-2 d(j-ntn-1)+(1-Rn-1)un p(j-ntn)
un p(j)=un d(j-ntn)
wherein u is0 dIs the ground seismic wave time course; u. ofi d、ui pThe seismic wave time courses of the lower interface and the upper interface of the formation i are respectively, wherein i is 1, 2. u. of0 d、ui d、ui pAre time series, j denotes a time sample number, j 0, 1.
In step 7, the seismic wave time travel u is calculated according to the following formula0 dCarrying out attenuation treatment according to time:
u(j)=u0 d(j)exp[-βjΔt] (9)
wherein the attenuation coefficient beta is-ln delta/tl(ii) a Delta is a calculation precision value for controlling the attenuation coefficient and satisfies 0<δ<0.01;j=0,1,......,nt-1。
In step 8, performing Discrete Fourier Transform (DFT) on the seismic time interval u (j) according to the following formula:
wherein i is an imaginary unit; k is 0, 1,.. cnt-1.
In step 9, whether the amplitude spectrum | fu (k) | has a maximum value is determined according to the following formula:
| Fu (k) | > | Fu (k-1) | and | Fu (k) | > | Fu (k +1) | (11)
|Fu(k)|>0.5max(|Fu(1)|,|Fu(2)|,......,|Fu(nT-1)|) (12)
If | Fu (k) | satisfies the above condition, the site excellent period is calculated by the following formula:
T0=0.5t0+kΔt (13)。
the invention has the beneficial effects that: the method adopts a simple algorithm with an excellent period to obtain an estimated value of the excellent period, and sets all calculation parameters except the field stratum information according to the value; taking the field soil as a seismic wave signal filter, calculating a ground seismic time course under the action of pulse waves, and performing Fourier transform on the ground seismic time course to obtain an amplitude spectrum; by judging the relation between the maximum value of the amplitude spectrum and the time (or frequency), a high-precision excellent period calculation value is obtained.
Drawings
FIG. 1 is a field stratum model in the method for simplifying calculation of excellent period of field based on stratum information;
FIG. 2 is a simplified calculation method of a distinguished period of a field based on stratum information, which is disclosed by embodiment 1, of an original ground earthquake time course;
FIG. 3 is a simplified calculation method of a distinguished period of a field based on stratum information, which is used for embodiment 1 of the invention, and the earthquake time course after attenuation processing is carried out;
FIG. 4 is a Fourier amplitude spectrum of seismic time course of embodiment 1 of the method for simplifying calculation of excellent period of field based on stratum information;
FIG. 5 is a simplified calculation method of excellent period of field based on stratum information, in accordance with embodiment 2, the time course of original ground earthquake;
FIG. 6 is a simplified calculation method of a site excellent period based on stratum information according to the invention, and the seismic time course after attenuation processing is performed in embodiment 2;
FIG. 7 is a Fourier amplitude spectrum of seismic time course of embodiment 2 of the method for simplifying calculation of excellent period of field based on formation information;
FIG. 8 is a simplified calculation method of excellent period of field based on stratum information in accordance with embodiment 3 of the present invention;
FIG. 9 is a simplified calculation method of a distinguished period of a field based on stratum information, which is shown in embodiment 3, and the seismic time course after attenuation processing;
FIG. 10 is a Fourier amplitude spectrum of seismic time course according to embodiment 3 of the method for simplifying calculation of excellent period of field based on formation information;
FIG. 11 is a simplified calculation method of excellent period of field based on stratum information, embodiment 4 of the invention is an original ground earthquake time course;
FIG. 12 is a simplified calculation method of a distinguished period of a field based on formation information, which is shown in FIG. 4, and the seismic time course after attenuation processing;
FIG. 13 is a Fourier amplitude spectrum of seismic time interval according to embodiment 4 of the method for simplifying calculation of excellent period of field based on formation information;
FIG. 14 is a simplified calculation method of excellent period of field based on stratum information, embodiment 5 of the invention is an original ground earthquake time course;
FIG. 15 is a simplified calculation method of a distinguished period of a field based on formation information, which is shown in FIG. 5, and shows a seismic time course after attenuation processing;
fig. 16 is a fourier amplitude spectrum of seismic time course of embodiment 5 of the method for simplifying calculation of the excellent period of the field based on the formation information.
Detailed Description
The present invention will be described in detail below with reference to the accompanying drawings and specific embodiments.
The invention relates to a simplified calculation method of a site excellent period based on stratum information, which is implemented according to the following steps:
Site excellent period estimation value t0The formula is as follows:
t0=4(h1/v1+h2/v2+……+hn/vn) (1)。
earthquakeDuration of wave tlThe formula is as follows:
tl≥10 t0 (2)
the time sampling interval Δ t is calculated by the following equation (3):
Δt<0.02min(h1/v1,h2/v2,……,hn/vn) (3)。
sampling point number nt of ground seismic wave holding time and sampling point number nt of seismic wave propagation time in stratum iiThe sampling point number nT of the Fourier cycle spectrum of the ground seismic wave is calculated by the following formula:
nt=Int(tl/Δt) (4)
nti=Int(hi/vi/Δt) (5)
nT=Int(t0/Δt) (6)
wherein Int is an integer function, i is 1, 2.
Reflection coefficient RiThe formula is as follows:
Ri=(ρi+1vi+1-ρivi)/(ρi+1vi+1+ρivi) Wherein, i is 1, 2, a. (7)
Seismic time-path u0 dThe formula is as follows:
u0 d(j)=-u1 p(j-nt1)
u1 d(j)=(1+R1)u0 d(j-nt1)-R1u2 p(j-nt2)
u1 p(j)=R1u0 d(j-nt1)+(1-R1)u2 p(j-nt2)
u2 d(j)=(1+R2)u1 d(j-nt2)-R2u3 p(j-nt3)
u2 p(j)=R2u1 d(j-nt2)+(1-R2)u3 p(j-nt3)
......
ui d(j)=(1+Ri)ui-1 d(j-nti)-Riui+1 p(j-nti+1)
ui p(j)=Riui-1 d(j-nti)+(1-Ri)ui+1 p(j-nti+1)
......
un-1 d(j)=(1+Rn-1)un-2 d(j-ntn-1)-Rn-1un p(j-ntn)
un-1 p(j)=Rn-1un-2 d(j-ntn-1)+(1-Rn-1)un p(j-ntn)
un p(j)=un d(j-ntn)
wherein u is0 dIs the ground seismic wave time course; u. ofi d、ui pThe seismic wave time courses of the lower interface and the upper interface of the formation i are respectively, wherein i is 1, 2. u. of0 d、ui d、ui pAre time series, j represents a time sample point serial number, j is 0, 1,.. multidot.nt-1,j-0 means that the time is equal to zero.
the seismic wave time course u is processed according to the following formula0 dCarrying out attenuation treatment according to time:
u(j)=u0 d(j)exp[-βjΔt] (9)
wherein the attenuation coefficient beta is-ln delta/tl(ii) a Delta is a calculation precision value for controlling the attenuation coefficient and satisfies 0<δ<0.01;j=0,1,......,nt-1。
performing Discrete Fourier Transform (DFT) on the seismic wave time interval u (j) according to the following formula:
wherein i is an imaginary unit; k is 0, 1,.. cnt-1.
Step 9, judging whether the amplitude spectrum | Fu (k) | has a maximum value, if not, expanding nT by 20%, and returning to the step 8; if yes, determining the site excellent period T0;
Whether the amplitude spectrum | fu (k) | has a maximum value is judged according to the following formula:
| Fu (k) | > | Fu (k-1) | and | Fu (k) | > | Fu (k +1) | (11)
|Fu(k)|>0.5max(|Fu(1)|,|Fu(2)|,......,|Fu(nT-1)|) (12)
If | Fu (k) | satisfies the above condition, the site excellent period is calculated by the following formula:
T0=0.5t0+kΔt (13)。
example 1
n=3
ρ1=2041 v1=200 h1=4
ρ2=2143 v2=300 h2=4
ρ3=2041 v3=200 h3=12
t0=4×(4/200+4/300+12/200)=0.3733
tl=10t0=10×0.3733=3.7333
Δt=0.0001<0.02min(4/200,4/300,12/200)=0.01333
nt=Int(tl/Δt)=Int(3.733/0.0001)=37333
nt1=Int(hi/vi/Δt)=Int(4/200/0.0001)=200
nt2=Int(hi/vi/Δt)=Int(4/300/0.0001)=133
nt3=Int(hi/vi/Δt)=Int(12/200/0.0001)=600
nT=Int(t0/Δt)=Int(0.3733/0.0001)=3733
R1=(2143×300-2041×200)/(2143×300+2041×200)=0.2233
R2=(2041×200-2143×300)/(2041×200+2143×300)=-0.2233
R3=1
Let nt, nt1,nt2,nt3And R1,R2,R3Substituting the correlated parameters into a time recurrence equation to obtain a ground seismic time course u0 d(j) J is 0, 1, a.
taking δ to be 0.005, the attenuation coefficient is calculated:
β=-lnδ/tl=-ln0.005/3.7333=1.4192
according to the formula u0 d(j)exp[-βjΔt]And (3) carrying out attenuation treatment to obtain the earthquake time course u (j), wherein j is 0, 1. (as shown in FIG. 3)
discrete Fourier Transform (DFT) is performed on the seismic time interval u to obtain Fourier spectrum Fu (k), k being 0, 1.
Step 9, judging whether the amplitude spectrum | Fu (k) | has a maximum value, if not, expanding nT by 20%, and returning to the step 8; if yes, determining the site excellent period T0;
In descending order, the maxima of the amplitude spectrum | fu (k) | are found, in 3731, 3730, 3729. According to calculation, k is 2069, and the site excellent period can be calculated as follows:
T0=0.5t0+kΔt=0.5×0.3733+2069×0.0001=0.3935
the excellent cycle estimates, excellent cycle calculations and the main calculation parameters are detailed in table 2.
Example 2
The formation parameters of a 2-layer soil field of 20m are shown in Table 3. The original ground seismic time interval, the seismic time interval after attenuation processing and the fourier amplitude spectrum thereof calculated according to the procedure of example 1 are shown in fig. 5, 6 and 7, and the estimated excellent period value, the calculated excellent period value and the main calculation parameters are detailed in table 4. As can be seen from the table, although the excellent cycle estimates, the main calculation parameters, are the same as in example 2, the excellent cycle calculations are not the same, 0.3935s and 0.4027s, respectively.
Example 3
The formation parameters of a certain 5-layer soil field are 10.7m, and are shown in Table 5. The original ground seismic time interval, the seismic time interval after attenuation processing and the fourier amplitude spectrum thereof calculated according to the procedure of example 1 are shown in fig. 8, 9 and 10, and the estimated excellent period value, the calculated excellent period value and the main calculation parameters are detailed in table 6. As can be seen, the calculated values for the excellent cycle are close to the estimated values, 0.2092s and 0.2345, respectively.
Example 4
The soil thickness of a certain site is 126m, the number of strata is more, the total number of the strata is 12, and the parameters of each stratum are shown in the table 7. The original ground seismic time interval, the seismic time interval after attenuation processing and the fourier amplitude spectrum thereof calculated according to the procedure of example 1 are shown in fig. 11, 12 and 13, and the estimated excellent period value, the calculated excellent period value and the main calculation parameters are detailed in table 8. The ground earthquake time course is complex (see fig. 11), and the excellent period reaches 1.3195 s.
Example 5
The soil thickness of a certain site is 250m, the stratum can reach up to 20 layers, and the parameters of each soil layer are shown in a table 9. The original ground seismic time interval, the seismic time interval after attenuation processing and the fourier amplitude spectrum thereof calculated according to the procedure of example 1 are shown in fig. 14, 15 and 16, and the estimated excellent period value, the calculated excellent period value and the main calculation parameters are shown in table 10. Due to the fact that the field is layered and thick, the ground earthquake time course is very complex (see fig. 14), the excellent period is large, and the value reaches 2.0319 s.
The simplified calculation method of the site excellent period based on the stratum information has the advantages that: the method adopts a simple algorithm with an excellent period to obtain an estimated value of the excellent period, and sets all calculation parameters except the field stratum information according to the value; taking the field soil as a seismic wave signal filter, calculating a ground seismic time course under the action of pulse waves, and performing Fourier transform on the ground seismic time course to obtain an amplitude spectrum; by judging the relation between the maximum value of the amplitude spectrum and the time (or frequency), a high-precision excellent period calculation value is obtained.
Table 1 example 1 formation parameters
Table 2 example 1 calculation of parameters
Table 3 example 2 formation parameters
Table 4 example 2 calculation of parameters
Table 5 example 3 formation parameters
Table 6 example 3 calculation of parameters
Table 7 example 4 formation parameters
Table 8 example 4 calculation of parameters
Table 9 example 5 formation parameters
Table 10 example 5 calculation of parameters
Claims (1)
1. A simplified calculation method for a site excellent period based on stratum information is characterized by comprising the following steps:
step 1, determining the number n of the site and the stratum density rho according to the investigation result of geotechnical engineeringiShear velocity viAnd the thickness h of the formationi(ii) a Wherein, i is 1, 2.. said., n;
step 2, calculating a site excellent period estimated value t0;
In the step 2, the site excellent period estimated value t0The formula is as follows:
t0=4(h1/v1+h2/v2+……+hn/vn) (1);
step 3, calculating the seismic wave duration tlAnd a time sampling interval Δ t;
in the step 3, the seismic wave duration tlThe formula is as follows:
tl≥10t0 (2)
the time sampling interval Δ t is calculated by the following equation (3):
Δt<0.02min(h1/v1,h2/v2,……,hn/vn) (3);
step 4, calculating the sampling point number nt of the ground seismic wave holding time and the sampling point number nt of the seismic wave propagation time in the stratum iiThe number nT of sampling points of the Fourier periodic spectrum of the ground seismic wave;
in the step 4, the number of sampling points nt of the ground seismic wave holding time and the number of sampling points nt of the seismic wave propagation time in the stratum iiThe sampling point number nT of the Fourier cycle spectrum of the ground seismic wave is calculated by the following formula:
nt=Int(tl/Δt) (4)
nti=Int(hi/vi/Δt) (5)
nT=Int(t0/Δt) (6)
wherein Int is an integer function, i is 1, 2.
Step 5, calculating the reflection coefficient R of the seismic waves at the boundary interface of the stratum i and the stratum i +1i;
In the step 5, the reflection coefficient RiThe formula is as follows:
Ri=(ρi+1vi+1-ρivi)/(ρi+1vi+1+ρivi) Wherein, i is 1, 2, a. (7)
Rn1, wherein i ═ n; (8);
step 6, calculating the time course u of the ground seismic wave0 d;
In step 6, seismic wave time course u0 dThe formula is as follows:
u0 d(j)=-u1 p(j-nt1)
u1 d(j)=(1+R1)u0 d(j-nt1)-R1u2 p(j-nt2)
u1 p(j)=R1u0 d(j-nt1)+(1-R1)u2 p(j-nt2)
u2 d(j)=(1+R2)u1 d(j-nt2)-R2u3 p(j-nt3)
u2 p(j)=R2u1 d(j-nt2)+(1-R2)u3 p(j-nt3)
......
ui d(j)=(1+Ri)ui-1 d(j-nti)-Riui+1 p(j-nti+1)
ui p(j)=Riui-1 d(j-nti)+(1-Ri)ui+1 p(j-nti+1)
......
un-1 d(j)=(1+Rn-1)un-2 d(j-ntn-1)-Rn-1un p(j-ntn)
un-1 p(j)=Rn-1un-2 d(j-ntn-1)+(1-Rn-1)un p(j-ntn)
un p(j)=un d(j-ntn)
wherein u is0 dIs the ground seismic wave time course; u. ofi d、ui pThe seismic wave time courses of the lower interface and the upper interface of the formation i are respectively, wherein i is 1, 2. u. of0 d、ui d、ui pAre time series, j represents a time sample serial number, j is 0, 1,.. said, nt-1, j is 0 represents that the time is equal to zero;
step 7, seismic wave time course u0 dCarrying out attenuation treatment according to time;
in the step 7, the seismic wave time travel u is processed according to the following formula0 dCarrying out attenuation treatment according to time:
u(j)=u0 d(j)exp[-βjΔt] (9)
wherein the attenuation coefficient beta is-ln delta/tl(ii) a Delta is a calculation precision value for controlling the attenuation coefficient and satisfies 0<δ<0.01;j=0,1,......,nt-1;
Step 8, performing discrete Fourier transform on the seismic wave time course u (j) to obtain Fu (k);
in the step 8, the seismic wave time interval u (j) is subjected to discrete fourier transform according to the following formula:
wherein i is an imaginary unit; k is 0, 1,.. cnt-1;
step 9, judging whether the amplitude spectrum | Fu (k) | has a maximum value, if not, expanding nT by 20%, and returning to the step 8; if yes, determining the site excellent period T0;
In step 9, it is determined whether the amplitude spectrum | fu (k) | has a maximum value according to the following formula:
| Fu (k) | > | Fu (k-1) | and | Fu (k) | > | Fu (k +1) | (11)
|Fu(k)|>0.5max(|Fu(1)|,|Fu(2)|,......,|Fu(nT-1)|) (12)
If | Fu (k) | satisfies the above condition, the site excellent period is calculated by the following formula:
T0=0.5t0+kΔt (13)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910903742.8A CN110609323B (en) | 2019-09-24 | 2019-09-24 | Simplified calculation method for site excellent period based on stratum information |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910903742.8A CN110609323B (en) | 2019-09-24 | 2019-09-24 | Simplified calculation method for site excellent period based on stratum information |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110609323A CN110609323A (en) | 2019-12-24 |
CN110609323B true CN110609323B (en) | 2022-02-22 |
Family
ID=68892038
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910903742.8A Active CN110609323B (en) | 2019-09-24 | 2019-09-24 | Simplified calculation method for site excellent period based on stratum information |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110609323B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112444875B (en) * | 2020-10-28 | 2023-10-03 | 西安理工大学 | Method for obtaining accurate solution of site excellent period |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4751655A (en) * | 1986-06-11 | 1988-06-14 | Amoco Corporation | Method of reconstituting seismic data |
CN108828661A (en) * | 2018-03-06 | 2018-11-16 | 西安理工大学 | Method based on seismic pulse response spectra measurement site predominant period |
CN108957537A (en) * | 2018-06-20 | 2018-12-07 | 西安理工大学 | Tilt the calculation method of site predominant period under ground seismic wave function |
-
2019
- 2019-09-24 CN CN201910903742.8A patent/CN110609323B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4751655A (en) * | 1986-06-11 | 1988-06-14 | Amoco Corporation | Method of reconstituting seismic data |
CN108828661A (en) * | 2018-03-06 | 2018-11-16 | 西安理工大学 | Method based on seismic pulse response spectra measurement site predominant period |
CN108957537A (en) * | 2018-06-20 | 2018-12-07 | 西安理工大学 | Tilt the calculation method of site predominant period under ground seismic wave function |
Non-Patent Citations (2)
Title |
---|
场地土卓越周期变化的试验分析;王家全等;《广西大学学报(自然科学版)》;20100820;第35卷(第04期);第684-687、692页 * |
成层土场地地震效应的时程算法研究;范留明等;《岩土力学》;20090910;第30卷(第09期);第2564-2568页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110609323A (en) | 2019-12-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107967240B (en) | Earthquake-resistant time-course acceleration optimization algorithm based on artificial earthquake motion generation process | |
Vaid et al. | Effect of static shear on resistance to liquefaction | |
Gourdin et al. | Flow stress of OFE copper at strain rates from 10− 3 to 104s− 1: grain-size effects and comparison to the mechanical threshold stress model | |
CN107505654A (en) | Full waveform inversion method based on earthquake record integration | |
CN110657935B (en) | Seismic wave acceleration integral processing method and system | |
CN110609323B (en) | Simplified calculation method for site excellent period based on stratum information | |
CN114442153B (en) | Near fault earthquake motion fitting method | |
CN106886047A (en) | A kind of method of receiver function and gravity Inversion CRUSTAL THICKNESS and ripple ratio | |
CN113640881B (en) | Multi-offset-distance two-dimensional transverse high-resolution transient surface wave detection method | |
CN108828661B (en) | Method based on seismic pulse response spectra measurement site predominant period | |
CN109187744B (en) | Shear wave velocity estimation method based on static cone penetration test | |
CN111045076A (en) | Multi-mode Rayleigh wave frequency dispersion curve parallel joint inversion method | |
CN116046307A (en) | Identification method suitable for earthquake damage mode of vibrating table with tunnel slope | |
CN103728661B (en) | A kind of high accuracy inverse Q filtering Processing Seismic Data | |
CN115235712A (en) | Multi-anchor-point anchor cable anti-slide pile damage deformation identification method suitable for vibration table test | |
CN110320400B (en) | Voltage flicker envelope parameter extraction method for quasi-synchronous sampling and improved energy operator | |
CN108957537A (en) | Tilt the calculation method of site predominant period under ground seismic wave function | |
Iyengar et al. | Generation of spectrum compatible accelerograms | |
CN114527506A (en) | Method for simplifying artificial generation of fault seismic oscillation | |
CN109541304B (en) | Power grid higher order weak amplitude harmonic detection method based on six-term minimum sidelobe window interpolation | |
CN105445724B (en) | Single hydrophone free field passive ranging method | |
CN110456413A (en) | The method for prejudging rock burst type | |
Arulanandan et al. | Level ground soil-liquefaction analysis using in situ properties: II | |
CN117761787A (en) | High-resolution processing method for defocusing of brine mud array induction logging data | |
CN113552633B (en) | Elastic wave frequency dispersion pressing method for optimizing difference coefficient and longitudinal and transverse wave separation FCT |
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 |