309 results about "Wavenumber" patented technology

The invention discloses a road cavity volume calculation method and device based on a three-dimensional FK migration algorithm, equipment and a medium, and relates to the technical field of road cavity detection, and the method comprises the steps: obtaining three-dimensional ground penetrating radar data, and carrying out the preprocessing of the data, and obtaining target radar data; converting the target radar data to a frequency wavenumber domain to obtain frequency wavenumber domain radar data; establishing a road hierarchical velocity model, and performing hierarchical wave field extrapolation and double weight correction based on the frequency wave number domain radar data and the road hierarchical velocity model to obtain corrected wave field data; performing sparse constraint optimization imaging processing on the corrected wave field data to obtain a target imaging result; and carrying out boundary extraction processing on the target imaging result to obtain a three-dimensional cavity boundary point set, and carrying out voxelization processing based on the three-dimensional cavity boundary point set to obtain a road cavity volume. According to the invention, high-precision restoration of the shape of the road cavity and automatic extraction of the volume can be realized in three-dimensional GPR imaging.

The invention discloses a multi-scale sea surface electromagnetic scattering calculation method and device based on marine environment elements and a storage medium, and belongs to the technical field of marine detection and physical remote sensing. The method comprises the following steps: calculating a tension wave sea spectrum according to an Elfouhail sea spectrum in combination with a Longuet-Higgins directional spectrum; a JONSWAP sea spectrum is combined with a Longuet-Higgins directional spectrum to calculate a gravity wave sea spectrum; superposing the tension wave on the gravity wave according to the slope by using a linear superposition method, and calculating a self-adaptive cut-off wave number; the contribution of small-scale waves to scattering is calculated by adopting an integral equation method, mirror scattering of large-scale waves is counted by adopting a Kirchhoff approximation method, and statistics and fusion are carried out on each scale wave according to the self-adaptive cut-off wave number, so that an electromagnetic scattering value is obtained. According to the method, a sea surface geometric calculation method based on real marine environment element driving is adopted, and accurate calculation is performed aiming at sea surface scattering characteristics under different sea conditions in combination with an integral equation method-Kirchhoff approximation electromagnetic scattering calculation model of an adaptive cut-off wave number, so that the calculation precision of scattering intensity is improved.

The invention provides a high-frequency electromagnetic excitation cable detection image reconstruction method and system, and belongs to the technical field of underground facility multi-field coupling detection and imaging, and the method comprises the steps: applying high-frequency electromagnetic excitation to an underground medium; receiving an electromagnetic signal transmitted through an underground medium, synchronously collecting a sound signal generated by electromagnetic excitation, and preprocessing the electromagnetic signal and the sound signal; performing structure imaging on the electromagnetic signal by adopting a frequency wave number migration imaging algorithm to obtain medium structure distribution; performing electrical parameter imaging on the electromagnetic signal by adopting a Born approximate inverse scattering imaging algorithm to obtain medium electromagnetic characteristic parameter distribution; sound source distribution reconstruction is carried out on the sound signals through a time reversal algorithm, and medium sound source response characteristic distribution is obtained; and inputting imaging results of the three algorithms into a pre-constructed deep reinforcement learning model for dynamic fusion to generate an underground medium comprehensive parameter reconstruction image.

The invention relates to the technical field of spectrum detection, in particular to a Raman spectrum-based rice oil detection method, which comprises the following steps: extracting main and auxiliary peaks to construct a peak group structure, performing sliding comparison to output a conformity region, identifying jump to remove an abnormal signal, calibrating a peak position to generate a correction structure, and performing homing numbering to complete a detection scheme. According to the method, the feature matching rule is constructed by combining the peak position spacing, the intensity ratio and the wavenumber index, the recognition precision of the multi-peak structure in the spectrum is enhanced, non-continuous response signals are recognized and eliminated through derivative change trend and proportion deviation constraint, and the accuracy of stable response extraction is improved. Dynamic calibration and position adjustment of peak positions are realized by utilizing an offset trend sequence, the consistency of a peak group structure in an overall direction and the continuity of a local structure are ensured, a numbering system and structure attribution division is completed after homing sorting, and the integrity of feature expression in a complex sample and the orderliness and reliability of spectral data processing are improved.

