121 results about "Wave equation" patented technology

The invention provides a frequency dispersion curve calculation method based on a decoupling type Legendre-Gauss-Labatto spectral element method, and belongs to the technical field of discrete solution of guided wave propagation in a thin plate, and the method comprises the following steps: S1, constructing a wave equation which comprises a Cartesian space coordinate system based on an orthotropic plate, and combining a constitutive equation, a geometrical relationship under the coordinate system and guided wave displacement to obtain a wave dispersion curve; substituting into a control equation when the volume force is zero to obtain an improved three-dimensional elastic wave equation; s2, a differential matrix is constructed, wherein discretization and decoupling are carried out on a wave equation on the basis of a Lagrange interpolation basis function of a Legendre-Gauss-Labatto node; s3, explicitly embedding the stress free boundary condition into the solving process through matrix line transformation, and solving to obtain a frequency dispersion curve, the derivation process is simplified, the boundary and the interior of the material are decoupled, and the physical significance is clearer. And meanwhile, the high-precision advantage of the spectrum method is maintained, and the calculation rate is remarkably improved.

A seismic energy estimation method based on generalized spherical wave propagation, relating to the field of earthquake engineering. The method comprises: determining initial seismic total energy and seismic power on the basis of seismic source characteristics and a seismic magnitude; determining a generalized spherical wave equation, and calculating an average energy flow density of propagated waves on the basis of the wave equation of spherical waves; using an earth surface spherical crown, which is projected on the zone of a seismic source surface, as a main impact zone of seismic spherical waves on an earth surface, the main impact zone being equivalent to a square and performing grid division, and calculating the hypocentral distance between the center of each grid and the center of a seismic source so as to form a hypocentral distance matrix; and expressing the generalized spherical wave energy flow density of each grid as the sum of the energy flow density of longitudinal waves and the energy flow density of transverse waves, performing grid subdivision, estimating the average amplitude of each subdivided grid mass unit within a seismic action period, and evaluating the seismic action energy of a site. The present invention estimates the action of seismic energy on the earth surface by defining the generalized spherical wave propagation of seismic source energy, thereby more accurately evaluating the distribution condition of seismic energy on the earth surface.

The invention discloses an elastic reverse time migration imaging method based on deviatoric stress characterization, and relates to the technical field of exploration geophysics. Deviatoric stress is introduced into a conventional elastic wave equation, and a speed-deviatoric stress elastic wave equation describing stress propagation characteristics is obtained through derivation; based on Hooke's law, P wave contribution in total partial strain is removed through partial strain decomposition, decoupling of P-S wave stress components is achieved, and correctness of amplitude, phase and physical unit of a decoupling wave field is guaranteed; and then a joint imaging condition with the average stress as a seismic source wave field and the multi-stress component as a receiving wave field is constructed, PP waves are subjected to normalization superposition through the average stress and the deviatoric stress, and PS waves are subjected to shear stress imaging. According to the method, the sensitivity to heterogeneity and local stress variation of underground media can be remarkably improved, P-S wave crosstalk and imaging noise are effectively eliminated, the imaging precision and resolution of a complex geologic structure are greatly improved, and a reliable technical means is provided for high-precision seismic exploration of complex oil and gas reservoirs.

The invention relates to a stable Q compensation reverse time migration method for seismic exploration data, a medium and equipment, and belongs to the technical field of geophysical exploration data processing. According to the method, based on a fractional order Laplace viscous acoustic wave equation, a physically guided time-varying regularization compensation operator is constructed in a wavenumber domain, the energy amplitude of a previous time step length is introduced as a reference, an energy growth monitoring mechanism is established, and the energy growth is monitored. And designing a three-level adaptive constraint system: local energy growth constraint, high-frequency adaptive attenuation and boundary region stabilization processing, and realizing adjustment of excessive compensation and adaptive suppression of noise in the wave field continuation process. According to the method, amplitude attenuation and phase frequency dispersion caused by stratum absorption are effectively recovered, meanwhile, the stability of numerical calculation can be guaranteed, the precision and resolution of deep structure imaging are remarkably improved, and good numerical stability and anti-noise capacity are achieved.

