283 results about "Wave equation" patented technology

According to the micro-seismic source intelligent positioning method and the model training method and system, physical laws such as a seismic wave equation are integrated, physical loss and data loss are comprehensively considered, the adaptability of the model to a complex geological environment is improved, micro-seismic source positioning challenges in the complex geological environment can be effectively dealt with, and the positioning accuracy of the micro-seismic source in the complex geological environment is improved. Noise immunity and positioning precision are improved, and more reliable technical support is provided for safety monitoring of underground engineering.

The invention discloses a microseismic source positioning method, device and system, and a storage medium. The method comprises the following steps: dividing a microseismic data sample data set into a training set and a test set; according to the training set, constructing a microseism inversion subnet containing a Swin Transform encoder, and according to the training set, constructing a microseism inversion subnet containing a Swin Transform encoder; the micro-seismic forward modeling subnet realizes seismic wave field continuation by inputting a micro-seismic source position and a speed model and utilizing a recurrent neural network structure and a convolution operator, and establishes a forward modeling subnet based on a wave equation; constructing an inversion-forward modeling closed-loop neural network according to the micro-seismic inversion subnet and the forward modeling subnet based on the wave equation; and inputting the test set into the inversion-forward closed-loop neural network to carry out micro-seismic source positioning. By adopting the technical scheme of the invention, the limitations on physical constraint, feature modeling and anti-noise capability in the prior art are overcome.

The invention provides a multiple longitudinal and transverse wave joint detection method based on a seabed crawler, which belongs to the technical field of geophysical exploration, and constructs an original longitudinal and transverse wave signal matrix by transmitting multi-frequency longitudinal and transverse wave signals and receiving stratum reflection signals. Fourier transform is utilized to extract frequency characteristics, physical analysis is carried out in combination with a wave equation, and a longitudinal and transverse wave interference matrix is constructed and weighted fusion is carried out. And constructing an error compensation matrix according to the crawl device attitude data, and performing interference compensation on the joint interference matrix. And performing time-frequency analysis through wavelet transform, extracting specific time-frequency features, inputting the specific time-frequency features into a pre-trained stratum classification prediction model, realizing anti-interference prediction on the type and thickness of the seabed stratum structure, finally drawing a seabed three-dimensional geological profile map according to a prediction result, and marking a key geological structure and a resource distribution region. The technical problem that accurate recognition and resource distribution positioning of the seabed stratum structure are affected due to much external interference in the seabed geological exploration process is solved.

The invention discloses a multi-parameter optimization method and system for a curved surface acoustic wave device. The method comprises the following steps: constructing theoretical models of the surface acoustic wave device under different strain conditions; solving the coupling wave equation to obtain the propagation velocity of the surface acoustic wave in the piezoelectric material under different strain conditions; acquiring at least one structural parameter of the curved surface acoustic wave device associated with the propagation speed; and performing global optimization on the at least one structure parameter according to a preset global topological optimization network to obtain a target structure parameter corresponding to the at least one structure parameter. A global topological optimization network is adopted, gradient back propagation is utilized, structural parameters of a device are intelligently designed, a quality factor value of the device is generated through finite element simulation calculation, and finally the quality factor distribution condition of the generated device is output. The problem that accurate real-time measurement of multiple parameters of the bent surface of the bearing cannot be achieved in the prior art can be solved.

The invention discloses a three-dimensional intelligent seismic velocity modeling method based on a two-dimensional self-attention U-shaped network and low-rank adaptation, and belongs to the technical field of seismic services of geophysical exploration. The problems that in the prior art, a traditional three-dimensional intelligent seismic velocity modeling method is low in efficiency and insufficient in model generalization ability are solved. The method comprises the following steps: generating a three-dimensional hybrid geological velocity model, obtaining multi-shot seismic record data based on a wave equation, obtaining a three-dimensional initial velocity model through Gaussian filtering, and obtaining an imaging result by adopting reverse time migration; carrying out dimension reduction processing on the three-dimensional data, and constructing a multi-dimensional pre-training sample set; carrying out pre-training by utilizing the training set to obtain an optimized LAUNet framework; and performing fine tuning by adopting a dual-stage collaborative optimization mechanism, and obtaining a fine-tuned LAUNet framework through a LoRA * method. The method effectively improves the model inversion efficiency, and can be applied to the establishment of a three-dimensional intelligent seismic velocity model.

