750 results about "Reflection coefficient" patented technology

The invention discloses a reconfigurable intelligent metasurface beam tracking device, and belongs to the technical field of wireless communication. Specifically, the reconfigurable intelligent metasurface module integrates a large number of metasurface units, a tunable element is embedded to achieve reconfiguration of electromagnetic characteristics, a receiving end of the communication receiving and transmitting module receives reflected signals of the metasurfaces and feeds back the quality of the received signals in the moving process, and the positioning sensing module senses the position of the receiving end in real time. And the codebook calculation module calculates a phase regulation codebook based on the signal quality and the position feedback information, converts the phase regulation codebook into a control voltage value, transmits a corresponding digital control signal to the voltage regulation module, and outputs multiple paths of analog voltages, so that the reflection coefficient of the metasurface array is regulated and controlled, and real-time dynamic wave beam tracking is realized. According to the invention, the coverage range, the transmission rate and the service quality of the communication network are obviously improved.

The invention provides a composite board detection method and system, and relates to the technical field of board detection, and the method comprises the steps: carrying out the ultrasonic scanning of a composite board, dividing the composite board into a plurality of acoustic characteristic layers, building a sound wave propagation model based on a recursive transfer matrix method, and calculating the reflection coefficient and transmission coefficient of each layer interface. A deep learning algorithm is combined with physical constraints to dynamically correct acoustic parameters, the corrected parameters are substituted into a model reconstruction signal, and defects are accurately positioned. And performing grid division on the defect position to analyze stress field distribution, and determining a stress concentration position. And fine scanning is carried out on a stress concentration area by adjusting probe parameters, and high-resolution data is obtained to update defect distribution information. And finally, outputting detection data including defect positions, stress concentration distribution and hazard levels, and realizing accurate detection and evaluation of the defects of the composite board.

The invention discloses a cable defect detection method and system based on double-path coupling frequency domain reflection, and the method comprises the steps: enabling a signal generator to generate two groups of correlated sine frequency sweeping signals, namely a test signal S1 and a reference signal S2; s1 is injected into the head end of the cable, and a reflected signal R1 is collected; respectively extracting the spectrum characteristics of the amplitudes and the phase differences of the R1, the S1 and the S2 by adopting a double-path signal comparison device, carrying out comprehensive judgment, and determining the reflection spectrum characteristic of the R1; the reflection coefficient spectrum is optimized through a window function, and signal features are obtained; converting the optimized frequency spectrum into a time domain reflection waveform, and extracting reflection peak time information; and calculating defect positions by combining electromagnetic wave propagation velocity. According to the cable defect detection method and system based on double-path coupling frequency domain reflection, high-precision detection of a long-distance cable is realized; and two paths of coupling signals are introduced for comparison, so that the accuracy of phase measurement is improved, and an efficient nondestructive detection scheme is provided for power grid operation and maintenance and the like.

A high-resolution wave-impedance inversion method and apparatus based on spectral fusion. The method comprises: under the constraints of horizon data and a spatial variation function, performing Kriging spatial interpolation processing on a low frequency band of well logging data to obtain a low-frequency-band wave-impedance data body zL; according to seismic wavelets and a mid-frequency model of Cauchy norm-constrained sparse spike inversion, performing Cauchy norm-constrained sparse spike inversion on seismic data to obtain a reflection coefficient, and performing trace integration to obtain a mid-frequency-band relative wave-impedance data body zM; under the constraints of the horizon data and the spatial variation function, performing sequential Gaussian stochastic simulation on a high frequency band of the well logging data to obtain a high-frequency-band wave-impedance data body zH; performing Fourier transform on zL, zM and zH to obtain frequency domain data bodies ZL, ZM and ZH; performing spectral fusion to obtain a full-frequency-band wave-impedance data body ZMGI; and performing inverse Fourier transform to obtain a high-resolution wave-impedance inversion result zMGI. High-resolution wave-impedance inversion is therefore effectively performed.

The invention belongs to the field of ultrasonic oil film thickness measurement, and discloses a minimum oil film thickness calibration method considering different friction pair shape structures, which utilizes an ultrasonic oil film thickness calibration platform capable of designing a corresponding calibration test piece according to different structure shapes of a measured friction pair, and aims at the shape structure characteristics of different friction pairs. The minimum oil film thicknesses of friction pairs with different shapes and structures are accurately controlled, ultrasonic reflection signals under the different minimum oil film thicknesses are collected and subjected to frequency domain analysis processing, and on the basis, polynomial fitting is conducted on the minimum oil film thicknesses and reflection coefficients in combination with a film thickness calculation model; and the ultrasonic minimum oil film thickness calibration under a special friction pair structure is realized.

