4059 results about "Earth quake" patented technology

The invention provides a tunnel disaster source multi-attribute joint inversion method and system based on physical property constraint, and relates to the technical field of tunnel engineering and joint inversion, and the method comprises the steps: constructing a tunnel face three-dimensional geologic model, applying controlled disturbance to surrounding rock physical property parameters, analyzing the change characteristics of earthquake and electromagnetic observation response, and obtaining a tunnel face three-dimensional geologic model; establishing a quantitative mapping relation between the surrounding rock physical property parameter disturbance quantity and the multi-field response signal variable quantity; constructing a multi-dimensional sensitivity evaluation index, and generating a joint sensitivity result; the surrounding rock units are identified and divided, differentiated partition rock physical relationship constraints are constructed, and a physical multi-attribute inversion target of joint inversion is determined based on the partition rock physical relationship constraints; and on the basis of a physical multi-attribute inversion target, a multi-source joint inversion target function is constructed in combination with the joint sensitivity result, joint solving and iterative updating are carried out, and a tunnel face advanced inversion result is output. According to the invention, calamity source fine identification and risk quantitative evaluation in front of the tunnel face during the tunnel construction period are realized.

Disclosed in the present invention is an anti-seismic simulation method for an existing building by considering a real damage state. The method comprises: determining a quantified component damage index on the basis of an observed component damage state; determining a material damage index of a damaged component on the basis of the component damage index; determining key parameters of a constitutive law for a damaged material on the basis of the material damage index; establishing a finite element model of a damaged building on the basis of the constitutive law for the material of the damaged component; and assessing the seismic risk and anti-seismic performance of the damaged building on the basis of the finite element model. In the present invention, a constitutive law for a damaged material is incorporated into a fiber model in OpenSEES, so that rapid modeling of a damaged building can be realized, thereby effectively reproducing the real damage of an existing building; and by using a seismic vulnerability analysis method, rapid assessment of the seismic risk and anti-seismic capability of the existing building can be realized, thereby further providing support for formulating reasonable reinforcement and reconstruction strategies by related departments.

The invention discloses a coal seam passive source seismic survey system and method based on deep learning, and relates to the technical field of coal mining safety monitoring and early warning, and the method comprises the following steps: S1, collecting the passive source seismic data of a coal seam area; s2, deep fusion features are extracted; s3, constructing a semi-supervised self-attention prediction network model; s4, extracting coal seam structure prediction parameters; s5, the corrected coal seam structure parameters are extracted; s6, constructing a fast seismic wave forward model, and calculating and outputting an initial model of the coal seam underground structure; s7, constructing a Transform optimization network, and extracting a depth feature vector; and S8, constructing a coal seam structure inversion model, and training iteration by using the improved batch depth Bayesian optimization algorithm. The method overcomes the limitation of manual analysis and slow response in a traditional coal seam seismic survey method, and provides an efficient and accurate solution for safe production of a coal mine.

The invention relates to the technical field of computational mechanics and geological disaster simulation, and discloses a full-process numerical simulation method suitable for simulating continuous medium damage to discontinuous medium movement. The invention provides a novel chain-type disaster simulation method for solving the problems that calculation efficiency and large deformation precision are difficult to consider at the same time and the evolution process of a slope from a continuum to a fragmented body cannot be dynamically reflected in the prior art. The core of the method is that a dynamic criterion is set based on a unit real-time damage variable, a failed finite element (FEM) unit is adaptively converted into smoothed particle hydrodynamics (SPH) particles, and the physical state of the particles is accurately mapped; and a virtual particle coupling algorithm is adopted to process a two-domain interface, so that bidirectional transmission of mechanical parameters is realized. The method is executed circularly, the whole process of'continuous deformation-fragmentation flow 'of the side slope is simulated dynamically, high-precision integrated simulation of a'seismic source-propagation-response-movement' disaster chain in the same frame is achieved, and an efficient tool is provided for disaster risk assessment and prevention and control.