The invention discloses a real-time wave number space resampling method and system for SS-OCT, and relates to the technical field of signal processing, and the method comprises the steps: synchronously collecting a K clock signal and an OCT interference signal, and carrying out the filtering preprocessing of the signals; carrying out Hilbert transformation on the preprocessed K clock signal, constructing an analysis signal, extracting an instantaneous phase, and forming an instantaneous phase sequence; using a preset ideal equal wave number interval phase sequence as a target phase sequence, matching the target phase sequence with the instantaneous phase sequence, calculating a corresponding position index of each target phase in the instantaneous phase sequence, and obtaining a remapping sampling point index; and processing the preprocessed OCT interference signal by using an SINC interpolation method according to the remapping sampling point index to obtain an OCT signal which is subjected to resampling and has equal wave number intervals. According to the invention, the SINC interpolation algorithm is utilized to complete resampling in real time and with high precision in the data acquisition stage, and the requirement of high-speed real-time imaging is met.

The invention relates to the technical field of signal processing, in particular to a pipeline guided wave signal processing method and system. The method comprises the following steps: acquiring a multi-channel ultrasonic guided wave signal through a sensor array circumferentially arranged around a pipeline, and converting the multi-channel ultrasonic guided wave signal into a frequency domain response matrix; calculating a circumferential amplitude distortion degree based on the power spectrum amplitude distribution of each channel; constructing a modal coupling matrix by combining the circumferential amplitude distortion degree and the theoretical wave number of the ideal pipeline; correcting the standard orthogonal modal matrix to generate a modal projection matrix; and decomposing the frequency domain response matrix by using the projection matrix, and separating a structure echo component and a defect signal component. According to the method, the influence of geometric defects such as ovality deviation is evaluated more accurately, energy exchange between modeling modes is accurately established, and orthogonality hypothesis is corrected; modal energy leakage is inhibited, and the decomposition accuracy is improved; signals are effectively separated, coherent noise interference is solved, and the accuracy and reliability of pipeline defect detection are improved.

The invention discloses a two-stage Raman hyperspectral imaging method, which relates to the technical field of spectral signal analysis and comprises the following steps: S0, preprocessing Raman hyperspectral data to remove abnormal values; the method comprises the following steps: S1, taking a wave number where a target characteristic peak is located as a target wave number, and based on a neighborhood of the target wave number in Raman hyperspectral data, performing spectrum denoising on the Raman hyperspectral data by using an adaptive low-rank matrix approximation algorithm to obtain the Raman hyperspectral data after spectrum denoising; and S2, using an improved BM3D method based on a rotating block to carry out image denoising on a Raman image after Raman imaging is carried out on a target wave number in the Raman hyperspectral data after spectrum denoising. According to the method, the neighborhood of the target wave number is used as the target area, the adaptive low-rank matrix approximation algorithm is used for spectrum denoising, then the improved BM3D method based on the rotating block is used for image denoising, Raman hyperspectral fast and efficient imaging is achieved, and the method has the advantages of being good in denoising effect, high in speed and light in weight.

The invention discloses a Fourier finite difference deconvolution imaging method based on inclination angle adaptive noise suppression, and relates to the technical field of geophysical exploration of petroleum. The method specifically comprises the following steps: designing a scattering point speed model to solve a partial point spread function of an underground space; solving a partial point spread function by a speed disturbance method; solving point spread functions of all positions of the underground space by carrying out Gaussian filtering technology and inverse distance weighted interpolation on part of the point spread functions; using a plane wave destruction method to solve the inclination angle attribute of the work area, and constructing an inclination angle self-adaptive diagonal diffusion Fourier finite difference migration method based on the inclination angle attribute to solve a high-quality initial image suitable for an inversion system; and performing multi-dimensional deconvolution in the wave number domain to obtain a final migration imaging result. According to the method, the point spread function is successfully introduced into Fourier finite difference imaging, and the resolution and deep amplitude of a conventional Fourier finite difference imaging method can be effectively improved.

