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

The present invention relates to the field of data visualization technology, and in particular to a method and system for rapid visualization of earthquake emergency information. The method comprises the following steps: acquiring real-time earthquake monitoring data; gridding the population density distribution and key facility locations in the earthquake zone based on the real-time earthquake monitoring data to generate a sensitivity matrix for the earthquake zone; dynamically simulating the seismic wave propagation process based on the real-time earthquake monitoring data and the seismic zone sensitivity matrix to generate wave propagation time series data; identifying epicenter topographic and geomorphological vector data from the real-time earthquake monitoring data, and separately evaluating the topographic relief and geological stability to obtain earthquake topographic relief index data and earthquake geological stability score data. The present invention utilizes data visualization technology to achieve rapid and accurate assessment and visualization of earthquake damage distribution, thereby improving the efficiency of earthquake emergency response and the scientific nature of disaster relief decision-making.

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.