The invention discloses a steel plate internal damage detection method based on guided waves and PINNs, and the method comprises the following steps: firstly establishing a rectangular plane coordinate system, selecting a guided wave application position and a guided wave receiving point on the surface of a steel plate, applying guided waves to the steel plate, and selecting M time points of free vibration of the steel plate, the method comprises the following steps: carrying out guided wave receiving on a steel plate, calculating out-of-plane displacements of the steel plate at all guided wave receiving points, constructing a PINNs model and a DNN model, training the PINNs model and the DNN model at the same time by using the out-of-plane displacements, coordinates and time of the steel plate at all guided wave receiving points, detecting the steel plate by using the PINNs model and the DNN model, and positioning the damage position of the steel plate. According to the method, the wave equation is approximately solved through the physical information neural network, a complex finite difference method formula is avoided, so that the calculation cost is reduced, compared with a three-dimensional wave equation, the adopted physical equation based on the plate-shell theory is simpler, the calculation cost is reduced, the accuracy of positioning the internal damage of the steel plate is improved, and the accuracy of positioning the internal damage of the steel plate is improved. The internal damage position of the steel plate can be positioned.

The invention discloses a method for generating a reflection angle domain point spread function based on an LSTM network, and belongs to the technical field of geophysical exploration, and the method comprises the following steps: S1, selecting a reference point, and constructing a label through ray tracing and wave equation forward modeling; s2, training an LSTM (Long Short Term Memory) network as a ray tracing agent model of the whole space of an imaging domain; s3, predicting a travel time gradient and an amplitude factor through a ray tracing agent model; s4, combining the seismic source and the geophone in pairs, and calculating a reflection angle through the travel time gradient; s5, solving an illumination wave number vector, and solving a wave number domain point spread function in combination with the illumination intensity; s6, performing inverse Fourier transform on the wavenumber domain point spread function to obtain a spatial domain point spread function; and S7, solving a deconvolution problem to optimize the angle gather. According to the method, the calculation efficiency of solving the Green function through seismic wave ray tracing can be improved, the imaging response described by the point spread function is eliminated from the imaging result, and the balance of the resolution of the imaging result and the illumination intensity can be improved.

The invention discloses a least square reverse time migration method and device for a three-dimensional anisotropic medium, and the method comprises the steps: building a forward equation of the three-dimensional anisotropic medium for a seismic block to be subjected to reverse time migration, and obtaining a background wave field under the three-dimensional anisotropic medium through regional forward modeling; determining a three-dimensional anisotropic medium adjoint wave field by using a pre-established wave equation of the adjoint wave field; according to the adjoint wave field and the background wave field, determining an initial gradient by using a pre-established least square reverse time migration cross-correlation equation, and obtaining a migration result according to the initial gradient; and updating an imaging result through a conjugate gradient method according to the offset result to obtain a reflected wave record.

The invention discloses a calculation method and a calculation device of an earthquake propagation direction. The calculation method comprises the steps that a wave equation of seismic time-space domain sound waves is acquired, and wave field information is acquired according to the wave equation; calculating a wave field component according to the wave field information; and substituting the wave field component and the wave field information into a pre-constructed seismic wave field propagation direction vector expression. According to the technical scheme, the situation that deviation occurs in the seismic wave propagation direction analysis can be effectively reduced, and the calculation precision is improved.

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

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

The present application belongs to the technical field of oil and gas exploration and seismic data processing, and discloses a seismic data interpolation method and system based on CBAM-Res2Unet network. The method solves the problems of limited receptive field caused by fixed convolution kernel in traditional ResUNet network and non-selective transmission of low-level features in skip connection. The method comprises: generating seismic records through wave equation forward modeling and constructing a training set containing random missing data; designing CBAM-Res2Unet network, wherein the encoder and the decoder adopt Res2Net modules containing compression excitation modules, and the multi-scale feature extraction capability is enhanced through grouped convolution and hierarchical residual connection; introducing CBAM modules in the skip connection to adaptively weight the features by using channel and spatial attention mechanisms; and adding a spectral normalization layer in the deep network to stabilize the training. The system comprises data acquisition, network generation, training and interpolation modules. Through multi-scale feature fusion, attention mechanism and stability optimization, the present application realizes high-precision interpolation of complex missing seismic data.

