CN109446587A - A kind of extensive Landslide Hazards Fast numerical simulation method and system - Google Patents
A kind of extensive Landslide Hazards Fast numerical simulation method and system Download PDFInfo
- Publication number
- CN109446587A CN109446587A CN201811166665.4A CN201811166665A CN109446587A CN 109446587 A CN109446587 A CN 109446587A CN 201811166665 A CN201811166665 A CN 201811166665A CN 109446587 A CN109446587 A CN 109446587A
- Authority
- CN
- China
- Prior art keywords
- grid
- zoning
- calculation
- roughening
- node
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 59
- 238000004088 simulation Methods 0.000 title claims abstract description 36
- 239000000463 material Substances 0.000 claims abstract description 27
- 238000004364 calculation method Methods 0.000 claims abstract description 24
- 230000008569 process Effects 0.000 claims abstract description 18
- 238000007670 refining Methods 0.000 claims abstract 2
- 238000013316 zoning Methods 0.000 claims description 37
- 238000007788 roughening Methods 0.000 claims description 29
- 230000008859 change Effects 0.000 claims description 5
- 230000006870 function Effects 0.000 claims description 5
- 238000004513 sizing Methods 0.000 claims description 4
- 238000013461 design Methods 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 2
- 230000033001 locomotion Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000005713 exacerbation Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000000116 mitigating effect Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000004540 process dynamic Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种大规模滑坡地质灾害快速数值仿真模拟方法及系统，包括如下步骤：获取滑坡地质灾害地形数据、物源区范围、高程以及物理力学参数；划分大规模滑坡地质灾害计算区域；对粗细化网格条件进行判断，并对计算区域网格粗化、网格细化和网格结构化；利用数值模拟软件对地表灾害动力过程进行数值模拟计算；输出n个计算区域的计算结果，并将n个计算区域的计算结果合并成一个计算区域输出。该方法通过对计算区域进行区域划分、网格粗细化和结构化网格，将有物源、需要计算的区域进行精细化计算，没有物源、不需要计算的区域网格粗化，不仅实现对大规模的计算网格进行计算，而且提高了计算速率，实现了对大规模滑坡地质灾害的快速模拟计算。
The invention discloses a rapid numerical simulation method and system for large-scale landslide geological disasters, comprising the following steps: acquiring landslide geological disaster terrain data, source area range, elevation and physical and mechanical parameters; dividing large-scale landslide geological disaster calculation areas; Judging the conditions of coarse and fine grids, and coarsening, refining and structuring the grids in the calculation area; using numerical simulation software to numerically simulate the dynamic process of surface disasters; outputting the calculation results of n calculation areas , and combine the calculation results of n calculation areas into one calculation area output. By dividing the calculation area, meshing coarsening and structuring the grid, the method performs refined calculation on the area with material source and needs calculation. The large-scale calculation grid is calculated, and the calculation speed is improved, and the rapid simulation calculation of large-scale landslide geological disasters is realized.
Description
Technical field
The present invention relates to geological disaster dynamic process numerical simulation calculating fields, more particularly to a kind of extensive landslide geology
Disaster Fast numerical simulation method and system.
Background technique
Geological disaster refers in the landslide evolution process of the earth, the disaster geology thing formed by various geologic processes
Part.China is disaster one of the countries with the most serious ... in the world, and since the nineties, the annual disaster in China causes thousands of people
Death causes direct economic loss to be up to tens billion of members to hundreds billion of members, the 3%-6% of suitable gross national product (GNP).With
Increasing for calamity kind, the increase of calamity rate, the exacerbation of calamity degree, harm seriousness as environmental problem, become restrict China
One of an important factor for sustainable development.And the geological disaster of General Promotion China is wanted to prevent and reduce natural disasters disaster relief ability, it not only needs to know
How road disaster is done after occurring, with greater need for sentencing knowledge before disaster occurs and predict the motion mode of disaster, track and endanger model
It encloses.The kinetics mechanism and kinetic characteristic and its complexity of geological disaster, numerical simulation and qualitative assessment based on dynamic process exist
Still lack good solution in the world, is current international frontier nature problem in science and the key areas urgently broken through.
Continuous media is based primarily upon, based on discrete and between the two three to mountain region disaster dynamic process numerical simulation at this stage
Kind method.Various method for numerical simulation have its unique function and solve the ability of problems, and take in integrated application
Obtained a large amount of fruitful research achievements.It is all earth's surface object under the force of gravity in view of geological disasters such as landslide, mud-rock flows
Matter movement, they all have a common trait, i.e., the characteristic dimension of its direction of motion is much larger than the feature in its vertical direction
Highly.But how to accomplish to imitate the dynamic process of extensive Landslide Hazards (especially disaster chain) and the quick of kinetic characteristic
True simulation is domestic or even problem to be solved in the world.Its Major Difficulties is embodied in the following aspects: 1. regions are big, net
Lattice quantity is more, and calculation scale is larger；2. computation rate is slow, it is long to calculate the time.And after disaster generation, the assessment of the condition of a disaster is carved
Do not allow to delay, how quickly to accomplish that even to large scale disasters preceding quick recognition and prediction calamity occur for the assessment to extensive landslide disaster
Harmful motion mode, track and damaging range provides timely data support for the correlation department of preventing and reducing natural disasters, to reduce very
It is still current urgent need to resolve to disaster is avoided the problem that loss caused by the mankind and Human dried bloodstains.
It is established with structured grid to the fast of extensive Landslide Hazards dynamic process currently, there is no and drawn based on region again
The report of fast Numerical Simulation method and system.
Summary of the invention
In order to overcome above-mentioned deficiency, the purpose of the present invention is by extensive grid carry out region division, to material resource ratio
The different region of example uses different grades of thick tessellated mesh, pays close attention to region refinement, and other area grid roughening provide one
Kind is drawn again based on region and the extensive Landslide Hazards Fast numerical simulation method and system of structured grid, overcomes
It is big to extensive Landslide Hazards numerical simulation time domain in the prior art, number of grid is more, computation rate is slow actually asks
Topic.
To achieve the goals above, the present invention adopts the following technical scheme:
A kind of extensive Landslide Hazards Fast numerical simulation method, includes the following steps:
Step 1: obtaining Landslide Hazards terrain data, source area range, elevation and physical and mechanical parameter；
Step 2: dividing extensive Landslide Hazards zoning；
Step 3: thick tessellated mesh condition being judged, and to the roughening of zoning grid, mesh refinement and grid knot
Structure；
Step 4: numerical simulation calculation being carried out to earth's surface disaster dynamic process using numerical simulation software；
Step 5: the calculated result of n zoning of output, and the calculated result of n zoning is merged into a meter
Calculate region output.
Preferably, the method for Landslide Hazards terrain data is obtained in the step 1 specifically: acquire using unmanned plane
Terrain data, and dem data is generated with 3 d modeling software.
It is furthermore preferred that taking fixed point to fly according to Landslide Hazards range when the acquisition terrain data using unmanned plane
Row, S-ing or the flight of the back-shaped line of flight, the 3 d modeling software are generated using smart3D, Pix4D or Photoscan
Dem data.
Preferably, the method that Landslide Hazards source area range and elevation acquisition are obtained in the step 1 specifically: benefit
The unmanned plane or remote sensing altitude data that front and back occurs with landslide disaster, are subtracted each other by grid, obtain source area range and elevation.
Preferably, the method for Landslide Hazards physical and mechanical parameter is obtained in the step 1 specifically: pass through correlative
Reason experiment of machanics obtains the parameters such as disaster volume density, coefficient of friction by way of tentative calculation.
Preferably, the method for extensive Landslide Hazards zoning is divided in the step 2 specifically: using simply
Be divided by and complementation, obtained extensive Landslide Hazards dem data is divided into n*n region, then will before great Qu
Calculating parameter involved in domain is assigned to n*n region, and the equivalence after guaranteeing region division, wherein n is according to different calculating
Precision and computer capacity carry out value.
Preferably, the method judged in the step 3 thick tessellated mesh condition specifically includes: according to calculating pair
As carrying out grid thickness to different zonings, using material resource as Rule of judgment, i.e., no material resource does not need the area calculated
Grid roughening in domain has the calculative area grid refinement of material resource, the grade that simultaneous selection is roughened, refines, such as zoning
Material resource amount is more than 50% grid constantly to refine to two-stage, zoning have material resource but material resource amount less than 50% by mesh refinement
Level-one；Or grid is roughened using Time step as Rule of judgment, such as initial time step, it is carried out with calculating, material resource
Amount increases, subsequent time step grid roughening；Different thickness conditions is used according to oneself calculating needs and computation rate etc. simultaneously.
Preferably, the method that area grid refinement is calculated in the step 3 specifically includes:
(1) refinement levels are set, and every level-one refinement includes to grid node number, sizing grid and grid virtual number of nodes
Refinement, implement such as formula 1:
(nx) '=2* (nx) -1
(ny) '=2* (ny) -1
(dx) '=0.5* (dx)
(dy) '=0.5* (dy)
(mbc) '=2* (mbc)
Formula 1
Wherein, nx, (nx) ' indicate to divide the direction forward and backward calculating grid x number of nodes；Ny, (ny) ' indicate to divide forward and backward meter
Calculate net (dy) ' the direction lattice y number of nodes；Dx, (dx) ' indicate to divide the direction forward and backward calculating grid x element length；Dy indicates to divide
The direction forward and backward calculating grid y element length；Mbc, (mbc) ' indicate to divide forward and backward calculating grid virtual number of nodes；
(2) to cell node assignment after refinement: respectively including carrying out assignment with the point of refinement front unit conode, in x and y
Direction newly refines node of the cell node between two cell nodes and carries out assignment, and the new node that refines is among rectangular element
Node carries out assignment；
(3) boundary condition is assigned to the region after thin draw, boundary condition includes symmetrical boundary condition, open boundaryconditions, wall
Boundary condition.
Preferably, the method that area grid roughening is calculated in the step 3 specifically includes:
(1) setting roughening grade, every level-one roughening includes grid node number, sizing grid and grid virtual number of nodes
Roughening, it is that odd number or even number save the last one if present node number is even number that every level-one roughening, which is divided into present node number,
Point is used as dummy node；
When number of nodes is odd number, every level-one roughening is such as formula 2:
(nx) '=((nx)+1)/2
(ny) '=((ny)+1)/2
(dx) '=2* (dx)
(dy) '=2* (dy)
(mbc) '=(mbc)/2
Formula 2
Wherein, nx, (nx) ' indicate to divide the direction forward and backward calculating grid x number of nodes；Ny, (ny) ' indicate to divide forward and backward meter
Calculate the direction grid y number of nodes；Dx, (dx) ' indicate to divide the direction forward and backward calculating grid x element length；Dy, (dy) ' indicate to divide
The direction forward and backward calculating grid y element length；Mbc, (mbc) ' indicate to divide forward and backward calculating grid virtual number of nodes；
(2) to cell node assignment after roughening: mainly carrying out assignment to the point of new roughening unit and former unit conode, examine
Considering number of nodes may be for even number, and assignment procedure need to judge to be roughened whether posterior nodal point is more than roughening prosthomere points in advance, then into
Row assignment；
(3) boundary condition is assigned to the region after roughening, boundary condition mainly includes symmetrical boundary condition, opens perimeter strip
Part, wall boundary condition.
Preferably, the method that area grid structuring is calculated in the step 3 specifically includes:
Which target area the dummy node unit of inquiry current region is located in；
The dummy node unit of current region is inquired in which unit of target area；
Utilize the design parameter of the dummy node unit of 5 shape function interpolation calculation current region of formula 3- formula:
Formula 3
N1=β1(x-1)(y-1)
N2=β2(x+1)(y-1)
N3=β3(x+1)(y+1)
N4=β4(x-1)(y+1)
Formula 4
Formula 5
Wherein, u indicates the relevant parameter of the dummy node unit of current region, such as height h, the momentum hu of the direction x speed,
The momentum hv etc. of the direction y speed.
Preferably, the step 4 specifically includes: the step 4 specifically includes: the numerical simulation software uses
Massflow software carries out numerical simulation meter to earth's surface disaster dynamic process according to the result of step 1, step 2 and step 3
It calculates.
Preferably, the step 5 specifically includes: the calculated result of n zoning of output, as a result with dem or
The calculated result of n zoning, is merged into according to different-format using Arcgis or tecplot by the output of tecplot format
One region is exported.
The present invention also provides a kind of extensive Landslide Hazards Fast numerical analogue simulation systems, including memory and place
Manage device, in which:
The memory is for storing program instruction；
The processor is for running described program instruction, to execute following steps:
Obtain geological disaster terrain data, source area range and elevation h and physical and mechanical parameter；
Divide extensive Landslide Hazards zoning；
The roughening of zoning grid, mesh refinement and network；
Numerical simulation calculation is carried out to earth's surface disaster dynamic process using simulation softward Massflow；
The calculated result of n zoning is exported, and the calculated result of n zoning is merged into a zoning
Output.
Compared with the prior art, the invention has the following advantages:
(1) extensive grid computing, improve computation rate: the present invention is right by carrying out region division to extensive grid
The different region of material resource ratio uses different grades of thick tessellated mesh, pays close attention to region refinement, and other area grids are roughened,
Due to many geological disasters, its motion path often only accounts for its slope surface landform sub-fraction, and above-mentioned strategy had both been avoided that extensive
Large-scale area is narrowed down to the zonule paid close attention to by the calculating in some non-material resource regions in region, guarantees computational accuracy,
Calculation amount can be saved again, improve computation rate.
(2) dynamic adjusts the thickness of zoning grid: the present invention not only needs and calculates speed to calculate according to oneself
Rate etc. uses different thickness conditions, more it is essential that Landslide Hazards are a Dynamic Evolutions, with landslide ground
The dynamic evolution of matter disaster, grid carries out dynamic thickness according to variations such as region material resource amounts, thus change of the grid with calculating
Change and changes, thus the region that more accurate calculating needs emphasis to calculate, such as in a super large basin, such as bridge pier,
The thin portions structure such as debris dam and drainage groove needs fine grid blocks that can just portray local feature, therefore structured grid and with calculating
Process grid divides with the obvious advantage again.
(3) structured grid and grid node interaction: the present invention is to realize the feature modeling in region and intending for boundary node
It closes, need to grid node structuring, traditional structured grid generation technique usually be utilized indirect method, be carried out by other grids
It generates, not anxious low efficiency, and there are certain error, the present invention is by by the dummy node unit adjacent domain of current region
Calculate node replace, it is theoretical by the shape function in finite element, the exchange of node is realized by " inquiry-exchange ", that is, is inquired
To where dummy node target area, after object element, the dummy node list of current region is solved using shape function interpolation
The design parameter of member.For this method relative to conventional method, computational efficiency is higher, while reducing the error after interpolation to the greatest extent.
Detailed description of the invention
It, below will be to use required in embodiment in order to illustrate more clearly of the technical solution of embodiment of the present invention
Attached drawing be briefly described, it should be understood that the following drawings illustrates only certain embodiments of the present invention, therefore is not to be seen as
It is the restriction to range, it for those of ordinary skill in the art, without creative efforts, can be with root
Other relevant attached drawings are obtained according to these attached drawings.
Fig. 1 is the result figure after landslide morphology data processing in Lishui provided in an embodiment of the present invention；
Fig. 2 is the result figure after being calculated in the embodiment of the present invention by numerical simulation software.
Specific embodiment
To keep the purposes, technical schemes and advantages of embodiment of the present invention clearer, implement below in conjunction with the present invention
The technical solution in embodiment of the present invention is clearly and completely described in attached drawing in mode, it is clear that described reality
The mode of applying is some embodiments of the invention, rather than whole embodiments.Based on the embodiment in the present invention, ability
Domain those of ordinary skill every other embodiment obtained without creative efforts, belongs to the present invention
The range of protection.
Therefore, requirement is not intended to limit to the detailed description of the embodiments of the present invention provided in the accompanying drawings below to protect
The scope of the present invention of shield, but it is merely representative of selected embodiment of the invention.Based on the embodiment in the present invention, ability
Domain those of ordinary skill every other embodiment obtained without creative efforts, belongs to the present invention
The range of protection.
Embodiment
By taking Lishui is come down as an example, a kind of extensive Landslide Hazards Fast numerical simulation method, including it is following several
A step:
1, (1) obtains image data using fixed point flight using unmanned plane by investigating Lishui landslide areas, by
Smart3D generates the Lishui DEM landslide image and DSM image.As shown in Figure 1, wherein (a) is Lishui landslide morphology digital elevation
(DEM) image (b) is Lishui landslide areas numerical cutting tool (DSM) image.
(2) it chooses wherein one section of landslide and carries out quantitative simulation calculating, while the DEM using Lishui landslide before sliding and after sliding,
Raster symbol-base is carried out using Arcgis, to obtain the landslide source area range and elevation.
Specific formula such as formula 6:
q2=(- 0.00247* (t-450)2+500.4)/6
H=0.4885*q2 0.6
Formula 6
Wherein, q2The momentum in the direction x speed is represented, t represents current time step, and h represents mud-rock flow source area elevation.
(3) measuring density of material by corresponding physical and mechanical experiment is 1000g/cm3, using coulomb model, wherein in
Poly- power is 10000, and Basement surge 0.35 does not consider interstitial hydraulic pressure.
2, more than 30 ten thousand steres are measured by Lishui landslide side, larger, therefore the region division is carried out at 25 zonules
It calculates.
3, to zoning grid thickness and structured grid, to zoning material resource amount be more than 50% can be by grid
Constantly refinement two-stage, zoning have material resource but material resource amount less than 50% by mesh refinement level-one, no material resource and do not need
The region of calculating is roughened two-stage.The thick tessellated mesh of dynamic simultaneously, i.e. guarantee grid can be flowed according to region material resource amount etc. change into
Mobile state thickness changes with the variation of calculating.
4, the above results are subjected to numerical simulation to earth's surface disaster dynamic process using numerical simulation softwares such as massflow
It calculates.
Wherein, massflow (i.e. surficial process dynamic numerical simulation software) is a earth's surface dynamic process numerical simulation
Software is that the exploitation of research team, the Chinese Academy of Sciences is completed, and pacifies a of Technology Co., Ltd.'s sale of preventing and reducing natural disasters by Chengdu mountainous region ring
Commercial software.
It is primary every mono- Time step output of 2s, the calculated result of 25 zonings of output output utilizes
Tecplot is merged into a zoning output, and some numerical results after being calculated by numerical simulation software are as shown in Figure 2.
Above-described embodiment, by dividing to zoning, the region different to material resource ratio is using different grades of
Thick tessellated mesh pays close attention to region refinement, other area grid roughening.A length of 12s when entire calculating, and utilize traditional number
It is 25s or so that value simulation (not having to region division and network), which calculates total duration, it is seen that has been obviously shortened calculating duration, has mentioned
High computation rate.Grid flows variation progress dynamic thicknessization adjustment according to region material resource amount simultaneously, pays close attention to material resource amount
More region, improves computational accuracy, can simulate extensive Landslide Hazards dynamical evolution process rapidly, thus to take precautions against natural calamities
Mitigation department provides faster more accurate data and supports.
The foregoing is merely the preferred embodiment of the present invention, are not intended to restrict the invention, for this field
For technical staff, the invention may be variously modified and varied.All within the spirits and principles of the present invention, made any
Modification, equivalent replacement, improvement etc., should all be included in the protection scope of the present invention.
Claims (9)
1. a kind of extensive Landslide Hazards Fast numerical simulation method, which comprises the steps of:
Step 1: obtaining Landslide Hazards terrain data, source area range, elevation and physical and mechanical parameter；
Step 2: dividing extensive Landslide Hazards zoning；
Step 3: thick tessellated mesh condition being judged, and to the roughening of zoning grid, mesh refinement and network；
Step 4: numerical simulation calculation being carried out to earth's surface disaster dynamic process using numerical simulation software；
Step 5: the calculated result of n zoning of output, and the calculated result of n zoning is merged into a calculating area
Domain output.
2. a kind of extensive geological disaster Fast numerical simulation method as described in claim 1, which is characterized in that described
The method of Landslide Hazards terrain data is obtained in step 1 specifically: acquire terrain data using unmanned plane, and built with three-dimensional
Mould Software Create dem data.
3. a kind of extensive Landslide Hazards Fast numerical simulation method as described in claim 1, which is characterized in that
The method that Landslide Hazards source area range and elevation acquisition are obtained in the step 1 specifically: front and back occurs using disaster
Unmanned plane or remote sensing altitude data, subtracted each other by grid, obtain source area range and elevation.
4. a kind of extensive Landslide Hazards Fast numerical simulation method as described in claim 1, which is characterized in that
The method of landslide disaster physical and mechanical parameter is obtained in the step 1 specifically: by related physical experiment of machanics or pass through examination
The mode of calculation obtains landslide disaster volume density, coefficient of friction parameter etc..
5. a kind of extensive Landslide Hazards Fast numerical analogue simulation side as described in claim 1-4 any claim
Method, which is characterized in that the method for extensive Landslide Hazards zoning is divided in the step 2 specifically: using simply
Be divided by and complementation, obtained extensive Landslide Hazards dem data is divided into n*n region, then will before great Qu
Calculating parameter involved in domain is assigned to n*n region, and the equivalence after guaranteeing region division, wherein n is according to different calculating
Precision and computer capacity carry out value.
6. a kind of extensive Landslide Hazards Fast numerical analogue simulation side as described in claim 1-4 any claim
Method, which is characterized in that
The method judged in the step 3 thick tessellated mesh condition specifically includes:
Grid thickness is carried out to different zonings according to computing object, using material resource as Rule of judgment, i.e., no material resource,
Do not need the area grid calculated roughening, there is a material resource calculative area grid refinement, simultaneous selection roughening, refinement etc.
Grade, such as zoning material resource amount are more than 50% grid constantly to refine to two-stage, and there is material resource in zoning but material resource amount is less than
50% by mesh refinement level-one；
The method that area grid refinement is calculated in the step 3 specifically includes:
(1) refinement levels are set, and every level-one refinement includes to the thin of grid node number, sizing grid and grid virtual number of nodes
Change, specific implementation are as follows:
(nx) '=2* (nx) -1
(ny) '=2* (ny) -1
(dx) '=0.5* (dx)
(dy) '=0.5* (dy)
(mbc) '=2* (mbc)
Formula 1
Wherein, nx, (nx) ' indicate to divide the direction forward and backward calculating grid x number of nodes；Ny, (ny) ' indicate to divide forward and backward calculating net
(dy) ' the direction lattice y number of nodes；Dx, (dx) ' indicate to divide the direction forward and backward calculating grid x element length；Dy indicates to divide forward and backward
Calculate the direction grid y element length；Mbc, (mbc) ' indicate to divide forward and backward calculating grid virtual number of nodes；
(2) to cell node assignment after refinement: respectively including carrying out assignment with the point of refinement front unit conode, in the direction x and y
Node of the new refinement cell node between two cell nodes carries out assignment, the new node for refining node among rectangular element
Carry out assignment；
(3) boundary condition is assigned to the region after thin draw, boundary condition mainly includes symmetrical boundary condition, open boundaryconditions, wall
Boundary condition；
The method that area grid roughening is calculated in the step 3 specifically includes:
(1) setting roughening grade, every level-one roughening includes the thick of grid node number, sizing grid and grid virtual number of nodes
Change, it is odd number or even number that every level-one roughening, which is divided into present node number, if present node number is even number, by the last one node
As dummy node；
When number of nodes is odd number, every level-one roughening is such as formula 2:
(nx) '=((nx)+1)/2
(ny) '=((ny)+1)/2
(dx) '=2* (dx)
(dy) '=2* (dy)
(mbc) '=(mbc)/2
Formula 2
Wherein, nx, (nx) ' indicate to divide the direction forward and backward calculating grid x number of nodes；Ny, (ny) ' indicate to divide forward and backward calculating net
The direction lattice y number of nodes；Dx, (dx) ' indicate to divide the direction forward and backward calculating grid x element length；Before dy, (dy) ' expression division,
The direction grid y element length is calculated afterwards；Mbc, (mbc) ' indicate to divide forward and backward calculating grid virtual number of nodes；
(2) to cell node assignment after roughening: mainly carrying out assignment to the point of new roughening unit and former unit conode, it is contemplated that
Number of nodes may be even number, and assignment procedure need to judge to be roughened whether posterior nodal point is more than roughening prosthomere points in advance, then be assigned
Value；
(3) boundary condition is assigned to the region after roughening, boundary condition mainly includes symmetrical boundary condition, open boundaryconditions, wall
Boundary condition；
The method that area grid structuring is calculated in the step 3 specifically includes:
(1) which target area the dummy node unit for inquiring current region is located in；
(2) the dummy node unit of current region is inquired in which unit of target area；
(3) design parameter of the dummy node unit of 5 shape function interpolation calculation current region of formula 3- formula is utilized:
N1=β1(x-1)(y-1)
N2=β2(x+1)(y-1)
N3=β3(x+1)(y+1)
N4=β4(x-1)(y+1)
Formula 4
Wherein, u indicates the relevant parameter of the dummy node unit of current region.
7. a kind of extensive Landslide Hazards Fast numerical analogue simulation side as described in claim 1-4 any claim
Method, which is characterized in that the step 4 specifically includes: the numerical simulation software use massflow software, according to step 1,
The result of step 2 and step 3 carries out numerical simulation calculation to earth's surface disaster dynamic process.
8. a kind of extensive Landslide Hazards Fast numerical analogue simulation side as described in claim 1-4 any claim
Method, which is characterized in that the step 5 specifically includes: the calculated result of n zoning of output, as a result with dem or
The calculated result of n zoning, is merged into according to different-format using Arcgis or tecplot by the output of tecplot format
One region is exported.
9. a kind of extensive Landslide Hazards Fast numerical analogue simulation system, which is characterized in that including memory and processing
Device, in which:
The memory is for storing program instruction；
The processor is for running described program instruction, to execute following steps:
Obtain Landslide Hazards terrain data, source area range, elevation and physical and mechanical parameter；
Divide extensive Landslide Hazards zoning；
Thick tessellated mesh condition is judged, and to the roughening of zoning grid, mesh refinement and network；
Numerical simulation calculation is carried out to earth's surface disaster dynamic process using numerical simulation software；
The calculated result of n zoning is exported, and it is defeated that the calculated result of n zoning is merged into a zoning
Out.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811166665.4A CN109446587A (en) | 2018-09-30 | 2018-09-30 | A kind of extensive Landslide Hazards Fast numerical simulation method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811166665.4A CN109446587A (en) | 2018-09-30 | 2018-09-30 | A kind of extensive Landslide Hazards Fast numerical simulation method and system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109446587A true CN109446587A (en) | 2019-03-08 |
Family
ID=65544514
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811166665.4A Pending CN109446587A (en) | 2018-09-30 | 2018-09-30 | A kind of extensive Landslide Hazards Fast numerical simulation method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109446587A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020221284A1 (en) * | 2019-04-29 | 2020-11-05 | 黄河勘测规划设计研究院有限公司 | Unmanned aerial vehicle monitoring method and system for basin-wide flood scene |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009074000A1 (en) * | 2007-12-11 | 2009-06-18 | Shengxie Xiao | Dynamic tracking method of depot geologic hazard evolvement process |
| CN106683184A (en) * | 2017-01-04 | 2017-05-17 | 朱军 | Mud-rock flow disaster process rapid simulation and visualization analysis method in network environment |
| CN107506566A (en) * | 2017-10-16 | 2017-12-22 | 中国科学院、水利部成都山地灾害与环境研究所 | A kind of new dynamics of debris flow Numerical Analysis methods and system |
| CN107860367A (en) * | 2017-11-03 | 2018-03-30 | 广州地理研究所 | A kind of Group-occurring landslides volume rapid extracting method based on low latitude unmanned aerial vehicle remote sensing |
| CN108010278A (en) * | 2017-12-25 | 2018-05-08 | 中国科学院、水利部成都山地灾害与环境研究所 | Mud-stone flow disaster dangerous situation dynamic early-warning method, the classified Monitoring that becomes more meticulous method for early warning |
-
2018
- 2018-09-30 CN CN201811166665.4A patent/CN109446587A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009074000A1 (en) * | 2007-12-11 | 2009-06-18 | Shengxie Xiao | Dynamic tracking method of depot geologic hazard evolvement process |
| CN106683184A (en) * | 2017-01-04 | 2017-05-17 | 朱军 | Mud-rock flow disaster process rapid simulation and visualization analysis method in network environment |
| CN107506566A (en) * | 2017-10-16 | 2017-12-22 | 中国科学院、水利部成都山地灾害与环境研究所 | A kind of new dynamics of debris flow Numerical Analysis methods and system |
| CN107860367A (en) * | 2017-11-03 | 2018-03-30 | 广州地理研究所 | A kind of Group-occurring landslides volume rapid extracting method based on low latitude unmanned aerial vehicle remote sensing |
| CN108010278A (en) * | 2017-12-25 | 2018-05-08 | 中国科学院、水利部成都山地灾害与环境研究所 | Mud-stone flow disaster dangerous situation dynamic early-warning method, the classified Monitoring that becomes more meticulous method for early warning |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020221284A1 (en) * | 2019-04-29 | 2020-11-05 | 黄河勘测规划设计研究院有限公司 | Unmanned aerial vehicle monitoring method and system for basin-wide flood scene |
| GB2590192A (en) * | 2019-04-29 | 2021-06-23 | Yellow River Engineering Consulting Co Ltd | Unmanned aerial vehicle monitoring method and system for basin-wide flood scene |
| GB2590192B (en) * | 2019-04-29 | 2023-09-27 | Yellow River Engineering Consulting Co Ltd | Unmanned aerial vehicle (UAV)-based monitoring method and system for basin-wide flood |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Xia et al. | A full-scale fluvial flood modelling framework based on a high-performance integrated hydrodynamic modelling system (HiPIMS) | |
| Kim et al. | Coupled modeling of hydrologic and hydrodynamic processes including overland and channel flow | |
| Kim et al. | Mesh type tradeoffs in 2D hydrodynamic modeling of flooding with a Godunov-based flow solver | |
| CN110750864B (en) | Dam break flood demonstration early warning method and platform based on BIM technology | |
| Begnudelli et al. | Conservative wetting and drying methodology for quadrilateral grid finite-volume models | |
| Kuiry et al. | Coupled 1D–Quasi-2D flood inundation model with unstructured grids | |
| Mizukami et al. | mizuRoute version 1: a river network routing tool for a continental domain water resources applications | |
| CN107451372B (en) | A kind of flood of a mountain area numerical simulation method that kinematic wave is combined with dynamic wave | |
| Wu et al. | An implicit 2-D shallow water flow model on unstructured quadtree rectangular mesh | |
| Avolio et al. | SCIDDICA-SS 3: a new version of cellular automata model for simulating fast moving landslides | |
| Zhao et al. | Novel variable reconstruction and friction term discretisation schemes for hydrodynamic modelling of overland flow and surface water flooding | |
| CN109190300A (en) | A kind of extensive Geological Hazards of debris Fast numerical simulation method and system | |
| Stecca et al. | Development of a numerical model for braided river morphology and vegetation evolution with application to the Lower Waitaki River (Aotearoa–New Zealand) | |
| Qu et al. | Numerical investigation of breaking solitary wave runup on permeable sloped beach using a nonhydrostatic model | |
| Krámer et al. | Solution-adaptivity in modelling complex shallow flows | |
| Jiang et al. | A dynamic bidirectional coupled surface flow model for flood inundation simulation | |
| CN109446587A (en) | A kind of extensive Landslide Hazards Fast numerical simulation method and system | |
| Unterweger et al. | Integration of FullSWOF2D and PeanoClaw: Adaptivity and local time-stepping for complex overland flows | |
| Christelis et al. | Employing surrogate modelling for the calibration of a 2D flood simulation model | |
| Roy et al. | A study on hydrodynamic and morphological behavior of Padma river using Delft3d model | |
| Vashist et al. | Minimisation of overestimation of river flows in 1D-hydrodynamic modeling | |
| Hraniciuc et al. | Creating Flood Hazard Maps Using 2d Hydraulic Models | |
| Cartier et al. | Longshore sediment transport measurements on sandy macrotidal beaches compared with sediment transport formulae | |
| Reis et al. | Wave overtopping of a porous structure: numerical and physical modeling | |
| CN106408661A (en) | Adaptive mixed interpolation method based on geological surface local complexity |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190308 |
|
| WD01 | Invention patent application deemed withdrawn after publication |