The invention discloses a fault co-seismic sliding surface inversion method and system based on Beidou and artificial intelligence. The method comprises the steps of collecting Beidou data and seismic waveform data of a target area in real time, and performing preprocessing to extract co-seismic displacement and waveform features; fusing the Beidou deformation features and the seismic wave features to form a multi-modal feature vector; a fault prediction model is constructed based on a graph neural network, forward modeling and historical earthquake example data training are utilized, and multi-modal features are input to obtain preliminary fault sliding distribution prediction; rapid optimization under physical constraint is carried out through the elastic dislocation model, and a final fault sliding model conforming to the geophysical law is obtained; finally, uncertainty quantification is carried out, and a visual product of the dynamic process including sliding distribution, seismic moments and fault deformation is generated. According to the system, full-process automatic processing is achieved through cooperation of all the units, the inversion speed, precision and physical credibility are improved, and reliable support is provided for earthquake emergency response and disaster assessment.

The invention provides a three-dimensional visualization analysis method for a real scene of an ancient earthquake exploring groove, which belongs to the technical field of seismic geological survey and comprises a data acquisition module used for acquiring point cloud data of target exploring grooves in different regions. And the model building module is used for carrying out scene format conversion on the point cloud data of the target exploring grooves in different regions, and carrying out exploring groove three-dimensional model building by utilizing the converted SLPK file. And the visualization module is used for carrying out live-action three-dimensional display on the three-dimensional models of the target probe grooves in different regions through a three-dimensional scene Web server based on the SLPK files in different regions, and rendering node data of different detail levels of the probe grooves according to user requirements. The ancient earthquake event can be better understood, and a reliable three-dimensional digital file and a new shared analysis normal form are provided for the ancient earthquake.

The invention belongs to the technical field of seismic exploration, and discloses a robust full-waveform inversion method, system and device based on generative modeling, and a medium, and the method comprises the steps: obtaining seismic observation data; the method comprises the following steps: starting from random Gaussian noise, establishing a model space comprising a plurality of initial velocity models through an unconditional score model, and determining a global optimal initial velocity model of the model space by adopting a strategy search algorithm; de-noising is carried out on the global optimal initial velocity model based on back diffusion, seismological observation data are introduced as observation constraints, the velocity model is updated by calculating the mismatch gradient of forward modeling data of the current velocity model and the seismological observation data, and an implicit condition sample is obtained; performing forward diffusion processing on the implicit condition sample to obtain a speed model for a subsequent annealing time step; and when the annealing process reaches a preset condition, outputting a final speed model. According to the method, the robustness and accuracy of full-waveform inversion are improved, and the bottleneck problem of traditional full-waveform inversion is solved.

The invention discloses a three-dimensional geological model dynamic updating system and method based on multi-modal data fusion, and belongs to the technical field of geological engineering and information. According to the method, multi-source data of geological radar, earthquake, drilling, remote sensing and the like are integrated, and after unified formatting and quality optimization preprocessing, geological features of a stratum interface, a fault zone and the like are automatically extracted by using a deep learning network; multi-source feature data fusion is realized through space-time registration and attention; constructing a three-dimensional geologic model by adopting a TIN and Octree mixed structure based on a fusion result; dynamic updating of the model is realized through incremental detection and a local reconstruction mechanism, and the boundary of the model is optimized by adopting a radial basis function interpolation. According to the method, the problems of single data, updating lag and the like of a traditional method are solved, the modeling precision and efficiency are remarkably improved, and reliable support is provided for geological engineering decision making.

A 3D inversion method of DC method using seismic data constraints is disclosed, comprising processing of seismic data and inversion processing of resistivity data; by utilizing seismic wave interval velocity information to constrain the high-density electrical inversion, which can play a complementary role to a certain extent, it is beneficial to supplement the resolution of resistivity inversion with 3D seismic image information, reflect a boundary of subsurface electrical anomalous body, and further characterizes the boundary characteristics of anomalous body through comprehensive geological-geophysical analysis. Seismic data constraints are applied to 3D resistivity inversion, performing numerical simulation and inversion of electrical data under undulating terrain using the finite volume method and Gauss-Newton inversion method. Furthermore, the seismic anisotropic velocity tensor constraint is introduced for inversion to mitigate the non-uniqueness inherent in single-method inversion, yielding a more accurate 3D subsurface electrical structure and improving the resolution of longitudinal boundaries.