The embodiment of the specification discloses a self-adaptive variable grid forward modeling method for middle-deep dry hot rock mass. p , a transverse wave velocity V s and a medium density p model, pre-setting forward modeling parameters; performing self-adaptive variable grid processing on the velocity model, determining a mapping relationship between coordinate systems, and determining a wave equation under variable grid based on the mapping relationship; performing self-adaptive variable grid processing on the velocity model, determining a mapping relationship between coordinate systems, obtaining the underground medium longitudinal wave velocity, transverse wave velocity and density after variable grid, determining a wave equation under variable grid based on the mapping relationship, and generating final regular grid horizontal direction longitudinal wave field, horizontal direction transverse wave field, vertical direction longitudinal wave field and vertical direction transverse wave field information by using the mapping relationship reversely according to the wave field information in each direction at a specific moment after the variable grid, so that the memory occupation is small and the calculation speed is faster.

The invention belongs to the technical field of remote detection acoustic logging migration imaging, and particularly relates to a high-precision iteration reverse time migration method and system for remote detection acoustic logging. The method comprises the following steps: collecting remote detection acoustic logging data, and processing the remote detection acoustic logging data to obtain reflected wave actual measurement data; performing remote detection reverse time migration processing on the reflected wave actual measurement data to obtain an initial imaging result; establishing a wave equation based on the initial imaging result, and solving the wave equation to obtain reflected wave prediction data; obtaining residual data according to the reflected wave prediction data and the reflected wave actual measurement data; and carrying out wave field back propagation by taking the residual data as a geophone source, constructing an imaging gradient and obtaining an imaging result. According to the method, the reverse time migration is used for replacing Kirchhoff migration and F-K migration commonly used in remote exploration acoustic logging actual data processing, vertical speed change in a logging observation system can be better adapted, and low-dip-angle fracture-cavity migration homing is achieved.

The invention discloses a power line real-time identification and spatial positioning method based on depth vision assistance, and relates to the technical field of power line real-time identification, and the method comprises the steps: constructing a lead identification and positioning frame with a strong generalization capability through deep fusion of a physical law and a depth vision technology; a vibration propagation model based on a wave equation can accurately capture signal interruption characteristics at a cross point, and a frequency domain fusion mechanism significantly improves the linear structure recovery capability of a backlight scene; dynamic vibration phases are coded into a continuous radiation density function by adopting an implicit neural radiation field technology, and the dependence of traditional point cloud reconstruction on dense sampling is broken through; and a dynamic deformable convolution kernel is introduced to realize self-adaptive tracking of the wire form through a physical deformation model, so that extreme working conditions such as icing can be effectively handled.

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 belongs to the technical field of oil-gas exploration and seismic data processing, and discloses a seismic data interpolation method and system based on a CBAM-Res2Unet network. According to the method, the problems that the receptive field is limited and low-level characteristics in jump connection are not selectively transmitted due to a fixed convolution kernel in a traditional ResUNet network are solved. The method comprises the following steps: generating a seismic record through wave equation forward modeling, and constructing a training set containing random missing data; a CBAM-Res2Unet network is designed, an encoder and a decoder of the CBAM-Res2Unet network adopt Res2Net modules containing compression excitation modules, and the multi-scale feature extraction capability is enhanced through grouping convolution and layered residual connection; introducing a CBAM module in jump connection, and adaptively weighting features by using a channel and a space attention mechanism; a spectrum normalization layer is added to the deep layer of the network to stabilize training. The system comprises a data acquisition module, a network generation module, a training module and an interpolation module. According to the invention, through multi-scale feature fusion and attention mechanism and stability optimization, high-precision interpolation of complex missing seismic data is realized.

The invention relates to the field of exploration geophysics, and discloses a wave equation inversion method for constructing a middle-deep part background velocity field, which comprises the steps of S1, inputting observation gun record data and an initial velocity model, S2, performing forward modeling on the initial velocity field to calculate a forward propagation wave field, and performing back propagation on observation data to calculate a back propagation wave field; s3, calculating and simulating primary reflected wave data according to the migration imaging result and the initial velocity field back migration, solving a data residual error by using an optimal transmission distance, constructing an adjoint source, substituting the adjoint source into reflected wave inversion, calculating a gradient under the optimal transmission measure, and calculating a gradient under the optimal transmission measure; and S4, iteratively updating the background velocity field according to the calculated gradient, and performing full-waveform inversion by using the background velocity field information to obtain a final high-precision inversion result. Wave field residual errors are extracted more accurately based on an optimal transmission distance strategy, and background field gradient information more conforming to an underground structure is constructed.

