139 results about "Seismic wave propagation" patented technology

The invention relates to the technical field of coal mine water disaster prevention and control, and particularly discloses a hidden structure detection system for mine water disaster intelligent monitoring, early warning and prevention, which is used for solving the problems of low mine hidden structure detection precision, mine water disaster early warning lag and large seismic wave propagation path simulation error. Comprising a multi-source data comprehensive detection module, a real-time imaging module, a micro-seismic multi-scale intelligent analysis module, a seismic wave signal identification module, a regional difference analysis module and an intelligent data processing module. According to the invention, through fusion of multi-source data, real-time imaging, micro-seismic analysis, seismic wave identification and intelligent data processing technologies, precise detection of a mine hidden structure, high-precision prediction of mine earthquakes and water disasters and real-time risk early warning are realized, and timeliness and accuracy of mine water disaster prevention and control are improved.

The invention relates to the field of seismic state recognition, and provides a converter station vector seismic state recognition method and system based on distributed optical fibers, and the method comprises the steps: collecting multi-dimensional vibration signal data of a converter station, and obtaining seismic monitoring data through processing; performing phase demodulation processing, time difference analysis and space vector decomposition algorithm processing on the seismic monitoring data to obtain seismic wave propagation vector characteristic data, and performing characteristic integration and optimization to obtain a seismic wave comprehensive state characteristic data set; identifying and classifying the seismic wave comprehensive state feature data set to obtain a seismic wave form classification result; performing analysis based on the seismic wave form classification result to obtain a converter station seismic state evaluation result; and carrying out comparative analysis on the evaluation result and a threshold value, generating earthquake early warning information of the converter station, and providing a corresponding emergency response suggestion according to the earthquake early warning information. According to the invention, accurate identification and classification of seismic waves are realized, and the earthquake safety early warning capability of the converter station is improved.

The invention relates to the technical field of data visualization, in particular to an earthquake emergency information rapid visualization method and system. The method comprises the following steps: acquiring earthquake monitoring real-time data; based on earthquake monitoring real-time data grid division, earthquake region population density distribution and key facility positions are generated, and an earthquake region sensitivity matrix is generated; dynamically simulating a seismic wave propagation process based on the seismic monitoring real-time data and the seismic region sensitivity matrix, and generating fluctuation propagation time sequence data; and identifying epicenter topography and landform vector data of the earthquake monitoring real-time data, and respectively evaluating topographic relief and geological stability to obtain earthquake topographic relief index data and earthquake geological stability score data. According to the invention, rapid and accurate assessment and visualization of earthquake damage distribution are realized through a data visualization technology, so that the earthquake emergency response efficiency and the disaster rescue decision scientificity are improved.

The invention discloses a micro-seismic multichannel data processing method and system based on a DBSCAN algorithm, and belongs to the technical field of micro-seismic monitoring. According to the invention, through feature vector construction and optimization distance calculation, the clustering reliability is enhanced; a multi-dimensional feature vector Xi containing P-wave refinement arrival time ti, an amplitude Ai and space coordinates (xi, yi, zi) is constructed, distance calculation is performed by combining a normalized cosine distance with wave velocity constraint, a micro-seismic wave propagation rule is fully utilized, the discrimination capability of DBSCAN clustering on physical consistency is enhanced, and the robustness of a screening result is improved, so that noise interference can be effectively processed in a complex environment, and the accuracy and the robustness of the screening result are improved. And the stability and the reliability of a micro-seismic positioning result under different geological conditions are ensured.