The invention discloses a three-dimensional tomography method and system for seismic exploration data migration, and relates to the technical field of geological exploration, and the three-dimensional tomography method comprises the following steps: laying a layered heterogeneous sensing array of adaptive terrain coupling in an exploration target area, the layered heterogeneous sensing array comprises a shallow fiber grating sensing layer, a middle MEMS sensing layer and a deep distributed fiber sensing layer, sensors are all provided with terrain adaptive coupling assemblies, and sensing shafts all point to the radial direction of an underground medium. According to the seismic exploration three-dimensional tomography imaging method, an iterative imaging mechanism of real-time geological constraint feedback is constructed, the limitation of traditional static geological constraint one-time importing is broken through, the imaging process can dynamically fuse newly added local geological data in the exploration process, the imaging model and the actual geological conditions are continuously adapted and optimized, and the imaging efficiency is improved. The integrating degree of the imaging result and the actual geological condition is enhanced, and the high requirement of deep complex area exploration for the imaging accuracy is met.

The invention discloses a seismic oscillation generation method based on a DFEM model, and the method comprises the steps: obtaining multi-source data in a target simulation region, and obtaining fusion stratum wave velocity structure data; determining the minimum grid size of the DFEM model according to the fused stratum wave velocity structure data and a preset simulation frequency upper limit, and dividing grids and determining nodes according to the minimum grid size; constructing a seismic propagation related physical constraint set; the physical constraint is embedded into a primary training process of a DFEM model, a Green function of a fracture surface sub-source in a target simulation area is obtained, convolution operation is carried out on the Green function, expanded seismic oscillation samples under different scenes are generated, and the expanded seismic oscillation samples and full-physical simulation seismic oscillation samples form a data set; performing secondary training on the DFEM model based on the data set to obtain a seismic motion prediction model; and generating diversified vibration samples based on the seismic motion prediction model. The problems that traditional pure physical simulation seismic oscillation simulation is high in cost, long in time and the like are solved.

The invention relates to the technical field of earthquake monitoring, and discloses a comprehensive earthquake prediction method, which comprises the following steps of: acquiring a continuous micro-motion data sequence to reconstruct an empirical green function sequence, and dividing the empirical green function sequence into a compression phase subset and a stretching phase subset by utilizing a synchronously acquired theoretical earth tide stress phase; respectively extracting a causal branch waveform corresponding to the positive time axis and a non-causal branch waveform corresponding to the negative time axis, and determining a weight index for inhibiting the non-structural interference by calculating a reciprocity residual error between a causal wave velocity offset and a non-causal wave velocity offset; correcting a wave velocity offset mean value by using the weight index to obtain a stress sensitivity coefficient reflecting the stress state of the earth crust medium; according to the method, a physical criterion is established by using an elastic wave propagation reciprocity principle, decoupling of medium attribute change and noise source spatial drift is realized, and visual wave velocity offset induced by environmental interference is eliminated.

The invention relates to the technical field of seismic data analysis, in particular to a seismic station ground motion noise quantification method and related equipment based on a GPU (Graphics Processing Unit), and the method comprises the steps: firstly obtaining waveform record data of a target seismic station in a set time period; performing sliding window segmentation processing on the waveform recording data in parallel by utilizing a GPU (Graphics Processing Unit) to obtain a plurality of waveform recording sub-segments; preprocessing the data of the plurality of waveform recording sub-segments in parallel through a GPU (Graphics Processing Unit); finally, the power spectrum density of the ground motion noise in the preprocessed sub-segments is calculated in parallel through the GPU; and carrying out parallel statistics on a power spectrum density probability density function result based on the power spectrum density by utilizing a GPU (Graphics Processing Unit) to complete the quantitative calculation of the earth motion noise. According to the method, the data processing efficiency is improved based on the parallel computing capability of the GPU, hardware computing resources are fully called, and the earthquake noise quantification of the seismic station record can be efficiently realized. The whole method is based on a CUDA programming model and a C language, cross-platform transplantation is easy, and the large-scale seismic network monitoring capability can be effectively improved.

