Seismic oscillation simulation method and device based on seismic source-propagation medium-sedimentary layer

Abstract

The invention provides a seismic oscillation simulation method and device based on a seismic source-propagation medium-sedimentary layer. The method comprises the steps that a large area model and a local area model for seismic oscillation simulation are determined; the large area model comprises a local area model; simulating a seismic source fracture process of the large-area model through the finite fault model; simulating a propagation process of seismic fluctuation of the large-area model in a first propagation medium through the viscoelastic constitutive model; the seismic fluctuation is generated by a seismic source fracture process; determining the position of the local area model corresponding to the large area model, obtaining the seismic acceleration at the position, and calculating an equivalent boundary force according to the acceleration; inputting the equivalent boundary force into the local area model, and generating transmission seismic fluctuation corresponding to the equivalent boundary force in the local area model; and simulating the propagation process of the transmission seismic fluctuation of the local area model in the second propagation medium through the nonlinear model, and determining the seismic oscillation response in the propagation process.

Description

Technical Field

[0001] The present invention relates to the field of earthquake technical analysis technology, and in particular to a seismic motion simulation method based on a seismic source-propagation medium-sedimentary layer, a seismic motion simulation device based on a seismic source-propagation medium-sedimentary layer, an electronic device, and a non-transient computer-readable storage medium. Background Art

[0002] Numerical simulation methods are an important means to analyze the laws of seismic motion propagation. However, current seismic motion simulation methods face two major difficulties: one is the contradiction between the large size of large-area simulation models and the coarse grid division and the demand for refined modeling of local sites, which makes it impossible to meet the needs of seismic safety analysis of buildings on the site; the other is the problem of selecting constitutive models, that is, the need to integrate large-area viscoelastic models and local site nonlinear models to accurately simulate the strong vibration response of the entire process from the seismogenic fault to the propagation medium and then to the sedimentary layer. Summary of the invention

[0003] The present invention aims to solve at least one of the technical problems in the related art to a certain extent. To this end, the first purpose of the present invention is to propose a seismic motion simulation method based on earthquake source-propagation medium-sedimentary layer, which can accurately simulate the whole process from earthquake source rupture to medium propagation and then to the local site response of the sedimentary layer, and realize the wide-band strong vibration simulation of the whole process of earthquake source-propagation medium-sedimentary layer by connecting the large-area viscoelastic model and the local site nonlinear model through equivalent boundary force, using different grid sizes and constitutive models, and obtaining accurate earthquake motion response law.

[0004] The second object of the present invention is to provide a seismic simulation device based on earthquake source-propagation medium-sedimentary layer.

[0005] A third objective of the present invention is to provide an electronic device.

[0006] A fourth object of the present invention is to provide a computer-readable storage medium.

[0007] To achieve the above-mentioned purpose, the first embodiment of the present invention proposes a method for simulating seismic motion based on source-propagation medium-sedimentary layer, the method comprising determining a large-area model and a local-area model for seismic motion simulation; wherein the large-area model comprises the local-area model; simulating the source rupture process of the large-area model through a finite fault model; simulating the propagation process of the seismic wave of the large-area model in the first propagation medium through a viscoelastic constitutive model; wherein the seismic wave is generated by the source rupture process, and the first propagation medium is the bedrock of the large-area model; determining the position of the local-area model corresponding to the large-area model, obtaining the seismic acceleration at the position, and calculating the equivalent boundary force according to the acceleration; inputting the equivalent boundary force into the local-area model, generating a transmission seismic wave corresponding to the equivalent boundary force in the local-area model; simulating the propagation process of the transmission seismic wave of the local-area model in the second propagation medium through a nonlinear model, and determining the seismic response during the propagation process; wherein the second propagation medium is the sedimentary layer of the local-area model.

[0008] In addition, the earthquake motion simulation method based on earthquake source-propagation medium-sedimentary layer according to the above embodiment of the present invention may also have the following additional technical features:

[0009] As an optional embodiment, the large area of ​​the large area model, the first grid size of the large area model, the first propagation medium parameter of the large area model and the first simulation frequency for the large area model are determined, the large area model is constructed according to the large area, the first grid size and the first propagation medium parameter, and the source rupture process of the large area model is simulated through a finite fault model according to the first simulation frequency; wherein the first propagation medium parameter is the medium parameter of the bedrock.

[0010] As an optional embodiment, the local area area of ​​the local area model, the second grid size of the local area model, the second propagation medium parameter of the local area model and the second simulation frequency for the local area model are determined, the local area model is constructed according to the local area area, the second grid size and the second propagation medium parameter, and the propagation process of the transmitted seismic fluctuations of the local area model in the second propagation medium is simulated through a nonlinear model according to the second simulation frequency; wherein the second propagation medium parameter is the medium parameter of the sedimentary layer, the second grid size is smaller than the first grid size, and the local area area is smaller than the large area.

[0011] As an optional embodiment, the first analog frequency and the second analog frequency are not less than 10 Hz.

[0012] As an optional embodiment, the first grid size is smaller than the wavelength of the seismic fluctuations of the large regional model.

[0013] As an optional embodiment, the second grid size is smaller than one tenth of the wavelength of the transmitted seismic wave of the local area model.

[0014] As an optional embodiment, the viscoelastic constitutive model and the nonlinear model each use an MPI-based parallel computing method during the simulation process.

[0015] According to the earthquake source-propagation medium-sedimentary layer earthquake motion simulation method based on the embodiment of the present invention, a large area model and a local area model for earthquake motion simulation are determined; wherein the large area model includes the local area model; the earthquake source rupture process of the large area model is simulated by a finite fault model; the propagation process of the earthquake wave of the large area model in the first propagation medium is simulated by a viscoelastic constitutive model; wherein the earthquake wave is generated by the earthquake source rupture process, and the first propagation medium is the bedrock of the large area model; the position of the local area model corresponding to the large area model is determined, the earthquake motion acceleration at the position is obtained, and the equivalent boundary force is calculated according to the acceleration; the equivalent boundary force is input into the local area model, and the transmission earthquake wave corresponding to the equivalent boundary force is generated in the local area model; the transmission process of the transmission earthquake wave of the local area model in the second propagation medium is simulated by a nonlinear model, and the earthquake motion response during the transmission process is determined; wherein the second propagation medium is the sedimentary layer of the local area model. The present invention can accurately simulate the entire process from earthquake source rupture to medium propagation and then to local site response of sedimentary layer. By connecting the large-area viscoelastic model and the local site nonlinear model through equivalent boundary force and adopting different grid sizes and constitutive models, it can realize wide-band strong vibration simulation of the entire process of earthquake source-propagation medium-sedimentary layer and obtain accurate earthquake motion response law.