The invention discloses a high-precision three-dimensional first-motion wave travel time chromatography velocity inversion method and device, and relates to the field of seismic exploration. A multi-template windward difference algorithm and an undulating surface narrow-band continuation strategy are adopted to determine three-dimensional global forward first-motion wave propagation travel time; ray tracing is carried out according to the three-dimensional global forward first-motion wave propagation travel time, and a first-motion wave propagation path and a Frechet differential quotient matrix are obtained; establishing a three-dimensional initial inversion velocity field constrained by micro-logging information; constructing a tomographic inversion equation based on the apparent slowness constraint and the regularization constraint according to the travel time residual error and the Frechet differential quotient matrix; and solving the tomographic inversion equation to obtain a tomographic inversion slowness solution, and converting the slowness solution into a velocity solution. The problem of non-uniform travel time delay of seismic wave propagation of a complex near-surface structure can be solved, a more accurate superficial layer velocity structure can be obtained, and the quality of seismic data processing and imaging inversion can be improved.

The invention provides an ocean seismic signal processing method based on distributed optical fiber vibration sensing, and belongs to the technical field of seismic signal processing. The method comprises the following steps: acquiring a marine seismic signal through a distributed optical fiber vibration sensor; performing data preprocessing on the collected seismic signals to construct a seismic oscillogram data set; according to the preprocessed data, feature extraction is carried out on the preprocessed seismic signals through a U-KAN network to obtain feature data; identifying and classifying the feature data to obtain a P wave and an S wave; and calculating the time difference between the P wave and the S wave according to the identified P wave and S wave, and estimating the seismic source distance and magnitude according to the time difference and the seismic wave propagation speed. The earthquake signal processing automation degree is improved, the earthquake wave arrival time is accurately determined, the marine earthquake monitoring response speed and accuracy are improved, and a scientific basis is provided for prevention and reduction of marine earthquake disasters.

The invention provides a seismic source positioning method and system based on three-dimensional seismic reflection data and electronic equipment, and relates to the technical field of seismic exploration, and the method comprises the steps: carrying out the initial positioning of a seismic source according to the three-dimensional geological reflection data of a target research region, and obtaining a first seismic source position; extracting a reflected wave time sequence array mapped to the detection point array from the three-dimensional geological reflection data; taking the first seismic source position as an optimization starting point, and performing correlation analysis on the three-dimensional geological reflection data and the reflected wave time sequence array to obtain seismic source space-time parameters; pre-constructing a three-dimensional seismic wave velocity model; and carrying out inversion iteration by taking the seismic source space-time parameters as initial conditions and taking the three-dimensional seismic wave velocity model as iteration constraints, and determining the position of a target seismic source. The technical problem of low positioning precision under a complex geological structure due to insufficient consideration of seismic wave propagation characteristics and reflection characteristics in the prior art is solved, and the technical effect of improving the seismic source positioning precision is achieved.

The invention discloses a seismic wave propagation numerical simulation method and a seismic wave propagation numerical simulation system suitable for a tunnel peripheral frozen soil layer, and relates to the field of seismic wave simulation propagation. According to the method, a layered model is constructed to describe physical property distribution of the frozen soil layer at different depths, and a viscoelastic model is combined to integrate the layered model into a finite element model of the tunnel peripheral frozen soil layer; dividing the finite element model by adopting an initialized uniform grid, then preliminarily simulating a seismic wave propagation process, and by comparing the spatial change rate of the wave velocity, identifying a region with a remarkably changed wave velocity gradient, and extracting the wave velocity change amplitude and directivity data to determine the boundary and shape of the region; and planning the density distribution of the updated grid in combination with the geological parameters of the gradient change region. According to the method, the boundary and internal characteristics of the wave velocity gradient region are identified, the grid shape and density are planned according to the spatial change rate and distribution characteristics of the gradient, it is ensured that the calculation region can reflect the multi-scale characteristics of the frozen soil layer, and the seismic wave propagation path and the numerical calculation of energy distribution can be optimized.

