88 results about "Seismic wave propagation" patented technology

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 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.

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 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 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 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.

This invention relates to the field of seismic data processing technology, specifically disclosing a method for removing strong reflection shielding in hidden river channels based on adaptive hybrid L0-L1 norm. The method includes: first, establishing an optimization objective function based on Bayes' theorem using hybrid L0-L1 norm; second, dynamically adjusting the L0-L1 norm weights using an adaptive weight function driven by the seismic signal, combined with reflection coefficient amplitude and residual information, enhancing sparsity constraints in strong reflection zones and reducing constraint strength in weak reflection zones; then, constructing a convex upper bound for the objective function using a minimization framework, and solving it iteratively in stages using an accelerated rapid iterative threshold shrinkage algorithm, while incorporating prior knowledge of seismic wave propagation laws and river channel deposition patterns to ensure the geological rationality of the solution. This invention solves the technical problems of traditional sparse processing methods, such as fixed parameters, lack of geological constraints, and inability to simultaneously address strong reflection suppression and weak signal protection, significantly improving the separation accuracy of strong reflections and the recovery rate of weak signals, effectively overcoming the strong reflection shielding effect.

The invention discloses a method and system for calculating a seismic oscillation amplification coefficient of a river valley field under SV wave incidence, and belongs to the technical field of rock-soil seismic engineering, and the method comprises the steps: building a model; analyzing an infinitesimal body; establishing a free field; establishing a scattering field; substituting boundary conditions; solving a matrix equation; and calculating an amplification coefficient. According to the method, a wave field caused by an irregular terrain of the earth surface is constructed by an interval decomposition strategy taking a terrain cuspidal point as a demarcation point, a method for calculating the seismic oscillation amplification coefficient of the river valley field under SV wave incidence is provided, and the influence of factors such as river valley geometric parameters, material parameters, the frequency and incident angle of incident waves and the like on the river valley seismic oscillation amplification coefficient is revealed; a theoretical basis is provided for aseismic design of major projects of river valleys; the calculation method is a semi-analytical method, reveals the physical nature of irregular terrains on seismic wave propagation and scattering, improves the calculation precision and efficiency, and can be used as a reference for verification of a numerical method.

The application discloses a low-frequency broadband seismic metamaterial filled with a five-mode dot array core and a design method thereof. The seismic metamaterial is composed of a five-mode dot array core, a concrete hollow pipe and an outermost soil layer. The five-mode dot array core is composed of a plurality of hexagonal crystal system units arranged periodically. The hexagonal crystal system unit is connected by an axial direction and is obtained by space structure conversion. The hexagonal crystal system unit is arranged periodically to form a three-dimensional dot array core structure composed of a plurality of hexagonal crystal system units. The low-frequency broadband seismic metamaterial filled with the five-mode dot array core can form a low-frequency broadband in a target frequency band. The seismic metamaterial is arranged in a vibration isolation area underground. When a seismic Rayleigh wave propagates to the area, the wave energy is significantly weakened under the action of a forbidden band, so that effective isolation and attenuation of the Rayleigh wave can be realized under the condition of a limited space.

This invention relates to the field of numerical simulation and computational technology for seismic wave propagation, and discloses a method and apparatus for generating seismic ground motions based on target spectra and Green's functions. It aims to address the difficulty in achieving a balance between physical consistency, computational efficiency, and spectral compatibility in existing methods. The scheme mainly includes: preprocessing environmental noise records and extracting empirical Green's functions through unsupervised clustering; constructing a spatially continuous propagation tensor by performing depth correction and gradient interpolation based on surface wave eigenfunctions; discretizing the kinematically finite fault model into sub-source units and generating long-period seismic ground motions including path effects using propagation tensor convolution; optimizing source parameters through Bayesian inversion using the GMPE target spectrum as a constraint; and finally generating high-frequency components and fusing them with the long-period waveform to obtain a broadband seismic ground motion time history. This invention achieves a unification of physical propagation mechanisms and statistical spectral constraints, significantly improving computational efficiency while ensuring spectral compatibility, and is particularly suitable for basin areas.

