CN108492236A - Multiple current Tsunami disaster appraisal procedure based on Monte Carlo stochastic simulation - Google Patents
Multiple current Tsunami disaster appraisal procedure based on Monte Carlo stochastic simulation Download PDFInfo
- Publication number
- CN108492236A CN108492236A CN201810124747.6A CN201810124747A CN108492236A CN 108492236 A CN108492236 A CN 108492236A CN 201810124747 A CN201810124747 A CN 201810124747A CN 108492236 A CN108492236 A CN 108492236A
- Authority
- CN
- China
- Prior art keywords
- tsunami
- earthquake
- seismic
- wave amplitude
- source
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 62
- 238000002948 stochastic simulation Methods 0.000 title claims abstract description 20
- 230000035939 shock Effects 0.000 claims abstract description 13
- 238000004088 simulation Methods 0.000 claims abstract description 11
- 238000011156 evaluation Methods 0.000 claims abstract description 8
- 238000004458 analytical method Methods 0.000 claims abstract description 6
- 230000006870 function Effects 0.000 claims description 19
- 238000007476 Maximum Likelihood Methods 0.000 claims description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 9
- 230000004907 flux Effects 0.000 claims description 6
- 238000003325 tomography Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 5
- 230000007246 mechanism Effects 0.000 claims description 5
- 238000010008 shearing Methods 0.000 claims description 5
- 238000006073 displacement reaction Methods 0.000 claims description 4
- 240000007594 Oryza sativa Species 0.000 claims description 3
- 235000007164 Oryza sativa Nutrition 0.000 claims description 3
- 239000003653 coastal water Substances 0.000 claims description 3
- 235000009566 rice Nutrition 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims 1
- 238000010276 construction Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 10
- 230000008569 process Effects 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000012502 risk assessment Methods 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 230000001066 destructive effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 240000006829 Ficus sundaica Species 0.000 description 1
- 241001364096 Pachycephalidae Species 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000007654 immersion Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000031068 symbiosis, encompassing mutualism through parasitism Effects 0.000 description 1
- 238000013076 uncertainty analysis Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
- G06Q50/265—Personal security, identity or safety
-
- 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
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Marketing (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- Computer Security & Cryptography (AREA)
- Human Resources & Organizations (AREA)
- Educational Administration (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geophysics And Detection Of Objects (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The multiple current Tsunami disaster appraisal procedure based on Monte Carlo stochastic simulation that the invention discloses a kind of, method include:S1 establishes potential tsunami source region seismicity parameters logic tree using a variety of focal shock parameter methods of estimation;S2 establishes tsunami unit source Green's function database using linear tsunami numerical simulation;S3 generates Stochastic earthquake event set according to above-mentioned logic tree, using Monte Carlo stochastic simulation;S4 is simulated using random slippage according to tsunami unit source Green's function database and Stochastic earthquake event set and is generated tsunami wave amplitude collection;S5 is according to tsunami wave amplitude collection, the analysis of uncertainty to multiple current Tsunami disaster distribution results statistics and result.The above method can solve the problems, such as that risk evaluation result is higher in the prior art and can not provide its correspondence probability of happening, a variety of uncertainties are fused in final assessment result, increase the confidence level of result, operational efficiency is improved simultaneously, policymaker is facilitated targetedly to make prevent and reduce natural disasters deployment and physical construction planning.
Description
Technical field
The invention belongs to Tsunami disaster assessment technologies more particularly to a kind of based on the multiple current of Monte Carlo stochastic simulation
Tsunami disaster appraisal procedure.
Background technology
Since 21st century, multiple earthquake Tsunami disaster has occurred in the whole world, and direct economic loss is more than 260,000,000,000
Dollar.Wherein, the Japanese 9.0 grades of seismic sea waves of the 9.2 grades of seismic sea waves of Indonesia Sumatra in 2004 and east in 2011 all give locality
Resident living and economic development cause destructive strike.It can be seen that tsunami has become the coastal residence in the threat whole world
One of the natural calamity of people's security of the lives and property most serious.
The core for carrying out Tsunami disaster risk assessment is to determine the tsunami danger intensity of a certain assessment area, specifically
Be to determine maximum tsunami wave amplitude, tsunami that the assessment area can suffer from climb, depth of immersion etc..
Currently, generally determining tsunami danger using certainty assessment technology route in the world.Certainty tsunami risk
Appraisal procedure is mainly inferred to most destructive seismic sea wave source according to historical events, utilizes its biography of tsunami model study
It broadcasts, the process climbed and flooded, provides that impact evaluation region is worst to flood scene, and assess tsunami risk accordingly.It should
The advantages of method is simple and practicable, can be obtained conclusion, and conclusion form letter often through the simulation of one or several scenes
It is single intuitive.But there is also certain deficiencies for this method:It is difficult to provide the probability of happening of these " the worst " so-called scenes, Huo Zheping
Estimate the possibility (return period) that different degrees of Tsunami disaster occurs in region.
In addition, the compound conservative of multi-parameter makes the scene probability of happening of these " the worst " minimum, even can not possibly
's.Fig. 1 is the schematic diagram for the tsunami risk that insider assesses marginal basins area using Manila trench " the worst " scene.
Although can intuitively find out which region is influenced maximum by Tsunami disaster, which can not provide the reproduction of the scene
The probability that phase, i.e. different regions are influenced by Tsunami disaster.
The tsunami risk of assessment area different reoccurrence how is reasonably provided as a result, to ask as what current needs solved
Topic.
Invention content
For the problems of the prior art, the present invention provides a kind of multiple current tsunami based on Monte Carlo stochastic simulation
Disaster Assessment method, this method can solve risk evaluation result in the prior art it is higher and can not provide its correspond to probability of happening
The problem of.
The present invention provides a kind of multiple current Tsunami disaster appraisal procedure based on Monte Carlo stochastic simulation, including:
S1, potential tsunami source region seismicity parameters logic tree is established using a variety of focal shock parameter methods of estimation;
S2, tsunami unit source Green's function database is established using linear tsunami numerical simulation;
S3, the potential tsunami source region seismicity parameters logic tree according to foundation are given birth to using Monte Carlo stochastic simulation
At Stochastic earthquake event set;
S4, according to tsunami unit source Green's function database and the Stochastic earthquake event set, using random slippage mould
It is quasi- to generate tsunami wave amplitude collection;
S5, according to the tsunami wave amplitude collection, to multiple current Tsunami disaster distribution results statistics and to multiple current
The analysis of uncertainty of Tsunami disaster distribution results.
Optionally, the step S1 includes:
Using in a variety of focal shock parameter methods of estimation maximum likelihood fit TGR relationships and seismic moment conservation theorem establish
The logic tree of the seismicity parameters of potential tsunami source region.
Optionally, using the magnitude parameter of TGR Relation acquisition tsunami source regions, including:
The TGR relationships are:
Wherein, M is seismic moment, MtIt is the complete Critical earthquake square of earthquake catalogue, McIt is corner seismic moment, β is TGR relationships
Slope, reflect the seismicity in region;F (M) is that seismic moment is more than MtEarthquake incidence;
The approximation relation of seismic moment M and earthquake magnitude m is M=101.5m+C(2);
The unit of M is ox rice, and C is constant;
Optimal parameter M is calculated with maximum likelihood functioncAnd β, maximum likelihood function are:
Wherein, N is seismic moment MiMore than Critical earthquake square MtEarthquake number;
Using the focal shock parameter method of estimation of seismic moment conservation theorem, including:
It is expressed as according to the seismic moment of geology parameter estimation:M=χ μ WLv (4);
Wherein, χ is the earthquake coefficient of coup, and μ is rigidity modulus of shearing, and W is the underriding width of earthquake zone, and L is the length of tomography
Degree, v are the Mean Speed of plate;
On the other hand, according to TGR relationships, seismic moment M is:
Wherein,
atIt is more than M for seismic momenttSeismic events year occurrence rate, Γ is gamma function, C1For constant;
In conjunction with formula (4) and formula (5), work as Mc> > MtWhen:Obtain corner seismic moment Mc:
Obtaining corner seismic moment McAfterwards, corner earthquake magnitude is inferred to according to formula (2);
The TGR relationships and seismic moment conservation theorem obtain multiple corner earthquake magnitude values, combining global using different parameters
Subduction zone TGR relationships establish potential tsunami source region seismicity parameters logic tree.
Optionally, before executing step S2, the method further includes:
S2a, tsunami source region is divided into the tsunami unit source that several equal in magnitude and slippages are 1 meter;Unit source
Tomography geometric parameter pass through historical earthquake focal mechanism solution and related geology parameter evaluation method obtains.
Optionally, the S3 includes:
According to the tsunami source parameter in the potential tsunami source region seismicity parameters logic tree, with linear tsunami number
Tsunami wave amplitude of the value each tsunami source of simulation calculating in assessment area;
Wherein, linear tsunami numerical model is expressed as:
In formula, η is the Free Surface displacement relative to mean sea level;For latitude;ψ is latitude;R is earth radius;H
For total depth of water;P is the flux along longitude unit width;Q is the flux along latitude unit width;F is Coriolis force coefficient;G attaches most importance to
Power acceleration;
Since tsunami progression is unsatisfactory for linear relationship in offshore, the tsunami wave amplitude being calculated needs to change by Green's rule
The tsunami wave amplitude for calculating offshore, is expressed as:
Wherein AcFor the tsunami wave amplitude of offshore output point, HcFor the offshore output point depth of water, A0For the tsunami wave of offshore output point
Width, H0For the littoral output point depth of water.
Optionally, the S3 includes:
Using Monte Carlo stochastic simulation based on the logic tree established in S1, generated at random within the scope of tsunami source region more
Seismic events collection in the following certain period of time of set;
The seismic events that each independent seismic events is concentrated random distribution in tsunami source region, and each sea
Earthquake magnitude-frequency relation of seismic events in howl source region meets wherein one in TGR logic trees.
Optionally, the S4 includes:
For each seismic events in Stochastic earthquake event set, using earthquake magnitude-, rupture range formula is determining and delimits
Calculate the unit source needed for tsunami wave amplitude;Earthquake magnitude-rupture range empirical equation is expressed as:
mw=4.868+1.392log10(L) (13)
mw=4.441+0.846log10(A) (14)
M in formulawFor magnitude, L is rupture length, and A is Strain energy, and the earthquake is selected according to rupture length and area
Event calculates required tsunami unit source;
The seismic events are generated at random every using the earthquake random breakage analogue technique based on von Karman functions
Rupture slippage on a unit source;In conjunction with tsunami unit source database calculate each event assessment area maximum tsunami wave
Width is generating maximum tsunami wave amplitude event set.
Optionally, the S5 includes:
For some seismic events collection, it is assumed that there is N number of seismic sea wave source region to contribute the tsunami risk of bank point, wherein
N-th of potential complications tsunami source region influences the outcross probability Pn (H >=h) of coastal waters bank site h tsunami wave amplitudes T, then the seas site h
Whistler waves width T must outcross probability be expressed as:
Output point is released by (15) formula and reaches specified tsunami wave amplitude HhReturn period ThFor:
For a certain specific Disaster Assessment output point O, can be transferred through in specific return period T, N number of earthquake catalogue
(16) formula calculates the maximum tsunami wave amplitude h on the aspecti, then the maximum tsunami wave amplitude for output point O in return period T be:
In addition, by maximum amplitude h caused by different earthquake catalogueiClassified statistic obtain different focal shock parameters pair
As a result influence degree.
In order to achieve the above objectives, the present invention also provides a kind of multiple current Tsunami disaster based on Monte Carlo stochastic simulation
Apparatus for evaluating, including memory, processor, bus and storage are on a memory and the computer journey that can run on a processor
Sequence, the processor are realized when executing described program such as the step of above method any one.
In order to achieve the above objectives, a kind of computer storage media, is stored thereon with computer program, and described program is handled
It is realized such as the step of the above method any one when device executes.
The device have the advantages that as follows:
The first, due to existing earthquake catalogue curtailment (longest was only more than 100 years), sample number is less, simple
Corner earthquake magnitude (the m in Fig. 4 therein can be caused using historical earthquake data fitting TGR relationshipsc, main function is to TGR relationships
High earthquake magnitude part is modified, and more than the return period of the seismic events of this earthquake magnitude, exponentially form increases rapidly in TGR relationships) partially
It is low, and then underestimate the risk of tsunami.In the present invention the potential tsunami source region TGR logics of relations are established in conjunction with seismic moment superposition principle
Tree, the law can pass through the Seismic annual occurrence rate of a period of time and regional structure parameter (underriding rate, the earthquake coupling of subduction zone
Close the factor, modulus of shearing etc.) calculate the corner earthquake magnitude in TGR relationships, corner earthquake magnitude is thus given in logic tree
A variety of possibility so that assessment result is more accurate.
The second, the certainty that any perfect evaluation system should all have uncertainty analysis system, however apply at present
Tsunami methods of risk assessment is can not provide result probabilistic.Monte carlo algorithm (Meng Teka is utilized in the present invention
Lip river stochastic simulation), simulate the seismic events of magnanimity.The position of each seismic events, earthquake magnitude or even rupture process are totally full
All it is random on the basis of the certain rule of foot.This just provides enough samples for analysis of uncertainty, also can will not really
It is qualitative to be fused in assessment result.
Third, during tsunami numerical simulation, tsunami initial displacement field model is typically the bullet that is proposed by Okada
Property dislocation FAULT MODEL calculate.The slippage of the model interrupting layer determines the sea level relief volume of initial fields, logarithm mould
Quasi- result is affected.However, in current tsunami risk assessment application, people are usual in the slippage of computed tomography
It assumes that it is uniformly distributed, and is calculated according to earthquake magnitude empirical equation.With the development of earthquake finite fault inversion technique, people
Find that the slippage of actual seismic rupture is unevenly distributed, and the influence to region tsunami numerical simulation result is very big.
Present invention application tsunami unit source Green's function database makes to calculate maximum tsunami wave amplitude, this method caused by each seismic events
Tsunami wave amplitude caused by the seismic events of magnanimity random breakage process must be simulated to be possibly realized, can both increase the accurate of result
Property, while uncertainty of the random breakage process to result can also be assessed.
It is higher and can not provide result pair to solve tsunami assessment result in the prior art for method using the present invention as a result,
A variety of uncertainties are fused in final assessment result by the problem of answering probability of happening, increase the confidence level of result, simultaneously
Operational efficiency is improved, policymaker is facilitated targetedly to make prevent and reduce natural disasters deployment and physical construction planning.
Description of the drawings
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technology description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Some embodiments of invention without having to pay creative labor, may be used also for those of ordinary skill in the art
To obtain other attached drawings according to these attached drawings.
Fig. 1 is the schematic diagram for the tsunami risk that Manila trench " the worst " scene assesses marginal basins area;
Fig. 2 a and Fig. 2 b are respectively method flow schematic diagram provided in an embodiment of the present invention;
Fig. 3 is distributed for Manila trench historical earthquake;
Fig. 4 is Manila trench TGR relation schematic diagrams;
Fig. 5 is Manila trench TGR parameter logistics trees;
Fig. 6 is that Manila trench tsunami unit source divides schematic diagram;
Fig. 7 be the marginal basins regional return period be 200 years maximum tsunami wave amplitude distribution schematic diagram;
Fig. 8 be the marginal basins regional return period be 500 years maximum tsunami wave amplitude distribution schematic diagram;
Fig. 9 be the marginal basins regional return period be 2000 maximum tsunami wave amplitude distribution schematic diagram;
Figure 10 (a) to Figure 10 (d) is respectively that the schematic diagram of marginal basins urban tsunami wave height curve (is followed successively by height
Hero, Ku Limaao, Hong Kong and Wenchang).
Specific implementation mode
In order to preferably explain the present invention, in order to understand, below in conjunction with the accompanying drawings, by specific implementation mode, this is sent out
It is bright to be used to assess marginal basins tsunami risk.
In the following description, by multiple and different aspects of the description present invention, however, for common skill in the art
For art personnel, the present invention can be implemented just with some or all structures or flow of the present invention.In order to explain
Definition for, specific number, configuration and sequence are elaborated, however, it will be apparent that these specific details the case where
Under can also implement the present invention.It in other cases, will no longer for some well-known features in order not to obscure the present invention
It is described in detail.
As shown in Figure 2 a and 2 b, Fig. 2 a and Fig. 2 b respectively illustrate one embodiment of the invention offer based on Monte Carlo
The flow chart of the multiple current Tsunami disaster appraisal procedure of stochastic simulation, the method for the present embodiment include the following steps:
S1, potential tsunami source region seismicity parameters logic tree is established using a variety of focal shock parameter methods of estimation.
For example, can be Manila sea according to geological structure theory, the main subduction zone of marginal basins in the present embodiment
Ditch, while being also the main potential tsunami source in the region.Historical earthquake in regional extent is as shown in figure 3, can in the present embodiment
Using the TGR relationships in the maximum likelihood fit region, the TGR relationships are:
Wherein, M is seismic moment, MtIt is the complete Critical earthquake square of earthquake catalogue, McIt is corner seismic moment, β is TGR relationships
Slope, reflect the seismicity in region;
The approximation relation of seismic moment M and earthquake magnitude m is M=101.5m+C(2);
The unit of M is ox rice, and C is constant.
Optimal parameter M is calculated with maximum likelihood methodcAnd β, maximum likelihood function are:
The TGR parameters in the region, wherein β=0.45, corner earthquake magnitude m are found out using maximum likelihood methodcValue be 7.5, it is bent
Line is as shown in Figure 4.Corner earthquake magnitude in the relationship is significantly lower than the average value of global subduction zone, mainly due to earthquake catalogue mistake
Caused by short, it is therefore desirable to be modified by other methods.
Using the focal shock parameter method of estimation of seismic moment conservation theorem, it may include:
It is represented by according to the seismic moment of geology parameter estimation:M=χ μ WLv (4);
Wherein χ is the earthquake coefficient of coup, and according to previous studies, the geology coefficient of coup in the region is between 0.1-0.3;μ
For rigidity modulus of shearing, value 3.0-4.9Gpa;W is the underriding width of earthquake zone, value 98km;L is the length of tomography, is taken
Value is 1100km;V is the Mean Speed of plate, value 90mm/yr;
On the other hand, according to TGR relationships, seismic moment M can be evaluated whether for:
Wherein,
atIt is more than M for seismic momenttSeismic events year occurrence rate, Γ is gamma function, C1For constant;
In conjunction with formula (4) and formula (5), work as Mc> > MtWhen:
Corner seismic moment M can be obtainedc:
Obtaining corner seismic moment McAfterwards, corner earthquake magnitude is inferred to according to formula (2).Bringing different parameters into can obtain
Modified four corner earthquake magnitudes value 7.9,8.1,8.2 and 8.4, β value can select the TGR fitting results 0.45 in Manila region and complete
The average result 0.66 of ball subduction zone.Thus Manila trench TGR logics of relations tree (as shown in Figure 5) is established.
S2, tsunami unit source Green's function database is established using linear tsunami numerical simulation.
For example, before executing S2, it is 1 meter that tsunami source region can be divided into several equal in magnitude and slippages
Tsunami unit source;The tomography geometric parameter in unit source passes through historical earthquake focal mechanism solution and related geology parameter evaluation method
It obtains.
For example, dividing Manila region in conjunction with existing geologic parameter near Manila trench and history focal mechanism solution
Tsunami unit source (as shown in Figure 6), each a length of 100km in unit source, width 50km, slippage be 1 meter.The tomography in unit source
Geometric parameter can be obtained by historical earthquake focal mechanism solution and relevant geologic parameter evaluation method.
Based on tsunami source parameter (being obtained in step S1) above, the tsunami numerical model in marginal basins region is established,
Terrain resolution is 2 points, chooses 355 output points in bank, is being commented with each tsunami source of linear tsunami numerical simulation calculation
The tsunami wave amplitude in region is estimated, wherein linear tsunami numerical model can be expressed as:
In formula, η is the Free Surface displacement relative to mean sea level;For latitude;ψ is latitude;R is earth radius;H
For total depth of water;P is the flux along longitude unit width;Q is the flux along latitude unit width;F is Coriolis force coefficient;G attaches most importance to
Power acceleration.
Since tsunami progression is unsatisfactory for linear relationship in offshore, the tsunami wave amplitude being calculated needs to change by Green's rule
The tsunami wave amplitude for calculating offshore, is represented by:
Wherein AcFor the tsunami wave amplitude of offshore output point, HcFor the offshore output point depth of water, A0For the tsunami wave of offshore output point
Width, H0For the littoral output point depth of water.
S3, the potential tsunami source region seismicity parameters logic tree according to foundation are given birth to using Monte Carlo stochastic simulation
At Stochastic earthquake event set.
Using Monte Carlo stochastic simulation (Monte carlo algorithm) based on the logic tree established in S1, in tsunami source region model
Enclose the seismic events collection in the following certain period of time of interior random generation mostly set.Symbiosis is at earthquake catalogue 100, each earthquake catalogue
Time is 100,000 year, amounts to seismic events 166096, largest magnitude 9.05.
S4, according to tsunami unit source Green's function database and the Stochastic earthquake event set, using random slippage mould
It is quasi- to generate tsunami wave amplitude collection.
For each seismic events in generated Stochastic earthquake event set in step S3, using earthquake magnitude-rupture model
It encloses formula determination and delimit the unit source calculated needed for tsunami wave amplitude.Earthquake magnitude-rupture range empirical equation is represented by:
mw=4.868+1.392log10(L) (13)
mw=4.441+0.846log10(A) (14)
M in formulawFor magnitude, L is rupture length, and A is Strain energy, and choosing can be determined according to rupture length and area
Determine seismic events and calculates required tsunami unit source.
Above-mentioned seismic events are given birth at random each using the earthquake random breakage analogue technique based on von Karman functions
Rupture slippage on unit source.In conjunction with tsunami unit source database calculate each event assessment area maximum tsunami wave amplitude
To generate maximum tsunami wave amplitude event set.
S5, according to the tsunami wave amplitude collection, to multiple current Tsunami disaster distribution results statistics and to multiple current
The analysis of uncertainty of Tsunami disaster distribution results.
For example, for some seismic events collection, it is assumed that there is N number of seismic sea wave source region to have tribute to the tsunami risk of bank point
It offers, wherein n-th of potential complications tsunami source region influences the outcross probability Pn (H >=h) of coastal waters bank site h tsunami wave amplitudes T, then
Site h tsunami wave amplitudes T must outcross probability be expressed as:
Output point can be released by (15) formula and reach specified tsunami wave amplitude HhReturn period ThFor:
For a certain specific Disaster Assessment output point O, can be transferred through in specific return period T, N number of earthquake catalogue
(16) formula calculates the maximum tsunami wave amplitude h on the aspecti, then the maximum tsunami wave amplitude for output point O in return period T be:
Maximum tsunami wave amplitude distribution in 200 years, 500 years and 2000 in marginal basins region is calculated according to above-mentioned formula
As shown in Fig. 7, Fig. 8 or Fig. 9, statistics obtains Kaohsiung, Ku Limaao, the tsunami maximum amplitude curve in Hong Kong and the city of Wenchang four
As shown in Figure 10, Figure 10 (a) to Figure 10 (d) is respectively the schematic diagram of marginal basins urban tsunami wave height curve, such as according to
Secondary is Kaohsiung, Ku Limaao, the tsunami wave height curve in Hong Kong and Wenchang city.
In the present embodiment the potential tsunami source region TGR logics of relations are established in conjunction with seismic moment superposition principle in the present invention
Tree, the law can pass through the Seismic annual occurrence rate of a period of time and regional structure parameter (underriding rate, the earthquake coupling of subduction zone
Close the factor, modulus of shearing etc.) calculate the corner earthquake magnitude in TGR relationships, corner earthquake magnitude is thus given in logic tree
A variety of possibility so that assessment result is more accurate.Using Monte carlo algorithm, the seismic events of magnanimity are simulated.Each earthquake
The position of event, earthquake magnitude or even rupture process are all random on the basis of totally meeting certain rule.This is just not true
Qualitative analysis provides enough samples, and also uncertainty can be fused in assessment result.Using tsunami unit source Green's letter
Number library makes the ground of simulation magnanimity random breakage process to calculate maximum tsunami wave amplitude, this method caused by each seismic events
Tsunami wave amplitude caused by shake event is possibly realized, and can both increase the accuracy of result, while can also assess random breakage
Uncertainty of the process to result.
Method using the present invention gives marginal basins area multiple current tsunami risk and urban as a result,
Tsunami amplitude curve.
Another aspect according to the ... of the embodiment of the present invention, the present invention also provides a kind of based on the more of Monte Carlo stochastic simulation
Return period Tsunami disaster apparatus for evaluating, the device include:Memory, processor, bus and storage are on a memory and can be by
The computer program that processor executes, the processor execute the method for realizing above-mentioned Fig. 2 a and Fig. 2 b when described program.For example,
Potential tsunami source region seismicity parameters logic tree is established using a variety of focal shock parameter methods of estimation;For example, answering
With in a variety of focal shock parameter methods of estimation maximum likelihood fit TGR relationships and seismic moment conservation theorem establish potential tsunami source
The logic tree of the seismicity parameters in area.
Tsunami unit source Green's function database is established using linear tsunami numerical simulation;
According to the potential tsunami source region seismicity parameters logic tree of foundation, using Monte Carlo stochastic simulation generate with
Machine seismic events collection;
According to tsunami unit source Green's function database and the Stochastic earthquake event set, simulates and give birth to using random slippage
At tsunami wave amplitude collection;
According to the tsunami wave amplitude collection, to multiple current Tsunami disaster distribution results statistics and to multiple current tsunami
The analysis of uncertainty of disaster distribution results.
In addition, the embodiment of the present invention also provides a kind of computer storage media, it is stored thereon with computer program, the journey
It is realized such as the step of the above method any one when sequence is executed by processor.
It should be clear that the invention is not limited in specific configuration described above and shown in figure and processing.
For brevity, it is omitted here the detailed description to known method.In the above-described embodiments, several tools have been described and illustrated
The step of body, is as example.But procedure of the invention is not limited to described and illustrated specific steps, this field
Technical staff can make various changes, modification and addition after the spirit for understanding the present invention, or suitable between changing the step
Sequence.
It should also be noted that, the exemplary embodiment referred in the present invention, is retouched based on a series of step or device
State certain methods or system.But the present invention is not limited to the sequence of above-mentioned steps, that is to say, that can be according in embodiment
The sequence referred to executes step, may also be distinct from that the sequence in embodiment or several steps are performed simultaneously.
Finally it should be noted that:Above-described embodiments are merely to illustrate the technical scheme, rather than to it
Limitation;Although the present invention is described in detail referring to the foregoing embodiments, it will be understood by those of ordinary skill in the art that:
It can still modify to the technical solution recorded in previous embodiment, or to which part or all technical features into
Row equivalent replacement;And these modifications or substitutions, it does not separate the essence of the corresponding technical solution various embodiments of the present invention technical side
The range of case.
Claims (10)
1. a kind of multiple current Tsunami disaster appraisal procedure based on Monte Carlo stochastic simulation, which is characterized in that including:
S1, potential tsunami source region seismicity parameters logic tree is established using a variety of focal shock parameter methods of estimation;
S2, tsunami unit source Green's function database is established using linear tsunami numerical simulation;
S3, the potential tsunami source region seismicity parameters logic tree according to foundation, using Monte Carlo stochastic simulation generate with
Machine seismic events collection;
S4, according to tsunami unit source Green's function database and the Stochastic earthquake event set, simulate and give birth to using random slippage
At tsunami wave amplitude collection;
S5, according to the tsunami wave amplitude collection, to multiple current Tsunami disaster distribution results statistics and to multiple current tsunami
The analysis of uncertainty of disaster distribution results.
2. according to the method described in claim 1, it is characterized in that, the step S1 includes:
Using in a variety of focal shock parameter methods of estimation maximum likelihood fit TGR relationships and seismic moment conservation theorem establish it is potential
The logic tree of the seismicity parameters of tsunami source region.
3. according to the method described in claim 2, it is characterized in that,
Using the magnitude parameter of TGR Relation acquisition tsunami source regions, including:
The TGR relationships are:
Wherein, M is seismic moment, MtIt is the complete Critical earthquake square of earthquake catalogue, McIt is corner seismic moment, β is the oblique of TGR relationships
Rate reflects the seismicity in region;F (M) is that seismic moment is more than MtEarthquake incidence;
The approximation relation of seismic moment M and earthquake magnitude m is M=101.5m+C(2);
The unit of M is ox rice, and C is constant;
Optimal parameter M is calculated with maximum likelihood functioncAnd β, maximum likelihood function are:
Wherein, N is seismic moment MiMore than Critical earthquake square MtEarthquake number;
Using the focal shock parameter method of estimation of seismic moment conservation theorem, including:
It is expressed as according to the seismic moment of geology parameter estimation:M=χ μ WLv (4);
Wherein, χ is the earthquake coefficient of coup, and μ is rigidity modulus of shearing, and W is the underriding width of earthquake zone, and L is the length of tomography, v
For the Mean Speed of plate;
On the other hand, according to TGR relationships, seismic moment M is:
Wherein,
atIt is more than M for seismic momenttSeismic events year occurrence rate, Γ is gamma function, C1For constant;
In conjunction with formula (4) and formula (5), work as Mc> > MtWhen:Obtain corner seismic moment Mc:
Obtaining corner seismic moment McAfterwards, corner earthquake magnitude is inferred to according to formula (2);
The TGR relationships and seismic moment conservation theorem obtain multiple corner earthquake magnitude values using different parameters, and combining global dives
Band TGR relationships, establish potential tsunami source region seismicity parameters logic tree.
4. according to the method described in claim 3, it is characterized in that, before executing step S2, the method further includes:
S2a, tsunami source region is divided into the tsunami unit source that several equal in magnitude and slippages are 1 meter;Break in unit source
Layer geometric parameter is obtained by historical earthquake focal mechanism solution and related geology parameter evaluation method.
5. according to the method described in claim 3, it is characterized in that, the S3 includes:
According to the tsunami source parameter in the potential tsunami source region seismicity parameters logic tree, with linear tsunami Numerical-Mode
Tsunami wave amplitude of the quasi- each tsunami source of calculating in assessment area;
Wherein, linear tsunami numerical model is expressed as:
In formula, η is the Free Surface displacement relative to mean sea level;For latitude;ψ is latitude;R is earth radius;H is total
The depth of water;P is the flux along longitude unit width;Q is the flux along latitude unit width;F is Coriolis force coefficient;G adds for gravity
Speed;
Since tsunami progression is unsatisfactory for linear relationship in offshore, the tsunami wave amplitude being calculated needs to be converted to by Green's rule
The tsunami wave amplitude of offshore, is expressed as:
Wherein AcFor the tsunami wave amplitude of offshore output point, HcFor the offshore output point depth of water, A0For the tsunami wave amplitude of offshore output point, H0
For the littoral output point depth of water.
6. according to the method described in claim 5, it is characterized in that, the S3 includes:
Using Monte Carlo stochastic simulation based on the logic tree established in S1, more sets are generated at random not within the scope of tsunami source region
Carry out the seismic events collection in certain period of time;
The seismic events that each independent seismic events is concentrated random distribution in tsunami source region, and each tsunami source
Earthquake magnitude-frequency relation of seismic events in area meets wherein one in TGR logic trees.
7. according to the method described in claim 6, it is characterized in that, the S4 includes:
For each seismic events in Stochastic earthquake event set, using earthquake magnitude-, rupture range formula is determining and delimits calculating
Unit source needed for tsunami wave amplitude;Earthquake magnitude-rupture range empirical equation is expressed as:
mw=4.868+1.392log10(L) (13);
mw=4.441+0.846log10(A) (14);
M in formulawFor magnitude, L is rupture length, and A is Strain energy, and the seismic events are selected according to rupture length and area
Tsunami unit source needed for calculating;
The seismic events are generated at random in each list using the earthquake random breakage analogue technique based on von Karman functions
Rupture slippage on the source of position;Each event is calculated in conjunction with tsunami unit source database to use in the maximum tsunami wave amplitude of assessment area
To generate maximum tsunami wave amplitude event set.
8. the method according to the description of claim 7 is characterized in that the S5 includes:
For some seismic events collection, it is assumed that there is N number of seismic sea wave source region to contribute the tsunami risk of bank point, wherein n-th
A potential complications tsunami source region influences the outcross probability Pn (H >=h) of coastal waters bank site h tsunami wave amplitudes T, then site h tsunamis
Wave amplitude T must outcross probability be expressed as:
Output point is released by (15) formula and reaches specified tsunami wave amplitude HhReturn period ThFor:
For a certain specific Disaster Assessment output point O, (16) formula is can be transferred through in specific return period T, N number of earthquake catalogue
Calculate the maximum tsunami wave amplitude h on the aspecti, then the maximum tsunami wave amplitude for output point O in return period T be:
In addition, by maximum amplitude h caused by different earthquake catalogueiClassified statistic obtain different focal shock parameters to result
Influence degree.
9. a kind of multiple current Tsunami disaster apparatus for evaluating based on Monte Carlo stochastic simulation, which is characterized in that including storage
Device, processor, bus and storage on a memory and the computer program that can run on a processor, the processor execution
It is realized such as the step of claim 1-8 any one when described program.
10. a kind of computer storage media, is stored thereon with computer program, it is characterised in that:Described program is held by processor
It is realized such as the step of claim 1-8 any one when row.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810124747.6A CN108492236B (en) | 2018-02-07 | 2018-02-07 | Multiple current Tsunami disaster appraisal procedure based on Monte Carlo stochastic simulation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810124747.6A CN108492236B (en) | 2018-02-07 | 2018-02-07 | Multiple current Tsunami disaster appraisal procedure based on Monte Carlo stochastic simulation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108492236A true CN108492236A (en) | 2018-09-04 |
CN108492236B CN108492236B (en) | 2019-06-14 |
Family
ID=63344743
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810124747.6A Expired - Fee Related CN108492236B (en) | 2018-02-07 | 2018-02-07 | Multiple current Tsunami disaster appraisal procedure based on Monte Carlo stochastic simulation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108492236B (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR102061500B1 (en) | 2019-10-30 | 2020-01-02 | (주)해양정보기술 | System predicting tsunami based on artificial neural network |
CN110727026A (en) * | 2019-10-10 | 2020-01-24 | 太原理工大学 | Inelastic displacement ratio spectrum model based on earthquake-tsunami continuous action |
CN110751399A (en) * | 2019-10-22 | 2020-02-04 | 中国地震局地球物理研究所 | Earthquake danger parallel analysis method, device and equipment |
CN111382908A (en) * | 2020-03-11 | 2020-07-07 | 中国地震局地球物理研究所 | Earthquake random event set simulation method considering large earthquake time correlation |
CN113268852A (en) * | 2021-04-14 | 2021-08-17 | 西南交通大学 | Monte Carlo simulation-based earthquake landslide probability risk analysis method |
CN114091756A (en) * | 2021-11-23 | 2022-02-25 | 国家海洋环境预报中心 | Township tsunami risk assessment method based on Thiessen polygon |
CN114114393A (en) * | 2021-10-11 | 2022-03-01 | 中国地震局地球物理研究所 | Method for evaluating historical seismic intensity characteristics |
CN115495941A (en) * | 2022-11-18 | 2022-12-20 | 海南浙江大学研究院 | Response surface-based rapid forecasting method for submarine landslide tsunami |
CN115542397A (en) * | 2022-10-08 | 2022-12-30 | 国家海洋环境预报中心 | Tsunami danger analysis method based on SVM |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104615847A (en) * | 2014-12-24 | 2015-05-13 | 中交天津港湾工程研究院有限公司 | Tsunami dangerousness forecasting method for the South China Sea based on probabilistic method |
CN106324674A (en) * | 2016-08-23 | 2017-01-11 | 中国石油大学(华东) | Shale gas TOC pre-stack seismic inversion prediction method |
-
2018
- 2018-02-07 CN CN201810124747.6A patent/CN108492236B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104615847A (en) * | 2014-12-24 | 2015-05-13 | 中交天津港湾工程研究院有限公司 | Tsunami dangerousness forecasting method for the South China Sea based on probabilistic method |
CN106324674A (en) * | 2016-08-23 | 2017-01-11 | 中国石油大学(华东) | Shale gas TOC pre-stack seismic inversion prediction method |
Non-Patent Citations (4)
Title |
---|
YAN Y.KAGAN: "Seismic moment distribution revisited: II. Moment conservation principle", 《GEOPHYS.J.INT》 * |
YUFANGRONG等: "俯冲带地震的震级极限", 《世界地震译丛》 * |
杨智博: "中国地震海啸危险性分析", 《中国优秀硕士学位论文全文数据库基础科学辑》 * |
王君成等: "基于海啸浮标数据反演的海啸数值预报方法的初步研究", 《海洋预报》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110727026A (en) * | 2019-10-10 | 2020-01-24 | 太原理工大学 | Inelastic displacement ratio spectrum model based on earthquake-tsunami continuous action |
CN110727026B (en) * | 2019-10-10 | 2021-05-28 | 太原理工大学 | Inelastic displacement ratio spectrum model based on earthquake-tsunami continuous action |
CN110751399B (en) * | 2019-10-22 | 2023-02-24 | 中国地震局地球物理研究所 | Earthquake danger parallel analysis method, device and equipment |
CN110751399A (en) * | 2019-10-22 | 2020-02-04 | 中国地震局地球物理研究所 | Earthquake danger parallel analysis method, device and equipment |
KR102061500B1 (en) | 2019-10-30 | 2020-01-02 | (주)해양정보기술 | System predicting tsunami based on artificial neural network |
CN111382908A (en) * | 2020-03-11 | 2020-07-07 | 中国地震局地球物理研究所 | Earthquake random event set simulation method considering large earthquake time correlation |
CN111382908B (en) * | 2020-03-11 | 2023-08-22 | 中国地震局地球物理研究所 | Earthquake random event set simulation method considering correlation of large earthquake time |
CN113268852B (en) * | 2021-04-14 | 2022-02-22 | 西南交通大学 | Monte Carlo simulation-based earthquake landslide probability risk analysis method |
CN113268852A (en) * | 2021-04-14 | 2021-08-17 | 西南交通大学 | Monte Carlo simulation-based earthquake landslide probability risk analysis method |
CN114114393A (en) * | 2021-10-11 | 2022-03-01 | 中国地震局地球物理研究所 | Method for evaluating historical seismic intensity characteristics |
CN114091756B (en) * | 2021-11-23 | 2022-07-05 | 国家海洋环境预报中心 | Township tsunami risk assessment method based on Thiessen polygon |
CN114091756A (en) * | 2021-11-23 | 2022-02-25 | 国家海洋环境预报中心 | Township tsunami risk assessment method based on Thiessen polygon |
CN115542397A (en) * | 2022-10-08 | 2022-12-30 | 国家海洋环境预报中心 | Tsunami danger analysis method based on SVM |
CN115495941A (en) * | 2022-11-18 | 2022-12-20 | 海南浙江大学研究院 | Response surface-based rapid forecasting method for submarine landslide tsunami |
CN115495941B (en) * | 2022-11-18 | 2023-03-28 | 海南浙江大学研究院 | Response surface-based rapid forecasting method for submarine landslide tsunami |
Also Published As
Publication number | Publication date |
---|---|
CN108492236B (en) | 2019-06-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108492236B (en) | Multiple current Tsunami disaster appraisal procedure based on Monte Carlo stochastic simulation | |
Takagi et al. | Track analysis, simulation, and field survey of the 2013 Typhoon Haiyan storm surge | |
Charvet et al. | A multivariate generalized linear tsunami fragility model for Kesennuma City based on maximum flow depths, velocities and debris impact, with evaluation of predictive accuracy | |
Glimsdal et al. | A new approximate method for quantifying tsunami maximum inundation height probability | |
Park et al. | Probabilistic seismic and tsunami hazard analysis conditioned on a megathrust rupture of the Cascadia subduction zone | |
Chen et al. | Extratropical storm inundation testbed: Intermodel comparisons in Scituate, Massachusetts | |
Tehranirad et al. | Far-field tsunami impact in the North Atlantic basin from large scale flank collapses of the Cumbre Vieja Volcano, La Palma | |
CN104615847A (en) | Tsunami dangerousness forecasting method for the South China Sea based on probabilistic method | |
Park et al. | Probabilistic Tsunami Hazard Assessment (PTHA) for resilience assessment of a coastal community | |
Knighton et al. | A proposed probabilistic seismic tsunami hazard analysis methodology | |
Sampurno et al. | Integrated hydrodynamic and machine learning models for compound flooding prediction in a data-scarce estuarine delta | |
CN115439029A (en) | Tsunami disaster key defense area determination method and system | |
El-Hussain et al. | Probabilistic and deterministic estimates of near-field tsunami hazards in northeast Oman | |
Williamson et al. | A source clustering approach for efficient inundation modeling and regional scale probabilistic tsunami hazard assessment | |
Annunziato | The Tsunami Assessment Modelling System by the Joint Research Centre | |
Liang et al. | Simulation of tsunami propagation using adaptive cartesian grids | |
Fukui et al. | Subgrid-scale modeling of tsunami inundation in coastal urban areas | |
Gonzalez et al. | Quantification of uncertainty in probabilistic storm surge models: literature review | |
Wang et al. | Deep Reinforcement Learning-Based Rescue Resource Distribution Scheduling of Storm Surge Inundation Emergency Logistics | |
Le Gal et al. | Reef influence quantification in light of the 1771 Meiwa tsunami | |
Hammond et al. | The development of a flood damage assessment tool for urban areas | |
de Lange et al. | Sea-Level Change | |
Roberts et al. | Simulation of tsunami and flash floods | |
El-Barmelgy et al. | Tsunami risk assessment a mitigation planning tool | |
Josiah et al. | Simulation of tsunami inundation in east coast of Sri Lanka |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190614 |
|
CF01 | Termination of patent right due to non-payment of annual fee |