The invention provides an anti-noise interference frequency parameter estimation method for high-precision optical coherence ranging, which comprises the following steps: acquiring a spectral interference signal of an optical low-coherence interference system and a reference light spectral signal output only by a reference arm, and converting the spectral interference signal into interference signals with wave number space arranged at equal intervals, carrying out spectral density amplitude normalization to obtain a standard interference signal so as to eliminate interference of light source spectrum fluctuation and non-uniform sampling on frequency estimation; hanning window modulation is applied to a standard interference signal, spectrum leakage is suppressed, spectrum analysis is performed on the modulated signal, a peak value and an adjacent secondary peak value of a spectrum main lobe are extracted, a frequency correction amount is calculated based on a peak value and secondary peak value amplitude ratio relation derived by a Hanning window spectrum mode function, and an accurate frequency of the interference signal is obtained; and calculating the displacement of the object to be measured according to the product of the accurate frequency and the spectral resolution of the optical low-coherence interference system in combination with the refractive index of the light in the measurement medium.

The invention relates to a multi-machine calibration method and system for a Raman spectrometer, and a storage medium, and the method comprises the steps: outputting an adjustable substrate voltage, and setting an initial signal amplification gain; controlling a CCD sensor to collect first Raman spectrum data of a standard sample with a known standard wave number; establishing a mapping model from a pixel domain to a wave number domain based on a plurality of known standard wave numbers and second pixel point positions of the characteristic peaks corresponding to the known standard wave numbers in the first Raman spectrum data to obtain calibration parameters of the Raman spectrometer; verifying the accuracy of the mapping model by verifying the characteristic peak, and storing calibration parameters after the verification is passed; controlling the plurality of Raman spectrometers to respectively detect the same sample to be detected, acquiring second Raman spectrum data of the sample to be detected, extracting a third pixel point position corresponding to a characteristic peak of the sample to be detected, calling the calibration parameter, and converting the calibration parameter into an actually measured wave number; and feeding back and adjusting the signal amplification gain of each instrument based on a signal intensity measurement value corresponding to the second Raman spectrum data by each Raman spectrometer.

The invention provides a stone analysis method based on attenuated total reflection Fourier transform infrared spectroscopy. The stone analysis method comprises the following steps: S1, placing and fixing a stone sample to be detected; s2, acquiring an infrared light signal and converting the infrared light signal into original spectral data; s3, processing the original spectrum data to obtain a corrected spectrum; s4, taking the correction spectrum as input, obtaining a preliminary prediction category by using a partial least square discriminant analysis algorithm, and calculating a feature importance index of each wavenumber at the same time; s5, according to the feature importance index, using a random forest algorithm to obtain a secondary prediction category; and combining the accuracy rate of the partial least square discriminant analysis model and the accuracy rate of the random forest model to obtain a final prediction category. According to the method, an attenuated total reflection Fourier transform infrared spectroscopy (ATR-FTIR) technology is taken as a core, and an innovative thought of'pretreatment-free sample treatment + multi-algorithm collaborative analysis' is combined, so that a complete stone analysis process from sample fixation to clinical report output is constructed.

The invention relates to the technical field of crossing of geophysical exploration and artificial intelligence, and discloses a three-dimensional VSP data reconstruction method based on a three-dimensional blind sidewalk time frequency-space multi-domain interpolation network model, and the method comprises the steps: obtaining a to-be-reconstructed three-dimensional VSP data volume containing an irregular missing sidewalk; inputting the three-dimensional VSP data volume into a trained three-dimensional space-time-frequency mixed domain convolutional network model; performing down-sampling on the input data through an encoder to extract initial spatial-temporal characteristics; inputting the initial spatiotemporal features into a residual spatiotemporal-frequency wavenumber domain joint convolution feature fusion module to obtain fusion features; and inputting the fusion features into a decoder, and reconstructing a complete three-dimensional VSP data volume through up-sampling operation. According to the invention, through joint learning of ground seismic data (providing a macrostructure) and illuminance data (providing a physical propagation path), a reconstruction result with a reasonable geologic structure and extremely high fidelity of wave field features can be generated.