The invention discloses a surface shallow creep detection method and system based on a time reversal technology, relates to the technical field of geological engineering, and solves the problem that the surface shallow creep detection precision is relatively low. The method comprises the following steps: acquiring multi-source physical field data of an earth surface shallow layer of a target area; establishing a simulation model of the surface shallow layer based on the multi-source physical field data; sending a first signal to the surface shallow layer through an antenna, and receiving a returned second signal; performing time reversal on the second signal to obtain a third signal; based on the simulation model, simulating a signal propagation path of a third signal sent by the antenna to the surface shallow layer to obtain target time inversion data; and comparing the target time inversion data with the historical time inversion data to obtain a creep abnormal region. Creep detection is carried out by adopting the time reversal technology, the creep abnormal area can be accurately positioned by utilizing the time reversibility of the wave equation and combining the self-consistent focusing characteristic of the time reversal signal at the real source point, and the creep detection precision is improved.

The invention relates to the technical field of numerical simulation, and discloses a seismic wave forward numerical simulation method, device and system and a storage medium, and the method comprises the steps: determining a first objective function based on a traditional regular grid finite difference format three-dimensional scalar wave equation; based on the first objective function, establishing a second objective function of a three-dimensional effective vector wave number range; a regularization parameter is introduced into the second objective function, an optimal solution is determined by solving a regularization problem, and an optimized finite difference coefficient is obtained; and substituting the optimized finite difference coefficient into the second objective function to obtain the maximum value of the effective vector wave number range. According to the method, the obtained three-dimensional space is established on the basis of the minimum norm, the domain second-order partial derivative rule grid finite difference coefficient is obtained, and error accumulation in numerical simulation of the deep target can be better reduced.

The invention relates to the technical field of coal seam inversion and imaging, in particular to a reflection channel wave imaging method and system considering the layering attenuation characteristic of a thick coal seam, and the method comprises the steps: constructing a coal seam layering geologic model, setting the initial wave velocity and initial quality factor of each layer, and forming an initial parameter field; setting an artificial seismic source and a plurality of detection points, exciting a seismic source signal at the artificial seismic source, recording seismic source information, and collecting a seismic wave response signal of each detection point; performing forward modeling by using a fractional order viscoelastic wave equation considering a spatial variable quality factor; carrying out full waveform inversion to obtain an inversion model; and on the basis of the inversion model, a level set function and a self-adaptive grid division strategy are adopted to express a barrier boundary in the coal seam to execute pre-stack wave field separation reverse time migration imaging, and a coal seam structure image is generated. The method can significantly improve the imaging resolution and calculation efficiency of the coal seam structure, and achieves the high-precision recognition of the boundary of the obstacle.

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 relates to the technical field of coal seam stress prediction, in particular to a coal seam three-dimensional stress prediction and visualization method, system and device and a medium, and the method comprises the steps: setting an artificial seismic source and a detection point, exciting a seismic source signal, and collecting a seismic wave response signal and a stress value; carrying out wave field forward modeling simulation through a time fractional order viscoelastic wave equation to obtain a forward simulation wave field, and carrying out full waveform inversion to output three-dimensional wave velocity data; performing reverse time migration imaging on the seismic wave response signal to obtain a three-dimensional coal seam structure mask pattern; performing slicing processing on the three-dimensional wave velocity data and the three-dimensional coal seam structure mask pattern along the same spatial dimension, then combining the three-dimensional wave velocity data and the three-dimensional coal seam structure mask pattern into an input tensor, inputting the input tensor into a convolutional neural network, and outputting a two-dimensional stress prediction slice; and splicing the two-dimensional stress prediction slices to construct a three-dimensional stress body, generating three-dimensional stress point cloud data, and further generating a three-dimensional stress grid model and a structure-stress integrated visual three-dimensional cloud picture. According to the invention, high-precision three-dimensional stress prediction and real-time early warning of the coal seam can be realized.

The invention discloses an intelligent exploration data analysis system and method, and relates to the technical field of geological exploration, and the system comprises a test module, an analysis module and a prediction module. According to the method, geological electromagnetic and stratum original data are collected in a multi-dimensional mode, forward modeling and inversion are carried out on electromagnetic data, the relation between underground conductivity and an electromagnetic field can be accurately simulated, data errors are effectively reduced, the electromagnetic data processing accuracy is improved, a geological unit calculates a seismic wave field based on wave equation numerical simulation, and the electromagnetic data processing accuracy is improved. The prediction wave field and the observation wave field are substituted into an objective function, the stratum velocity and density are accurately inverted, the geological data processing precision is guaranteed, the joint unit organically combines an electromagnetic inversion objective function, an earthquake full-waveform inversion objective function and a parameter correlation constraint term through a joint optimization formula, multi-source data depth fusion is achieved, and the geological data processing precision is guaranteed. The limitation of single data processing is avoided, and the underground geological condition is reflected more comprehensively.