[0016] To achieve the above-mentioned purpose, the second embodiment of the present invention proposes a seismic motion simulation device based on source-propagation medium-sedimentary layer, the device comprising a first determination module, configured to determine a large-area model and a local-area model for seismic motion simulation; wherein the large-area model comprises a local-area model; a first simulation module, configured to simulate the source rupture process of the large-area model through a finite fault model; a second simulation module, configured to simulate the propagation process of the seismic wave of the large-area model in the first propagation medium through a viscoelastic constitutive model; wherein the seismic wave is generated by the source rupture process, and the first propagation medium is the bedrock of the large-area model; a second determination module, configured to determine the position of the local-area model corresponding to the large-area model, obtain the seismic acceleration at the position, and calculate the equivalent boundary force according to the acceleration; a transmission module, configured to input the equivalent boundary force into the local-area model, and generate a transmission seismic wave corresponding to the equivalent boundary force in the local-area model; a third simulation module, configured to simulate the propagation process of the transmission seismic wave of the local-area model in the second propagation medium through a nonlinear model, and determine the seismic response during the propagation process; wherein the second propagation medium is the sedimentary layer of the local-area model.

[0017] According to an embodiment of the present invention, a seismic motion simulation device based on a source-propagation medium-sedimentary layer is provided, by determining a large-area model and a local-area model for seismic motion simulation; wherein the large-area model includes the local-area model; the source rupture process of the large-area model is simulated by a finite fault model; the propagation process of the seismic fluctuations of the large-area model in the first propagation medium is simulated by a viscoelastic constitutive model; wherein the seismic fluctuations are generated by the source rupture process, and the first propagation medium is the bedrock of the large-area model; the position of the local-area model corresponding to the large-area model is determined, the seismic acceleration at the position is obtained, and the equivalent boundary force is calculated according to the acceleration; the equivalent boundary force is input into the local-area model, and a transmission seismic fluctuation corresponding to the equivalent boundary force is generated in the local-area model; the transmission process of the transmission seismic fluctuations of the local-area model in the second propagation medium is simulated by a nonlinear model, and the seismic motion response during the transmission process is determined; wherein the second propagation medium is the sedimentary layer of the local-area model. The present invention can accurately simulate the entire process from earthquake source rupture to medium propagation and then to local site response of sedimentary layer. By connecting the large-area viscoelastic model and the local site nonlinear model through equivalent boundary force and adopting different grid sizes and constitutive models, it can realize wide-band strong vibration simulation of the entire process of earthquake source-propagation medium-sedimentary layer and obtain accurate earthquake motion response law.

[0018] To achieve the above-mentioned purpose, the third aspect of the present invention proposes an electronic device, including a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that when the processor executes the program, the seismic motion simulation method based on source-propagation medium-sedimentation layer as described above is implemented.

[0019] In order to achieve the above-mentioned purpose, the fourth aspect of the present invention proposes a computer-readable storage medium, which stores computer instructions, and the computer instructions are used to enable the computer to execute the above-mentioned seismic motion simulation method based on source-propagation medium-sedimentation layer.

[0020] Additional aspects and advantages of the present invention will be given in part in the following description and in part will be obvious from the following description, or will be learned through practice of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0021] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only for the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0022] Figure 1 A schematic flow chart of a method for simulating earthquake motion based on a source-propagation medium-sedimentary layer provided in an embodiment of the present invention.

[0023] Figure 2 A schematic diagram of the framework of a method for simulating earthquake motion based on a source-propagation medium-sedimentary layer provided in an embodiment of the present invention.

[0024] Figure 3 A schematic diagram of a seismic simulation device based on a seismic source-propagation medium-sedimentary layer provided in an embodiment of the present invention.

[0025] Figure 4 A schematic diagram of the structure of an electronic device provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0026] In order to make the objectives, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with specific embodiments and with reference to the accompanying drawings.

[0027] It should be noted that, unless otherwise defined, the technical terms or scientific terms used in the present invention should be understood by people with ordinary skills in the field to which the present invention belongs. The "first", "second" and similar words used in the present invention do not indicate any order, quantity or importance, but are only used to distinguish different components. "Include" or "comprise" and similar words mean that the elements or objects appearing before the word include the elements or objects listed after the word and their equivalents, without excluding other elements or objects. "Connect" or "connected" and similar words are not limited to physical or mechanical connections, but may include electrical connections, whether direct or indirect. "Up", "down", "left", "right" and the like are only used to indicate relative positional relationships. When the absolute position of the described object changes, the relative positional relationship may also change accordingly.

[0028] In the field of earthquake science and engineering, understanding the propagation laws of earthquake motion is the basis for assessing earthquake risks, designing earthquake-resistant structures, and formulating emergency response strategies. Currently, there are three main methods to understand the propagation laws of earthquake motion: measured data, analytical solutions, and numerical simulations.

[0029] Measured data: Although the seismic monitoring network is expanding globally, the number of measured stations is still limited, and the data of large and rare earthquakes that can be recorded is extremely scarce. This makes it difficult to fully and deeply understand the response laws of ground motion under different terrain and geological conditions. In particular, under complex terrain and geological conditions, the lack of measured data limits the in-depth understanding of the propagation characteristics of ground motion.

[0030] Analytical solution: The analytical solution method relies on a highly simplified topographic and geological model, and solves the propagation problem of seismic motion through mathematical formulas and theoretical derivations. However, the topographic and geological conditions in real life are extremely complex, and the geological conditions of the sedimentary layer are highly variable. This severely limits the analytical solution method in practical applications and cannot accurately reflect the propagation law of seismic motion under complex topographic and geological conditions.

[0031] Numerical simulation: In view of the limitations of measured data and analytical solution methods, numerical simulation methods have gradually become an important means of analyzing the propagation laws of seismic motion. Numerical simulation methods can simulate the propagation process of seismic waves under complex terrain and geological conditions, including source rupture, propagation of seismic waves in the propagation medium, and interaction between seismic waves and sedimentary layers. Through numerical simulation, we can have a deeper understanding of the propagation characteristics of seismic motion and provide a scientific basis for earthquake risk assessment, seismic structure design, and the formulation of emergency response strategies.