The application provides a warning method and system for the tipping moment of a thermoacoustic system under Lévy noise, and relates to the technical field of thermoacoustic system stability analysis and early warning control, which comprises the following steps: based on the physical model of the thermoacoustic system, a non-uniform wave equation describing the pressure and velocity fluctuation of the system is constructed; the Galerkin modal expansion and Taylor expansion truncation are performed on the equation, and a random differential equation model is established by combining the time-varying control parameters and Lévy noise; the random differential equation model is solved by the Monte Carlo simulation and the random fourth-order Runge-Kutta numerical method, and the amplitude sequence of the thermoacoustic system changing with time is obtained; based on the probability distribution of the amplitude sequence, the change of the system dynamics characteristics is quantified by calculating the Shannon entropy, and the sequence of the Shannon entropy changing with time is obtained; the warning moment when the Shannon entropy sequence first passes through the preset threshold is determined, and the warning domain is generated according to the time difference between the warning moment and the actual tipping moment of the thermoacoustic system, so that the warning for the tipping moment of the thermoacoustic system is realized.

The invention discloses a seismic wave field simulation method and system in an attenuation medium. A new constant Q viscous sound wave equation is constructed to accurately meet the dispersion relation of the Kjartansson constant Q theory, and an amplitude attenuation item and a phase dispersion item are subjected to explicit separation, so that complete decoupling of attenuation and dispersion can be realized while the theoretically accurate constant Q dispersion relation is kept, and the frequency dispersion accuracy is improved. And solving a constant Q viscous acoustic wave equation through a generalized pseudo-spectral numerical simulation method based on space Fourier transform to obtain a seismic wave field with spatio-temporal variation, so that the precision of viscous acoustic wave simulation is remarkably improved.

This invention relates to the field of seismic imaging technology and discloses a method, apparatus, equipment, and medium for forward modeling of seismic acoustic waves in a viscoelastic-acoustic (VTI) medium. The method includes: obtaining a desired Q value; calculating the viscoelastic velocity of the seismic acoustic wave based on the desired Q value; fitting parameters of a generalized standard linear body model; constructing a wave equation for the seismic acoustic wave in the VTI medium based on the generalized standard linear body model; and solving the wave equation for the seismic acoustic wave in the VTI medium using the finite difference method and the viscoelastic velocity. This invention solves the problem that traditional seismic simulations cannot simultaneously consider seismic anisotropy and seismic attenuation, significantly improving the accuracy and applicability of seismic simulations, as well as the accuracy of seismic simulations in complex terrain, geology, and oil and gas reservoir areas. It provides a foundation for subsequent high-precision seismic velocity modeling and seismic migration imaging.

The invention belongs to the technical field of petroleum geophysical exploration, and relates to an overlapped meshless undulating surface seismic wave field numerical simulation method and system. The method comprises the following steps: establishing an undulating surface seismic wave speed model and an absorption attenuation Q model based on geological parameters of an exploration area; based on the undulating surface seismic wave speed model and the absorption attenuation Q model, adopting an overlapping meshless node subdivision method to obtain a node data set of an overlapping region; taking the node data set as input, expanding a radial basis function based on a Chebyshev polynomial and forming a difference matrix, and solving to obtain a node difference coefficient set; and taking the node data set and the node difference coefficient set as input, performing finite difference discretization on a Kelvin-Voigt model-based viscous-acoustic wave equation based on a radial basis function, implementing forward continuation of the wave field, and completing the seismic wave field numerical simulation of the overlapped meshless undulating surface. According to the method, the constant Q model is accurately described, amplitude-phase decoupling is realized, and the absorption attenuation effect of the medium can be effectively described.

A method and an apparatus for measuring ultrasonic viscoelasticity, a device and a medium are provided. The method includes: establishing an inhomogeneous viscoelastic wave equation; coding the inhomogeneous viscoelastic wave equation as a loss function, and constructing a physics-informed neural network, where the spatial-temporal neural network is configured to predict a stream function, the inhomogeneous viscoelastic wave equation is used to obtain vertical vibration velocity of particles based on the stream function outputted, the spatial neural network is configured to predict shear modulus and viscosity, and the loss function is used to calculate an aggregate loss; inputting the multi-frame vertical vibration velocity of particles of the object to be measured as supervisory data into the physics-informed neural network, to obtain a spatial distribution of the shear modulus and a spatial distribution of the viscosity of the object to be measured.