The invention discloses a random fluctuating wind field simulation method for power transmission tower structure dynamic response analysis. The method comprises the following steps: S1, determining a wave number-frequency spectrum model of a fluctuating wind field of the power transmission tower according to wind field parameters of a place where the power transmission tower is located; s2, based on the eta point set and a wave number-frequency spectrum power function, non-uniform discretization is carried out on a wave number-frequency domain, and wave number-frequency discrete points are obtained; and S3, simulating the fluctuating wind field in the one-dimensional vertical space of the power transmission tower by adopting a spectral expression method. According to the method, the characteristic of height uniformity of the eta point set is utilized, so that the precision of wind field simulation is ensured while fewer discrete wave number-frequency points are used, and the calculated amount of harmonic superposition is greatly reduced in the process of power transmission tower wind field spectrum expression simulation, so that the simulation efficiency of the wind field is improved.

The invention provides a frequency-wavenumber domain-based matrix array full-matrix capture three-dimensional ultrasonic imaging method and system, and belongs to the technical field of ultrasonic imaging. Five-dimensional time-space domain radio frequency data are collected through a two-dimensional full-sampling matrix array and are converted into a frequency-wave number domain through five-dimensional fast Fourier transform, angular frequency is mapped into axial wave number through Stolt interpolation, data resampling is achieved, then three-dimensional image frequency spectrum data are synthesized by defining spatial wave number coordinates of an image and conducting coherence stacking, and the image frequency spectrum data are obtained. And finally, reconstructing a high-quality three-dimensional ultrasonic image through three-dimensional inverse fast Fourier transform. Compared with a traditional delay summation beam forming method, the method has the advantages that the calculation efficiency is remarkably improved, the speed is increased by several times to dozens of times, the spatial resolution, the contrast noise ratio and the spot signal-to-noise ratio are remarkably improved, and the method is suitable for high-resolution real-time three-dimensional and four-dimensional ultrasonic imaging application.

The invention relates to a Lamb wave frequency dispersion characteristic sparse reconstruction method based on non-convex optimization, and belongs to the technical field of high-end equipment perception and test.The method comprises the steps that excitation is conducted on a U-shaped suspension arm, and time domain sampling signals at the junction of a U-shaped suspension arm structure are obtained; the time domain sampling signals are preprocessed; constructing a non-convex sparse frequency dispersion curve estimation model based on the preprocessed signal; and converting the non-convex sparse frequency dispersion curve estimation model into an optimization problem, searching a convex preserving condition of a target function, selecting a forward and backward splitting algorithm for solving, obtaining frequency-wave number distribution, and further extracting a frequency dispersion curve. The method solves the problems that a traditional convex regularization method is limited in induction sparse capability and underestimates signal amplitude. Through a non-convex optimization method, the purpose of accurate estimation of the frequency dispersion curve based on fewer measurement points is achieved.

The invention provides a method for processing and imaging sound wave remote detection data in a well, and belongs to the field of oil-gas exploration. The method comprises the following steps of: firstly, evaluating full wave train data by adopting an energy statistical method, eliminating data with problems, and recovering damaged data by utilizing a POCS interpolation method; then, completing reflection wave extraction based on polarization characteristics of direct waves and reflection waves, and performing optimization through multi-stage filtering processing and amplitude restoration; and then, optimizing a speed extraction method in dynamic correction by using a competition-cooperation particle swarm algorithm, and providing an accurate stacking speed for common sound source gather stacking. And finally, utilizing a frequency-wavenumber domain migration method to realize well-side fracture cavity imaging. According to the method, through an innovative data processing algorithm and flow optimization, imaging of the well-side fracture-cavity body is achieved, and a reliable technical means is provided for oil-gas exploration of a complex reservoir.