The invention relates to a BIM (Building Information Modeling)-based building anti-seismic information determination method and device and a computer. The method comprises the following steps: acquiring regional geological data, regional historical seismic data and a building structure information model corresponding to a target building; according to the regional geological data and the regional historical seismic data, performing seismic wave propagation path simulation on a building site corresponding to the target building to obtain building seismic wave field information; performing seismic excitation coupling simulation according to the building seismic wave field information and the building structure information model to obtain component vibration response information; according to the component vibration response information, performing damage prediction analysis on each component in the target building to obtain damage information of each predicted component; and according to the predicted member damage information and the building structure information model, analyzing the structure safety of the target building to obtain building anti-seismic information corresponding to the target building. By means of the method, refined and personalized safety design can be achieved for specific buildings and complex sites.

The invention discloses a non-natural earthquake phase identification method based on image processing, and relates to the technical field of image processing, and the method comprises the following steps: capturing a continuous signal with the amplitude changing with time in a seismic wave propagation process in real time based on a sensor distributed by an earthquake monitoring station network, forming an original amplitude sequence, and carrying out the calculation of the original amplitude sequence; collected continuous signals are converted into three-channel images, and the image brightness is automatically adjusted according to the seismic wave data dynamic range. According to the method, a seismic wave multi-dimensional analysis system is constructed through a three-channel construction strategy, the first channel clearly shows the time domain difference between natural seismic progressive fluctuation and non-natural seismic pulse type sudden change, and the second channel accurately captures microsecond-level sudden change which is easy to neglect in traditional time domain analysis; the third channel clearly distinguishes natural earthquake low-frequency energy concentration and non-natural earthquake intermediate-frequency characteristic spectral lines, and the adaptive brightness adjustment technology aims at the intensity difference of multi-source earthquake signals, improves the distinguishing degree of non-natural and natural earthquake edges, and builds a recognition foundation firmly.

The invention discloses a dynamic polarization window seismic waveform separation method based on a time-frequency ridge line, which utilizes a time-frequency ridge line analysis technology and combines a dynamic polarization window calculation method to accurately calculate the polarization characteristic of a seismic signal. According to the method, the polarization information in the seismic wave propagation process is captured by analyzing the dynamic change of the seismic wave in the time-frequency domain, so that the seismic phase is distinguished, and the purpose of waveform separation is achieved. And more accurate data support and analysis means are provided for the application fields of seismic exploration, seismic early warning, seismic disaster assessment and the like.

The invention discloses a thermal-fluid-solid medium seismic wave simulation method, device, program and medium, and relates to the technical field of oil and gas resource geophysical exploration, and the method comprises the following steps: constructing a first-order thermal-fluid-solid coupling wave equation in an equation coupling form; constructing a coupling wave equation staggered network finite difference discrete format; analyzing numerical stability conditions; constructing a complete matching layer absorption boundary; and performing forward modeling simulation of seismic waves in the three-layer heat-fluid-solid coupling model according to the model parameters and the seismic source parameters. The method breaks through the assumption that the classical seismic wave propagation theory is completely elastic under the isothermal or adiabatic deformation condition, considers the influence of temperature change on the propagation of seismic waves in fluid and solid media, and improves the description precision of the propagation process of the seismic waves in the heat-fluid-solid media. And the capability of inverting the underground temperature field based on seismic wave data in the later period can be effectively improved.

The invention belongs to the field of seismic exploration of oil and gas resources, and provides a numerical simulation method for a seismic dispersion attenuation wave field containing fluid viscoelastic rock. According to the method, a reflection coefficient calculation method for accurately describing the energy conversion relation of seismic waves at the viscoelastic rock interface is established, the reflection time difference caused by seismic wave speed dispersion can be reflected, and the dispersion attenuation characteristics of the seismic waves propagating in viscoelastic stratum rocks are completely described; the method can be used for numerical simulation of seismic wave field propagation characteristics of various different complex rocks, and an intuitive and efficient implementation way is provided for researching seismic wave propagation rules of complex viscoelastic media.