The invention relates to the field of seismic data reconstruction, in particular to a seismic data low-rank reconstruction method fusing double-domain transformation, and the method comprises the steps: obtaining original incomplete seismic data, and initializing to-be-reconstructed seismic data and other model variables and parameters under an alternating direction multiplier method frame; performing fractional order gradient transformation on seismic data to be reconstructed along different directions, and mapping to obtain fractional order gradient domain data; performing space-time Hankel tensor expansion transformation on each piece of fractional order gradient domain data, and converting the data into a space-time Hankel matrix; performing low-rank constraint on each space-time Hankel matrix by adopting a Schatten norm to obtain a seismic data low-rank reconstruction model fused with double-domain transformation; and carrying out iterative solution on the model by using an ADMM framework until the model is stable, and obtaining a final seismic data reconstruction result. The method provided by the invention can effectively reconstruct the seismic data with the missing trace, and has a reconstruction effect with strong robustness and relatively good precision for incomplete seismic data with different missing degrees.

The invention discloses an earthquake forecasting and tsunami early warning method and system, and the method comprises the steps: monitoring the earth electromagnetic field data in real time, and recognizing and extracting a regional transient electromagnetic signal which accords with a preset abnormal feature as an earthquake precursor signal; taking the position of the precursor signal as a reference, and predicting the earthquake generation time, the magnitude and the potential influence area based on an earth crust stress wave propagation model; when the early warning condition is met, earthquake early warning information containing risk avoiding guidance is generated and issued, and if the prediction area relates to the coastal area and the earthquake magnitude reaches a threshold value, tsunami early warning information is synchronously generated and issued. The corresponding system comprises a geomagnetic monitoring module, a forecast analysis module and an early warning issuing module. According to the method, pre-earthquake prediction based on the physical precursor is realized, the prediction time is remarkably prolonged, the problems that earthquake and tsunami early warning information is separated and published lags are solved, and the disaster response capability is improved.

The present application belongs to the technical field of geophysical well logging data processing, and particularly relates to a multi-scale dual-branch neural network well logging curve reconstruction method based on seismic constraint, which comprises reading in well-seismic and interpreted horizon data, target interval sublayering and facies belt division based on three-dimensional seismic interpretation, pyramid multi-scale training set sample construction by fusing multi-curve and well flanking seismic trace information, residual-Transformer dual-branch neural network training model establishment taking a composite loss function as the core, three-stage progressive training strategy and reconstruction of target well missing section and result evaluation. The method still maintains high stability and robustness under the condition of complex heterogeneous reservoirs; on the basis of seismic data constraint, geological sublayering and facies belt division and other interpretation information are fused, so that the reconstruction result is more in line with geological laws; and the obtained reconstructed curve provides reliable support for re-interpretation of old data and full use of incomplete well logging data.

The invention discloses an earthquake disaster loss dynamic assessment method and system based on multi-source disaster situation data, and relates to the technical field of earthquake disaster loss assessment, and the method comprises the following steps: S1, obtaining the multi-source original data of an earthquake disaster, and carrying out the preprocessing, the multi-source original data comprises seismic oscillation parameter data and multi-source disaster data; and S2, based on the preprocessed multi-source data, constructing a coupling driving model of the earthquake influence field and the multi-source disaster situation data, and generating a regional disaster damage thermodynamic map. According to the earthquake disaster loss dynamic evaluation method and system based on the multi-source disaster situation data, the earthquake multi-source disaster situation data is acquired, the coupling driving model of the earthquake influence field and the multi-source disaster situation data is constructed, and the influence of various factors on the earthquake disaster loss is comprehensively considered, so that the limitation that a traditional method depends on a single data source is overcome; the real situation of earthquake disaster loss can be reflected more comprehensively and accurately, and a more reliable basis is provided for rescue decision making.