The invention relates to the technical field of elastic waves, and discloses an elastic wave full-waveform inversion method, which comprises the steps of inputting observed seismic data into a preset neural network encoder, and extracting to obtain high-dimensional feature representation; inputting the high-dimensional feature representation into a preset neural network decoder to obtain a gradient update value; obtaining an elastic parameter model based on the gradient update value and the initial elastic parameter model; obtaining Lame parameters based on the longitudinal wave velocity, the transverse wave velocity and the medium density; based on a Lame parameter and an elastic wave equation, obtaining synthetic seismic data; obtaining an objective function based on the observed seismic data and the synthesized seismic data; using the target function to iteratively update the preset neural network weight until the loss function reaches a first preset threshold value or the number of iterations reaches a second preset threshold value, and obtaining a target gradient update value; and based on the target gradient update value and the target elastic parameter model, obtaining a final target elastic parameter model so as to relieve discomfort caused by simultaneous inversion of multiple elastic parameters.

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 meshbelt structure simulation method with heterogeneous left and right selvedges based on parameterized segmented modeling, which comprises the following steps of: deconstructing a meshbelt into three heterogeneous functional areas of left overlock, middle main body and right connection, and independently controlling by adopting a Bezier curve, a composite wave equation and a catenary fusion model respectively; the synchronous accurate expression of the closeness of the left side overlock coil, the warp and weft interweaving stability of the main body and the mechanical rationality of the right side U-shaped connection is realized, and the problem of geometric distortion caused by unified modeling in the prior art is solved. According to the method, a cross-scale coupling equation of tension (T), a friction coefficient (mu) and linear density (rho) is introduced, and a closed-loop mapping relation of technological parameters, geometrical morphology and mechanical response is established. For example, when the tension of the braid is increased, the curvature radius of the U-shaped connecting part is automatically reconstructed according to a formula, so that parameter regulation and control are upgraded from experience trial and error to intelligent optimization of physical driving.

The invention relates to a high-order staggered grid viscous sound medium reverse time migration tunnel advanced detection method and system, and belongs to the field of tunnel engineering geophysical prospecting. According to the method, a viscoelastic anisotropic TTI medium anisotropy wave equation expressed by a fractional Laplace operator is deduced, and a viscoelastic anisotropic TTI medium first-order speed-stress elastic wave equation is solved by adopting a high-order staggered grid finite difference method. According to the invention, the wave field information at any moment is obtained, so that the propagation process of seismic waves in an underground medium is simulated more accurately. In order to solve the problem of reflection possibly caused by artificial boundaries, a PML (perfectly matched layer) absorption boundary processing method of a new attenuation function is adopted. According to the method, seismic waves can be effectively absorbed, reflection is avoided, and therefore simulation stability and accuracy are guaranteed. Generally speaking, the invention provides an accurate and stable tunnel advanced detection forward modeling method, and a reverse time migration algorithm is carried out in an attenuation TTI medium.

The invention discloses a generalized finite difference scalar wave unfavorable geology frequency dispersion analysis method and related equipment, relates to the field of unfavorable geology analysis, and mainly aims to solve the problem that the conventional unfavorable geology analysis is relatively low in flexibility and accuracy. The method comprises the following steps: generating field nodes of a target area; performing field variable approximation on the field nodes to determine a discretized wave equation; and performing frequency dispersion analysis and wave field simulation based on the discretized wave equation to obtain an output result of unfavorable geology. The method is used for the unfavorable geological frequency dispersion analysis process of the generalized finite difference scalar waves.

The present invention provides an optical property modeling method and device, and an optical parameter measurement method, wherein the optical property modeling method comprises: obtaining characteristic parameters of a periodic medium, and constructing a first coupled wave equation group based on the characteristic parameters and a strict coupled wave analysis method; the periodic medium is symmetrical about the incident plane of the light beam; obtaining a linear combination of two orthogonal components of an electric field or two orthogonal components of a magnetic field, and constructing a second coupled wave equation group based on the first coupled wave equation group and a linear combination of the two orthogonal components of the electric field or the two orthogonal components of the magnetic field; simplifying the second coupled wave equation group according to the symmetry of the periodic medium; solving the simplified second coupled wave equation group to obtain a theoretical spectrum. The present invention improves the efficiency of optical modeling to obtain a theoretical spectrum.