This invention relates to a method for extracting guided wavenumbers in variable thickness structures, belonging to the field of nondestructive testing technology. It utilizes an excitation signal to generate guided wave signals on the surface of the object under test. Based on a set window function and a set window length, the two-dimensional wavefield signal obtained by sampling the guided wave signal is divided into multiple local wavefield signals. A short-space Fourier transform is then performed on each local wavefield signal to obtain the amplitude-frequency-wavenumber array corresponding to each local space. A frequency range is then selected by summing the amplitude-frequency-wavenumber arrays of each local space, using the center frequency of the excitation signal as the midpoint. The weighted average wavenumber of each frequency point within the selected frequency range is calculated using the amplitude corresponding to the wavenumber sequence at each frequency point as the weight. The mean of the sum of the weighted average wavenumbers at each frequency point is then used as the wavenumber corresponding to the center position of each local signal space, and these values ​​are sequentially spliced ​​to obtain the wavenumber curve of the guided wave signal. This method solves the problem in existing technologies where the obtained guided wavenumbers cannot accurately reflect the thickness changes of variable thickness structures.

The application discloses a space-time-frequency wave number domain synchronous extrusion superimposed dense channel reservoir characterization method, which comprises the following steps: S1, inputting original seismic data s(t, x) to be analyzed, wherein x is the number of channels, and t is time; S2, obtaining short-time Fourier transform V2(t, x, omega, k) of the signal by using a Gaussian window function g(t, x), wherein omega is frequency, and k is wave number; S3, calculating an instantaneous frequency and an instantaneous wave number estimation of the signal according to S2; S4, defining an extrusion operator TDSTO according to S3; S5, constructing a space-time-frequency wave number domain synchronous extrusion transform characterization method TDSST(t, omega, x, k) according to the extrusion principle; S6, fixing the wave number by maximizing the energy of TDSST, so as to obtain a TDSST time-space-frequency domain three-dimensional data body TDSST(t, x, omega, k(t, x, omega)); and S7, performing feature extraction of frequency f0 in S6, so as to obtain a single-frequency profile to characterize the morphological structure characteristics of a complex superimposed dense channel sandstone reservoir. The application can provide a space-time reservoir morphological characterization method considering the lateral characteristics of a complex superimposed dense channel sandstone.

A video device (100) is provided with: a plurality of transmitters (101) for transmitting fluctuations to a measurement region; a plurality of receivers (102) that receive fluctuating scattered waves from the measurement region; and an information processing circuit (103) that uses the measurement data of the scattered waves to image an object in the measurement region, and that derives a scattered field function using the measurement data and the velocity vector of the object. Deriving a video function that is determined using the amount that is output from the scattered field function by inputting the position of the object to be video into the scattered field function, and video the object in the measurement region using the video function; the information processing circuit (103) changes the number of scattered waves reflected by the scattered field function due to the Doppler effect corresponding to the velocity vector.

The application provides a feature reflection wave travel time inversion method, device, equipment and storage medium. The method comprises the following steps: performing extended domain migration imaging on observed seismic shot gather data to obtain an extended domain common imaging point gather; selecting a target stratum from a zero offset profile of the extended domain common imaging point gather based on the extended domain common imaging point gather; selecting a common imaging point gather of the target stratum from the extended domain common imaging point gather based on the target stratum; performing de-migration forward simulation on the common imaging point gather to obtain feature reflection wave data of the target stratum; establishing a travel time inversion framework based on the feature reflection wave data to obtain simulated reflection wave data; calculating reflection wave travel time difference based on the feature reflection wave data and the simulated reflection wave data; calculating a gradient of travel time inversion based on the travel time difference; and updating a velocity model based on the gradient of travel time inversion. The method can restore various wave number components of the velocity model, thereby improving inversion efficiency.

The invention discloses a fluid filling pipeline guided wave modal tracking method based on an optimal transmission distance. The method comprises the following steps: firstly, calculating a characteristic solution by using a semi-analytical finite element method; secondly, introducing an optimal transmission distance (EMD) in an optimal transmission theory as modal similarity measurement, measuring modal similarity by calculating the minimum physical work of kinetic energy distribution rearrangement, and combining a section geometric geodesic distance to enhance physical interpretability; meanwhile, a continuous stabilization strategy based on wavenumber prediction is introduced to suppress mismatching; and finally, realizing global optimal modal matching by using a Hungary algorithm. According to the method, modal exchange errors can be effectively eliminated, tracking robustness is kept in a complex frequency dispersion area, and meanwhile, visual graphic interpretation of modal evolution can be provided.