The application provides a seismic motion simulation method based on a stochastic finite fault model. The method comprises the following steps: determining an equivalent magnitude and an equivalent epicentral distance from different exceedance probability levels of a target area, determining geometric parameters of a fault model according to the equivalent magnitude, wherein the geometric parameters comprise a fault burial depth, a fault plane size and a number of sub-faults; dividing the fault plane into a plurality of sub-faults according to the number of sub-faults; calculating rupture propagation time and seismic wave propagation time of each sub-fault to an observation point by using a bilateral rupture mode; generating a seismic motion time history of each sub-fault by using the stochastic finite fault model; introducing a frequency attenuation factor and a path attenuation factor to correct the seismic motion time history; and superimposing corrected seismic motion time histories corresponding to all sub-faults to generate a seismic motion field of the target area. The technical scheme provided by the application realizes fine simulation of a fault rupture process based on a multi-probability scenario, and improves the physical rationality and spatial resolution of target area seismic motion field prediction.

The invention discloses a rock physical forward modeling method applicable to formation pore fluid of carbon dioxide, and relates to the technical field of petroleum and natural gas engineering. The method comprises the steps of firstly obtaining a carbon dioxide physical property constraint condition based on a state equation of carbon dioxide under a reservoir temperature and pressure condition, and then based on the constraint condition and according to logging information of a target well, core analysis and a laboratory test result, carrying out formation evaluation to calculate reservoir parameters. Then, based on reservoir parameters, the sandstone content is obtained through calculation, a rock physical volume model of a target well is constructed, and finally, based on the model, corresponding rock physical forward modeling is conducted according to different carbon dioxide saturations by means of the equivalent medium theory to generate stratum elastic attribute parameters. And the change relation between the corresponding saturation and the stratum elastic attribute parameters is quantitatively represented, so that the CO2 quantitative monitoring precision and reliability can be remarkably improved by comprehensively considering the special physical properties of CO2, the fluid replacement mechanism and the seismic wave propagation characteristics.

The invention relates to a ground motion field simulation method based on a river valley site terrain effect, which is mainly applied to the field of slope anti-seismic analysis of water-power engineering. The method comprises the following steps: firstly, aiming at topographic features of a symmetric V-shaped river valley, constructing a planar model by combining a cylindrical SH wave theory, obtaining a seismic wave propagation rule and scattering characteristics by solving a wave field structure, and verifying a topographic amplification effect by utilizing actually measured data; and carrying out seismic oscillation execution foundation terrain simulation, flood discharge tunnel outlet slope simulation and landslide accumulation body simulation in three stages. According to the method, three types of simulation results are integrated, and the most unfavorable seismic oscillation input scheme for engineering anti-seismic analysis is determined, wherein seismic oscillation peak acceleration data of the wave-facing side during horizontal incidence are adopted by a flood discharge tunnel outlet slope and a landslide accumulation body. By means of the method, efficient coupling simulation of the river valley terrain effect and the non-uniform seismic oscillation field is achieved, and theoretical support is provided for aseismic design of complex site slopes.

The invention discloses a construction method, system, equipment and medium of an oil and gas reservoir velocity-temperature model, and particularly relates to the technical field of model construction, and the technical key points are as follows: obtaining environmental thermodynamic parameters of an underground oil and gas reservoir medium and sound pressure components of seismic waves; constructing a thermal-acoustic coupling medium seismic wave propagation model based on an elastic theory and an acoustic approximation theory; based on a thermal-acoustic coupling medium seismic wave propagation model, forward modeling is carried out by using a staggered grid finite difference method, and a simulated sound pressure component of a seismic wave is obtained through simulation; acquiring an actual seismic sound pressure component, and constructing a target function by minimizing an error between the actual seismic sound pressure component and the simulated sound pressure component; implementing full-waveform inversion by using a gradient descent method, and solving an objective function to obtain an optimal model parameter; and performing parameter updating on the thermal-acoustic coupling medium seismic wave propagation model by using the optimal model parameters to obtain an updated oil and gas reservoir velocity-temperature model.