The invention discloses a method and system for measuring the ultrasonic proportionality coefficient of the thickness of a lubricating film of a sliding bearing, and the method comprises the steps: firstly, obtaining the comprehensive response coefficient of sound waves in a transition combination layer of a substrate and a lining through calculating the proportionality coefficient of a reflection echo signal of a bearing bush in an air medium and a reflection echo signal of a substrate-lining interface; characterizing the propagation rule of sound waves in a base body and lining combination layer; secondly, a proportionality coefficient is constructed based on a lubricating film reflection echo signal and a matrix-lining interface reflection echo signal, a lubricating film reflection coefficient is calculated in combination with a comprehensive response coefficient of sound waves in a matrix-lining transition combination layer, and finally the thickness of the lubricating film is solved through a resonance model method or a composite model method. According to the invention, the influence of the change of the binding force between the sensor and the substrate and the unstable voltage of ultrasonic pulse emission equipment on the reference signal can be avoided, the robustness is high, and the long-term reliable monitoring of the thickness of the lubricating film of the sliding bearing is facilitated.

The invention provides a blade lightning protection channel online test method and system based on a pulse current method, and belongs to the technical field of blade lightning protection detection.The method comprises the steps that a pulse waveform is generated through a pulse current generator to simulate lightning impulse, and the lightning impulse is injected into the surface of a lightning arrester of a blade lightning protection system; measuring channel currents at different positions of the surface of the lightning arrester in the length direction by using a distributed sensor; a high-voltage differential probe is used for collecting and measuring the waveform of an incident wave and the waveform of a reflected wave of an injection point; detecting blade surface potential distribution by using a capacitive coupling sensor; comparing the waveform of the incident wave, the waveform of the reflected wave and the theoretical waveform, determining a time domain reflection coefficient and a frequency domain impedance spectrum, and determining a fault type; after the fault type is determined, the fault point is positioned based on the channel current and the blade surface potential distribution, the lightning impulse is simulated through the pulse waveform, the lightning protection channel is further tested, shutdown is not needed, the testing efficiency of the blade lightning protection channel is effectively improved, and the power generation efficiency is ensured.

Disclosed are systems, apparatuses, processes, and computer-readable media for wireless communications. For example, an example of a process includes transmitting, to a reconfigurable intelligent surface (RIS), a plurality of reference signal beams in a plurality of directions. The process may further include receiving, from the RIS, a plurality of reflection reference signal beams in the plurality of directions based on a plurality of backtracking reflection coefficients. The process may include measuring a signal strength of each reflection reference signal beam in each direction of the plurality of directions for each backtracking reflection coefficient of the plurality of backtracking reflection coefficients. The process may further include determining, basing on measuring the signal strength, a selected beam direction from the plurality of directions and a selected backtracking reflection coefficient of the plurality of backtracking reflection coefficients.

A method and apparatus for detecting a local defect of a cable, and a storage medium. The method comprises: acquiring a reflected signal of a testing signal at a head end of a cable to be subjected to detection; generating a first reflection coefficient of the head end of said cable, and on the basis of a real part of the first reflection coefficient, determining whether said cable has a defect, and if so, determining the location of the defect on the basis of the real part of the first reflection coefficient, and acquiring a frequency component and a first reflection intensity at the location of the defect, and a corresponding correction angle value; and determining the polarity of the defect, and generating a first-reflection intensity cosine function, calculating the difference between the first-reflection intensity cosine function and the real part of the first reflection coefficient to generate a second reflection coefficient, wherein the first reflection coefficient is a reflection coefficient when said cable has no defect, and determining the type of the defect on the basis of the ratio of a second reflection intensity regarding the second reflection coefficient to the first reflection intensity. Thus, cable defect locating and defect type identification and classification are realized, thereby improving the cable fault detection efficiency.

Systems, apparatus, processes, and computer-readable media for wireless communication are disclosed. For example, according to aspects described herein, a reconfigurable smart surface (RIS) may receive a configuration message including an indication to reduce a gain of at least a portion of a reflected beam. The RIS may generate reflection coefficients for a metadata element of the RIS based on the configuration message. The RIS may configure the metadata element based on the reflection coefficient. For example, based on configuring the metadata element, the RIS may generate or generate the reflected beam (e.g., a sensing reflected beam or an interfering reflected beam) by reflecting a sensing signal received from a sensing transmitter device.

