CN113987888A - Bedrock incident-downlink wave number value separation method based on array observation - Google Patents
Bedrock incident-downlink wave number value separation method based on array observation Download PDFInfo
- Publication number
- CN113987888A CN113987888A CN202111355758.3A CN202111355758A CN113987888A CN 113987888 A CN113987888 A CN 113987888A CN 202111355758 A CN202111355758 A CN 202111355758A CN 113987888 A CN113987888 A CN 113987888A
- Authority
- CN
- China
- Prior art keywords
- wave
- incident
- observation
- array
- bedrock
- 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
- 238000000926 separation method Methods 0.000 title claims abstract description 17
- 239000011159 matrix material Substances 0.000 claims abstract description 16
- 238000004364 calculation method Methods 0.000 claims abstract description 15
- 239000000463 material Substances 0.000 claims description 8
- 230000035484 reaction time Effects 0.000 claims description 7
- 239000002689 soil Substances 0.000 claims description 6
- 239000011435 rock Substances 0.000 claims description 5
- 238000000547 structure data Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 abstract description 11
- 238000000034 method Methods 0.000 abstract description 9
- 238000010586 diagram Methods 0.000 description 9
- 238000004088 simulation Methods 0.000 description 6
- 238000005452 bending Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 230000002411 adverse Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
Images
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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a matrix incident-downlink wave number value separation method based on array observation, which comprises the steps of establishing a two-dimensional stratified field finite element model according to wave speed structure information of a field where an array is located, inputting array matrix observation waves into grid nodes at the bottom, extracting a counter force time course RF (t) from middle grid nodes at the bottom after submitting operation, dividing the counter force time course by the grid node width L to obtain a node counter force time course F (t) with unit length, and converting the station matrix observation waves a (t) into an equivalent node force time course F (t)1And (t) calculating an equivalent node force time course f (t) of the bedrock incident wave, substituting the equivalent node force time course of the bedrock incident wave into an equivalent node force calculation formula, and performing inverse solution calculation to obtain a bedrock incident wave u (t). The method can separate the incident wave and the down-going wave reflected by the earth surface in the observation record of the vertical array bedrock, obtain the accurate input data for the earthquake-resistant analysis of the underground structure, and realize the accurate analysis of the earthquake field reaction.
Description
Technical Field
The invention relates to a matrix incident-downlink wave number value separation method based on array observation, which can be used for separating incident waves contained in a field matrix observation record of a vertical array of a global strong earthquake observation platform network and belongs to the technical field of earthquake engineering.
Background
The complex field reaction caused by the near-surface propagation of the vibration is always an important subject of seismic engineering. Researches show that the propagation process of seismic motion on the near-surface comprises incident waves transmitted from bedrock to the surface and downlink waves reflected to the bedrock through the surface, which causes certain difficulties in the analysis of the propagation process of the seismic on the near-surface and the seismic analysis of underground structures. The accurate input of the earthquake motion of the bedrock is very important, the existing simulation analysis underground structure earthquake resistance mostly adopts a mode of combining equivalent node force and a viscoelastic boundary as an earthquake motion input mode and a bottom boundary condition, bedrock waves recorded by a typical station field close to a researched engineering field are selected as input waves, and the wave band is superposed waves of incident waves and reflected waves of the bedrock of the station observation field, so that the superposed waves are directly converted into estimated earthquake motion input energy with overhigh equivalent node force, and the adverse effect is generated on the analysis of damage of the underground structure in the earthquake.
Based on the fact that the existence of the down waves causes that the bedrock observation records are directly converted into seismic force input is inaccurate, and a more reasonable calculation result can be obtained only by converting bedrock incident waves into the seismic force as the input, so that the research of the incident wave separation method has important value for analyzing the underground structure seismic design.
Disclosure of Invention
In order to solve the defects of the prior art, the invention provides a bedrock incident-downlink wave number value separation method based on array observation, which can separate incident waves and downlink waves reflected by the earth surface in the observation records of the bedrock of a vertical array by applying finite element numerical simulation software ABAQUS and the analysis algorithm provided by the invention, obtain accurate input data for performing earthquake-resistant analysis of an underground structure and realize accurate reaction analysis of an earthquake field.
The technical scheme adopted by the invention for solving the technical problem is as follows: the utility model provides a bedrock incident-downlink wave number value separation method based on array observation, which comprises the following steps:
s1, establishing a two-dimensional layered site finite element model according to wave velocity structure data, provided by a seismic observation table network, of S waves and P waves of a site soil layer where a vertical array is located, changing along the depth of the soil layer, inputting an array bed rock observation wave a (t) recorded by an array bed rock at the bottom of the two-dimensional layered site finite element model, establishing a dynamic analysis step, and submitting the dynamic analysis step to matrix iterative operation;
s2, extracting a mesh node reaction time course RF (t) at the middle node at the bottom of the calculated two-dimensional layered field finite element model, and dividing the mesh node reaction time course RF (t) by the mesh node width L to obtain a node reaction time course F (t) of unit length;
s3, converting the array bedrock observation wave a (t) into an equivalent node force time course F by adopting an equivalent node force calculation formula1(t)；
S4, calculating the equivalent node force time course f (t) of the bedrock incident wave by the following formula:
f(t)＝(F(t)+F1(t))/2
and S5, substituting the equivalent node force time course of the bedrock incident wave into an equivalent node force calculation formula for inverse solution calculation to obtain an incident velocity time course u (t).
In step S1, a two-dimensional layered site finite element model is established by using ABAQUS/CAE, and both the lateral boundary and the bottom boundary of the two-dimensional layered site finite element model adopt vertical constraints to control the mesh nodes of the two-dimensional layered site finite element model to move only in the horizontal direction.
In step S3, the following equivalent node force calculation formula is used to convert the array bedrock observation wave path a (t) into an equivalent node force time path F1(t)：
F1(t)＝AbVsρu0(t)
Wherein A isbDenotes the grid size, p denotes the material density, VsDenotes the material shear wave velocity u0(t) represents the observation speed of the bedrock of the station, and is calculated by the following formula:
step S5 specifically uses the following formula to calculate the incident velocity time interval u (t):
u(t)＝f(t)/(AbρVs)。
the invention has the beneficial effects based on the technical scheme that: the existing simulation analysis of underground structure earthquake resistance mostly adopts a mode of combining equivalent node force and viscoelastic boundary as an earthquake motion input form and bottom boundary condition, selects basement wave recorded by a typical station field close to a researched engineering field as input wave, but the wave band is superposition wave of basement incident wave and reflected wave of the station observation field, directly converts the superposition wave into energy of estimation earthquake motion input with overhigh equivalent node force, and generates adverse effect on the evaluation of damage of the underground structure in the earthquake. The matrix incident-downlink wave number value separation method based on array observation provided by the invention can separate incident waves in vertical array matrix observation records from the downlink waves reflected by the earth surface, obtain accurate input data for underground structure earthquake resistance analysis, realize accurate reaction analysis of earthquake fields, and has important engineering significance for underground structure earthquake resistance simulation.
Drawings
Fig. 1 is a schematic flow chart of a matrix incident-downlink wave number value separation method based on array observation provided by the invention.
Figure 2 is a schematic diagram of a vertical array bedrock and surface accelerometer arrangement.
FIG. 3 is a schematic diagram of a two-dimensional layered yard finite element model.
Fig. 4 is a schematic structural diagram of the S wave and P wave speed of the soil layer of the field where the array is located.
Fig. 5 is a schematic diagram of incident waves and observed waves under an ideal two-dimensional field model.
FIG. 6 is a schematic diagram of an equivalent node force and a unit node counter force time course.
Fig. 7 is a schematic diagram of the time course of the incident wave calculated by the separation of the incident-downlink wave number values from the theoretical solution.
Fig. 8 is a schematic diagram of seismic simulation analysis of a two-dimensional field tunnel cross section.
FIG. 9 is a schematic diagram of the effect of separating the downgoing wave on the bending moment envelope of the tunnel cross section.
Detailed Description
The invention is further illustrated by the following figures and examples.
The invention provides a bedrock incident-downlink wave number value separation method based on array observation, which comprises the following steps with reference to fig. 1:
s1, arranging the seismic observation station net as shown in figure 2, providing wave velocity structure data of S waves and P waves of a soil layer of a vertical station array as shown in figure 4, and establishing a two-dimensional layered site finite element model shown in figure 3 in ABAQUS/CAE according to the data. The ABAQUS/CAE is a preprocessing module of ABAQUS software and can directly establish a finite element model. The height of the two-dimensional layered field finite element model is the drilling depth of a vertical array, square grids with the side length of L are used for grid dispersion in ABAQUS, grid nodes of the lateral boundary and the bottom boundary of the two-dimensional layered field finite element model both adopt vertical constraint conditions provided by ABAQUS, the lateral boundary and the bottom boundary grid nodes of the two-dimensional layered field finite element model are controlled to move only along the horizontal direction, the acceleration time course a (t) recorded by array bed rock is input to the grid nodes at the bottom of the model, and dynamic analysis step submission matrix iterative operation is established.
S2, extracting a reaction force time course RF (t) at the middle node at the bottom of the calculated two-dimensional layered field finite element model, and dividing the node reaction force time course RF (t) by the grid node width L to obtain a node reaction force time course F (t) with unit length.
S3, converting the array bedrock observation wave path a (t) into an equivalent node force time path F by adopting the following equivalent node force calculation formula1(t)：
F1(t)＝AbVsρu0(t)
Wherein A isbDenotes the grid size, pDenotes the density of the material, VsDenotes the material shear wave velocity u0(t) represents the observation speed of the bedrock of the station, and is calculated by the following formula:
s4, calculating the equivalent node force time course f (t) of the bedrock incident wave by the following formula:
f(t)＝(F(t)+F1(t))/2
s5, substituting the equivalent node force time course of the bedrock incident wave into an equivalent node force calculation formula to perform inverse solution calculation to obtain the bedrock incident wave, and specifically calculating by adopting the following formula to obtain an incident velocity time course u (t):
u(t)＝f(t)/(AbρVs)。
the feasibility can be verified through an ideal two-dimensional field model:
a finite field with the height of 50m and the length of 2m is established in the ABAQUS/CAE to serve as a research area, and middle nodes at the top and the bottom of an ideal two-dimensional field finite element model serve as monitoring points A and B of the model. The material parameters were as follows: density rho 1000kg/m3Elastic modulus E is 24MPa, Poisson ratio V is 0.2, and material shear wave velocity VsAnd (3) performing grid discretization on the ideal two-dimensional site finite element model by adopting a four-node plane strain unit (CPE4) at 100m/s, wherein the unit size is equal to Δ x and equal to Δ y and equal to 1 m. Applying a unit pulse wave of the following formula to the bottom boundary in the X direction:
under an ideal two-dimensional field model, the obtained pulse wave is used as an input wave, converted into an equivalent node force adapted to a viscoelastic boundary and used as the input wave to the model bedrock, and the time course of the observation wave at the bottom of the model is shown in fig. 5.
Regarding the model bottom observation wave of FIG. 5 as an array bedrock observation record, extracting a middle node reaction time course RF (t) at the bottom of the model, wherein the node reaction is a unit node reaction F (t) because the size of a model unit is 1 m; the time courses of the bedrock observation wave, F1(t), F (t) and F1(t) are converted by applying the formula (1) to the equivalent node force time course as shown in fig. 6, and it can be seen that the first peaks of F (t) and F1(t) are identical and the second peaks are opposite.
Unit node reaction force F (t) at bottom of model and equivalent node force time course F recorded by bedrock observation1The first peak of (t) is the same, and the second peak is exactly opposite, so F (t) ═ F (t) + F is used1(t))/2, the equivalent node force f (t) corresponding to the incident wave can be calculated, and u (t) is obtained by inverse solution of an equivalent node force calculation formula, wherein u (t) is the bedrock wave time course, and the node force f (t) and the bedrock wave u (t) are shown in fig. 7. As can be seen from fig. 7, the time course of incident wave displacement obtained by applying the method is basically consistent with the theoretical solution of fig. 5, and the feasibility and effectiveness of the method are proved.
The method provided by the invention is applied to the anti-seismic simulation analysis of the cross section of a tunnel in a two-dimensional field, the matrix incident wave obtained by separating the matrix incident wave value and the downlink wave value in the array observation field is used as the matrix seismic input, and the influence of the separated and undiseparated downlink wave on the anti-seismic analysis of the tunnel is compared:
fig. 8 is a schematic diagram of seismic simulation analysis of a two-dimensional field tunnel cross section, in which acceleration time courses of separated and non-separated traveling waves are applied to the bottom boundary of a numerical model of the two-dimensional field tunnel cross section as input waves respectively, and a bending moment envelope curve of the tunnel cross section is calculated. Fig. 9 shows the effect of separating and not separating the down-traveling wave on the bending moment envelope curve of the cross section of the tunnel, and it can be seen from fig. 9 that the incident wave obtained by applying the method of the present invention is used as the seismic input of the bedrock, and the calculated bending moment envelope curve of the cross section of the tunnel is smaller compared with the traditional input method, which shows that the energy of the seismic input is estimated in consideration of the over-high input method of the down-traveling wave separation.
According to the matrix incident-downlink wave number value separation method based on array observation provided by the invention, incident waves in vertical array matrix observation records and downlink waves reflected by the earth surface can be separated, accurate input data for underground structure seismic analysis is obtained, and accurate reaction analysis of a seismic field can be realized.
Claims (4)
1. An incident-downgoing wave field separation method based on array observation is characterized by comprising the following steps:
s1, establishing a two-dimensional layered site finite element model according to wave velocity structure data, provided by a seismic observation table network, of S waves and P waves of a site soil layer where a vertical array is located, changing along the depth of the soil layer, inputting an array bed rock observation wave a (t) recorded by an array bed rock at the bottom of the two-dimensional layered site finite element model, establishing a dynamic analysis step, and submitting the dynamic analysis step to matrix iterative operation;
s2, extracting a mesh node reaction time course RF (t) at the middle node at the bottom of the calculated two-dimensional layered field finite element model, and dividing the mesh node reaction time course RF (t) by the mesh node width L to obtain a node reaction time course F (t) of unit length;
s3, converting the array bedrock observation wave a (t) into an equivalent node force time course F by adopting an equivalent node force calculation formula1(t)；
S4, calculating the equivalent node force time course f (t) of the bedrock incident wave by the following formula:
f(t)＝(F(t)+F1(t))/2
and S5, substituting the equivalent node force time course of the bedrock incident wave into an equivalent node force calculation formula for inverse solution calculation to obtain an incident velocity time course u (t).
2. The array observation-based incident-downgoing wavefield separation method of claim 1, wherein: in step S1, a two-dimensional layered site finite element model is established by using ABAQUS/CAE, and both the lateral boundary and the bottom boundary of the two-dimensional layered site finite element model adopt vertical constraints to control the mesh nodes of the two-dimensional layered site finite element model to move only in the horizontal direction.
3. The array observation-based incident-downgoing wavefield separation method of claim 1, wherein: in step S3, the following equivalent node force calculation formula is used to convert the array bedrock observation wave path a (t) into an equivalent node force time path F1(t)：
F1(t)＝AbVsρu0(t)
Wherein A isbDenotes the grid size, p denotes the material density, VsDenotes the material shear wave velocity u0(t) represents the observation speed of the bedrock of the station, and is calculated by the following formula:
4. the array observation-based incident-downgoing wavefield separation method of claim 1, wherein: step S5 specifically uses the following formula to calculate the incident velocity time interval u (t):
u(t)＝f(t)/(AbρVs)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111355758.3A CN113987888B (en) | 2021-11-16 | 2021-11-16 | Bedrock incidence-downlink wave numerical value separation method based on array observation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111355758.3A CN113987888B (en) | 2021-11-16 | 2021-11-16 | Bedrock incidence-downlink wave numerical value separation method based on array observation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113987888A true CN113987888A (en) | 2022-01-28 |
| CN113987888B CN113987888B (en) | 2024-07-09 |
Family
ID=79748814
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111355758.3A Active CN113987888B (en) | 2021-11-16 | 2021-11-16 | Bedrock incidence-downlink wave numerical value separation method based on array observation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113987888B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108845351A (en) * | 2018-06-26 | 2018-11-20 | 中国石油大学(华东) | A kind of VSP seismic data converted wave full waveform inversion method |
| US20190250291A1 (en) * | 2018-02-09 | 2019-08-15 | China University Of Geosciences, Beijing | Method and system for acquiring probability of slope failure and destabilization caused by earthquake |
| US20200124755A1 (en) * | 2018-10-22 | 2020-04-23 | Pgs Geophysical As | Correction of source motion effects in seismic data recorded in a marine survey using a moving source |
-
2021
- 2021-11-16 CN CN202111355758.3A patent/CN113987888B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20190250291A1 (en) * | 2018-02-09 | 2019-08-15 | China University Of Geosciences, Beijing | Method and system for acquiring probability of slope failure and destabilization caused by earthquake |
| CN108845351A (en) * | 2018-06-26 | 2018-11-20 | 中国石油大学(华东) | A kind of VSP seismic data converted wave full waveform inversion method |
| US20200124755A1 (en) * | 2018-10-22 | 2020-04-23 | Pgs Geophysical As | Correction of source motion effects in seismic data recorded in a marine survey using a moving source |
Non-Patent Citations (1)
| Title |
|---|
| 张雨霆;肖明;张志国;: "大型地下洞室群地震响应分析的动力子模型法", 岩石力学与工程学报, no. 2, 15 September 2011 (2011-09-15) * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113987888B (en) | 2024-07-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108090268B (en) | Integrated addition method of seismic time-course wave under viscoelastic boundary | |
| CN112287578B (en) | Slope stability calculation method under end slope exploitation condition | |
| CN103413185B (en) | A kind of rocker arm of coal mining machine vibration transducer optimization placement method | |
| Liu et al. | A coupled mathematical model for accumulation of wave-induced pore water pressure and its application | |
| Gao et al. | Strength reduction model for jointed rock masses and peridynamics simulation of uniaxial compression testing | |
| CN111737895A (en) | Method for dynamically evaluating stability of roof of underground goaf of strip mine | |
| CN104992013A (en) | Mountain tunnel anti-seismic analysis method | |
| CN107203652A (en) | The analogy method that becomes more meticulous of underground structure floating centrifuge test in earthquake liquefaction | |
| CN113987888B (en) | Bedrock incidence-downlink wave numerical value separation method based on array observation | |
| CN109063324A (en) | Finite element crustal stress analogue technique method based on Corner-point Grids | |
| CN108427849A (en) | It is a kind of fall story frame structure earthquake proof performance analysis method | |
| CN112364543A (en) | Seepage simulation method based on soft plastic loess tunnel double-row type cluster well precipitation model | |
| CN107066744B (en) | Earthquake-resistant analysis method for subway station in water-rich stratum | |
| CN114660659B (en) | Longitudinal and transverse wave decoupling high-precision microseism positioning method | |
| CN116432498A (en) | Earthquake motion input implementation method suitable for slope power analysis | |
| CN116822267A (en) | Efficient parallel computing simulation method considering local discontinuous domain based on OPenSeesMP | |
| Kato et al. | Fully integrated 3D analysis on site-city interaction in an urban transport hub | |
| CN115392032A (en) | GIS-MPM seamless integrated dynamic three-dimensional geological model construction method | |
| Shafiqu et al. | Numerical analysis of a pile-soil system under earthquake loading | |
| Kamijo et al. | Seismic tests of a pile-supported structure in liquefiable sand using large-scale blast excitation | |
| CN106646605A (en) | Complex fracture zone rock formation curvature optimization calculation method | |
| KR100592219B1 (en) | An Estimation Method of Shear Wave Velocities and Damping Ratios in Soil Layers Using Forced Vibration Test Data | |
| Ye et al. | A generalized framework of two-way coupled numerical model for fluid-structure-seabed interaction (FSSI): Explicit algorithm | |
| Yang et al. | ALEM-FEM-BEM coupling for response analysis of pile group in laminated fluid-saturated geomaterials induced by tunnel excavation | |
| Elgamal | Nonlinear modeling of large-scale ground-foundation-structure seismic response |
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 |