The present application belongs to the technical field of oil and gas seismic exploration, and particularly relates to a pre-stack depth migration imaging method based on generalized diffraction stacking. The present application uses one-way wave equation to calculate the Green's function required for imaging, thereby avoiding the defects of two-way wave equation, making the imaging result not containing the wave field component useless for imaging, ensuring good imaging effect, and reducing the calculation amount, accelerating the operation speed, and improving the imaging efficiency. Moreover, the Green's function is solved in the frequency domain, and there is no sampling rate dispersion problem in the time domain, so that the calculation mode of one-way wave equation can adapt to larger spatial grid interval. Moreover, the relative positions of the imaging point, the source point and the receiver point are judged through angle calculation, so as to perform regional decomposition on the velocity model, and only the necessary region needs to be recursively calculated when solving the Green's function, without recursively calculating in the whole velocity model space. On the basis of ensuring the recursive effect, the useless calculation region is removed, and the imaging efficiency is further improved.

This disclosure pertains to the field of exploration geophysics and presents a method and system for combined active and passive source reverse time migration imaging of metallic minerals. The method includes: performing time difference correction on active and passive source seismic data from the same subsurface structure; performing wavefield separation on the time-aligned active and passive source seismic data to reconstruct the P-wave and S-wave fields of the active source and the passive source; performing active and passive source wavefield extrapolation based on the two-way wave equation; obtaining active and passive source PS-converted wave imaging results using cross-correlation imaging conditions; and performing weighted superposition and fusion to obtain the combined active and passive source imaging results. This disclosure leverages the high resolution advantage of active sources combined with the deep detection capability of passive sources, replacing the limitations of single-source detection and achieving high-precision deep reverse time migration imaging.

The embodiment of the specification discloses a geothermal rock mass elastic wave reverse time migration method based on adaptive variable grid. By inputting the migration velocity model and seismic record of the geothermal rock mass, and setting the imaging parameters, the seismic wave is generated according to the migration velocity model, the seismic record and the imaging parameters; the migration velocity model is adaptively processed by variable grid, the mapping relationship between the coordinate systems is determined, and the wave equation under the variable grid is determined based on the mapping relationship; according to the wave equation, the reverse time continuation of the wave field of the receiver point wave field is carried out in the variable grid along the time reverse propagation direction, and the reverse continuation wave field is generated; the cross-correlation imaging condition is used for the forward continuation wave field and the reverse continuation wave field, the reverse time migration profile under the variable grid is generated, and the imaging result under the variable grid is obtained by stacking; the imaging result is reversely processed by the variable grid according to the mapping relationship, and the reverse time migration imaging result in the conventional coordinate system is obtained; thereby the memory occupation is reduced, and the calculation is fast.

The invention discloses a microseism synchronous inversion method based on multi-task learning, and belongs to the field of microseism inversion, and the method comprises the steps: firstly constructing a ternary data set corresponding to seismic data, a speed model and seismic source position probability distribution; then a double-branch U-Net network model is designed, a shared encoder of the double-branch U-Net network model is integrated with a multi-scale convolution module to extract multi-scale features, and double decoders realize feature decoupling and reinforcement through a task-adaptive guided attention gate; meanwhile, based on wave equation Frechet derivative analysis, a physically-driven dynamic weight adjustment strategy is provided, and optimization of focus positioning and speed modeling tasks is adaptively balanced in the training process. According to the method, the seismic source position and the underground velocity model can be synchronously and precisely inverted from the microseismic data, the problems of parameter coupling, multiplicity, unstable training and the like in a traditional method are effectively solved, and the inversion precision and robustness are remarkably improved.

The invention provides a foundation pile low-strain reflection wave method test method, and belongs to the technical field of pile foundation detection. The method comprises the following operation steps: S1, building a column-bearing platform-pile soil system three-dimensional dynamic model under a rectangular coordinate system according to an elastic theory; s2, a pile, a bearing platform, a column, soil around the pile and soil at the bottom of the pile are regarded as isotropic linear elastomers, physical power is not considered, and system dynamic response is solved based on an elastic wave equation; s3, analyzing the influence of different excitation force positions and different vibration pickup positions on the speed response curve; s4, verifying the correctness of the theoretical model by comparing the numerical calculation curve with the actual measurement curve, and applying the theoretical model to engineering practice to assess the completeness of the foundation pile. Through three-dimensional theoretical innovation, test technology optimization and engineering verification, the accuracy, reliability and practicability of low-strain detection of the column-bearing platform-pile system are comprehensively improved; an effective technical solution is provided for foundation pile integrity evaluation under complex conditions, and the method has important theoretical value and engineering application prospects.