The invention relates to a method for rapidly calculating the RCS of a coated unmanned aerial vehicle under a complex background, and belongs to the field of stealth target electromagnetic scattering calculation. Comprising the following steps: determining reflection coefficients of a single-layer wave-absorbing material with a PEC substrate at different incident angles according to an equivalent transmission line theory; calculating induced current generated by incident electromagnetic waves on the surface of the target coated with the single-layer wave-absorbing material according to the reflection coefficient; based on the induced current, utilizing physical optics and physical diffraction theories to calculate scattered fields of all the illuminated surface elements, and superposing all the calculated scattered fields; the scattering field of each surface element on the sea surface model at a fixed angle is obtained through vector superposition of the scattering fields obtained through an EPO'field and a TSM; and finally, superposing all the scattered fields, and calculating the total field RCS of the complex sea background and unmanned aerial vehicle composite model based on the superposed total scattered field. According to the method, the RCS of the coating unmanned aerial vehicle under the complex background can be quickly calculated.

The invention discloses an omni-directional reconfigurable metasurface-assisted array radar transmitted beam forming method for non-line-of-sight target perception, and aims at an omni-directional reconfigurable metasurface-assisted array radar to maximize the main lobe gain of a transmission region and a reflection region as a target function. Constructing an optimization model by taking the array emission energy and the transmission coefficient and the reflection coefficient of the omnidirectional reconfigurable metasurface as constraints, and performing equivalent transformation on the optimization model; constructing a scaling form augmented Lagrange function of the optimization model after the equivalent transformation, decomposing the optimization model after the equivalent transformation into a plurality of sub-problems based on an ADMM criterion, and solving each sub-problem by adopting non-convex quadratic constraint quadratic programming and positive semidefinite relaxation technologies, and obtaining an optimized emission weight and a transmission coefficient and a reflection coefficient of the omni-directional reconfigurable metasurface through iteration. The STARS-assisted array radar provided by the invention realizes full space coverage, and has better performance in the aspect of main lobe gain compared with a standard RIS-assisted array radar.

The invention belongs to the technical field of advanced detection, and provides a shield tunnel optical fiber earthquake forecasting method and system, and the method comprises the steps: obtaining the carrying position and number of seismic wave excitation seismic source equipment based on a shield tunneling machine structure type, installing a pneumatic hammering seismic source, calculating the dynamic correlation between the air pressure and the hammering force, and determining the shield air pressure range. Determining the relative spatial position of the seismic source and the optical fiber sensor, and establishing an observation mode; according to the relative spatial position and the observation mode, arranging an optical fiber sensor in a spatial range behind the shield tail of the shield tunnel, controlling a pneumatic hammering seismic source to hammer surrounding rock to generate seismic waves to obtain geological forecast data, preprocessing the obtained geological forecast data to obtain geological wave field data, and calculating the geological wave field data according to the geological wave field data. And calculating rock mass mechanical parameters, a reflection coefficient and a stress gradient of the shield excavation surface, obtaining a three-dimensional image, carrying out feature analysis, obtaining a disaster-causing structure position, a disaster-causing structure scale and an underground water occurrence condition, and completing prediction of the shield tunnel optical fiber earthquake.

The invention provides a crack pore medium multi-domain seismic inversion method and device, and relates to the technical field of geological prospecting, and the method employs a reflection coefficient calculation mode based on an accurate Graebner equation, and employs a Cauchy prior constraint to combine a target function constructed by a time domain and a frequency domain. The method realizes simultaneous inversion of three parameter profiles of porosity, crack weakness and water saturation, an accurate Graebner equation can accurately describe reflection coefficient characteristics of a crack medium, time domain inversion can enhance noise immunity, and frequency domain inversion can improve resolution, so that the accuracy, the noise immunity and the resolution of an inversion result are considered, and the method is suitable for large-scale popularization and application. And high-precision and high-resolution estimation of the porosity, the crack weakness and the water saturation is realized.