The invention discloses a stratigraphic water content layering inversion method based on a shallow seismic full waveform, and relates to the technical field of geophysical exploration data processing, and the method comprises the following steps: 1, under the same survey line condition, synchronously obtaining a seismic wave full waveform, complex impedance and an induced polarization frequency spectrum; 2, constructing a layered initial model containing a layer boundary line and an in-layer water content initial value; based on the layered initial model, executing dual-domain quantile alignment interface keeping iteration until convergence, and obtaining a final layer boundary line and an in-layer water content field; the double-domain quantile alignment interface maintaining iteration comprises the following steps: updating an in-layer moisture content field and a layer boundary line according to a seismic wave full waveform; and performing total variation regularization processing on the updated in-layer water content field to maintain the transition characteristic of the layer boundary line. According to the method, an efficient and reliable solution with actual guidance value is provided for mineral resource exploration, geological structure detection, groundwater survey and engineering geology evaluation.

The invention relates to the field of environmental analysis, and provides a gas hydrate three-phase coexistence region modeling method and system. The method comprises the following steps: acquiring three-dimensional seismic data and original logging data, and preprocessing the original logging data to obtain corrected logging data; respectively performing logging interpretation and structural interpretation according to the corrected logging data and the three-dimensional seismic data to obtain interpretation results; converting a structure interpretation result in the interpretation results into a three-dimensional geological structure model, and performing grid division on the research area by adopting a corner grid system to obtain a structure grid model; and performing attribute modeling on the constructed grid model based on variation function analysis, and predicting porosity and saturation in a three-dimensional space through a Kriging interpolation algorithm to obtain a three-phase coexistence region model. According to the method, the influence of complex geological conditions on the accuracy of the model is reduced, and the precision of the three-phase coexistence region model is remarkably improved.

The invention relates to the technical field of oil and gas geophysical exploration, and particularly discloses a tight sandstone gas reservoir fluid mobility gradient attribute prediction method and system.The method comprises the steps that a reflection coefficient expression containing a fluid mobility item is constructed, and a formula for calculating fluid mobility through frequency dispersion attributes is deduced; constructing a fluid mobility gradient attribute related to the frequency and the incident angle of the seismic wave; the method comprises the following steps: establishing a double-pore rock physical model, generating a synthetic seismic record, applying frequency variation AVO inversion to the synthetic seismic record, and extracting a bulk modulus dispersion attribute and a shear modulus dispersion attribute; and calculating the fluid fluidity changing along with the frequency and the incident angle, and calculating the gradient attribute of the fluid fluidity based on a gradient attribute calculation formula and the fluid fluidity. According to the method, seismic wave frequency dispersion and attenuation response caused by permeability and gas saturation are fully utilized, the seismic gradient attribute of the fluid fluidity and the calculation method are constructed, and high-precision prediction of the fluid fluidity by applying seismic signal elasticity and frequency dispersion rock physical information is realized.

The utility model relates to a sliding shock insulation modular building which comprises a module group forming a building main body. The building foundation is used for bearing the module group; the sliding layer is arranged between the module group and the building foundation and comprises an upper sliding interface and a lower sliding interface which are fixedly connected with the bottom of the module group and the top of the building foundation respectively, and a sliding medium arranged between the upper sliding interface and the lower sliding interface; and the sliding starting control part is connected with the module group and the building foundation and is used for controlling the sliding starting state. Compared with the prior art, the earthquake response of an upper modular building can be reduced, non-structural members in the modules can be protected, the module units are lighter in design, the connection between the adjacent modules is simpler, and the connecting structure adapts to the characteristics of various layouts and stress; the anti-seismic performance of the modular building is improved through simple sliding seismic isolation measures which are low in manufacturing cost and convenient to construct.