The invention provides a compressed air energy storage chamber optimization design method, and relates to the technical field of energy storage chamber optimization design, and the method comprises the steps: collecting surrounding rock stress data, seismic wave propagation speed, sound wave attenuation coefficient and other information, building a multi-dimensional stress distribution model, and precisely recognizing a stress concentration region and a change trend thereof; through iterative simulation and optimization of the process, parameters such as elastic modulus, wave velocity ratio and fracture density are comprehensively considered, the basic shape and the supporting structure of the chamber are optimized, and the damage risk caused by stress concentration is remarkably reduced; according to the method, the safety and stability of the chamber are improved, the service life is prolonged, the bearing capacity is improved, and the large-scale energy storage application requirement is met.

Methods and systems are disclosed. The method may include obtaining an anisotropic velocity model for a subterranean region of interest, where the anisotropic velocity model represents a propagation velocity of seismic waves discretized on a first grid of nodes representing the subterranean region of interest and initiating an initial anisotropic traveltime for each node of a second grid representing the subterranean region of interest, where the anisotropic traveltime comprises a seismic traveltime from a source location. The method may further include forming a computational system, comprising a discretization of a traveltime equation for the anisotropic velocity model, and determining an updated anisotropic traveltime for each node on the second grid based, at least in part, on a parallel fast sweeping Cuthill-McKee ordering solution to the computational system and the initial anisotropic traveltime for each node.

The invention discloses a robust full-waveform inversion method of a first-arrival constraint physical embedded recurrent neural network, which comprises the following steps of: solving a wave equation through forward propagation of the physical embedded recurrent neural network, generating a forward shot gather record, screening out a seismic trace with high matching degree with observation data by using first-arrival information, and participating in loss calculation to obtain a seismic trace with high matching degree with observation data; and along with the proceeding of the inversion process, the range of seismic traces participating in the inversion is gradually expanded, it is ensured that the inversion result gradually tends to a globally optimal solution, back propagation is carried out by using space-time field gradient information recorded by forward propagation of the recurrent neural network, the speed parameter is further optimized, and the inversion precision is improved. According to the method, the physical process of seismic wave propagation is embedded into the recurrent neural network, the first arrival constraint and the progressive inversion strategy are introduced, the dependence of the FWI on the initial model is remarkably reduced, and the underground velocity field structure can still be accurately reconstructed even if the initial model and a real model have remarkable deviation.

The invention discloses a seismic wave velocity change analysis method under a hydraulic fracturing condition, and belongs to the technical field of seismic wave propagation analysis. The method comprises the following steps: S1, acquiring data by using a plurality of seismic stations and removing transient signals by using a long and short window method; s2, calculating a cross-correlation function between the two stations; s3, extracting a tail wave signal and estimating relative velocity change by using mobile window cross-spectrum analysis; s4, extracting a group velocity frequency dispersion curve based on a cross-correlation function, and constructing a two-dimensional velocity model before and after fracturing; and S5, generating a speed image through inversion, and calculating a difference image to reveal speed change characteristics. According to the method, the tail wave sensitivity and the moving window cross-spectrum analysis are combined, the speed change and the anisotropy characteristics of the speed change are accurately quantified, the method is suitable for the complex geological environment, and a new interpretation framework is provided for the influence of hydraulic fracturing on the underground structure.

The invention discloses a seismic wave line analysis processing method based on a deep learning technology, and the method comprises the following steps: 1, inputting seismic data, and carrying out the normalization processing, so as to eliminate the magnitude difference between different seismic events, and make preparations for the input of a subsequent deep learning model; step 2, feature extraction and analysis are performed on the normalized seismic data by using a deep learning model for training a large number of seismic waveform data sets so as to realize classification and identification of different types of seismic events; step 3, associating the extracted features with a known earthquake directory event to determine detailed information of an earthquake, such as an epicentral position, an earthquake magnitude and the like; and step 4, based on the output of the deep learning model, classifying and identifying seismic events, and the method solves the problems that a traditional seismic waveform analysis method depends on expert experience and a signal processing technology, but the methods are often low in efficiency and difficult to process a large-scale data set.

