CN113552628A - Earthquake surge height calculation method - Google Patents
Earthquake surge height calculation method Download PDFInfo
- Publication number
- CN113552628A CN113552628A CN202110817856.8A CN202110817856A CN113552628A CN 113552628 A CN113552628 A CN 113552628A CN 202110817856 A CN202110817856 A CN 202110817856A CN 113552628 A CN113552628 A CN 113552628A
- Authority
- CN
- China
- Prior art keywords
- earthquake
- calculating
- unit
- seismic
- surge
- 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
- 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. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Abstract
The invention discloses a method for calculating earthquake surge height, comprising the following steps of S1, acquiring an earthquake wave time-course curve; s2, obtaining a first-order dominant frequency P of the seismic waves according to the seismic wave time-course curve, wherein the unit is Hz; the unit of the seismic wave peak acceleration a is m/s; triggering the maximum earth surface motion speed v by the earthquake, wherein the unit is m/s; s3, calculating a dimensionless quantity G of the peak acceleration; the calculation process is the ratio of the seismic wave peak acceleration a to the gravity acceleration g; s4, obtaining the average water depth h of the reservoir to be predicted0In the unit of m; s5, calculating the Froude number F of the seismic waves, wherein the calculation process is as follows:
(ii) a S6, calculating the earthquake surge height of the reservoir to be predicted caused by the earthquakehIn the unit m. The earthquake surge height with higher accuracy is obtained by fitting by comprehensively considering the earthquake wave peak acceleration, the earthquake wave dominant frequency, the maximum earth surface movement speed and the initial water depth, so that accurate prediction is provided for the collapse of the dam body and the dammed lake in the earthquake, and the disaster prediction and prevention are improvedTreating capability and reducing loss.
Description
Technical Field
The invention relates to the field of disaster prediction and prevention, in particular to a method for calculating earthquake surge height.
Background
Earthquake is one of the main triggering factors for the collapse of the barrage lake and the earth dam reservoir.
In the past, people mainly pay attention to the influence of earthquakes on the stability of reservoir dams, but the image data of the 5.12 Wenchuan earthquake shows that the reservoir surge phenomenon can be triggered by the earthquake, so that the damage of the earthquake surge cannot be ignored in the investigation of reservoir and lake burst risk sources.
In Wenchuan earthquake, a video that the reservoir water surge beats the bank of the reservoir is shot by a monitoring device beside a purple plateau reservoir, and the video shows that the surge height triggered by the earthquake can reach more than 3 m.
The surge of this magnitude can completely cause the burst for naturally formed barrage lakes and earth dam reservoirs.
However, in the current disaster prediction and prevention of dams and barrage lakes, the earthquake surge height is usually estimated only according to the designed earthquake intensity, the water depth in front of the dam or the dam height roughly, and the problems that the estimated value range is large, the estimation precision is not high, and the influence of each earthquake element on the earthquake surge height is not comprehensively considered exist.
For example, the surging height estimation range is generally uniformly determined to be 0.5-1.5m in China according to the design seismic intensity and the water depth before the dam, and the seismic surging height is estimated according to 1% of the dam height in Japan.
The above methods are all that the earthquake surge height is roughly estimated through a single dimension, and the earthquake surge condition cannot be accurately estimated, so that no effective guidance is provided for disaster prevention and control possibly caused by earthquake surge.
Disclosure of Invention
The invention aims to provide a method for calculating the earthquake surge height, which can accurately predict the earthquake surge height and improve the prevention, judgment and coping capability of disasters caused by earthquake surge.
In order to achieve the purpose, the invention adopts the following technical scheme:
the method for calculating the earthquake surge height comprises the following steps:
s1, acquiring a seismic wave time-travel curve;
s2, obtaining a first-order dominant frequency P of the seismic waves according to the seismic wave time-course curve, wherein the unit is Hz; the unit of the acceleration a of the peak value of the seismic wave is m/s2(ii) a Triggering the maximum earth surface motion speed v by the earthquake, wherein the unit is m/s;
s3, calculating a dimensionless quantity G of the peak acceleration; the calculation process is the ratio of the seismic wave peak acceleration a to the gravity acceleration g;
s4, obtaining the average water depth h of the reservoir to be predicted0In the unit of m;
s5, calculating the Froude number F of the seismic waves, wherein the calculation process is as follows:
s6, calculating the earthquake surge height of the reservoir to be predicted caused by the earthquakehIn the unit of m; the calculation process is as follows:
preferably, in S2, a first-order dominant frequency P of the seismic waves is obtained by performing Fast Fourier Transform (FFT) on the seismic wave time-course curve; seismic wave peak acceleration a; the earthquake triggers a maximum surface motion velocity v.
The method has the advantages that the seismic wave peak acceleration, the seismic wave dominant frequency, the maximum earth surface movement speed and the initial water depth are comprehensively considered, the seismic surge height with higher accuracy is obtained through fitting, accurate prediction is provided for the collapse of the dam body and the dammed lake in the earthquake, the disaster prediction and prevention capability is improved, and the loss is reduced.
Drawings
FIG. 1 is a flow chart of a method for calculating the height of a seismic surge according to the invention.
FIG. 2 is a graph of the comparison of the values measured in the test of example 2 with the values calculated by the method of the invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be described clearly and completely below, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all embodiments.
All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Example 1:
as shown in FIG. 1, the method for calculating the earthquake surge height comprises the following specific steps:
s1, acquiring a seismic wave time-travel curve;
s2, performing Fast Fourier Transform (FFT) on the seismic wave time-course curve to obtain a first-order main frequency P of the seismic waves, wherein the unit is Hz; the unit of the seismic wave peak acceleration a is m/s; triggering the maximum earth surface motion speed v by the earthquake, wherein the unit is m/s;
s3, calculating a dimensionless quantity G of the peak acceleration; the calculation process is the ratio of the seismic wave peak acceleration a to the gravity acceleration g;
s4, obtaining the average water depth h of the reservoir to be predicted0In the unit of m;
s5, calculating the Froude number F of the seismic waves, wherein the calculation process is as follows:
s6, calculating the earthquake surge height of the reservoir to be predicted caused by the earthquakehIn the unit of m; the calculation process is as follows:
example 2:
the invention discloses a laboratory simulation verification method of a seismic surge height calculation method, which comprises the following steps:
sending the seismic waves to a one-dimensional electro-hydraulic servo drive type vibration table through a computer, wherein the vibration frequency of the vibration table is 0.3 Hz-15 Hz, the maximum vibration acceleration is 1.2g, the acceleration repetition precision is +/-3%, and the maximum amplitude is +/-100 mm; the size of the water tank of the vibration table is 4m multiplied by 2m, and the load is 25 t.
Meanwhile, the seismic waveforms used for the experiments include the Wolongbo, the Mexican-Ribo and the Tongmai-wave.
The waveform of the seismic waves was corrected before the start of each experiment.
And when the precision of the real vibration waveform fed back by the vibration table is more than 95% of the seismic wave output by the computer, the real vibration waveform of the vibration table is considered to be the seismic wave output by the computer.
During formal experiments, the computer sends seismic waves to the vibrating table, and simultaneously, the wave height acquisition system is used for recording the seismic surge height value generated by the water tank of the vibrating table, and each experiment is carried out twice to reduce errors.
The experimental data are collated into table 1 and compared with the predicted value of the seismic surge height calculation method of the invention, as shown in fig. 2, it can be seen that 50 groups of the 65 groups of experimental data have an error within 15%.
Therefore, the method can obtain more accurate landslide surge prediction height, and can greatly improve the applicability of active disaster prevention.
TABLE 1
Maximum wave height (cm) of earthquake surge under various test conditions
Example 3:
case verification of the earthquake surge height calculation method is achieved.
At 15 hours, 40 minutes and 42 seconds, 4 months and 4 days 2010, earthquake surge occurred in a swimming pool of mercapali, caused by a 7.2 earthquake in california, usa.
And analyzing the swimming pool earthquake surge videos of the earthquake in the current time, which are monitored and shot beside the swimming pool by Premiere Pro video processing software, frame by frame, and estimating the earthquake surge heights at different moments by taking the table and the chair on the shore as reference objects to obtain a time-course curve of the earthquake surge heights.
From the time course curve of the seismic surge height, it can be seen that the seismic surge height of the pool caused by this earthquake reaches a maximum height of about 0.6m at 23s of the video.
The method for calculating the earthquake surge height predicts that the earthquake surge height of the swimming pool can reach 0.57 m according to earthquake motion data of an earthquake monitoring point which is close to the mercantile city and the water depth of the swimming pool obtained by video analysis, and the calculation result is basically consistent with the actual situation.
Therefore, the method and the device can improve the accuracy of earthquake surge height prediction, improve the disaster prediction prevention and control capability and reduce loss.
Claims (2)
1. A method for calculating earthquake surge height is characterized by comprising the following steps: the method comprises the following steps:
s1, acquiring a seismic wave time-travel curve;
s2, obtaining a first-order dominant frequency P of the seismic waves according to the seismic wave time-course curve, wherein the unit is Hz; the unit of the acceleration a of the peak value of the seismic wave is m/s2(ii) a Triggering the maximum earth surface motion speed v by the earthquake, wherein the unit is m/s;
s3, calculating a dimensionless quantity G of the peak acceleration; the calculation process is that the seismic wave peak acceleration a is divided by the gravity acceleration g;
s4, obtaining the average water depth h of the reservoir to be predicted0In the unit of m;
s5, calculating the Froude number F of the seismic waves, wherein the calculation process is as follows:
s6, calculating the earthquake surge height of the reservoir to be predicted caused by the earthquakehIn the unit of m; the calculation process is as follows:
2. the method of calculating seismic surge height of claim 1, wherein: in S2, performing Fast Fourier Transform (FFT) on the seismic wave time-course curve to obtain a first-order main frequency P of the seismic wave; seismic wave peak acceleration a; the earthquake triggers a maximum surface motion velocity v.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110817856.8A CN113552628B (en) | 2021-07-20 | 2021-07-20 | Method for calculating height of earthquake wave |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110817856.8A CN113552628B (en) | 2021-07-20 | 2021-07-20 | Method for calculating height of earthquake wave |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113552628A true CN113552628A (en) | 2021-10-26 |
| CN113552628B CN113552628B (en) | 2023-08-15 |
Family
ID=78132238
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110817856.8A Active CN113552628B (en) | 2021-07-20 | 2021-07-20 | Method for calculating height of earthquake wave |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113552628B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1330771A (en) * | 1998-12-18 | 2002-01-09 | 施鲁博格控股有限公司 | Method for determination of local wave heights and acoustic sensor in marine seismic signals |
| CN107016185A (en) * | 2017-03-31 | 2017-08-04 | 中国地质科学院探矿工艺研究所 | Calculation method for prediction of peak flow of collapse flood of tillite lake |
| KR101864307B1 (en) * | 2017-08-10 | 2018-06-05 | 한국지질자원연구원 | Method of seismic survey data processing for detecting sub-surface structure and swell effect correction using gradient analysis |
| KR102003466B1 (en) * | 2019-02-08 | 2019-07-24 | 한국지질자원연구원 | Method for swell effect correction of offshore 3d seismic survey data at shallow tratum and marine 3d seismic survey mehod using the same |
| US20190250291A1 (en) * | 2018-02-09 | 2019-08-15 | China University Of Geosciences, Beijing | Method and system for acquiring probability of slope failure and destabilization caused by earthquake |
| CN111797552A (en) * | 2020-06-12 | 2020-10-20 | 中国海洋大学 | Numerical data simulation method for undulating sea surface seismic wave field based on sea wave spectrum |
| RU2738589C1 (en) * | 2020-03-25 | 2020-12-14 | Владимир Васильевич Чернявец | Method for determining tsunami hazard |
-
2021
- 2021-07-20 CN CN202110817856.8A patent/CN113552628B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1330771A (en) * | 1998-12-18 | 2002-01-09 | 施鲁博格控股有限公司 | Method for determination of local wave heights and acoustic sensor in marine seismic signals |
| CN107016185A (en) * | 2017-03-31 | 2017-08-04 | 中国地质科学院探矿工艺研究所 | Calculation method for prediction of peak flow of collapse flood of tillite lake |
| KR101864307B1 (en) * | 2017-08-10 | 2018-06-05 | 한국지질자원연구원 | Method of seismic survey data processing for detecting sub-surface structure and swell effect correction using gradient analysis |
| US20190250291A1 (en) * | 2018-02-09 | 2019-08-15 | China University Of Geosciences, Beijing | Method and system for acquiring probability of slope failure and destabilization caused by earthquake |
| KR102003466B1 (en) * | 2019-02-08 | 2019-07-24 | 한국지질자원연구원 | Method for swell effect correction of offshore 3d seismic survey data at shallow tratum and marine 3d seismic survey mehod using the same |
| RU2738589C1 (en) * | 2020-03-25 | 2020-12-14 | Владимир Васильевич Чернявец | Method for determining tsunami hazard |
| CN111797552A (en) * | 2020-06-12 | 2020-10-20 | 中国海洋大学 | Numerical data simulation method for undulating sea surface seismic wave field based on sea wave spectrum |
Non-Patent Citations (1)
| Title |
|---|
| 严佳宏等: "基于振动台造波模拟实验的地震涌浪频谱与共振现象研究", 地球物理学进展, no. 06 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113552628B (en) | 2023-08-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Cuomo et al. | Wave-in-deck loads on coastal bridges and the role of air | |
| Alsina et al. | Sediment suspension events in the inner surf and swash zone. Measurements in large-scale and high-energy wave conditions | |
| Hoechner et al. | Probabilistic tsunami hazard assessment for the Makran region with focus on maximum magnitude assumption | |
| Vestbøstad et al. | Column slamming loads on a TLP from steep and breaking waves | |
| CN113552628A (en) | Earthquake surge height calculation method | |
| Tehranirad et al. | Tsunami benchmark results for nonhydrostatic wave model NHWAVE (Version 1.1) | |
| Cecioni et al. | Real-time inversion of tsunamis generated by landslides | |
| KR101219352B1 (en) | Dynamic calculation range type numerical simulation method of tsunami based on flux | |
| Hu et al. | A study of the Coulomb stress and seismicity rate changes induced by the 2008 Mw 7.9 Wenchuan earthquake, SW China | |
| Hughes et al. | Continuity of instantaneous wave overtopping discharge with application to stream power concepts | |
| Martinelli et al. | Wave loads on exposed jetties: description of large scale experiments and preliminary results | |
| Gomez et al. | Numerical validation of an effective slender fault source solution for past tsunami scenarios | |
| Grilli et al. | Wave impact forces on mixed breakwaters | |
| Tajima et al. | Impact of dynamic morphology change on storm surge disaster risks along the Meghna Estuary | |
| KR101273977B1 (en) | Dynamic calculation range type numerical simulation method of tsunami | |
| Sim et al. | An experimental study of tsunami amplification by a coastal cliff | |
| Mikami et al. | Flow geometry of overflowing tsunamis around coastal dykes | |
| Ruol et al. | Experimental investigation on caisson breakwater sliding | |
| CN112379423A (en) | Method and device for determining moment-magnitude of microseism event | |
| Van Dang et al. | Physical model comparison of gray and green mitigation alternatives for flooding and wave force reduction in an idealized urban coastal environment | |
| Guannel et al. | Influence of velocity moments on sand bar movement during CROSSTEX | |
| JP6314065B2 (en) | Vibration response evaluation system, vibration response evaluation method, and vibration response evaluation program | |
| Majewski et al. | Laboratory studies of the impact of non-breaking waves on a horizontal emerged deck | |
| Moon et al. | Experimental Modelling of Tsunami-Like Waves Approaching the Gentle Coast of Malaysia | |
| Cao et al. | An experimental study of wave runup: cylinder fixed in waves versus cylinder surging in still water |
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 |