CN109190300A - A kind of extensive Geological Hazards of debris Fast numerical simulation method and system - Google Patents
A kind of extensive Geological Hazards of debris Fast numerical simulation method and system Download PDFInfo
- Publication number
- CN109190300A CN109190300A CN201811166384.9A CN201811166384A CN109190300A CN 109190300 A CN109190300 A CN 109190300A CN 201811166384 A CN201811166384 A CN 201811166384A CN 109190300 A CN109190300 A CN 109190300A
- Authority
- CN
- China
- Prior art keywords
- grid
- debris
- zoning
- mud
- geological hazards
- 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 60
- 238000004088 simulation Methods 0.000 title claims abstract description 38
- 238000013316 zoning Methods 0.000 claims abstract description 44
- 238000007788 roughening Methods 0.000 claims abstract description 31
- 239000000463 material Substances 0.000 claims abstract description 28
- 230000008569 process Effects 0.000 claims abstract description 19
- 238000004364 calculation method Methods 0.000 claims abstract description 16
- 239000011435 rock Substances 0.000 claims description 32
- 230000006870 function Effects 0.000 claims description 5
- 238000004513 sizing Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 238000013461 design Methods 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 claims description 3
- 239000004575 stone Substances 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims description 2
- 238000007670 refining Methods 0.000 claims 1
- 230000033001 locomotion Effects 0.000 description 7
- 230000008901 benefit Effects 0.000 description 4
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 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
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004540 process dynamic Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 230000009897 systematic effect Effects 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
-
- 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- General Business, Economics & Management (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of extensive Geological Hazards of debris Fast numerical simulation method and system, include the following steps: to obtain Geological Hazards of debris terrain data, source area range, elevation and physical and mechanical parameter；Divide extensive Geological Hazards of debris 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 the calculated result of n zoning is merged into a zoning output；This method is by carrying out region division, grid thickness and structured grid to zoning, to there are material resource, calculative region to carry out fining calculating, there is no material resource, need the area grid calculated roughening, it not only realizes and the large-scale grid that calculates is calculated, and computation rate is improved, realize the Fast simulation calculating to extensive Geological Hazards of debris.
Description
Technical field
The present invention relates to geological disaster dynamic process numerical simulation calculating field more particularly to a kind of extensive mud-rock flow
Matter 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.
At this stage to Geological Hazards of debris dynamic process numerical simulation be based primarily upon continuous media, based on it is discrete and between
Three kinds of methods of the two.Various method for numerical simulation have its unique function and solve the ability of problems, and answer in synthesis
A large amount of fruitful research achievements are achieved in.It is all earth's surface object under the force of gravity in view of geological disasters such as 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 the quick of the dynamic process of extensive Geological Hazards of debris (especially disaster chain) and kinetic characteristic
Analogue 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,
Number of grid is more, and calculation scale is larger；2. computation rate is slow, it is long to calculate the time.And the assessment after disaster generation, to the condition of a disaster
It is very urgent, how quickly to accomplish the assessment to extensive mud-stone flow disaster or even to quick recognition before large scale disasters generation and in advance
Motion mode, track and the damaging range of disaster are surveyed, timely data support is provided for the correlation department of preventing and reducing natural disasters, to subtract
It is few even to avoid the problem that disaster to loss caused by the mankind and Human dried bloodstains is still current urgent need to resolve.
It is established with structured grid to extensive Geological Hazards of debris dynamic process currently, there is no and drawn based on region again
Fast numerical simulation method and systematic account.
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 Geological Hazards of debris Fast numerical simulation method and system of structured grid, gram
Take the reality big to extensive Geological Hazards of debris numerical simulation time domain in the prior art, number of grid is more, computation rate is slow
Border problem.
To achieve the goals above, the present invention adopts the following technical scheme:
A kind of extensive Geological Hazards of debris Fast numerical simulation method, includes the following steps:
Step 1: obtaining Geological Hazards of debris terrain data, source area range, elevation and physical and mechanical parameter；
Step 2: dividing extensive Geological Hazards of debris 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 Geological Hazards of debris terrain data is obtained in the step 1 specifically: adopt using unmanned plane
Collect terrain data, and generates dem data with 3 d modeling software.
It is furthermore preferred that taking fixed point according to Geological Hazards of debris range when the acquisition terrain data using unmanned plane
Flight, S-ing or the flight of the back-shaped line of flight, the 3 d modeling software are raw using smart3D, Pix4D or Photoscan
At dem data.
Preferably, the method that Geological Hazards of debris source area range and elevation acquisition are obtained in the step 1 specifically:
It is averaged mud depth and mud-rock flow hydraulic gradient according to the flow velocity of the mud-rock flow of measurement, mud-rock flow section, using mud-rock flow
Section part ditch bed longitudinal slope is calculated the coarse rate coefficient of debris flow gully bed by 1 debris flow velocity formula of formula, and is calculated according to formula 2
Obtain mud-rock flow source area range and elevation:
VC=MC×HC2/3×iC1/2
Formula 1
Wherein, VCIndicate the flow velocity of mud-rock flow, MCIndicate debris flow gully bed roughness coefficien；HCIndicate that mud-rock flow section is average
Mud is deep (m)；iCIt indicates mud-rock flow hydraulic gradient (‰)；
H=k5/3*(q2)0.6
K=Mc*ic 1/2
q2=Q/L
Formula 2
Wherein, h indicates source area mud-rock flow height；K indicates mud-rock flow silt correction factor, passes through MC、iCIt is solved,
q2Indicate that mud-rock flow channel cross section in the speed in the direction x, is solved by flow Q, L；Q indicates horizontal by mud-rock flow channel
The flow of section；L indicates to pass through the length in mud-rock flow channel cross section.
Preferably, the method for Geological Hazards of debris physical and mechanical parameter is obtained in the step 1 specifically: pass through correlation
Physical and mechanical experiment obtains the parameters such as disaster volume density, coefficient of friction by way of tentative calculation.
Preferably, the method for extensive Geological Hazards of debris zoning is divided in the step 2 specifically: utilize letter
It is single to be divided by and complementation, obtained extensive Geological Hazards of debris dem data is divided into n*n region, then before general
Calculating parameter involved in big region is assigned to n*n region, and the equivalence after guaranteeing region division, wherein n is according to different
Computational accuracy 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 3:
(nx) '=2* (nx) -1
(ny) '=2* (ny) -1
(dx) '=0.5* (dx)
(dy) '=0.5* (dy)
(mbc) '=2* (mbc)
Formula 3
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 4:
(nx) '=((nx)+1)/2
(ny) '=((ny)+1)/2
(dx) '=2* (dx)
(dy) '=2* (dy)
(mbc) '=(mbc)/2
Formula 4
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 7 shape function interpolation calculation current region of formula 5- formula:
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 6
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 numerical simulation software uses 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.
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 Geological Hazards of debris Fast numerical analogue simulation system, including memory and
Processor, 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 Geological Hazards of debris 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 Geological Hazards of debris is a Dynamic Evolution, with mudstone
The dynamic evolution of geological disaster is flowed, grid carries out dynamic thickness according to variations such as region material resource amounts, so that grid is with calculating
Variation and change, thus the region that more accurate calculating needs emphasis to calculate, such as in a super large basin, such as bridge
The thin portions structure such as pier, debris dam and drainage groove needs fine grid blocks that can just portray local feature, therefore structured grid and with meter
Calculation 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 mud gulf ditch topographic data processing provided in an embodiment of the present invention；
Fig. 2 is in the embodiment of the present invention to region division and to division rear region zoning grid thickness and grid knot
The result figure of structure；
Fig. 3 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 the gully mud-rock flow of mud gulf as an example, a kind of extensive Geological Hazards of debris Fast numerical simulation method, including
Following steps:
1, (1) is by obtaining image data using the flight of s shape using unmanned plane, borrowing to mud gulf gully mud-rock flow Regional survey
Smart3D is helped to generate DEM mud gulf ditch image and DSM image.As shown in Figure 1, wherein (a) is mud gulf ditch landform digital elevation
(DEM) image (b) is mud gulf ditch numerical cutting tool (DSM) image.
(2) it chooses wherein one section of counter drain and carries out quantitative simulation calculating, while being obtained using debris flow velocity and flow formula
The mud-rock flow source area range and elevation h, specific formula such as formula 8:
q2=(- 0.00247* (t-450)2+500.4)/6
H=0.4885*q2 0.6
Formula 8
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 1825g/cm3, and cohesive force 10000, substrate is rubbed
Wiping coefficient is 0.35.
2, mud gulf Watershed area 10.3km2, it is larger, it can only be calculated by the region division at 36 zonules
Mud-rock flow movement path sections, other parts not calculate.
3, to zoning grid thickness and structured grid, to zoning material resource amount be more than 50% by grid not
Disconnected 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 to count
The region of calculation is roughened two-stage.The thick tessellated mesh of dynamic simultaneously, i.e. guarantee grid are moved according to variations such as region material resource amount flowings
State thickness changes with the variation of calculating.As shown in Fig. 2, wherein (a), (b), (c), (d) output respectively time step be 1s,
Result after 3s, 5s, 7s.
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 of 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 36 zonings of output output utilizes
Tecplot is merged into a zoning output.Some numerical results after being calculated by numerical simulation software are as shown in Figure 3.
Due to many geological disasters, its motion path often only accounts for its slope surface landform sub-fraction, and above-mentioned strategy both can guarantee
Precision can save calculation amount again, improve computation rate.In a super large basin, bridge pier, debris dam and drainage groove etc. are thin
Portion's structure needs fine grid blocks that can just portray local feature, therefore structured grid and divides advantage again with calculating process grid
Obviously.
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 30.2s when entire calculating, and utilize traditional
It is 40s or so that numerical simulation (not having to region division and network), which calculates total duration, sees and has been obviously shortened calculating duration, mentions
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, simulates extensive geological disaster dynamical evolution process rapidly, thus for the portion that prevents and reduces natural disasters
Door 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 Geological Hazards of debris Fast numerical simulation method, which comprises the steps of:
Step 1: obtaining Geological Hazards of debris terrain data, source area range, elevation and physical and mechanical parameter；
Step 2: dividing extensive Geological Hazards of debris 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 Geological Hazards of debris terrain data is obtained in step 1 specifically: acquire terrain data using unmanned plane, and with three-dimensional
Modeling software generates dem data.
3. a kind of extensive Geological Hazards of debris Fast numerical simulation method as described in claim 1, feature exist
In the method for obtaining Geological Hazards of debris source area range and elevation acquisition in the step 1 specifically:
It is averaged mud depth and mud-rock flow hydraulic gradient according to the flow velocity of the mud-rock flow of measurement, mud-rock flow section, using mudstone flow section
Locate ditch bed longitudinal slope, the coarse rate coefficient of debris flow gully bed is calculated by 1 debris flow velocity formula of formula, and mud is calculated to obtain according to formula 2
Rock glacier source area range and elevation:
VC=MC×HC2/3×iC1/2
Formula 1
Wherein, VCIndicate the flow velocity of mud-rock flow, MCIndicate debris flow gully bed roughness coefficien；HCIndicate that the mud-rock flow section mud that be averaged is deep
(m)；iCIt indicates mud-rock flow hydraulic gradient (‰)；
H=k5/3*(q2)0.6
K=Mc*ic 1/2
q2=Q/L
Formula 2
Wherein, h indicates source area mud-rock flow height；K indicates mud-rock flow silt correction factor, passes through MC、iCIt is solved, q2Table
Show that mud-rock flow channel cross section in the speed in the direction x, is solved by flow Q, L；Q indicates to pass through mud-rock flow channel cross section
Flow；L indicates to pass through the length in mud-rock flow channel cross section.
4. a kind of extensive Geological Hazards of debris Fast numerical simulation method as described in claim 1, feature exist
In the method for acquisition mud-stone flow disaster physical and mechanical parameter in the step 1 specifically: by related physical experiment of machanics or lead to
The mode for crossing tentative calculation obtains the parameters such as mud-stone flow disaster volume density, coefficient of friction.
5. a kind of extensive Geological Hazards of debris Fast numerical analogue simulation as described in claim 1-4 any claim
Method, which is characterized in that the method for extensive Geological Hazards of debris zoning is divided in the step 2 specifically: utilize
Simply it is divided by and complementation, obtained extensive Geological Hazards of debris dem data is divided into n*n region, then by it
Calculating parameter involved in preceding big region is assigned to n*n region, and the equivalence after guaranteeing region division, wherein n is according to difference
Computational accuracy and computer capacity carry out value.
6. a kind of extensive Geological Hazards of debris Fast numerical analogue simulation 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 3
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 4:
(nx) '=((nx)+1)/2
(ny) '=((ny)+1)/2
(dx) '=2* (dx)
(dy) '=2* (dy)
(mbc) '=(mbc)/2
Formula 4
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 7 shape function interpolation calculation current region of formula 5- 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 6
Wherein, u indicates the relevant parameter of the dummy node unit of current region.
7. a kind of extensive Geological Hazards of debris Fast numerical analogue simulation as described in claim 1-4 any claim
Method, which is characterized in that the step 4 specifically includes: the numerical simulation software uses 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 Geological Hazards of debris Fast numerical analogue simulation 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 Geological Hazards of debris Fast numerical analogue simulation system, which is characterized in that 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 Hazards of debris terrain data, source area range and elevation h and physical and mechanical parameter；
Divide extensive Geological Hazards of debris 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 simulation softward Massflow；
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 |
|---|---|---|---|
| CN201811166384.9A CN109190300A (en) | 2018-09-30 | 2018-09-30 | A kind of extensive Geological Hazards of debris Fast numerical simulation method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811166384.9A CN109190300A (en) | 2018-09-30 | 2018-09-30 | A kind of extensive Geological Hazards of debris Fast numerical simulation method and system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109190300A true CN109190300A (en) | 2019-01-11 |
Family
ID=64947208
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811166384.9A Pending CN109190300A (en) | 2018-09-30 | 2018-09-30 | A kind of extensive Geological Hazards of debris Fast numerical simulation method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109190300A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112733472A (en) * | 2021-01-11 | 2021-04-30 | 中国科学院、水利部成都山地灾害与环境研究所 | Debris flow slope surface source starting amount dynamic calculation method and system |
| CN113255197A (en) * | 2021-07-13 | 2021-08-13 | 西南交通大学 | Physical mud flow erosion and deposition distance simulation method based on discrete element method |
| CN117951965A (en) * | 2024-03-26 | 2024-04-30 | 四川省蜀通勘察基础工程有限责任公司 | Debris flow material source analysis method and system |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105160193A (en) * | 2015-09-21 | 2015-12-16 | 中国科学院、水利部成都山地灾害与环境研究所 | Dynamic process based debris flow dynamic risk analysis system and method |
| JP5875669B1 (en) * | 2014-12-24 | 2016-03-02 | エー・シー・エス株式会社 | Earth and sand disaster prediction system |
| 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 |
| 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 CN201811166384.9A patent/CN109190300A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5875669B1 (en) * | 2014-12-24 | 2016-03-02 | エー・シー・エス株式会社 | Earth and sand disaster prediction system |
| CN105160193A (en) * | 2015-09-21 | 2015-12-16 | 中国科学院、水利部成都山地灾害与环境研究所 | Dynamic process based debris flow dynamic risk analysis system and method |
| 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 |
| 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 (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112733472A (en) * | 2021-01-11 | 2021-04-30 | 中国科学院、水利部成都山地灾害与环境研究所 | Debris flow slope surface source starting amount dynamic calculation method and system |
| CN112733472B (en) * | 2021-01-11 | 2021-10-26 | 中国科学院、水利部成都山地灾害与环境研究所 | Debris flow slope surface source starting amount dynamic calculation method and system |
| CN113255197A (en) * | 2021-07-13 | 2021-08-13 | 西南交通大学 | Physical mud flow erosion and deposition distance simulation method based on discrete element method |
| CN117951965A (en) * | 2024-03-26 | 2024-04-30 | 四川省蜀通勘察基础工程有限责任公司 | Debris flow material source analysis method and system |
| CN117951965B (en) * | 2024-03-26 | 2024-05-28 | 四川省蜀通勘察基础工程有限责任公司 | Debris flow material source analysis method and system |
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 | |
| Debreu et al. | Two-way nesting in split-explicit ocean models: Algorithms, implementation and validation | |
| Moodie et al. | Modeling deltaic lobe‐building cycles and channel avulsions for the Yellow River Delta, China | |
| Maheswaran et al. | Regional scale groundwater modelling study for Ganga River basin | |
| Lin et al. | Developing river rehabilitation scenarios by integrating landscape and hydrodynamic modeling for the Ciliwung River in Jakarta, Indonesia | |
| CN109190300A (en) | A kind of extensive Geological Hazards of debris Fast numerical simulation method and system | |
| Wei et al. | Long-term process-based morphodynamic modeling of the Pearl River Delta | |
| Zhao et al. | Novel variable reconstruction and friction term discretisation schemes for hydrodynamic modelling of overland flow and surface water flooding | |
| Liang et al. | New prospects for computational hydraulics by leveraging high-performance heterogeneous computing techniques | |
| Kumar | Toward a hydrologic modeling system | |
| Woodruff et al. | Subgrid corrections in finite-element modeling of storm-driven coastal flooding | |
| Sridharan et al. | Local-inertial shallow water model on unstructured triangular grids | |
| Krámer et al. | Solution-adaptivity in modelling complex shallow flows | |
| Desta et al. | Dam breach analysis: a case of Gidabo dam, Southern Ethiopia | |
| Zhang et al. | An examination of estuary stability in response to human interventions in the South Branch of the Yangtze (Changjiang) estuary, China | |
| Moritan et al. | Fluvial terrace mapping from ALOS-PALSAR data for archaeological prospection in the middle Uruguay River basin, Argentina | |
| CN115983158B (en) | Loose coupling method for groundwater model and two-dimensional hydrodynamic model | |
| Pender et al. | 12 Flood Inundation Modelling to Support Flood Risk Management | |
| Chaney et al. | Two-way coupling between the sub-grid land surface and river networks in Earth system models | |
| CN109446587A (en) | A kind of extensive Landslide Hazards Fast numerical simulation method and system | |
| Kalmbacher et al. | Effects of tides and currents on tsunami propagation in large rivers: Columbia River, United States | |
| Liu et al. | Modelling morphodynamic responses of a natural embayed beach to Typhoon Lekima encountering different tide types | |
| Overes et al. | The importance of time-varying, non-tidal currents in modelling in-situ sand wave dynamics | |
| Wu et al. | A three-dimensional model and its potential application to spring protection |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190111 |