The invention discloses a seismic reflection coefficient determination method and apparatus, and an electronic device. The method comprises the following steps: constructing a constraint relationship between seismic wavelets and L1 norms according to a convolution model of a target detection area; acquiring three-dimensional stratigraphic dip angle data of a target layer in the target monitoring area; according to the three-dimensional stratigraphic dip angle data and the constraint relation, the seismic reflection coefficient of the target stratum is determined. According to the technical scheme provided by the embodiment of the invention, the constraint relationship between the seismic wavelet and the L1 norm is constructed according to the convolution model of the target detection area, the three-dimensional stratigraphic dip angle data of the target stratum in the target monitoring area is acquired, and the seismic reflection coefficient of the target stratum is determined according to the three-dimensional stratigraphic dip angle data and the constraint relationship. Therefore, the seismic reflection coefficient is accurately determined, the fidelity of seismic data is ensured, the high-resolution result is finer, and a high-quality data basis is provided for reservoir prediction and seismic inversion.

The invention provides a crustal stress prediction method and device, electronic equipment and a storage medium, and relates to the technical field of seismic survey. The method comprises the following steps: obtaining observation seismic data of seismic wavelets in a strong VTI medium, and constructing a to-be-inverted parameter matrix based on a PP wave reflection coefficient corresponding to the strong VTI medium; constructing a posterior probability function obeyed by an inversion parameter matrix corresponding to the to-be-inverted parameter matrix based on a Bayesian inversion theory and observation seismic data, and determining a target functional based on a prior probability function and a likelihood function corresponding to the posterior probability function; determining medium density and each stiffness matrix coefficient based on an inversion parameter matrix solving result of the target functional, and determining a flexibility matrix of the strong VTI medium based on each stiffness matrix coefficient; and predicting the ground stress distribution of the target profile in the strong VTI medium based on the medium density and the positive strain matrix and the flexibility matrix corresponding to the strong VTI. Therefore, the crustal stress prediction accuracy under the strong VTI medium is improved.

The invention provides a photovoltaic power generation optimization method for a fishing-light complementary site and a related device. The method comprises the following steps: acquiring a photovoltaic panel sunlight angle, a water surface wind speed, a water surface wave height, a water depth, a photovoltaic panel core parameter, a reflector parameter and double-measuring-point actually-measured light intensity data; determining the initial receiving radiation quantity of the photovoltaic panel according to the sunlight angle of the photovoltaic panel, the core parameters of the photovoltaic panel, the parameters of the reflector and the actually measured light intensity data of the double measuring points; according to the water surface wind speed, the water surface wave height, the water depth, the double-measuring-point actually-measured light intensity data and the reflector parameters, the water surface correction reflection coefficient is determined; correcting the initial received radiation quantity of the photovoltaic panel by adopting the water surface correction reflection coefficient to obtain the actual received radiation quantity of the photovoltaic panel; and according to the actually received radiation quantity of the photovoltaic panel, the inclination angle of the photovoltaic panel is adjusted, and then the photovoltaic power generation process of the fishing-light complementary site is optimized. Excessive concentration of local radiation quantity caused by secondary reflection of sunlight through a water surface-reflector can be avoided, and generation of hot spots on the photovoltaic panel is inhibited.

The invention provides a multiple prediction method and device for a converted wave free interface, equipment and a medium, and belongs to the field of seismic data processing. The method comprises the following steps: acquiring a seismic signal, and extracting a longitudinal wave component downlink wave record of the seismic signal; determining a downlink wave imaging data volume according to the downlink wave record of the longitudinal wave component based on Kirchhoff offset; based on Kirchhoff demigration, determining a longitudinal wave component demigration seismic record according to the downlink wave imaging data volume, and taking the longitudinal wave component demigration seismic record as a free interface multiple prediction factor; determining a free interface factor based on the longitudinal wave component demigration seismic record and a free interface preset reflection coefficient; and constructing a free interface multiple prediction model based on the free interface factor and the free interface multiple prediction factor. According to the method, all free interface multiple components in converted wave data can be accurately predicted, so that the multiple elimination precision can be further improved.

A method includes defining a design domain for a meta-element configured to be attached to an edge of thin wall structure with an arbitrary boundary condition, and executing a topological optimization process on the design domain and providing a topology optimized shape for the meta-element. The topological optimization process includes an objective function that minimizes a reflection coefficient of flexural waves in the audible frequency range propagating towards and impinging the meta-element.