[0032] However, current methods for numerical simulation of earthquake motion still face many challenges. Especially in large-area simulations, due to the large range of earthquake-causing faults, large model size, and limited computing power, the grid division is often coarse. This results in the inability to accurately characterize the strong variability of the sedimentary layer when simulating the propagation of local earthquake motion, thus affecting the accuracy and reliability of the simulation results. Therefore, this patent aims to propose an innovative earthquake motion simulation method to solve the problems existing in the current numerical simulation methods and improve the accuracy and efficiency of earthquake motion simulation.

[0033] Specifically, one difficulty of the current earthquake motion simulation method is that the range of the earthquake fault is large, often tens of kilometers, and the simulation model for earthquake source rupture and medium propagation often needs to be greater than 2,500 square kilometers. The earthquake propagation simulation of the local site to be analyzed requires refined modeling to well characterize the strong variability of the local sedimentary layer. The simulation model of the large area is large in size and the grid division is rough, otherwise the computing power is difficult to meet the requirements. If the grid size is too large, the simulation frequency of the local sedimentary layer will be greatly reduced, which cannot meet the needs of the seismic safety analysis of the buildings on the site. Another difficulty in the strong vibration simulation of the whole process of earthquake source-propagation medium-sedimentary layer is the selection of the constitutive model. The earthquake fault and the medium propagation of the large area are often on the bedrock. The vibration simulation of the bedrock can be simulated by viscoelastic simulation, which has high computational efficiency while ensuring the simulation accuracy. The local sedimentary layer is often a soil layer with strong nonlinearity. The viscoelastic model cannot simulate the response of the soil layer under the action of strong earthquakes, and a nonlinear model needs to be used for simulation. Therefore, it is necessary to establish a strong vibration simulation method that integrates the large-area viscoelastic model and the local site nonlinear model.

[0034] The technical solution of the present invention is further described in detail below through specific embodiments.

[0035] refer to Figure 1 , which is a schematic flow chart of a method for simulating earthquake motion based on source-propagation medium-sedimentary layer provided in an embodiment of the present invention.

[0036] Step S101, determining a large-area model and a local-area model for earthquake simulation; wherein the large-area model includes the local-area model.

[0037] Specifically, a large regional model can be determined based on the research purpose and geographical scope, which covers a wide area that may be affected by an earthquake. Within the large regional model, one or more local regional models are determined based on specific research interests or concerns (such as sedimentary layers). These local regional models are usually located in areas where the seismic response may be more significant or complex.

[0038] As an optional embodiment, the large area of ​​the large area model, the first grid size of the large area model, the first propagation medium parameter of the large area model and the first simulation frequency for the large area model can be determined first, and the large area model can be constructed according to the large area, the first grid size and the first propagation medium parameter. The source rupture process of the large area model can be simulated through a finite fault model according to the first simulation frequency; wherein the first propagation medium parameter is the medium parameter of the bedrock.

[0039] In a specific implementation, a large area is determined according to the research purpose and the potential impact range of seismic activity. This area should be large enough to include all key areas that can be affected by earthquakes (for example, greater than 2,500 square kilometers, covering the main areas of potential earthquake impact). Considering the balance between computing power and simulation accuracy, the first grid size of the large area model is determined. The selection of the grid size should ensure that the simulation is completed within a reasonable time frame while maintaining sufficient accuracy to capture the propagation characteristics of seismic waves. The first propagation medium parameter describes the physical properties of the bedrock, such as density, elastic modulus, Poisson's ratio, etc. These parameters are crucial for simulating the propagation of seismic waves in the bedrock. In this embodiment, the first propagation medium parameter specifically refers to the medium parameter of the bedrock. According to the research purpose and the required simulation accuracy, a first simulation frequency is determined. This frequency will determine the choice of time step in the simulation process, thereby affecting the accuracy and computational efficiency of the simulation. Further, a large area model is constructed based on the determined large area area, the first grid size, and the first propagation medium parameter. This model will be used to simulate the propagation process of seismic waves in the bedrock. In the constructed large area model, a finite fault model is used to simulate the rupture process of the earthquake source. The finite fault model can take into account factors such as the geometry, rupture velocity and direction of the earthquake fault, thereby more realistically reflecting the initiation and propagation process of the earthquake. According to the determined first simulation frequency, the time step in the simulation process is set, and the finite fault model is started to simulate the rupture of the earthquake source. During the simulation process, the seismic waves will propagate in the bedrock and be affected by the physical properties of the bedrock.

[0040] It should be noted that when determining the area of ​​a large region, the potential impact range of seismic activity should be fully considered to ensure the accuracy of the simulation results. When selecting the first grid size, it is necessary to weigh the relationship between computing power and simulation accuracy. A grid size that is too large may lead to a decrease in simulation accuracy, while a grid size that is too small may increase the computational burden. When determining the first propagation medium parameters, it should be based on reliable geological surveys and laboratory test results to ensure the reliability of the simulation results. When setting the first simulation frequency, it is necessary to consider the balance between the required time resolution and computational efficiency. A frequency that is too high may result in a long calculation time, while a frequency that is too low may not capture some important characteristics of seismic waves.

[0041] As an optional embodiment, the local area area of ​​the local area model, the second grid size of the local area model, the second propagation medium parameter of the local area model and the second simulation frequency for the local area model can be determined, and the local area model can be constructed according to the local area area, the second grid size and the second propagation medium parameter. The propagation process of the transmitted seismic fluctuations of the local area model in the second propagation medium is simulated through a nonlinear model according to the second simulation frequency; wherein the second propagation medium parameter is the medium parameter of the sedimentary layer, the second grid size is smaller than the first grid size, and the local area area is smaller than the large area.