The application provides a silt geology building pile foundation positioning method, and belongs to the technical field of pile foundation construction. The silt geology building pile foundation positioning method comprises the following steps: determining a main target detection area and a subordinate target area of silt geology; controlling a seismic source to send seismic waves from the subordinate target area to the main area; obtaining a first first-arrival time of the seismic waves from any geophone; determining a first target velocity of seismic wave propagation at each position in the main target detection area and a velocity inversion graph by using a seismic wave tomography method according to the first first-arrival time; determining the range of a silt geology high stress area in the velocity inversion graph according to the first target velocity and the corresponding relationship between velocity and stress; comparing the values of the vibration sensors in each detection hole in the high stress area range to determine the best piling position. The application solves the problems that the existing silt geology building pile foundation positioning method cannot quickly determine the range of the high stress area and cannot determine the best piling position.

The invention discloses a high-speed rail seismic source seismic data wave field separation method and system using random Cadzow filtering, and belongs to the technical field of exploration geophysics. The method comprises the following steps: acquiring seismic data when a high-speed rail passes; estimating the seismic wave propagation speed through Radon transformation; performing time difference correction on the same-direction wave field event to obtain aligned data; blocking the aligned data and randomly rearranging the aligned data; fourier transform is carried out on the rearranged data, and a Hankel matrix is constructed; carrying out SVD (Singular Value Decomposition) decomposition and low-rank reconstruction on the Hankel matrix; reconstructing a frequency domain data block through an anti-Hankel operation; inverse Fourier transform is carried out, and a trace sequence is recovered; and extracting and averaging a central channel to obtain unidirectional wave field data. The system comprises modules with corresponding functions. The method can effectively separate the high-speed rail seismic source mixed wave field, improves the signal-to-noise ratio of data, and provides a reliable data basis for underground imaging.

The application discloses a kind of microseismic event precision positioning method and system based on knowledge enhanced model, comprising the following steps: collecting the waveform signal of downhole and surface multi-array geophone, obtaining spatial coordinates and synchronous clock information, generating propagation time field representing seismic wave propagation;Determine the spatial constraint range, set continuity and velocity condition to form knowledge constraint domain;Geometric energy model of propagation time distribution is established in the constraint domain, and the propagation time field is iteratively updated by variational solution, and converges to stable space-time manifold;Identify the point that meets the condition on the manifold as the location and time of microseismic event, and according to local residual and noise level, grid refinement and recalculation are carried out, and the reliability evaluation of event positioning result is completed.The application is based on knowledge enhanced propagation time field modeling and variational manifold solution, and realizes high-precision positioning and reliability evaluation of microseismic event in complex geological environment.

The invention belongs to the technical field of rock physics and multiphase flow coupling seismic wave propagation simulation, and discloses a wave velocity-attenuation whole-process calculation method based on an Lo two-phase fluid model. According to the method, thermodynamic characteristics of CO2-saline two-phase fluid are obtained by utilizing temperature and pressure conditions of a target reservoir, seepage characteristics are fitted according to experimental data, and a unified medium physical property layer is constructed by combining fluid parameters and rock skeleton parameters; after an elastic matrix, an inertia matrix and a viscosity matrix of a multiphase medium are calculated according to physical property parameters, phase velocity and attenuation characteristics of three types of longitudinal waves and one type of transverse waves are obtained through a Lo two-phase fluid model control equation. According to the method, a complete technical process from physical property construction, control equation solution to sensitivity judgment and mechanism analysis is formed, the method can be used for seismic monitoring, reservoir feature recognition and relevant parameter inversion in the CO2 geological storage process, and the method has good engineering applicability and popularization value.

The invention discloses a meshless seismic wave forward modeling method based on a difference numerical method and related equipment, and belongs to the field of exploration geophysics. According to the method, flexible mesh generation is carried out based on an undulating surface seismic wave speed model, a surface interface is completely fitted, and a high-precision difference coefficient set is obtained without manually selecting deformation parameters; by adding a filter related to a second-order time item in a two-dimensional scalar time domain acoustic wave equation, the robustness and accuracy of numerical simulation are effectively improved, so that a high-precision and stable seismic wave numerical simulation result is obtained, the precision and stability of seismic data processing can be effectively improved by utilizing the method, and the seismic data processing efficiency is improved. The method provides solid technical support for loess tableland seismic wave forward modeling, and has important practical significance for evaluating the influence of the earth surface effect on exploration data in the loess tableland, improving the accuracy and reliability of the exploration data, analyzing the physical properties and change rules of underground media and simulating the physical properties and seismic wave propagation conditions of underground structures.