The invention discloses an ARM (Advanced RISC Machines) architecture-oriented CGFDM (China General Microbiological Fused Deposition Modeling) seismic wave solver for hybrid precision and vectorization collaborative optimization and a simulation method, and aims to break through the bottleneck of memory bandwidth and calculation efficiency of ARM platform seismic simulation. According to the method, firstly, an elastic wave equation is reconstructed in a dimensionless mode, key physical quantities are made to adapt to the FP16 range, the mixed precision strategy of FP16 storage and FP32 calculation is adopted, and on the premise that precision is guaranteed, about 50% of memory consumption is reduced. And secondly, for ARM scalable vector expansion (SVE), the CGFDM differential template is mapped to a vector register, SIMD parallel computing is realized, and the throughput is remarkably improved. Finally, the FP16-SVE solver realizes data rearrangement and conversion between the FP16 and the FP32, the problems of alignment and efficiency in mixed precision operation are solved, and synchronous improvement of storage and calculation efficiency is realized. Tests show that the method has the advantages that while the precision is maintained, the memory is halved, nearly three times of speed-up ratio is obtained, and an efficient scheme is provided for large-scale earthquake simulation and real-time disaster assessment of an ARM platform.

The application discloses an ultrasonic CT waveform inversion method and device based on a dynamic aperture, and belongs to the technical field of ultrasonic image processing, and comprises the following steps: acquiring an observed sound pressure signal at a position of a receiving element of an ultrasonic CT data acquisition system; inputting a current sound velocity distribution into a wave equation to obtain a synthesized sound pressure signal at the position of the receiving element; for each transmitting element, extracting the synthesized sound pressure signals at the positions of M receiving elements away from the transmitting element; calculating the difference between the synthesized sound pressure signals at the positions of the M receiving elements corresponding to all the transmitting elements and the corresponding observed sound pressure signals; calculating the gradient of the current sound velocity distribution by using an adjoint state method according to the difference, and then obtaining an updated sound velocity distribution; re-substituting the updated sound velocity distribution into the wave equation for iteration, and simultaneously increasing the value of M until the iteration is stopped, and outputting the updated sound velocity distribution. The application can improve the accuracy of a reconstructed image under the condition of missing low-frequency information.

The invention discloses a three-phase current-based power transmission line ground fault phase identification method, system and device and a medium, and belongs to the technical field of ground fault identification, and the method comprises the steps: calculating a virtual space derivative of a zero-sequence current according to a time derivative of a three-phase current transient component and a line distributed capacitance parameter, and forming a space-time coupling tensor; performing spatial modal decomposition by adopting an intrinsic function of the zero-sequence traveling wave equation to obtain spatial modal coefficients of all orders; constructing a phase sequence leakage matrix by taking the spatial modal coefficient as a weight; performing singular value decomposition on the phase sequence leakage matrix, and extracting main singular vector feature parameters; transforming the composite mapping matrix to obtain an island vector; determining a suspected fault phase according to the projection coefficient energy ratio of the island vector; and reversely solving the distributed capacitance and calculating a theoretical value, and outputting a final fault phase when the deviation between the theoretical value and the measured value meets a condition. According to the method, accurate fault phase identification under high-resistance grounding and long-distance line working conditions is realized through virtual space derivative mapping and modal energy decomposition.

Systems and methods are disclosed. The method includes, for each of a plurality of spatial partial differential equations (sPDEs) within a wave equation, determining a linear system of equations using an approximate solution at a plurality of grid nodes. The linear system of equations includes an inverse matrix, a first vector, and a second vector. The method further includes, for each of the plurality of sPDEs, evaluating the inverse matrix by evaluating a first portion of the inverse matrix using a first deconvolution filter and evaluating a second portion of the inverse matrix using a second deconvolution filter. The method further still includes, for each of the plurality of sPDEs, evaluating the first vector using the evaluated inverse matrix and the second vector as well as determining the wavefield using the evaluated first vector for each of the plurality of sPDEs, a velocity model, and the wave equation.