The invention discloses an edge sea petroliferous basin paleotopography inversion method based on a machine learning algorithm, and belongs to the technical field of geological exploration and data processing. The method comprises the following steps: constructing a seismic profile data pool and carrying out image preprocessing; establishing an affine transformation relationship between pixel coordinates and physical coordinates by using an image digitization tool, and extracting a two-dimensional data sequence; establishing a mapping function of the measuring line distance and the longitude and latitude, and constructing a three-dimensional scatter data set; carrying out gridding interpolation by adopting a K nearest neighbor regression algorithm, and constructing a three-dimensional grid model; performing motion back-pushing and deformation-removing correction by using a plate reconstruction technology to obtain a paleotopography model; and generating a three-dimensional dynamic visualization result of paleotopography evolution through time interpolation. The method solves the problems that seismic profile images cannot be quantitatively analyzed, sparse data interpolation is difficult, and paleotopography inversion neglects structural deformation, and provides technical support for multiple fields such as oil-gas exploration.

The invention provides a TLF-GV signal correlation earthquake three-element joint probability prediction method and system, belongs to the technical field of earthquake prediction, and aims to solve the problems that traditional three-element prediction lacks joint fusion, evidence conflicts are difficult to process, and weights are fixed. Comprising the following steps: extracting physical features of TLF-GV signals, carrying out preorder three-element prediction results and model credibility parameters, splicing the physical features into feature vectors, inputting the feature vectors into a dynamic weight prediction model constructed by a shallow MLP, and outputting time, magnitude and position weights; dividing the three elements into 3 * 3 * 3 = 27'time-magnitude-position 'propositions, generating basic probability assignment (BPA) of each proposition based on weights, and quantifying evidence contradictions through a conflict coefficient K; and adopting a Dempster combination rule to obtain a joint BPA, calculating a mean value of a reliability function and a likelihood function of each proposition as a comprehensive probability, and extracting a maximum probability proposition as a coupling result. According to the method, the scientificity and practicability of strong earthquake short-term and temporary prediction are remarkably improved.

The invention discloses a seismic image super-resolution reconstruction method based on a double-branch multiple attention network, and belongs to the technical field of seismic exploration, and the method comprises the steps: constructing a data sample library for model training; adopting a super-resolution backbone network and a gradient branch auxiliary network to construct a double-branch multiple attention network model; and training the double-branch multiple attention network model by adopting the data sample library and the mixed loss function to obtain a network model for realizing seismic image super-resolution reconstruction and high-frequency texture information reconstruction, and performing reconstruction based on the trained network model to obtain a seismic image subjected to super-resolution reconstruction. According to the method, a gradient branch is introduced as an auxiliary network, so that the super-resolution reconstruction effect is enhanced. Meanwhile, a mixed loss function is adopted to train the network, and a loss function for constraining high-frequency information is particularly introduced, so that the common structural deformation problem in a super-resolution task is effectively solved, the generation of artifacts is prevented, and a new technical approach is provided for high-precision seismic interpretation.

One or more embodiments of the invention provide a Scholte wave frequency dispersion curve adaptive inversion method and a related device, and the method comprises the steps: obtaining and processing the original seismic data of a target seabed, and obtaining a Scholte wave frequency dispersion curve; constructing an initial shear wave velocity model through linear inversion based on a Schlte wave frequency dispersion curve; and performing fitting optimization on the Schlte wave frequency dispersion curve based on the initial shear wave velocity model to obtain a two-dimensional shear wave velocity profile. According to the technical scheme, the problem that a traditional linear inversion method depends on an initial model is solved, the calculation efficiency is remarkably improved, the problem of local extremum is avoided, and an automatic high-precision inversion solution suitable for batch data processing is provided for submarine seismic exploration.

The invention provides a three-dimensional seismic data mixed noise suppression method based on an MSAT-Unet. The method comprises the following specific steps: constructing an MSAT-Unet network comprising a multi-scale expansion convolution residual module, a channel-space attention mechanism and a Transform convolution module; the encoder is improved into multi-scale expansion residual convolution, so that the receptive field is expanded, and the capability of capturing local details and global semantic information is enhanced; a channel-space attention mechanism is integrated behind the decoder, a direction sensitive context is extracted through multi-dimensional adaptive pooling, and details and edge recovery are enhanced; meanwhile, an improved decoder is a Transform convolution module, and the feature reconstruction and complex structure recovery capability is improved; the optimized model carries out training and reasoning on the three-dimensional seismic data, mixed noise can be effectively suppressed, a clear data basis is provided for subsequent interpretation, and the method is high in generalization and robustness and good in performance in the aspect of three-dimensional seismic data denoising.