[0042] In a specific implementation, a local area area is determined according to the research purpose and specific focus. This area should include areas where the seismic response can be more significant or complex, and is smaller than the area of ​​the large area model. Considering that the local area model requires higher simulation accuracy, a second grid size smaller than the first grid size of the large area model is determined. This helps to capture subtle changes in seismic waves in the sedimentary layer. The second propagation medium parameter describes the physical properties of the sedimentary layer, such as density, shear modulus, damping ratio, etc. These parameters are crucial for simulating the propagation of seismic waves in the sedimentary layer. In this embodiment, the second propagation medium parameter specifically refers to the medium parameter of the sedimentary layer. According to the research purpose and the required simulation accuracy, a second simulation frequency is determined. Since the local area model requires higher accuracy, the second simulation frequency can be higher than the first simulation frequency of the large area model. According to the determined local area area, the second grid size and the second propagation medium parameter, a local area model is constructed. This model will be used to simulate the propagation process of seismic waves in the sedimentary layer. Further, in the constructed local area model, a nonlinear model is used to simulate the propagation process of seismic waves in the sedimentary layer. The nonlinear model can take into account the complex physical characteristics and nonlinear behavior of the sedimentary layer, thereby more accurately describing the propagation and attenuation characteristics of seismic waves in the sedimentary layer. According to the determined second simulation frequency, the time step in the simulation process is set, and the nonlinear model is started to simulate the propagation of seismic waves in the sedimentary layer. During the simulation process, the seismic waves will be affected by the physical characteristics of the sedimentary layer and produce corresponding ground motion responses.

[0043] It should be noted that when determining the area of ​​the local area, the research purpose and specific concerns should be fully considered to ensure the pertinence and accuracy of the simulation results. When selecting the second grid size, it is necessary to weigh the relationship between simulation accuracy and computational burden. A grid size that is too small may lead to excessive computational burden, while a grid size that is too large may reduce simulation accuracy. When determining the parameters of the second propagation medium, it should be based on reliable geological surveys and laboratory test results to ensure the reliability of the simulation results. At the same time, considering the complexity and variability of the sedimentary layer, a more sophisticated parameterization method may be required to describe its physical properties. When setting the second simulation frequency, it is necessary to consider the balance between the required time resolution and computational efficiency. A frequency that is too high may result in a long calculation time, while a frequency that is too low may not capture some important characteristics of seismic waves.

[0044] As an optional embodiment, the first analog frequency and the second analog frequency are not less than 10 Hz.

[0045] As an optional embodiment, the first grid size is smaller than the wavelength of the seismic fluctuations of the large regional model.

[0046] As an optional embodiment, the second grid size is smaller than one tenth of the wavelength of the transmitted seismic wave of the local area model.

[0047] Specifically, the grid size is smaller than the wavelength of the seismic wave in the medium to ensure the accuracy and resolution of the simulation. The simulation frequency is higher than 10 Hz to capture the high-frequency components of the seismic wave, which is crucial for understanding the detailed characteristics of ground motion and the potential risk of earthquake disasters. Among them, when the spectral element method is used, the grid size is smaller than the wavelength of the seismic wave in the sedimentary layer to ensure the accuracy and resolution of the simulation. When the finite element method is used, the grid size is less than 1 / 10 of the wavelength to capture finer wave characteristics and nonlinear responses of the sedimentary layer.

[0048] Step S102, simulating the earthquake source rupture process of the large regional model by using a finite fault model.

[0049] The finite fault model is a physical model used to describe the rupture process of an earthquake fault. It takes into account factors such as the geometry of the fault, the rupture velocity, the rupture mode, and the slip distribution on the fault plane, and can more realistically reflect the initiation and expansion process of an earthquake.

[0050] Specifically, firstly, the key parameters such as the geometry (such as length, width, strike, etc.), position, and rupture velocity of the fault in the large regional model can be determined based on seismological research and geological survey data. Further, using these parameters, a finite fault model is constructed in the large regional model. The fault can usually be divided into multiple small blocks or grid points, and corresponding physical and mechanical parameters are assigned to each small block or grid point. Further, according to the principles of seismology and observational data, the boundary conditions and initial conditions of the model are set. These conditions may include the propagation direction, velocity, acceleration, etc. of seismic waves. Further, the model is run using appropriate numerical simulation methods (such as finite difference method, finite element method, or spectral element method, etc.) to simulate the rupture process of the earthquake fault. During the simulation process, key information such as the propagation path, amplitude, and frequency of the seismic wave is recorded and analyzed. Finally, the simulation results are analyzed in detail, including the waveform, amplitude, and frequency distribution characteristics of the seismic wave. At the same time, the accuracy and reliability of the simulation results need to be evaluated to verify the applicability of the model.

[0051] Step S103, simulating the propagation process of the seismic wave of the large-area model in the first propagation medium through the viscoelastic constitutive model; wherein the seismic wave is generated by the earthquake source rupture process, and the first propagation medium is the bedrock of the large-area model.

[0052] The viscoelastic constitutive model is a mathematical model that describes the material's behavior when subjected to external forces, showing both elastic and viscous properties. In seismology, the viscoelastic model is widely used to simulate the propagation of seismic waves in the strata because it can more realistically reflect the physical properties of the strata.

[0053] Specifically, seismic fluctuations are generated by the source rupture process. When an earthquake occurs, the fault in the crust suddenly ruptures, releasing huge energy and forming seismic waves. These seismic waves include longitudinal waves (P waves) and transverse waves (S waves), which propagate at different speeds in the crust. In the simulation, the seismic waves start from the source and first enter the bedrock part of the large regional model. The bedrock usually has relatively stable physical properties, such as high density, hardness and elastic modulus. Therefore, the propagation process of seismic waves in the bedrock is relatively simple, but the influence of viscoelastic effects on wave velocity and waveform still needs to be considered. In the simulation, the viscoelastic constitutive model needs to be constructed first. Usually, suitable viscoelastic elements (such as springs, dampers, etc.) and connection methods (such as series, parallel, etc.) can be selected to form a complex model that can reflect the physical properties of the formation. After the model is built, the corresponding parameters need to be set. These parameters include the density, elastic modulus, viscosity coefficient, etc. of the formation, which directly affect the propagation speed and waveform of seismic waves. In actual simulations, these parameters are usually determined based on geological survey data and geophysical exploration data. In order to simulate the propagation process of seismic waves in viscoelastic formations, appropriate numerical simulation methods are required. Numerical simulation methods include but are not limited to finite difference method, finite element method and spectral element method. These methods are able to handle complex geometric shapes and physical parameters and generate high-precision simulation results. After the simulation is completed, the results need to be analyzed in detail. This includes the extraction and analysis of the waveform, amplitude, frequency and other characteristics of the seismic waves, as well as the evaluation of wave velocity and attenuation characteristics. The accuracy and reliability of the model can be verified by comparison with actual observation data.