The invention relates to a spherical wave oscillation-slow nuclear high-order decoupling method and device for physical property inversion, and belongs to the field of geophysical detection and acoustic detection. The method comprises the following steps: obtaining and inputting observation data, and initializing seabed elastic parameters; determining an integral interval according to the observation data, decomposing the integral interval into a plurality of sub-intervals, and setting a decoupling order; calculating an oscillation nuclear local integral matrix of each order corresponding to the observation data incident angle; calculating the angle deviation guiding quantity of each order of the slowly-varying kernel; calculating a reflection coefficient vector of each order of spherical wave; calculating a spherical wave reflection coefficient of high-order decoupling; and solving a subsequent inversion problem by using the obtained high-order decoupled spherical wave reflection coefficient. According to the method, the calculation efficiency of the seabed spherical wave reflection coefficient in the inversion process is effectively improved; the decoupling order is improved, the calculation precision is improved, or the number of the subintervals can be reduced by improving the decoupling order under the condition that the calculation precision is certain.

The invention discloses a polyether-ether-ketone-based composite material with low reflection and high electromagnetic shielding effectiveness and a preparation method of the polyether-ether-ketone-based composite material, and belongs to the technical field of electromagnetic shielding materials. The problems that an existing electromagnetic shielding material is high in reflection coefficient, insufficient in shielding effectiveness, poor in high temperature resistance and corrosion resistance and difficult to adapt to severe environment service are solved. The composite material is composed of a porous conductive substrate and an asymmetric conductive layer. Wherein the porous conductive substrate is formed by blending a soluble ketimine type polyether-ether-ketone solution and a conductive filler suspension and then carrying out non-solvent induced phase separation and acidification reduction reaction, the electric conductivity is 2.8 * 10 <-4 >-0.1839 S / cm, and the porosity is 80%-90%; the asymmetric conductive layer is plated on the surface of one side of the porous conductive substrate, and the asymmetric conductive layer is made of nickel-phosphorus alloy. The composite material significantly reduces the electromagnetic wave reflection coefficient through an absorption-reflection-reabsorption mechanism, has high electromagnetic shielding effectiveness and excellent high temperature resistance and corrosion resistance, and is suitable for serving in a severe environment.

The invention discloses a method for estimating a radar scattering cross section of a coated medium target, relates to the field of radar application, is suitable for an electrically large coated medium target, and improves the calculation efficiency of the radar scattering cross section while ensuring the calculation precision. The method comprises the following steps: calculating a reflection coefficient of a target surface by adopting a stepping method based on target coating medium characteristics; respectively calculating induced electromagnetic currents of different surface elements after the target is subdivided; based on the induction electromagnetic currents of different surface elements, the scattering electric field of each surface element of the target is calculated, the scattering electric fields of the visible surface elements are summed to obtain the total scattering electric field of the target, and the radar cross section RCS is calculated according to the total scattering electric field. The method can be suitable for electromagnetic characteristic analysis of the electrically large-size coating target, and the calculation precision is improved while the calculation efficiency is ensured.

The invention provides a sea surface oil spill GNSS reflection signal simulation system and method, and the method comprises the steps: building an oil spill sea surface dielectric constant module, determining the effective dielectric constant of an oil spill sea surface, and obtaining a Fresnel reflection coefficient; establishing an oil spill sea surface rough spectrum module, obtaining a model of a sea wave spectrum of the oil spill sea surface according to a traditional model of the sea wave spectrum and the influence of oil spill on the sea wave spectrum, and obtaining a surface spectral density function and a surface root-mean-square height; establishing a microwave scattering module of the oil spill sea surface, and determining a microwave scattering coefficient according to the Fresnel reflection coefficient, the surface spectral density function and the surface root-mean-square height; a DDM simulation module is established, and a DDM waveform is obtained through simulation according to the GNSS signal parameters and the microwave scattering coefficients; and the DDM waveform is obtained through simulation of each module. According to the simulation system and method, a unified mathematical model is established to describe the relation between the oil film and the GNSS reflection signal, and accurate simulation of the influence of the oil film on the GNSS reflection signal is achieved.

The present invention discloses a method for obtaining the reflection coefficient of a floating breakwater, and belongs to the field of marine engineering technology. The present invention constructs a numerical simulation water tank, adopts a mass source wave generation method to generate waves, the wave height is H, and the period is T; obtains the spatial distribution curve of the free liquid surface height of the synthetic reflection wave; and calculates the wave height H of the synthetic reflection wave according to the spatial distribution curve of the free liquid surface height of the synthetic reflection wave. 合 ; Through the formula H 反 =#imgabs0#get H 反 The value of, where a=H 合 / 2, the radian value corresponding to the phase angle A is equal to the phase difference Φ / 2, and the phase difference Φ is obtained by formula 1. Without the use of measuring instruments such as wave height meters and Doppler current meters, the reflection coefficient of the floating embankment can be obtained by vector superposition analysis using the spatial distribution data of the free liquid surface, avoiding the previous method of extracting data through measuring instruments and solving with the help of complex formulas.