The invention discloses a seismic wave velocity model inversion method based on time-space coupling deep learning, and the method achieves a time-space coupling mechanism through a Fourier neural operator framework of a wavelet transform-attention mechanism, and supports the reconstruction of a spatial two-dimensional velocity model from a time sequence track of multiple seismic sources. Firstly, the local time-frequency representation of each seismic time-domain trajectory is extracted by adopting continuous wavelet transform, and then a wavelet time-frequency spectrogram is coded through a convolution feature extractor. Then, information from multiple sources is fused at each receiver point location through a multi-head attention mechanism, enhancing correlation between multiple sources. And finally, projecting the fused features into a two-dimensional Fourier neural operator to efficiently learn space mapping and reconstruct a bottom-layer velocity model based on sound velocity. According to the method, the time sequence characteristic and the spatial distribution characteristic of the seismic channel are fully utilized, the precision and robustness of an inversion result are improved, and a more reliable speed model is provided for geological structure interpretation and acoustic imaging.

The invention relates to a seabed elastic parameter inversion method, system and device under an earthquake spherical reflection wave oscillation-gradual change nuclear decoupling framework, and belongs to the field of geophysical detection and acoustic detection. The method comprises the following steps: (1) inputting observation data and initializing seabed elastic parameters; (2) determining an integral interval and decomposing the integral interval into a plurality of sub-intervals; (3) calculating an oscillation nuclear local integral matrix corresponding to the incident angle of the observation data; (4) calculating a slowly-varying kernel local mean vector; (5) calculating a spherical wave reflection coefficient of oscillation-gradual change kernel decoupling; and (6) constructing a target function by using the decoupled spherical wave reflection coefficient, and solving the target function. In the objective function solving process, only the steps (4) and (5) need to be repeatedly executed after the seabed elastic parameters are updated, the oscillation nuclear local integral matrix in the step (3) does not need to be repeatedly calculated, redundant calculation about oscillation nuclear integration is avoided, and the calculation efficiency is improved.

The invention relates to the technical field of seismic exploration, in particular to a seismic offset angle gather generation method and device, equipment and a medium. The method comprises the following steps: acquiring pre-stack seismic data of a work area after preprocessing; performing pre-stack time migration on the pre-stack seismic data, and performing grouping accumulation according to inclination angles in the pre-stack time migration process to generate an inclination angle gather; according to the inclination angle gather, picking up Fresnel zone boundaries of reflection events of all imaging points according to preset intervals; stable-phase offset is carried out by using the picked Fresnel zone boundary, and an offset angle gather of stable-phase control is generated according to grouping accumulation of incident angles in the offset process. The problems that an existing method for generating an angle gather through pre-stack depth migration is difficult in interval velocity modeling and high in calculation cost are solved.

The invention belongs to the technical field of underground engineering and computer-aided engineering analysis, and particularly relates to an anti-seismic safety collaborative analysis method for multi-modal data of proximity tunnel engineering. By integrating engineering geology, structural design, construction monitoring and multi-level seismic oscillation data, a three-dimensional numerical model fusing geology and structural information is constructed. Firstly, refined simulation of the whole construction process and the operation static load is carried out, and a reference mechanical state is obtained; carrying out multi-task condition earthquake dynamic time history analysis on the basis, dynamically capturing dynamic response of surrounding rock and a structure, carrying out collaborative interpretation on static and dynamic analysis results, emphatically analyzing spatio-temporal evolution and connectivity risks of a tunnel anchor and an adjacent highway tunnel surrounding rock plastic zone, and extracting a plurality of quantitative collaborative indexes; and an interactive analysis report integrating a three-dimensional visual scene, a multi-dimensional comparison chart and a rule base-based safety interpretation conclusion is automatically generated, and systematic decision support is provided for seismic design, construction control and safety evaluation of proximity projects.