[0054] Step S104, determining the position of the local area model corresponding to the large area model, obtaining the earthquake acceleration at the position, and calculating the equivalent boundary force according to the acceleration.

[0055] Specifically, a geographic information system (GIS) or other spatial analysis tools can be used to determine the exact position of the local area model in the large area model. If the large area model and the local area model use different coordinate systems, coordinate transformation is required to ensure the consistency of the position. Define the boundary of the local area model according to the research purpose and accuracy requirements. This can be based on geological characteristics, terrain changes or other relevant factors. Further, the ground motion acceleration data can be obtained from the earthquake monitoring network, earthquake simulation results or historical earthquake records. If there is no direct ground motion acceleration observation data at the location of the local area model, an interpolation method (such as Kriging interpolation, linear interpolation, etc.) can be used to estimate the ground motion acceleration at that location. According to the principles of seismology and observation data, the acquired ground motion acceleration data is corrected to ensure its accuracy and reliability. Finally, the equivalent boundary force refers to the boundary force that causes the local area model to produce the same response as the actual ground motion under the action of an earthquake. It can be calculated by methods such as dynamic equations or finite element analysis. The calculation formula of the equivalent boundary force may vary depending on the specific physical model and numerical method. In general, it can be calculated by combining the ground motion acceleration with the mass, stiffness and other parameters of the local area model. For example, in finite element analysis, the following formula can be used to calculate the equivalent node force: F = -Ma (where F is the equivalent node force, M is the mass matrix, and a is the acceleration vector). Use finite element analysis or other numerical methods to calculate the equivalent boundary force. Usually, a numerical model of the local area model can be established, and the seismic acceleration is applied as an input condition. The calculation results are verified to ensure their accuracy and reliability. This can be achieved by comparison with actual observation data, model testing, or other methods.

[0056] Step S105, inputting the equivalent boundary force into the local area model, and generating the transmission seismic fluctuation corresponding to the equivalent boundary force in the local area model.

[0057] Specifically, a numerical model of the local area model can be constructed according to the simulation requirements and research purposes. It can include finite element models, discrete element models or other types of numerical models. Ensure that the geometry, material properties and boundary conditions of the model are consistent with the actual engineering situation, mesh the local area model to ensure the accuracy and efficiency of the numerical calculation, and select the appropriate mesh density and shape according to the complexity of the model and the computing resources. Further, according to the equivalent boundary force calculated previously, determine the magnitude and direction of the force to be input into the local area model. Ensure that the application position of the equivalent boundary force matches the boundary conditions of the local area model. Apply the equivalent boundary force on the boundary of the local area model. This can be achieved by applying concentrated force or distributed force on the boundary nodes of the numerical model. Ensure that the application method of the equivalent boundary force meets the requirements of physical principles and numerical methods. Finally, use appropriate numerical methods (such as finite element method, discrete element method, etc.) to run the local area model to simulate the propagation process of seismic waves. During the simulation process, record and analyze key information such as the propagation path, amplitude, and frequency of seismic waves. Conduct detailed analysis of the simulation results, including the extraction and analysis of seismic wave waveforms, amplitudes, frequencies, etc. Verify the accuracy and reliability of the simulation results by comparing them with actual observation data or other reliable information.

[0058] Step S106, simulating the propagation process of the transmitted seismic wave of the local area model in the second propagation medium through a nonlinear model, and determining the seismic response during the propagation process; wherein the second propagation medium is the sedimentary layer of the local area model.