The invention provides a remote sensing water depth detection method, medium and system based on blue-green laser, and belongs to the technical field of water depth measurement, and the method comprises the steps: firstly, obtaining blue and green laser signals with known transmitting power, and collecting water surface reflection and water bottom reflection data under different incident angles; based on the data, a laser incident angle matrix is established, and a water surface reflection coefficient and a water surface transmission coefficient are calculated. Meanwhile, an atmospheric transmission matrix is established according to the meteorological parameters, and the transmission loss of the laser signal in the atmosphere is calculated. Then, a water surface fluctuation model is established, and water surface reflection and water bottom reflection signals are separated; and establishing a water depth calculation equation set by using the blue-green dual-wavelength signal intensity ratio matrix and combining with the optical characteristics of the water body. And finally, according to the propagation time difference and the equivalent propagation speed of the laser signal in water, solving a water depth equation set, and finally obtaining a water depth topographic map. The problem that in the prior art, a water surface reflection signal and a water bottom reflection signal cannot be separated, and the measurement precision is affected is solved.

The invention provides a seismic high-resolution inversion method and device. The method comprises the following steps: acquiring angle domain pre-stack seismic trace gather data, seismic horizon data, seismic wavelet data, initial low-frequency longitudinal wave impedance, initial low-frequency transverse wave impedance and initial low-frequency density; obtaining a reflection coefficient according to the angle domain pre-stack seismic trace set data, the seismic horizon data, the seismic wavelet data and a Fatti longitudinal wave reflection coefficient equation; establishing an inversion objective function according to the initial low-frequency longitudinal wave impedance, the initial low-frequency transverse wave impedance, the initial low-frequency density and the reflection coefficient, and obtaining an inversion equation through the inversion objective function; presetting an inversion iteration frequency and an inversion error, and carrying out iterative solution on the inversion equation to obtain intermediate-frequency longitudinal wave impedance, intermediate-frequency transverse wave impedance and intermediate-frequency density; and processing the intermediate-frequency longitudinal wave impedance, the intermediate-frequency transverse wave impedance and the intermediate-frequency density by using a preset first formula to obtain final high-resolution longitudinal wave impedance, transverse wave impedance and density. According to the invention, the precision and resolution of seismic inversion can be effectively improved.

The invention discloses a full-waveform inversion method for enhancing boundary perception. The full-waveform inversion method comprises the following steps: acquiring seismic data; constructing a seismic waveform inversion network model; the seismic waveform inversion network model comprises an encoder, a decoder and a jump connection; carrying out the size adjustment of the seismic data through an encoder, and then extracting the advanced features of the reserved low-frequency information, and obtaining a multi-scale coding feature map; inputting the multi-scale encoding feature maps extracted by the encoder into different levels of the decoder by utilizing jump connection; a decoder is used for obtaining distinguishing information of the spatial relation contained in the multi-scale coding feature map, and a prediction speed model is generated; multi-scale analysis is executed through a wavelet convolution layer in an encoder, and low-frequency features are enhanced; then, distinguishing information is extracted by utilizing the spatial relationship of the features through a spatial attention module in a decoder so as to improve the resolution; finally, a reflection coefficient is introduced into the loss function to emphasize punishment on the low-contrast boundary; therefore, the full-waveform inversion precision is improved.

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 a seismic wavelet and noise simultaneous estimation method and device, and the method comprises the steps: obtaining a first forward model of a reflection coefficient and seismic data, and a second forward model of seismic wavelets and seismic data, and building a probability graph model according to the relation between random variables in the first forward model and the second forward model; establishing a prior model based on the variables with the interdependency relationship being deterministic dependency in the probability graph model, and establishing a posterior model based on the variables with the interdependency relationship being probabilistic dependency in the probability graph model; logging data and seismological observation data are obtained, and simultaneous iteration calculation is performed on the seismic wavelets and the seismic noise through a Gibbs algorithm by using the logging data, the seismological observation data, the prior model and the posterior model to obtain a target seismic wavelet and noise sequence. According to the method, the coupling between the seismic wavelets and the seismic noise is considered, the seismic wavelets and the noise are iteratively solved at the same time through the Gibbs algorithm, and the seismic wavelet extraction precision is improved.