CN113552628B - Method for calculating height of earthquake wave - Google Patents
Method for calculating height of earthquake wave Download PDFInfo
- Publication number
- CN113552628B CN113552628B CN202110817856.8A CN202110817856A CN113552628B CN 113552628 B CN113552628 B CN 113552628B CN 202110817856 A CN202110817856 A CN 202110817856A CN 113552628 B CN113552628 B CN 113552628B
- Authority
- CN
- China
- Prior art keywords
- earthquake
- calculating
- seismic wave
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 230000001133 acceleration Effects 0.000 claims abstract description 22
- 238000004364 calculation method Methods 0.000 claims abstract description 13
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 12
- 230000001960 triggered effect Effects 0.000 claims abstract description 5
- 230000005484 gravity Effects 0.000 claims abstract description 4
- 230000002265 prevention Effects 0.000 abstract description 6
- 230000004888 barrier function Effects 0.000 abstract description 5
- 230000009182 swimming Effects 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 241000209140 Triticum Species 0.000 description 1
- 235000021307 Triticum Nutrition 0.000 description 1
- BBBFJLBPOGFECG-VJVYQDLKSA-N calcitonin Chemical compound N([C@H](C(=O)N[C@@H](CC(C)C)C(=O)NCC(=O)N[C@@H](CCCCN)C(=O)N[C@@H](CC(C)C)C(=O)N[C@@H](CO)C(=O)N[C@@H](CCC(N)=O)C(=O)N[C@@H](CCC(O)=O)C(=O)N[C@@H](CC(C)C)C(=O)N[C@@H](CC=1NC=NC=1)C(=O)N[C@@H](CCCCN)C(=O)N[C@@H](CC(C)C)C(=O)N[C@@H](CCC(N)=O)C(=O)N[C@@H]([C@@H](C)O)C(=O)N[C@@H](CC=1C=CC(O)=CC=1)C(=O)N1[C@@H](CCC1)C(=O)N[C@@H](CCCNC(N)=N)C(=O)N[C@@H]([C@@H](C)O)C(=O)N[C@@H](CC(N)=O)C(=O)N[C@@H]([C@@H](C)O)C(=O)NCC(=O)N[C@@H](CO)C(=O)NCC(=O)N[C@@H]([C@@H](C)O)C(=O)N1[C@@H](CCC1)C(N)=O)C(C)C)C(=O)[C@@H]1CSSC[C@H](N)C(=O)N[C@@H](CO)C(=O)N[C@@H](CC(N)=O)C(=O)N[C@@H](CC(C)C)C(=O)N[C@@H](CO)C(=O)N[C@@H]([C@@H](C)O)C(=O)N1 BBBFJLBPOGFECG-VJVYQDLKSA-N 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
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/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
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 method for calculating the height of a seismic surge, which comprises the following steps of S1, acquiring a seismic wave time course curve; s2, obtaining a first-order main frequency P of the seismic wave according to a seismic wave time-course curve, wherein the unit is Hz; the peak acceleration a of the seismic wave is expressed as m/s; the maximum earth surface movement velocity v triggered by earthquake is m/s; s3, calculating a dimensionless quantity G of the peak acceleration; the calculating process is the ratio of the peak acceleration a of the seismic wave to the gravity acceleration g; s4, obtaining the average water depth h of the reservoir to be predicted 0 The unit is m; s5, calculating the Froude number F of the seismic waves, wherein the calculation process is as follows:the method comprises the steps of carrying out a first treatment on the surface of the S6, calculating the earthquake surge height of the reservoir to be predicted caused by the earthquakehThe unit is m. The method comprehensively considers the peak acceleration of the earthquake waves, the main frequency of the earthquake waves, the maximum surface movement speed and the initial water depth, fits to obtain the earthquake surge height with higher accuracy, provides accurate prediction for the breaking of the dam body and the barrier lake in the earthquake, improves the disaster prediction prevention and control capability, and reduces the loss.
Description
Technical Field
The invention relates to the field of disaster prediction prevention and control, in particular to a method for calculating the height of earthquake waves.
Background
Earthquake is one of the main triggers for dam break of barrier lake and earth dam reservoirs.
In the past, people mainly pay attention to the influence of the earthquake on the stability of a reservoir dam body, but the image data of the '5.12' Wenchuan earthquake show that the earthquake can trigger the reservoir surge phenomenon, so that the damage of the earthquake surge cannot be ignored in the investigation of a reservoir-lake burst risk source.
In the Wenchuan earthquake, monitoring equipment beside a purple apron reservoir shoots videos of the pool water surge to strike the pool bank, and the videos indicate that the surge height triggered by the earthquake can reach more than 3 m.
The surge of this magnitude can completely trigger a burst for naturally occurring barrier lakes, earth dam reservoirs.
However, in the current disaster prediction control of the dam body and the barrier lake, the earthquake surge height is generally estimated only according to the designed earthquake intensity, the water depth in front of the dam or the dam body height, the estimated value range is large, the estimated precision is not high, and the problem that the influence of each element of the earthquake on the earthquake surge height is not comprehensively considered is solved.
For example, the estimated range of the surge height is generally uniformly determined to be 0.5-1.5m in China according to the designed seismic intensity and the depth of water in front of the dam, and the estimated surge height of the earthquake is estimated according to the dam height of 1% in Japan.
The method is characterized in that the seismic surge height is roughly estimated through a single dimension, and the situation of the seismic surge cannot be accurately estimated, so that the disaster prevention and control possibly caused by the seismic surge is not effectively guided.
Disclosure of Invention
The invention aims to provide a calculation method for the height of an earthquake surge, which can accurately predict the height of the earthquake surge and improve the prevention, judgment and coping capacity for disasters caused by the earthquake surge.
In order to achieve the above purpose, the invention adopts the following technical scheme:
the invention relates to a method for calculating the height of an earthquake surge, which comprises the following steps:
s1, acquiring a seismic wave time course curve;
s2, obtaining a first-order main frequency P of the seismic wave according to a seismic wave time-course curve, wherein the unit is Hz; peak acceleration of seismic wave a in m/s 2 The method comprises the steps of carrying out a first treatment on the surface of the The maximum earth surface movement velocity v triggered by earthquake is m/s;
s3, calculating a dimensionless quantity G of the peak acceleration; the calculating process is the ratio of the peak acceleration a of the seismic wave to the gravity acceleration g;
s4, obtaining the average water depth h of the reservoir to be predicted 0 The unit is 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 earthquakehThe unit is m; the calculation process is as follows:
。
preferably, in S2, the first order dominant frequency P of said seismic wave is obtained by performing a Fast Fourier Transform (FFT) on the seismic wave time-course curve; peak acceleration a of seismic wave; the earthquake triggers the maximum surface movement velocity v.
The method has the advantages that the peak acceleration of the earthquake waves, the main frequency of the earthquake waves, the maximum earth surface movement speed and the initial water depth are comprehensively considered, the earthquake wave height with higher accuracy is obtained by fitting, the accurate prediction is provided for the breaking of the dam body and the barrier lake in the earthquake, the disaster prediction prevention and control 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 test measurements of example 2 versus the calculated values of the method of the present invention.
Detailed Description
The following description of the technical solutions in the embodiments of the present invention will be clear and complete, and it is obvious that the described embodiments are only some embodiments of the present invention, but not all embodiments.
All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Example 1:
as shown in FIG. 1, the method for calculating the height of the earthquake surge comprises the following specific steps:
s1, acquiring a seismic wave time course curve;
s2, performing Fast Fourier Transform (FFT) on a seismic wave time course curve to obtain a first-order main frequency P of the seismic wave, wherein the unit is Hz; the peak acceleration a of the seismic wave is expressed as m/s; the maximum earth surface movement velocity v triggered by earthquake is m/s;
s3, calculating a dimensionless quantity G of the peak acceleration; the calculating process is the ratio of the peak acceleration a of the seismic wave to the gravity acceleration g;
s4, obtaining the average water depth h of the reservoir to be predicted 0 The unit is 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 earthquakehThe unit is m; the calculation process is as follows:
。
example 2:
laboratory simulation verification of the seismic surge height calculation method provided by the invention:
the method comprises the steps that earthquake waves are sent to a one-dimensional electrohydraulic servo driving type vibrating table through a computer, the vibrating frequency of the vibrating table is 0.3 Hz-15 Hz, the maximum vibrating acceleration is 1.2g, the acceleration repetition precision is +/-3%, and the maximum amplitude is +/-100 mm; the size of the vibrating table water tank is 4m multiplied by 2m, and the load is 25t.
Meanwhile, the seismic waveforms used in the experiment include a bedroom wave, a Mexicoli wave and a through wheat wave.
The waveform of the seismic wave is corrected before each experiment begins.
When the accuracy of the real vibration waveform fed back by the vibration table and the seismic wave output by the computer is more than 95%, the real vibration waveform of the vibration table is considered to be the seismic wave output by the computer.
During formal experiments, the computer transmits earthquake waves to the vibrating table, the wave height acquisition system is used for recording the height value of the earthquake waves generated by the water tank of the vibrating table, and each experiment is carried out twice, so that errors are reduced.
The experimental data are collated in table 1 and compared with the predicted value of the seismic surge height calculation method according to the present invention, and as shown in fig. 2, it can be seen that 50 sets of errors are within 15% in 65 sets of experimental data.
Therefore, the invention can obtain more accurate prediction of the landslide surge height, and can greatly improve the applicability of active disaster prevention.
TABLE 1
Maximum wave height (cm) of earthquake surge under various test working conditions
Example 3:
the invention relates to a case verification method for calculating the height of a seismic surge.
The time of 15 minutes and 42 seconds from 4.4.2010, the earthquake surge occurs in a swimming pool in the ink-west kari city caused by the 7.2 earthquake in california.
And analyzing the swimming pool seismic surge video frame by frame in the current earthquake shot by monitoring beside the swimming pool through Premiere Pro video processing software, and estimating the seismic surge heights at different moments by taking a bank table and a bank chair as a reference object to obtain a time course curve of the seismic surge heights.
From the time course of the seismic surge height, it can be seen that this earthquake causes the seismic surge height of the pool to reach a maximum height of about 0.6m at 23s of video.
According to the earthquake surge height calculation method, the earthquake surge height of the swimming pool can be predicted to be 0.57 m according to the earthquake vibration data of one earthquake monitoring point which is closer to the ink-jet Karil city and the water depth of the swimming pool obtained through video analysis, and the calculation result is basically consistent with the actual situation.
Therefore, the method can improve the accuracy of earthquake surge height prediction, improve disaster prediction control capability and reduce loss.
Claims (2)
1. A method for calculating the height of an earthquake surge is characterized by comprising the following steps: the method comprises the following steps:
s1, acquiring a seismic wave time course curve;
s2, obtaining a first-order main frequency P of the seismic wave according to a seismic wave time-course curve, wherein the unit is Hz; peak acceleration of seismic wave a in m/s 2 The method comprises the steps of carrying out a first treatment on the surface of the The maximum earth surface movement velocity v triggered by earthquake is m/s;
s3, calculating a dimensionless quantity G of the peak acceleration; the calculation process is that the peak acceleration a of the seismic wave is divided by the gravity acceleration g;
s4, obtaining the average water depth h of the reservoir to be predicted 0 The unit is 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 earthquakehThe unit is m; the calculation process is as follows:
。
2. the method of calculating a seismic surge height of claim 1, wherein: s2, obtaining a first-order main frequency P of the seismic wave by performing Fast Fourier Transform (FFT) on a seismic wave time-course curve; peak acceleration a of seismic wave; the earthquake triggers the maximum surface movement 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 CN113552628A (en) | 2021-10-26 |
| CN113552628B true 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 (6)
| 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 |
| 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 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108254782B (en) * | 2018-02-09 | 2019-11-05 | 中国地质大学（北京） | A kind of acquisition methods and system of side slope eaerthquake damage failure probability |
-
2021
- 2021-07-20 CN CN202110817856.8A patent/CN113552628B/en active Active
Patent Citations (6)
| 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 |
| 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 |
|---|
| 严佳宏等.基于振动台造波模拟实验的地震涌浪频谱与共振现象研究.地球物理学进展.2016,(第06期),全文. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113552628A (en) | 2021-10-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| McAllister et al. | Laboratory recreation of the Draupner wave and the role of breaking in crossing seas | |
| CN104266819A (en) | Device for simulating liquid sloshing generated under random wave action and wave generation method thereof | |
| D'Alessandro et al. | Wave–dune interaction and beach resilience in large-scale physical model tests | |
| Bolshakova et al. | The properties of co-seismic deformations of the ocean bottom as indicated by the slip-distribution data in tsunamigenic earthquake sources | |
| CN113552628B (en) | Method for calculating height of earthquake wave | |
| Zhao et al. | A numerical study of tsunami wave impact and run-up on coastal cliffs using a CIP-based model | |
| Vestbøstad et al. | Column slamming loads on a TLP from steep and breaking waves | |
| 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 | |
| St-Germain et al. | Numerical modeling of tsunami-induced hydrodynamic forces on onshore structures using SPH | |
| Gao et al. | Experimental investigation of breaking regular waves slamming on offshore wind jacket structure | |
| CN106910507B (en) | Detection and identification method and system | |
| Bullock et al. | Violent breaking-wave impacts. Part 4: A detailed analysis and comparison of field and 1: 4 scale measurements on sloping and vertical walls including the influence of air and scale effects | |
| CN104596735B (en) | Optimal arrangement method of bridge group | |
| Zhang et al. | Classification of breaking wave impact loads on a fixed surface-piercing square column with an overhanging deck | |
| CN112379423B (en) | Method and device for determining moment and vibration level of microseism event | |
| Ito et al. | Hybrid two-phase flow measurements in a narrow channel using neutron radiography and liquid film sensor | |
| Merchant et al. | A dynamic analysis of an elastic model of the human head | |
| Chopakatla | A CFD model for wave transformations and breaking in the surf zone | |
| Dassanayake et al. | Efficacy of analysis techniques in assessing broken wave loading on a cylinder upon a shoal | |
| CN117252013B (en) | Method for constructing broadband limited seismic source model | |
| Birknes-Berg et al. | Methods for establishing governing deck impact loads in irregular waves | |
| JP6314065B2 (en) | Vibration response evaluation system, vibration response evaluation method, and vibration response evaluation program | |
| CN116804768B (en) | Earthquake motion determination method for earthquake-resistant analysis of near-span fault structure | |
| CN111783542B (en) | Method and device for automatically extracting stress wave reflection period | |
| CN117571213B (en) | Submarine carbon dioxide sealing monitoring method, system, equipment and storage medium based on sound-image combination |
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 |