[0059] Specifically, a three-dimensional geological model of the local area can be established based on the geological structure, topography and other characteristics of the study area. The model should include sedimentary layers and other possible rock layers or geological structures. A physical model that can describe the nonlinear characteristics of the sedimentary layer is used, such as a constitutive model based on the nonlinear characteristics of the soil (such as the Martin-Seed-Davidenkov model). These models can consider the stress-strain relationship of the sedimentary layer under the action of seismic waves, thereby more accurately simulating the propagation process of seismic waves. Further, according to the historical seismic data or potential earthquake threats in the study area, appropriate seismic waves are selected as input. Seismic waves should include components of different frequencies and amplitudes to simulate the complexity of real earthquakes. The selected seismic waves are applied to the boundaries of the local area model in an appropriate manner, such as the bottom boundary or the side boundary. At the same time, factors such as the incident angle and polarization direction of the seismic wave should be considered. Further, numerical methods such as spectral element method, finite element method or finite difference method can be used, but not limited to, to discretize the nonlinear model and solve the propagation process of seismic waves in the sedimentary layer, which can handle complex nonlinear problems and provide high-precision simulation results. The parameters of the model are set according to the physical properties (such as density, elastic modulus, Poisson's ratio, etc.) and geometric properties (such as thickness, inclination, etc.) of the sedimentary layer. These parameters will directly affect the propagation process of seismic waves and the ground motion response. Finally, by analyzing the simulation results, the ground motion response of seismic waves in the process of propagation in the sedimentary layer is obtained, including parameters such as acceleration, velocity, and displacement. These parameters can reflect the impact of seismic waves on the sedimentary layer and the amplification or attenuation of seismic waves by the sedimentary layer. The impact of the physical and geometric properties of the sedimentary layer on the ground motion response can also be further analyzed. For example, factors such as the thickness, density, and elastic modulus of the sedimentary layer may affect the propagation speed and amplitude of the seismic wave. At the same time, factors such as the inclination and impedance ratio of the sedimentary layer may also have a significant impact on the ground motion response.

[0060] The simulation results can be applied to engineering seismic design to provide reasonable seismic parameters for buildings, bridges and other structures. These parameters can reflect the actual response of the study area under the action of earthquakes, thereby ensuring the safety and stability of the structure. The simulation results can also be used to assess the earthquake disaster risk in the study area. By analyzing the intensity and distribution characteristics of the seismic response, the degree of damage and impact range that may be caused by an earthquake can be predicted, providing a scientific basis for disaster prevention and mitigation.

[0061] As an optional embodiment, the viscoelastic constitutive model and the nonlinear model each use an MPI-based parallel computing method during the simulation process.

[0062] MPI is a standard interface and message passing protocol for distributed or parallel computing. It allows multiple processes to run in parallel on multiple computers or multi-core processors on a single computer and communicate over a network or shared memory. MPI provides a wealth of features such as process management, message passing, synchronization, and error handling, which can support the efficient execution of large-scale parallel computing tasks.

[0063] When simulating the propagation process of seismic waves in local area models, especially in complex geological structures containing sedimentary layers, the use of viscoelastic constitutive models and nonlinear models, combined with parallel computing methods based on MPI (Message Passing Interface), can significantly improve computing efficiency and simulation accuracy.

[0064] refer to Figure 2 , which is a schematic diagram of the framework of the earthquake motion simulation method based on source-propagation medium-sedimentary layer provided in an embodiment of the present invention.

[0065] Specifically, in step S201, after starting the earthquake wave simulation, a finite fault model is used to simulate the source rupture process, and the earthquake force is equivalent to a series of earthquake moment tensors distributed along time and space. The earthquake moment tensors are applied to the grid of the spectral element method, and the imbalance of force will produce earthquake waves.

[0066] In step S202, the generated seismic waves will propagate to a distant place along the propagation medium. When simulating the propagation process of seismic motion in the medium, a viscoelastic model is used to simulate the propagation medium, and the spectral element method is used for simulation. The grid size of the model is less than one wavelength. The parallel computing method based on MPI is used to improve the computing efficiency and achieve a simulation frequency of more than 10 Hz.

[0067] Step S203, when a calculation step is completed, extract the earthquake acceleration at the position of the local site model boundary corresponding to the large regional model, calculate the equivalent boundary force, and input the equivalent boundary force into the outer boundary of the local site model, which will generate earthquake fluctuations in the local site model.

[0068] Step S204, the seismic fluctuations generated in the local site model will propagate in the local sedimentary layer. A nonlinear model is used to simulate the seismic response of the sedimentary layer in the simulation, and the spectral element method or finite element software is used. When the spectral element method is used, the grid size is smaller than the wavelength, and when the finite element method is used, the grid size is smaller than 1 / 10 of the wavelength. The MPI-based parallel computing method is used to improve the computing efficiency and achieve a simulation frequency of more than 10 Hz.

[0069] Step S205, after obtaining the seismic response of the local site, a judgment is made according to the model calculation step requirements. If the simulation calculation step requirements have been met, the calculation is terminated, otherwise the calculation of steps S201 to S204 is repeated.

[0070] In summary, the earthquake motion simulation method based on earthquake source-propagation medium-sedimentary layer provided by the present invention carries out earthquake source rupture simulation and earthquake wave propagation process simulation in the propagation medium in a large area (greater than 2500 square kilometers), adopts viscoelastic constitutive model, and simulates using spectral element method software, and the simulation frequency is higher than 10Hz; in the local site (less than 20 square kilometers), the earthquake wave simulation in the sedimentary layer is carried out, and the nonlinear model is adopted, and the spectral element method or finite element software is used for simulation, and the simulation frequency is higher than 10Hz. The equivalent boundary force is used as a connection channel between the large-area viscoelastic model and the local site nonlinear model, and the accurate simulation of earthquake motion in the large area and the local site is achieved through the transmission of equivalent boundary force. Among them, the sedimentary layer adopts a nonlinear model, and the spectral element method or finite element software is used for simulation. When the spectral element method is used, the grid size is less than the wavelength, and when the finite element method is used, the grid size is less than 1 / 10 of the wavelength, and the simulation frequency is higher than 10Hz. The equivalent boundary force is used as a connection channel between the large-area viscoelastic model and the local site nonlinear model. Among them, the equivalent boundary force is calculated by the earthquake acceleration at the boundary of the local site model. The equivalent effect boundary is transferred between the large-area viscoelastic model and the local site nonlinear model, and is transferred from the large-area viscoelastic model to the local site nonlinear model as the transmitted seismic fluctuation of the input ground motion of the local site nonlinear model.

[0071] It should be noted that the equivalent boundary force transmitted by the present invention is dynamically stored, and only one equivalent boundary force is retained each time. The next calculation covers the result of the previous calculation, which greatly reduces the demand for storage capacity and improves the calculation speed.

[0072] As can be seen from the above, the earthquake motion simulation method based on source-propagation medium-sedimentary layer provided by the present invention is through determining a large area model and a local area model for earthquake motion simulation; wherein, the large area model includes the local area model; the earthquake source rupture process of the large area model is simulated by a finite fault model; the propagation process of the earthquake wave of the large area model in the first propagation medium is simulated by a viscoelastic constitutive model; wherein, the earthquake wave is generated by the earthquake source rupture process, and the first propagation medium is the bedrock of the large area model; the position of the local area model corresponding to the large area model is determined, the earthquake motion acceleration at the position is obtained, and the equivalent boundary force is calculated according to the acceleration; the equivalent boundary force is input into the local area model, and the transmission earthquake wave corresponding to the equivalent boundary force is generated in the local area model; the transmission process of the transmission earthquake wave of the local area model in the second propagation medium is simulated by a nonlinear model, and the earthquake motion response during the transmission process is determined; wherein, the second propagation medium is the sedimentary layer of the local area model. The present invention can accurately simulate the entire process from earthquake source rupture to medium propagation and then to local site response of sedimentary layer. By connecting the large-area viscoelastic model and the local site nonlinear model through equivalent boundary force and adopting different grid sizes and constitutive models, it can realize wide-band strong vibration simulation of the entire process of earthquake source-propagation medium-sedimentary layer and obtain accurate earthquake motion response law.

[0073] It should be noted that the method of the embodiment of the present invention can be performed by a single device, such as a computer or a server. The method of this embodiment can also be applied in a distributed scenario and completed by multiple devices cooperating with each other. In the case of such a distributed scenario, one of the multiple devices can only perform one or more steps in the method of the embodiment of the present invention, and the multiple devices will interact with each other to complete the described method.

[0074] It should be noted that some embodiments of the present invention are described above. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recorded in the claims can be performed in an order different from that in the above embodiments and still achieve the desired results. In addition, the processes depicted in the accompanying drawings do not necessarily require the specific order or continuous order shown to achieve the desired results. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0075] Based on the same inventive concept, corresponding to the method provided in any of the above embodiments, the present invention further provides a seismic simulation device based on source-propagation medium-sedimentary layer.