The invention relates to the technical field of earthquake engineering, geotechnical engineering and geological disaster assessment, in particular to an earthquake landslide risk assessment method and system, and the method comprises the following steps: obtaining regional topographic data, and constructing a topographic geometric model based on the regional topographic data; on the basis of the terrain geometric model, identifying a specific landform unit in which energy convergence may exist in the region; for the specific landform unit, simulating a propagation path of a seismic wave in the terrain geometric model, and determining the energy influence degree of the specific landform unit on the seismic wave according to the change of the propagation path; according to the method, the influence of complex terrains on the local energy of seismic waves can be effectively considered by constructing a terrain geometric model, identifying an energy convergence landform unit, simulating seismic wave propagation and correcting seismic oscillation parameters, so that the influence of the complex terrains on the local energy of the seismic waves can be effectively considered.

The invention discloses a pre-stack waveform inversion method based on a viscoelastic layered medium, which relates to the field of geophysical exploration of petroleum, and comprises the following steps: carrying out forward modeling based on a viscoelastic reflection matrix method, and accurately simulating the propagation characteristics of seismic waves in the viscoelastic layered medium, thereby effectively improving the inversion prediction precision. According to the method, the elastic characteristic and the attenuation characteristic of the stratified medium are considered at the same time in the seismic wave propagation process, the P-wave speed, the S-wave speed, the density, the P-wave quality factor and the S-wave quality factor can be inverted at the same time, the precision of the inversion result is improved through fine simulation of the underground reservoir medium, and the method has the advantages of being high in practicability and the like. The multi-parameter inversion method has the advantages that strong support is provided for reservoir prediction, a Bayesian inversion framework is adopted, a multi-parameter inversion process is optimized by fusing prior information, robustness and noise immunity of the inversion method are improved, and stability and accuracy of the multi-parameter inversion method in actual seismic data application are guaranteed.

The invention provides a sound wave remote exploration data reconstruction technology and device based on high-fidelity data mapping, and the technology comprises the steps: converting sound wave remote exploration data with missing channels or bad channels into a mapping domain through employing a high-resolution Radon operator, and carrying out the separation and removal of the data in the mapping domain if the data have noise; the mapping domain data are converted back to a time domain through a high-resolution inverse radon mapping operator, and reconstruction data corresponding to a missing track or a bad track are obtained; and performing energy balance processing on the reconstructed single-channel data and the original data, and outputting reconstructed sound wave remote detection data. Compared with the prior art, the fidelity data reconstruction technology based on logging sound wave remote exploration data mapping is good in stability, high in adaptability and good in fidelity. According to the method, based on the seismic wave propagation theory, the high-fidelity data mapping operator is utilized to reconstruct the sound wave remote exploration data, the reconstruction result complies with the geophysical law, and the reasonability is high.

The invention relates to the technical field of seismic engineering, and provides a method and system for solving a seismic wave propagation equation based on a physical information neural network, and the method comprises the steps: determining a space domain, a time domain, an initial condition and a boundary condition of a one-dimensional seismic wave equation; initial boundary value coordinate points are randomly extracted to serve as a training data set constraint network, and Latin hypercube sampling is utilized to generate space-time domain configuration points to meet partial differential equation constraint; a physical information neural network model is established, after space-time coordinates are input, equation residual errors are calculated through automatic differential, and initial boundary value loss and physical residual errors are combined to construct a composite loss function for training; and proposing a self-adaptive region sampling strategy based on residual errors, resampling a high residual error region according to an average residual error value of the configuration points in a training period, and iteratively optimizing network parameters after a newly added point and an initial configuration point are fused. According to the method, by dynamically enhancing the sampling density of a high-error region, the L2 relative error of a prediction solution is reduced by one order of magnitude, the network is effectively prevented from falling into local optimum, the convergence efficiency is remarkably improved, grid discretization or prior data support is not needed, and a high-precision meshless solution scheme is provided for a one-dimensional seismic wave equation.