The invention provides an elastic wave CT (Computed Tomography) signal denoising method for improving the inversion precision of karst cave exploration, which belongs to the technical field of engineering exploration and comprises the following steps of: performing normalization preprocessing on a received signal, suppressing random noise by adaptive filtering, converting a time domain signal into an F-K domain and identifying a surface wave, filtering a surface wave signal by a Gaussian band elimination filter, and converting the F-K domain signal back to the time domain. The method comprises the following steps of: highlighting a take-off characteristic of a weak P wave covered by a strong surface wave through normalization; using a minimum mean square error adaptive filter to accurately filter random noise; subsequently converting to an F-K domain, and smoothly attenuating the surface wave through a Gaussian band elimination filter by utilizing the frequency wave number distribution difference of the surface wave and the P wave; the method realizes layered removal of random noise and surface waves, finally outputs high-purity direct P-wave signals, greatly improves P-wave travel time pickup and wave velocity calculation precision, and remarkably improves inversion precision of an elastic wave CT exploration method in detection of geologic bodies such as karst caves and fracture zones.

The application discloses a filtering method in marine climate research and application thereof, and comprises the following steps: (1) Fourier transform is conducted on time or space data x to obtain results in a frequency domain or a wave number domain; (2) a function is used to conduct a truncation processing on the frequency domain or the wave number domain, and low-pass and high-pass filtering frequency domain results are obtained respectively; (3) energy of the low-pass and high-pass filtering results is calculated and analyzed; (4) inverse Fourier transform is conducted on the low-pass filtering frequency domain or wave number domain result, and a large-scale characteristic distribution in the time or space domain can be obtained; (5) when the filtering result meets a filtering quality index in the frequency or wave number domain obtained in the step (3) by 20%, the function parameter in the step (2) belongs to a range of function parameters suitable for scale division; and (6) the filtering result is brought into a marine climate system of the earth to improve accuracy. The application aims to provide a filtering method and application for judging and improving scale division accuracy in marine climate research.

The invention discloses a ship identification method and system based on dual physical model fusion verification, and relates to the technical field of optical fiber sensing and ocean monitoring. The method aims at solving the problems that ship noise in submarine optical cable DAS signals is prone to interference of environment noise such as low-speed guided waves, the false alarm rate is high, and accurate positioning cannot be achieved. The method comprises the steps that DAS phase data is acquired and preprocessed to generate a channel-time energy diagram; suppressing low-speed guided wave noise through frequency-wavenumber domain slowness scanning, detecting an energy ridge line on the purified energy diagram, and calculating an observation apparent velocity; acquiring AIS data to calculate and predict the apparent velocity; the relative error of the observed apparent velocity and the predicted apparent velocity is firstly matched, and then linear regression is carried out on the arrival time of the candidate ridge line and the AIS distance to verify and fit the root mean square error; the validated ridge line is confirmed as an effective ship event and a directory is generated. According to the technical scheme, through physical mechanism denoising and dual model verification, the accuracy and reliability of ship identification are effectively improved.

This application provides a method and system for locating low-frequency broadband pulse signals under unknown environmental information. The method includes: obtaining the broadband wavenumber spectrum of the pulse signal using a hypothetical dedispersion formula; obtaining preliminary estimates of each mode using a variational Bayesian method; searching for optimal waveguide invariants and reference sound velocities; re-estimating each mode; and finally, substituting the wavenumber and phase information of each mode into a distance estimation cost function to estimate the distance to the pulse source. The advantages of this application are: by using the dedispersion formula and the variational Bayesian method, mode separation of broadband pulse signals can be achieved even with a short horizontal array aperture; and the phase difference between different modes can be used to estimate the source distance. Through this method, the optimal waveguide invariants and reference sound velocities can be searched under conditions lacking environmental information, reducing the error in source distance estimation.

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.