[0076] refer to Figure 3 , which is a schematic diagram of a seismic simulation device based on a source-propagation medium-sedimentary layer provided in an embodiment of the present invention.

[0077] The device comprises:

[0078] The first determination module 301 is configured to determine a large-area model and a local-area model for earthquake simulation; wherein the large-area model includes the local-area model;

[0079] The first simulation module 302 is configured to simulate the earthquake source rupture process of the large regional model through a finite fault model;

[0080] The second simulation module 303 is configured to simulate the propagation process of the seismic wave of the large-area model in the first propagation medium through the viscoelastic constitutive model; wherein the seismic wave is generated by the earthquake source rupture process, and the first propagation medium is the bedrock of the large-area model;

[0081] The second determination module 304 is configured to determine the position of the local area model corresponding to the large area model, obtain the earthquake acceleration at the position, and calculate the equivalent boundary force according to the acceleration;

[0082] The transfer module 305 is configured to input the equivalent boundary force into the local area model, and generate a transfer seismic wave corresponding to the equivalent boundary force in the local area model;

[0083] The third simulation module 306 is configured to simulate the propagation process of the transmitted seismic wave of the local area model in the second propagation medium through a nonlinear model, and determine the seismic response during the propagation process; wherein the second propagation medium is the sedimentary layer of the local area model.

[0084] According to an embodiment of the present invention, a seismic motion simulation device based on a source-propagation medium-sedimentary layer is provided, by determining a large-area model and a local-area model for seismic motion simulation; wherein the large-area model includes the local-area model; the source rupture process of the large-area model is simulated by a finite fault model; the propagation process of the seismic fluctuations of the large-area model in the first propagation medium is simulated by a viscoelastic constitutive model; wherein the seismic fluctuations are generated by the source rupture process, and the first propagation medium is the bedrock of the large-area model; the position of the local-area model corresponding to the large-area model is determined, the seismic acceleration at the position is obtained, and the equivalent boundary force is calculated according to the acceleration; the equivalent boundary force is input into the local-area model, and a transmission seismic fluctuation corresponding to the equivalent boundary force is generated in the local-area model; the transmission process of the transmission seismic fluctuations of the local-area model in the second propagation medium is simulated by a nonlinear model, and the seismic motion response during the transmission process is determined; wherein the second propagation medium is the sedimentary layer of the local-area model. The present invention can accurately simulate the entire process from earthquake source rupture to medium propagation and then to local site response of sedimentary layer. By connecting the large-area viscoelastic model and the local site nonlinear model through equivalent boundary force and adopting different grid sizes and constitutive models, it can realize wide-band strong vibration simulation of the entire process of earthquake source-propagation medium-sedimentary layer and obtain accurate earthquake motion response law.

[0085] For the convenience of description, the above system is described as being divided into various modules according to their functions. Of course, when implementing the present invention, the functions of each module can be implemented in the same or multiple software and / or hardware.

[0086] The system of the above embodiment is used to implement the corresponding method in any of the above embodiments, and has the beneficial effects of the corresponding method embodiment, which will not be repeated here.

[0087] Based on the same inventive concept, corresponding to the method described in any of the above embodiments, the present invention also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the method described in any of the above embodiments when executing the program.

[0088] Figure 4 A more specific schematic diagram of the hardware structure of an electronic device provided in this embodiment is shown, and the device may include: a processor 410, a memory 420, an input / output interface 430, a communication interface 440, and a bus 450. The processor 410, the memory 420, the input / output interface 430, and the communication interface 440 are connected to each other through the bus 450 in the device.

[0089] The processor 410 can be implemented by a general-purpose CPU (Central Processing Unit), a microprocessor, an application-specific integrated circuit (ASIC), or one or more integrated circuits, and is used to execute relevant programs to implement the technical solutions provided in the embodiments of this specification.

[0090] The memory 420 may be implemented in the form of ROM (Read Only Memory), RAM (Random Access Memory), static storage device, dynamic storage device, etc. The memory 420 may store an operating system and other application programs. When the technical solutions provided in the embodiments of this specification are implemented by software or firmware, the relevant program codes are stored in the memory 420 and are called and executed by the processor 410.

[0091] The input / output interface 430 is used to connect the input / output module to realize information input and output. The input / output module can be configured in the device as a component (not shown in the figure), or it can be externally connected to the device to provide corresponding functions. The input device may include a keyboard, a mouse, a touch screen, a microphone, various sensors, etc., and the output device may include a display, a speaker, a vibrator, an indicator light, etc.

[0092] The communication interface 440 is used to connect a communication module (not shown) to realize communication interaction between the device and other devices. The communication module can realize communication through a wired mode (such as USB, network cable, etc.) or a wireless mode (such as mobile network, WIFI, Bluetooth, etc.).

[0093] The bus 450 comprises a pathway for transmitting information between the various components of the device (eg, the processor 410, the memory 420, the input / output interface 430, and the communication interface 440).

[0094] It should be noted that, although the above device only shows the processor 410, the memory 420, the input / output interface 430, the communication interface 440 and the bus 450, in the specific implementation process, the device may also include other components necessary for normal operation. In addition, it can be understood by those skilled in the art that the above device may also only include the components necessary for implementing the embodiments of the present specification, and does not necessarily include all the components shown in the figure.

[0095] The electronic device of the above embodiment is used to implement the corresponding method in any of the above embodiments, and has the beneficial effects of the corresponding method embodiment, which will not be described in detail here.

[0096] Based on the same inventive concept, corresponding to the method described in any of the above embodiments, the present invention also provides a non-transitory computer-readable storage medium, wherein the non-transitory computer-readable storage medium stores computer instructions, and the computer instructions are used to enable the computer to execute the method described in any of the above embodiments.

[0097] The above-mentioned non-transitory computer-readable storage medium can be any available medium or data storage device that can be accessed by a computer, including but not limited to magnetic storage (such as floppy disks, hard disks, magnetic tapes, magneto-optical disks (MO), etc.), optical storage (such as CD, DVD, BD, HVD, etc.), and semiconductor storage (such as ROM, EPROM, EEPROM, non-volatile memory (NAND FLASH), solid-state drive (SSD)), etc.

[0098] The computer instructions stored in the storage medium of the above embodiment are used to enable the computer to execute the method described in any embodiment of the above exemplary method part, and have the beneficial effects of the corresponding method embodiment, which will not be repeated here.