The invention discloses a ground penetrating radar reverse time migration imaging method, system and device based on local Nyquist attenuation compensation and a medium. The method comprises the steps that ground penetrating radar observation data are collected, and a speed model is constructed; calculating to obtain a local time window; executing forward propagation wave field simulation to obtain a sparse forward propagation wave field; executing back propagation wave field simulation according to the wave equation, and performing symbol inversion on a conductivity item in the wave equation to obtain a compensation back propagation wave field; performing cross-correlation calculation to obtain an imaging value; and generating a reverse time migration imaging profile. According to the invention, through cooperation of triple mechanisms of a local time window, Nyquist sparse sampling and attenuation compensation operation, the wave field storage demand and the calculation burden are reduced while the imaging precision is maintained; and the energy attenuation of electromagnetic waves in the conducting medium is compensated, the imaging energy of deep reflection signals is enhanced, the technical problems of high storage cost and insufficient deep imaging capability of a traditional method are solved, and a solution is provided for high-precision ground penetrating radar imaging under complex geological conditions.

The application discloses a method for pre-stack seismic velocity inversion of natural gas hydrate reservoir based on deep learning. According to the occurrence characteristics of natural gas hydrate and free gas under actual geological conditions, combined with the real seabed layered sedimentary area and the public ocean velocity model, a synthetic velocity model of thin layer associated structure containing high and low speed alternating hydrate and free gas is constructed; the model is forward by solving the wave equation to generate forward seismic records; the forward seismic records and the synthetic velocity model are registered, and the training, testing and verification sets are divided; the deep neural network is constructed, and the training environment and hyperparameters are configured; the perception loss is introduced, the mixed loss function is designed, the network training is multi-dimensionally constrained, the network is optimized and the optimal parameters are saved; the optimal parameters are called to batch inversion of seismic records, and high-precision velocity inversion and thin layer identification are realized. Through the above method, the inversion efficiency and reservoir characterization ability are improved, and effective technical means are provided for natural gas hydrate exploration.

The invention provides a geophysical exploration reverse time migration imaging calculation method. According to the method, a seismic wave field computational domain is divided into parallelogram space-time slices by introducing a cache precise time tilt strategy, and an index table is constructed to record the position, size and dependency relationship of each slice. During calculation, seismic wave field data are carried to a unified buffer area from a global memory in batches according to an index table, efficient parallel calculation is carried out by utilizing a matrix calculation unit, a wave equation is solved by adopting a finite difference format, forward propagation of a seismic source wave field and reverse propagation of a received wave field are completed, and finally an imaging result is written back to the global memory.

The application provides a wave equation full-wave Q tomography method, which comprises the following steps: performing FFT transformation on existing seismic records to determine the frequency of a Ricker wavelet; performing forward simulation using the frequency of the Ricker wavelet determined in step 1 to obtain seismic correlation data; converting the seismic correlation data to a local domain and performing FFT on the data in the local domain to obtain the frequency of a single seismic event; obtaining the peak frequency shift of each seismic event in the local domain according to the frequency of the single seismic event; obtaining a companion source according to the peak frequency shift; performing wave field continuation calculation according to the companion source to obtain a back-propagating wave field; obtaining a gradient using the back-propagating wave field; and performing iterative updating using a conjugate gradient method according to the gradient to finally obtain an updated Q model. The method provided by the application can be simulated to obtain a relatively accurate Q model.

The application discloses a deep-water shallow seismic data structure-constrained attenuation compensation velocity modeling method, which comprises the following steps: obtaining an initial velocity model by using a tomographic inversion method and performing reverse-time migration imaging to obtain an initial imaging result; and obtaining a seismic data dip angle field by using a dip angle prediction technology; a structure-constrained velocity modeling target functional is established by using the dip angle field and multi-frequency band seismic data; wave field continuation for Q compensation is performed by using a decoupled viscous wave equation, and the gradient of the functional is calculated to update the velocity model; a multi-scale multi-frequency band inversion strategy is used to solve a high-precision velocity inversion optimization problem to obtain a high-precision deep-water shallow velocity field. By applying attenuation compensation to full waveform inversion, a structure guide constraint is introduced into a full waveform inversion target function, and a velocity parameter model is constrained by using seismic imaging structure information, so that the precision and efficiency of inversion can be effectively improved.

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.