The invention belongs to the technical field of seismic exploration, and relates to a porous rock seismic wave propagation characteristic simulation method based on multi-mechanism coupling, which comprises the following steps: simulating the propagation of seismic waves in viscoelastic porous rock by adopting a seismic wave propagation model; wherein the modeling process of the seismic wave propagation model comprises the following steps: coupling a jet flow mechanism with a Biot theory framework; the method comprises the following steps: representing the inherent viscoelasticity of a solid skeleton by adopting a Zener model, and coupling with a jet flow mechanism and a Biot theoretical framework to obtain a total volume modulus and a total shear modulus; calculating a memory variable according to the total volume modulus and the total shear modulus; obtaining a time domain relaxation function containing a memory variable; substituting the time domain relaxation function containing the memory variables into the frequency domain stress-strain relation to obtain a time domain speed-stress equation set containing all the memory variables; and obtaining a seismic wave propagation model by combining the first-order differential equation of the memory variable. According to the invention, wave field characteristics such as reflection and diffraction in a complex heterogeneous medium can be accurately and comprehensively simulated.

The application belongs to the technical field of earthquake prevention and disaster reduction, and proposes a seismic isolation and damping system based on spatial arrangement of trees. The system comprises a density gradient tree arrangement area, a periodic tree array area, a target protection area and a buffer transition zone. The density gradient arrangement area is located at the front end of the system, and the tree density gradually increases along the direction of seismic wave propagation, which is used to realize impedance smooth transition and guide wave energy into the array area; the trees in the periodic array area are arranged at equal intervals to form a quasi-periodic elastic structure, which is used to produce Bragg scattering band gap effect, reflect and dissipate seismic waves of a specific frequency band; the target protection area is used to arrange buildings or infrastructure that need to be protected; the buffer transition zone is arranged between the array end and the target area, which is used to absorb residual wave energy and reduce echo amplification effect. The application constructs a green, low-cost and sustainable natural seismic isolation system, which has good seismic wave control effect and engineering adaptability.

The present application relates to the technical field of machine learning, in particular to a seismic wave propagation feature analysis and prediction system, the system comprising: a feature decoupling coding module, a conditional model generation module, a forward propagation prediction module, and a feature quantization analysis module.In the present application, the intrinsic response features representing the medium are separated from the actual multi-channel seismic data channel set, and based on the features, a series of possible velocity model samples are probabilistically constructed using a conditional generation model, and then all the velocity model samples are iteratively simulated for forward propagation, thereby obtaining a wave field state snapshot sequence set containing multiple possibilities.Through statistical calculation of the set, the expected propagation amplitude can be obtained and the predicted uncertainty can be quantified, and finally a comprehensive atlas of expected wave field propagation and uncertainty analysis is generated.

The application discloses a fractional-order viscous wave equation numerical simulation method, system and device, relates to the technical field of seismic wave propagation numerical simulation, and comprises the following steps: substituting a wave equation plane wave solution into a fractional-order Laplacian operator finite difference formula to obtain a generating function corresponding to the finite difference formula; replacing regular grid finite difference coefficients with optimized regular grid finite difference coefficients, performing inverse Fourier transform on the generating function to obtain optimized difference coefficients; based on the optimized difference coefficients of the finite difference formula and the fractional-order Laplacian operator finite difference formula, performing numerical dispersion on a preset fractional-order Laplacian viscous wave equation to obtain a fractional-order Laplacian viscous wave equation finite difference wave field recursion formula of seismic wave field data, and then simulating seismic response characteristics in an attenuation medium. The application can efficiently simulate numerical values of seismic response characteristics in the attenuation medium.