[0099] In addition, although the operations of the method of the present invention are described in a specific order in the accompanying drawings, this does not require or imply that the operations must be performed in this specific order, or that all the operations shown must be performed to achieve the desired results. On the contrary, the steps depicted in the flow chart can be performed in a different order. Additionally or alternatively, some steps can be omitted, multiple steps can be combined into one step, and / or one step can be decomposed into multiple steps.

[0100] It should be understood that the various parts of the present invention can be implemented by hardware, software, firmware or a combination thereof. In the above-mentioned embodiments, a plurality of steps or methods can be implemented by software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented by hardware, as in another embodiment, it can be implemented by any one of the following technologies known in the art or their combination: a discrete logic circuit having a logic gate circuit for implementing a logic function for a data signal, a dedicated integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

[0101] It should be noted that, unless otherwise defined, the technical terms or scientific terms used in the embodiments of the present invention should be understood by people with ordinary skills in the field to which the present invention belongs. The "first", "second" and similar words used in the embodiments of the present invention do not indicate any order, quantity or importance, but are only used to distinguish different components. "Include" or "comprise" and similar words mean that the elements or objects appearing before the word include the elements or objects listed after the word and their equivalents, without excluding other elements or objects. "Connect" or "connected" and similar words are not limited to physical or mechanical connections, but may include electrical connections, whether direct or indirect. "Up", "down", "left", "right" and the like are only used to indicate relative positional relationships. When the absolute position of the described object changes, the relative positional relationship may also change accordingly.

[0102] Although the spirit and principle of the present invention have been described with reference to several specific embodiments, it should be understood that the present invention is not limited to the disclosed specific embodiments, and the division of various aspects does not mean that the features in these aspects cannot be combined to benefit, and this division is only for the convenience of expression. The present invention is intended to cover various modifications and equivalent arrangements included in the spirit and scope of the attached claims. The scope of the attached claims conforms to the broadest interpretation, thereby including all such modifications and equivalent structures and functions.

Claims

1. A method for simulating earthquake motion based on earthquake source-propagation medium-sedimentary layer, characterized in that: include: Determine a large-area model and a local-area model for earthquake simulation; wherein the large-area model includes the local-area model; Simulating the earthquake source rupture process of the large regional model through a finite fault model; The propagation process of the seismic wave of the large-area model in the first propagation medium is simulated by a viscoelastic constitutive model; wherein the seismic wave is generated by the earthquake source rupture process, and the first propagation medium is the bedrock of the large-area model; Determine the position of the local area model corresponding to the large area model, obtain the earthquake acceleration at the position, and calculate the equivalent boundary force according to the acceleration; Inputting the equivalent boundary force into the local area model, and generating a transmission seismic fluctuation corresponding to the equivalent boundary force in the local area model; The propagation process of the transmitted seismic wave of the local area model in the second propagation medium is simulated by a nonlinear model, and the seismic response during the propagation process is determined; wherein the second propagation medium is the sedimentary layer of the local area model.

2. The earthquake motion simulation method based on source-propagation medium-sedimentary layer according to claim 1 is characterized in that: Also includes: Determine the large area of ​​the large area model, the first grid size of the large area model, the first propagation medium parameter of the large area model and the first simulation frequency for the large area model, construct the large area model according to the large area, the first grid size and the first propagation medium parameter, and simulate the source rupture process of the large area model through a finite fault model according to the first simulation frequency; wherein the first propagation medium parameter is the medium parameter of the bedrock.

3. The earthquake motion simulation method based on source-propagation medium-sedimentary layer according to claim 2 is characterized in that: Also includes: Determine the local area area of ​​the local area model, the second grid size of the local area model, the second propagation medium parameter of the local area model and the second simulation frequency for the local area model, construct the local area model according to the local area area, the second grid size and the second propagation medium parameter, and simulate the propagation process of the transmitted seismic waves of the local area model in the second propagation medium through a nonlinear model according to the second simulation frequency; wherein the second propagation medium parameter is the medium parameter of the sedimentary layer, the second grid size is smaller than the first grid size, and the local area area is smaller than the large area.

4. The earthquake motion simulation method based on source-propagation medium-sedimentary layer according to claim 3 is characterized in that: The first analog frequency and the second analog frequency are not less than 10 Hz.

5. The earthquake motion simulation method based on source-propagation medium-sedimentary layer according to claim 4 is characterized in that: The first grid size is smaller than a wavelength of seismic fluctuations of the large regional model.

6. The earthquake motion simulation method based on source-propagation medium-sedimentary layer according to claim 5 is characterized in that: The second grid size is smaller than one tenth of a wavelength of transmitted seismic fluctuations of the local area model.

7. The earthquake motion simulation method based on source-propagation medium-sedimentary layer according to claim 6 is characterized in that: The viscoelastic constitutive model and the nonlinear model each use a parallel computing method based on MPI during the simulation process.

8. A seismic simulation device based on source-propagation medium-sedimentary layer, characterized in that: include: A first determination module is configured to determine a large-area model and a local-area model for earthquake simulation; wherein the large-area model includes the local-area model; A first simulation module is configured to simulate the earthquake source rupture process of the large regional model through a finite fault model; A second simulation module is configured to simulate the propagation process of the seismic wave of the large-area model in the first propagation medium through a viscoelastic constitutive model; wherein the seismic wave is generated by the earthquake source rupture process, and the first propagation medium is the bedrock of the large-area model; A second determination module is configured to determine a position of the local area model corresponding to the large area model, obtain the earthquake acceleration at the position, and calculate an equivalent boundary force according to the acceleration; A transfer module, configured to input the equivalent boundary force into the local area model, and generate a transfer seismic fluctuation corresponding to the equivalent boundary force in the local area model; The third simulation module is configured to simulate the propagation process of the transmitted seismic waves of the local area model in the second propagation medium through a nonlinear model, and determine the seismic response during the propagation process; wherein the second propagation medium is the sedimentary layer of the local area model.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the program, the earthquake motion simulation method based on source-propagation medium-sediment layer as described in any one of claims 1 to 7 is implemented.

10. A non-transitory computer-readable storage medium storing computer instructions, characterized in that: The computer instructions are used to enable the computer to implement the earthquake motion simulation method based on source-propagation medium-sedimentary layer as described in any one of claims 1 to 7.

Images