This invention provides a method and system for three-dimensional high-frequency controllable source electromagnetic numerical simulation with dual-source excitation. The method constructs a single-source electromagnetic field grid model and a three-dimensional spatial domain total field equation, converting it into a quadratic field control equation by separating the primary and secondary fields. Fourier transform is used to transform the equations to the wavenumber domain, and shape function interpolation is performed iteratively to obtain the wavenumber domain numerical solution. A single-source forward excitation solution is obtained through inverse Fourier transform. A single-source reverse excitation solution is obtained by reversing the current direction. Finally, the forward and reverse solutions are vector-superimposed to obtain the dual-source excitation simulation result. This invention solves the technical problems of existing technologies, such as simulation distortion caused by neglecting displacement current in the 100kHz-10MHz frequency band, numerical calculation interference from field source singularities, and low spatial resolution of complex underground anomalies under single-source excitation. By simulating gradient field response characteristics through dual-source reverse superposition, the sensitivity and exploration resolution of underground target boundaries are significantly improved.

The present application relates to a kind of low wave number initial model establishment method, device and medium of depth domain seismic inversion, method includes: extracting the well logging data of target layer section, the logging attribute of the low wave number initial model to be built is calculated;According to well seismic calibration result, extract well seismic data of well logging section, and according to logging depth sampling interval is up-sampled, obtain well seismic data on well logging depth section;The amplitude attribute of well seismic data section is calculated;With the amplitude attribute of well seismic data section as independent variable, the logging attribute of the low wave number initial model to be built is function value, linear fitting is carried out;From the fitting result, the amplitude attribute of optimal linear correlation is selected, and the optimal slope parameter and optimal intercept parameter are obtained;The minimum amplitude attribute of corresponding depth domain seismic data body is calculated;Depth domain seismic data body amplitude attribute is converted into the attribute of low wave number initial model to be built;Low wave number initial model attribute to be built is carried out Gaussian filtering, and the final low wave number initial model is obtained.

The invention relates to the technical field of radar detection and application thereof, and provides a data processing method for ground penetrating radar data, which comprises the following steps: firstly, carrying out preprocessing and gain recovery on original radar data, then carrying out mixed denoising by adopting empirical mode decomposition and Wiener filtering, screening and reconstructing effective signals by utilizing cross correlation coefficients, and finally, carrying out data processing on the effective signals. The method comprises the following steps: firstly, constructing a filter in a frequency wave number domain to remove random noise and low-frequency interference, then constructing a fan-shaped filter in the frequency wave number domain, introducing transition band processing to suppress linear coherent interference, then carrying out wave velocity estimation based on an image entropy minimization criterion to obtain an optimal imaging wave velocity, and finally executing Stolt migration imaging by using the optimal wave velocity and outputting a result. According to the method, edge details of effective signals are effectively reserved through multi-scale decomposition, the problem of imaging defocusing is solved through self-adaptive wave velocity search, and the signal-to-noise ratio, the resolution ratio and the focusing precision of ground penetrating radar data in a complex and strong-interference environment are remarkably improved.

The invention discloses an OCT (Optical Coherence Tomography) broadband spectrum nonlinear interference signal wave number calibration and depth reconstruction method, which comprises the following steps of: acquiring a nonlinear sampling broadband spectrum interference signal of a single-layer reflector by utilizing frequency domain OCT, removing a direct current item and noise, performing continuous wavelet transform on the signal by taking Morlet complex wavelet as a generating function, detecting and extracting a wavelet ridge line through a wavelet modulus maximum value, and reconstructing the depth of the nonlinear sampling broadband spectrum interference signal. Obtaining an instantaneous phase and converting the instantaneous phase into wave number information, and constructing a nonlinear Fourier transform matrix by using the wave number information; and testing a to-be-tested target by using the frequency domain OCT to obtain a broadband spectrum interference signal, carrying out direct current item removal and noise treatment, and then carrying out calculation by using the constructed nonlinear Fourier transform matrix to obtain reflectivity information of the to-be-tested target at different depths. According to the method, the problem of nonlinear sampling of the broadband spectrum interference signal in the frequency domain OCT is solved, and high-precision wave number calibration and accurate recovery of nonlinear sampling interference signal depth information are realized.