The invention relates to an imaging gather Q value compensation processing method, device and equipment and a storage medium. The imaging gather Q value compensation processing method comprises the following steps: acquiring an imaging gather; performing windowing processing on the seismic signals in the imaging gather by using a Hamming window; grouping the imaging gathers according to offsets, and determining relative Q values of time window positions corresponding to different offset groups relative to the minimum offset group based on a mean spectrum ratio method; and based on stable inverse Q filtering, compensation is carried out through the relative Q value. Relative compensation is carried out on the gathers according to the spectrum consistency of the gathers, the defect that a seismic wave propagation rule is not considered in a traditional single-channel compensation algorithm is overcome, seismic data is compensated through the method, the in-phase property of the imaging gathers can be effectively improved, and due to the fact that the relative Q value is used, the signal-to-noise ratio cannot be remarkably reduced, and the imaging accuracy is improved. Finally, the imaging effect can be effectively improved, the imaging result is stable, the imaging resolution is high, and the compensation precision is relatively good.

The invention provides a tunnel micro-seismic blasting targeted protection control method and device, and relates to the technical field of tunnel blasting, and the method comprises the steps: obtaining the distance between a protected building and a tunnel, and building a blasting vibration wave propagation model based on the distance and an initial attenuation coefficient; calculating predicted vibration data of seismic waves generated by each blasting section in a preset blasting scheme, wherein the seismic waves are transmitted to the protected building; taking the minimum theoretical peak vibration velocity of the protected building as a target, and calculating by adopting a numerical optimization algorithm to obtain a non-equal-difference detonation delay sequence of each blasting section; when the detonating time interval of the adjacent detonating sections is not equal to the integral multiple of the inherent main earthquake period, detonating is carried out according to the non-equal-difference detonating delay sequence; and comparing the actually measured vibration data with the predicted vibration data, and adjusting the attenuation coefficient based on the comparison result, thereby effectively preventing the resonance effect, preventing the vibration energy from being amplified, and ensuring the construction safety.

This invention discloses an artificial intelligence-based emergency drill visualization generation system, belonging to the technical field of data visualization. It includes a marking module, a judgment module, and an association module. The marking module marks the first and second timestamps of buildings. The judgment module determines whether to perform a display or rendering operation. The association module performs association operations based on the judgment results. This invention achieves synchronous visualization mapping between the drill process and the evolution of the disaster situation through a dual timestamp comparison mechanism. By setting a three-level state judgment logic, it achieves automatic phased switching of the disaster response. Through the sequential construction and coverage update of P-wave diffusion rings and S-wave diffusion rings, it can dynamically display the spatiotemporal diffusion process of seismic wave propagation centered on the epicenter, intuitively presenting the spatial evolution trend of the disaster over time.

This invention discloses a modeling method for layered soft soil subgrades reinforced with pile foundations under simulated seismic loading, belonging to the field of seismic numerical simulation in geotechnical engineering. This method constructs a multi-field coupled numerical model based on a finite element platform, defines the contact mechanisms at the pile-soil and reinforced-subgrade interfaces, and completes the construction of in-situ stress equilibrium and initial pore water pressure fields. Static results are output through fill-consolidation coupling analysis. Distributed viscoelastic artificial boundaries are automatically constructed using Python scripts, and seismic waves are converted into equivalent nodal loads using wave theory. Multi-physics field connections between static and dynamic analysis boundaries are achieved, and core response parameters such as embankment settlement and pile-soil stress are output through implicit dynamic analysis. This invention considers the seismic wave propagation effect in layered soils, reduces spurious oscillations caused by the transformation of static and dynamic boundaries, and improves simulation accuracy and efficiency. It is suitable for seismic design and analysis of transportation infrastructure in soft soil areas.