A computer-implemented method for determining modeled seismic data of a subsurface region includes receiving observed seismic data from the subsurface region captured by one or more seismic receivers as one or more reflected seismic signals; determining, based on a velocity model of the subsurface region and the observed seismic data, a traveltime function associated with the subsurface region; migrating, using the traveltime function and a Kirchhoff migration operator, the observed seismic data to produce a migrated seismic image; and performing reverse-time demigration on the migrated seismic image, using the velocity model and a solved full wave-equation, to produce modeled seismic data.

The present application relates to a kind of back wave field synthesis processing method, device, computer equipment and storage medium, method includes the following steps: respectively obtaining source back wave data and source wavelet data;According to source back wave data and source wavelet data, determine sparse linear filter coefficient;Using the numerical value of non-zero sparse linear filter coefficient and delay, time delay weighted superposition of forward wave field, synthesis back wave field.The above-mentioned back wave field synthesis processing method, based on source back data is the convolution of source wavelet data and sparse linear filter, determine sparse linear filter coefficient by source back data and source wavelet data, then according to the principle of superposition, time delay weighted superposition of forward wave field, synthesis back wave field.No need to rely on back data to obtain back wave field by solving wave equation, power is greatly reduced, can improve processing efficiency, reduce processing time, can also reduce memory consumption, so as to reduce power cost.

A computer-implemented method for determining modeled seismic data of a subsurface region includes receiving observed seismic data from the subsurface region captured by one or more seismic receivers as one or more reflected seismic signals; determining, based on a velocity model of the subsurface region and the observed seismic data, a traveltime function associated with the subsurface region; migrating, using the traveltime function and a Kirchhoff migration operator, the observed seismic data to produce a migrated seismic image; and performing reverse-time demigration on the migrated seismic image, using the velocity model and a solved full wave-equation, to produce modeled seismic data.

The invention discloses a VSP velocity modeling method based on a Stein variational gradient descent method, and relates to the technical field of geophysical exploration, and the method comprises the steps: S1, generating a certain number of random particles as an initial velocity model based on Gaussian random field disturbance; s2, performing wave equation forward modeling simulation on the initial model to obtain synthetic seismic data; s3, calculating the optimal updating direction of the particles in the regeneration kernel Hilbert space according to the observed and synthesized seismic data; s4, updating the particles based on the optimal updating direction calculated in the step S3; and S5, repeating the steps S2 to S4 until the maximum number of iterations is reached, introducing a Stein variation strategy, replacing random sampling with optimized particles, providing technical support for VSP speed modeling, solving the problem of dependence of a traditional full-waveform inversion method on an initial model, improving the calculation efficiency of the inversion method, and quantifying the uncertainty of an inversion result at the same time.

This disclosure relates to the field of seismic imaging technology, and particularly to a method, apparatus, device, and medium for generating common imaging point gathers for wave equations. The method includes: identifying a source wavelet based on a selected source location; synthesizing the return wavefield of a single-channel seismic record corresponding to a geophone based on the forward propagation wavefield of the source wavelet and a pre-determined linear filter; synthesizing a single-channel seismic image corresponding to the single-channel seismic record using the return wavefield and the forward propagation wavefield; calculating the offset distance corresponding to each single-channel seismic image; classifying the single-channel seismic images according to the offset distance; and superimposing the classified single-channel seismic images to obtain an offset common imaging point gather; performing amplitude processing on the offset common imaging point gathers; and determining the amplitude-processed offset imaging point gathers as the common imaging point gathers for wave equations. This disclosure can improve the efficiency of common imaging point gather generation.

The present application relates to a kind of stable Q compensation reverse time migration method, medium and equipment for seismic exploration data, belong to geophysical exploration data processing technical field.The method is based on fractional order Laplace viscous wave wave equation, constructs physical guide time-varying regularization compensation operator in wave number domain, by introducing the energy amplitude of previous time step as reference benchmark, establishes energy growth monitoring mechanism, and designs three-level adaptive constraint system: local energy growth constraint, high-frequency adaptive attenuation and boundary region stabilization processing, in wave field continuation process, realize the adjustment of overcompensation and the adaptive suppression of noise.The method can effectively restore the amplitude attenuation and phase dispersion caused by stratum absorption while ensuring the stability of numerical calculation, significantly improve the accuracy and resolution of deep structure imaging, with good numerical stability and noise resistance.