The invention belongs to the technical field of landslide prediction, and relates to a shallow earthquake-stress joint inversion tunnel portal landslide prediction method, which comprises the following steps of: obtaining slope body first-arrival wave travel time data through a high-frequency seismic source array, constructing a geologic model containing a rock mass integrity coefficient matrix, identifying rock mass structure deterioration through elastic wave velocity distribution, and predicting the landslide at a tunnel portal through shallow earthquake-stress joint inversion. The method comprises the following steps: locking a potential unstable structure weak area in advance, then embedding a rock mass integrity coefficient matrix as priori knowledge into a stress field for solving, and by establishing a coupling relationship between seismic wave propagation and stress, inverting slope stress field distribution, capturing an early stress redistribution signal caused by structural degradation of a deep rock mass, and finally obtaining an initial stress redistribution signal of the deep rock mass. Early perception of mechanical instability precursor is realized; and then calculating strain energy density distribution of the potential slip surface based on the stress and strain data, and accurately judging the potential slip surface before the slip surface is not completely cut through by judging whether the strain energy density exceeds a threshold value, whether a continuous sheet is formed and whether the expansion rate exceeds the standard, thereby avoiding lagging caused by dependence on a displacement signal.

The application discloses a method for extracting and correcting the time difference of synchronous seismic data collection while digging, which comprises the following steps: 1, establishing a seismic data collection system while digging, reading all seismic wave data in the collection system and filtering; 2, performing velocity analysis on the filtered seismic wave to obtain the propagation velocity of the seismic direct wave while digging; 3, using the seismic wave propagation velocity original data to correct the direct wave time difference of each seismic trace, and eliminating the influence of the seismic wave propagation time difference; 4, calculating the model trace data of each collection substation through the data after the time difference correction; 5, determining the reference substation according to the model trace data of each collection substation and the data of each seismic trace in the corresponding collection substation; and 6, calculating the time difference of the corresponding collection substation according to the cross-correlation time difference between the model trace of the reference substation and the model trace of each collection substation, and correcting the original data.

The application discloses a two-dimensional distributed optical fiber sensing full waveform simulation method based on a discrete elastic lattice, and belongs to the technical field of seismic wave field simulation and distributed optical fiber sensing. The application comprises the following steps: receiving seismic wave data of a target geological medium layer; setting a two-dimensional calculation region according to the spatial distribution of the longitudinal wave velocity, the transverse wave velocity and the density parameters of the seismic wave, and taking the two-dimensional calculation region as a numerical simulation carrier of seismic wave propagation; discretizing the two-dimensional calculation region into a plurality of particle nodes, and connecting adjacent particle nodes through elasticity to form a discrete elastic lattice, so as to represent the mechanical properties of the stratum medium; and according to the density parameters, the longitudinal wave velocity, the transverse wave velocity or the equivalent elastic parameters of the stratum medium, the elastic connection between the particle nodes is endowed with corresponding elastic coefficients. The application provides a reliable numerical tool for DAS deployment design, data interpretation and seismological inversion, and can be widely applied to the fields of earthquake monitoring, oil and gas exploration, geological disaster warning and the like.

The application discloses an explicit finite element parallel computing method for seismic wave simulation, comprising the following steps: constructing a three-dimensional computing domain, generating an unstructured finite element grid based on real terrain and velocity structure data; performing domain decomposition on the computing domain, and dividing the grid into multiple sub-regions; distributing the tasks of the sub-regions to parallel computing nodes, discretely solving a three-dimensional elastic wave equation containing a damping term by using explicit time integration and a lumped mass matrix, synchronously calculating time sequences of wave fields of the sub-regions; integrating the results to obtain seismic wave field distribution, and outputting ground motion parameters of a target region. By constructing an integrated explicit finite element parallel computing framework of a seismic source-medium-terrain, efficient, stable and high-precision simulation of seismic wave propagation under complex geological conditions is realized, and the reliability of strong ground motion prediction is significantly improved.