846 results about "Geodesy" patented technology

The invention discloses an InSAR (Interferometric Synthetic Aperture Radar) and GNSS (Global Navigation Satellite System) robust adaptive fusion method for high-precision three-dimensional surface deformation monitoring, which belongs to the technical field of space geodetic survey and remote sensing, and comprises the following steps of: constructing a time sequence deformation field based on a small baseline set time sequence InSAR processing technology, and resampling GNSS observation data to a grid consistent with the InSAR by using Kriging spatial interpolation; a Helmert variance component estimation method is adopted to adaptively estimate variance components of InSAR and GNSS data, and reasonable weight fixing of a multi-source observation value is achieved; iterative reweighted least square estimation is introduced, outlier influences in various observation data are dynamically restrained through a Tukey double-weight function, robust optimization of a fusion model is achieved, and therefore a high-precision three-dimensional deformation field is reconstructed. The method comprises the steps of Helmert weight fixing, IRLS robust, variance component updating and a closed-loop feedback mechanism of weight matrix optimization.

The invention provides a three-dimensional flow velocity field inversion method and system based on coherence-assisted SAR adaptive offset tracking, and a storage medium, so as to improve the problem of decoherence caused by time heterogeneity in a synthetic aperture radar interferometry (InSAR) technology, and relates to the technical field of electric power engineering geological disaster monitoring. The method comprises the following steps: establishing a space-time baseline network between SAR images; obtaining coherence data after geocoding; selecting a pair of SAR images from the space-time baseline network as a reference unit, calculating an offset result through a pixel offset tracking technology for the reference unit, and determining an optimal matching window parameter of the reference unit; establishing a mapping function between coherence and optimal matching window parameters; according to the established mapping function, calculating an optimal matching window parameter of each remaining SAR image pair; and in combination with the line-of-sight average deformation rate fields of the ascending and descending rails, inverting a surface three-dimensional flow velocity deformation field of the target area through a geometric vector decomposition model.

The invention relates to the technical field of deep foundation pit monitoring, and discloses a deep foundation pit micro-deformation InSAR remote sensing method and system and a storage medium. The method comprises the following steps: synchronously acquiring data through a ground radar and an InSAR satellite, and obtaining a point target set and an original scattering characteristic sequence; obtaining a net deformation phase sequence without influence through registration residual correction, phase deduction and weighted fusion processing, and carrying out phase consistency quality control and unwrapping processing to obtain an accumulated deformation quantity sequence; constructing a deformation mechanism correlation characteristic matrix based on correlation analysis of the cumulative deformation quantity sequence and radar reflection intensity amplitude attenuation characteristics; and fusing the deformation mechanism correlation feature matrix and the real-time differential interference image group, extracting deformation rate and phase gradient change, and further generating a deep foundation pit local micro-deformation risk early warning index set. According to the invention, the tiny deformation of the surrounding area of the foundation pit can be effectively detected, potential risks are warned in advance, and the safety and monitoring efficiency in the construction process are improved.

The embodiment of the invention provides a method and a system for processing ocean multi-parameter wide-width imaging data of a Beidou interference altimeter, and relates to the technical field of GNSS-R. The method comprises the following steps: constructing a set of complete and strict data processing chain; through precise calculation of a GNSS-R direct signal and reflection signal delay inequality interference geometric model and correction of actual measurement delay inequality errors, systematic errors in the process of converting original signals into geophysical parameters are remarkably reduced, and the precision and reliability of parameters obtained through final inversion are improved. And on the basis of a self-adaptive data collection strategy of parameter space-time transformation characteristics and a spatial interpolation reconstruction process, discrete strip data is intelligently'woven 'into a continuously-covered wide-swath imaging type data product, so that the data format better meets the analysis requirements in the fields of oceanography, climatics and the like, and the intuition and usability of the product are improved.

The invention belongs to the technical field of geophysical exploration, and particularly relates to a method and a system for quantitatively measuring intrinsic parameters of bedrock fractures based on all-angle domain amplitude and phase joint inversion. The method comprises the following steps: acquiring geological radar reflection data of a target fracture at a plurality of incident angles; synchronously extracting data sequences of amplitude change along with angles and phase change along with angles; establishing a fracture reflection theoretical model; constructing a target function jointly constrained by the amplitude residual and the phase residual; and adopting a global optimization algorithm to minimize the target function, and synchronously inverting to obtain the fracture opening, the filler dielectric constant and the occurrence parameter. According to the invention, through the combined constraint of the amplitude and the phase, the inversion multiplicity is effectively inhibited; accurate system effect correction is adopted, so that the reliability of large-angle data is ensured; according to the method, multi-parameter synchronous inversion without prior hypothesis is realized, the inversion precision is small in relative error, the method is obviously superior to a traditional method, and the method has important engineering application value.

The invention discloses a load saturation prediction method and device for satellite communication, and the method is applied to a gateway station, and comprises the steps: constructing a first data set corresponding to each sampling time period, and carrying out the fitting of each first data set, and generating a corresponding reference distribution curve; constructing a second data set corresponding to the current time period, and fitting the second data set to generate a corresponding first actual distribution curve; pre-determining the channel capacity between each satellite and the corresponding coverage area, and broadcasting the first actual distribution curve, a first reference distribution curve corresponding to the current time period and the channel capacity to each satellite; determining a second reference distribution curve corresponding to the next time period, calculating a deformation parameter between the first reference distribution curve and the second reference distribution curve, and predicting a second actual distribution curve; and selecting a prediction point in the second actual distribution curve, and determining the load saturation of each satellite in the next time period according to the prediction point and the channel capacity.

The invention provides a building three-dimensional automatic modeling method based on oblique photogrammetry data, and belongs to the technical field of three-dimensional modeling. The method comprises the following steps: converting measurement data into three-dimensional point clouds, and obtaining non-ground point clouds in the three-dimensional point clouds by adopting a CSF algorithm; non-ground point clouds are clustered on the XY plane by adopting a DBSCAN algorithm to obtain independent point cloud clusters corresponding to the building, and two-dimensional vector contours of the independent point cloud clusters are obtained by adopting an Alpha Shape algorithm. The independent point cloud clusters are horizontally sliced along the Z axis, a plurality of cross section point sets are extracted in the Z-axis direction, two-dimensional convex hulls of the cross section point sets are calculated, the intersection-to-union ratio of the two-dimensional convex hulls and preset convex hulls is calculated, and when the intersection-to-union ratio exceeds a threshold value for the first time and is kept stable, the boundary height of the building is obtained. And stretching the two-dimensional vector contour along the Z axis according to the boundary height to obtain a first three-dimensional model of the building. According to the invention, manual intervention can be reduced, and the building modeling efficiency is improved.

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

The invention provides a method and a device for correcting a backscattering coefficient of a wide scanning synthetic aperture radar, and the method comprises the steps: determining a backscattering coefficient of each pixel in a to-be-processed image, and then determining a backscattering coefficient correction coefficient of the pixel through employing the data of a mode sea surface wind field, a geophysical model function and the like; and correcting the backscattering coefficient of the pixel by using the backscattering coefficient correction coefficient to generate a target backscattering coefficient corresponding to the pixel. According to the method and the device, the backscattering coefficients of the pixels are corrected, so that the amplitude band inconsistency between the SAR image strips in the wide-width scanning mode can be eliminated, the more uniform and higher-precision SAR image in the wide-width scanning mode can be obtained, and better data support is provided for subsequent processing.

The invention relates to a three-component high-dimensional seismic data reconstruction method and device based on vector tensor reduced rank, a medium and equipment, and the method comprises the steps: firstly, placing three-component high-dimensional seismic data at three imaginary part positions of a quaternion array in a time domain through employing a quaternion tool, and setting a real part of the quaternion array to be 0 value; and realizing vector joint representation of the quaternion-based three-component seismic data. Quaternion representation can protect the nonlinear orthogonal relation of three component seismic data in the space direction, and vector field structural characteristics of seismic waves are kept. And secondly, transforming a time domain quaternion array into a frequency domain by using quaternion Fourier transform, extracting a spatial high-dimensional quaternion array corresponding to each frequency domain slice, and performing missing data restoration and reconstruction on the high-dimensional quaternion array by using a tensor reduced-rank reconstruction technology based on a quaternion parallel matrix decomposition algorithm. According to the method, the memory occupation amount is remarkably reduced, the calculation amount is small, the calculation cost is low, and the reconstruction calculation efficiency is high.

The invention provides a laser radar attitude compensation wind measurement method and system, and relates to the technical field of laser radar wind measurement, and the method comprises the steps: obtaining the six-degree-of-freedom motion data of a carrier; driving the two-axis stabilized platform to compensate based on the three-axis rotation angle so as to stabilize the direction of the laser beam, and scanning to obtain initial radial wind speed data; then synchronizing the translation speed with the wind speed data by using a linear fitting algorithm to obtain a synchronous data set; a laser radar scanning mode is determined, and a corresponding attitude compensation strategy is selected according to the mode; inputting the synchronous data set and the strategy into a pre-training machine learning model to generate a compensation parameter, and processing the initial data through Kalman filtering based on the parameter to obtain target radial wind speed data; and finally, performing wind field inversion calculation based on the target data to generate a wind field measurement result. According to the invention, the wind field measurement precision of the laser radar on the motion platform under the complex sea condition is improved.

The invention provides a physical information long-short-term memory network site seismic response prediction method based on kinematics derivative constraint and data folding. The method comprises the following steps: collecting and preprocessing bedrock and ground surface seismic oscillation data; displacement, speed and acceleration supervision labels are constructed through surface acceleration integration; constructing a long-short-term memory network containing a data folding and unfolding module; physical constraint loss is constructed by utilizing kinematics derivative relations between displacement and speed and between speed and acceleration, and a network is trained by combining an adaptive weight strategy; and outputting earth surface displacement, speed and acceleration time history by adopting the trained network. According to the method, under the condition that constitutive parameters of a soil layer are not explicitly given, physical consistency constraint can be applied to the field seismic response prediction process, the long-time program train training burden is reduced through sequence folding processing, and the method is suitable for field seismic response prediction.

The invention relates to the technical field of disaster site surveying and mapping, and discloses a disaster site three-dimensional surveying and mapping system based on an unmanned aerial vehicle. According to the system, an unmanned aerial vehicle is controlled to carry a multispectral sensor to carry out omnibearing scanning on a disaster site, and earth surface deformation data and multiband reflectivity data are collected; carrying out time sequence alignment and noise filtering on the collected data, constructing an initial three-dimensional point cloud through dynamic voxelization processing, and fusing inertial measurement unit data to correct a spatial position; identifying earth surface fault zone edge and elevation mutation features from the point cloud, and dividing geological stability partitions by combining reflectivity data; a self-adaptive octree segmentation algorithm is adopted to carry out multi-level detail division on the stability partitions, the grid resolution is automatically adjusted according to the fault zone density, and a disaster site three-dimensional topological structure is generated; by comparing historical and current three-dimensional data, the earth surface displacement vector and the settlement gradient are calculated, a potential secondary disaster risk area is marked, and accurate assessment and risk early warning of a disaster site are achieved.

The invention discloses a seismic-electromagnetic joint inversion method based on bidirectional cross gradient constraint. The method comprises the following steps: firstly, establishing three physical field initial models; taking a seismic grid as a main grid, and mapping respective models of a direct current method and a transient electromagnetic method to the main grid through interpolation; constructing a joint inversion objective function containing data fitting terms of the three physical fields, model constraint terms and cross gradient constraint terms between every two of the three physical fields; decomposing the joint inversion objective function into three inversion objective functions which are coupled with one another, and respectively carrying out iterative updating on the seismic velocity model, the direct current method resistivity model and the transient electromagnetic method resistivity model by adopting a collaborative alternative iteration strategy; each iteration update of any model is constrained by the cross gradient exported by the last iteration result of the other two models, so that the structural consistency between the models is forcibly realized; and after completion, outputting a geologic structure model fusing the three-field information, thereby effectively improving the detection precision and interpretation reliability of the geologic structure under the complex geological condition.

The invention discloses a side slope three-dimensional deformation monitoring and early warning method and system based on a cascade single-line multipoint crack meter. The method and system are suitable for deformation monitoring and landslide state forecasting of a side slope area without macroscopic cracks. The method comprises the steps that a monitoring main line is planned in a potential slip area in the main slip direction, intelligent fusion nodes are arranged on the main line, multi-direction monitoring branch lines are led out, and a tree-shaped monitoring network is formed; two ends of the main line are anchored in a stable area, and tail ends of the branch lines are fixed on a slope; the node devices synchronously collect relative displacement and three-dimensional attitude data in the main line direction and the branch line direction, a three-dimensional deformation field is built in combination with an early warning model, and landslide state study and judgment and threat range prediction are achieved. The system achieves multi-point and multi-direction displacement monitoring through a tree-shaped structure, has the advantages of being wide in coverage, low in cost, flexible in arrangement and the like, and supports early-stage automatic early warning of slope disasters.

The invention relates to a subway track line intelligent optimization method based on a hierarchical archiving differential evolution algorithm. The method comprises the following steps: obtaining an actual topographic change curve through topographic relief characteristics of a railway line design plane so as to construct a standardized coordinate axis model; obtaining a terrain height change sequence by discretizing the terrain height change of the coordinate axis model; determining technical limitation conditions of the rail according to a railway line design specification; selecting a construction type according to a construction rule and terrain information; calculating to obtain the total construction cost; according to the terrain height change sequence and the technical limitation condition of the rail, establishing a Ground object and preprocessing the Ground object; and an HAPI-DE algorithm model is established, and the Ground object is sent into the algorithm model to be evolved so as to obtain an optimal solution railway route and the lowest cost. According to the method, the optimal railway longitudinal section design scheme which is lower in cost, more accurate and more reliable can be obtained.

The invention discloses a reflected wave travel time difference measurement method and system based on transverse constraint, and the method comprises the steps: carrying out the migration imaging of observed seismic data under a given speed model, obtaining an imaging result through the cross-correlation of a seismic source end wave field and a receiving end wave field, and obtaining the information of an underground reflection structure based on a multi-attribute Markov decision process. And then solving a background wave field and a primary reflection wave field in a given speed model and an underground reflection structure, and extracting a simulated seismic record at a detection point position. And calculating zero-offset travel time by using the depth of the reflection horizon, and executing transverse tracking in the reflection event direction by taking the travel time as a seed to obtain simulated travel time of each offset position. A local transverse constraint time window is constructed based on simulation travel time, cross-correlation operation is performed on observation data and simulation data in the time window, and a local extreme point is extracted. And finally, performing transverse tracking by taking the zero-offset travel time position as a starting point to obtain continuous and consistent reflected wave travel time difference, thereby realizing stable and accurate measurement of the reflected travel time.

A single-station positioning method based on Doppler frequency shifts and time difference information, and a system and an electronic device. The method comprises: acquiring an actual period of a pilot signal sent by an unmanned aerial vehicle (S1); on the basis of the actual period of the pilot signal sent by the unmanned aerial vehicle, determining two adjacent target signals sent by the unmanned aerial vehicle in a target area, and acquiring actual center frequencies of target signal information, an actual time difference of arrival between the two target signals at a receiving base station, actual angles of arrival, and actual Doppler frequency shifts of the unmanned aerial vehicle (S2); importing an actual range of the target area, the actual period of the pilot signal, the actual center frequencies of the target signals, the actual time difference and the actual angles of arrival into variables in a preset theoretical Doppler frequency shift model, so as to obtain a theoretical Doppler frequency shift value table, wherein the table comprises different theoretical Doppler frequency shift values corresponding to different distances from the unmanned aerial vehicle to the receiving base station (S3); performing error metric calculation processing on the actual Doppler frequency shifts and each theoretical Doppler frequency shift value, so as to obtain a theoretical Doppler frequency shift value corresponding to the minimum error metric (S4); and on the basis of the theoretical Doppler frequency shift value corresponding to the minimum error metric, determining the positions at which the unmanned aerial vehicle sends the target signals (S5). The method can ensure the accuracy of single-station positioning in a complex scenario.

The invention discloses an iron tower PS pixel point simulation and differential interference deformation monitoring method and system based on three-dimensional point cloud data, and the method comprises the steps: obtaining the point cloud data of an iron tower body structure under a radar coordinate system through a radar distance Doppler positioning method through the three-dimensional point cloud data of the iron tower body structure; identifying pixels containing the iron tower point cloud data in the SAR image as iron tower PS pixel points; constructing a radar scattering weight model based on elevation information of the iron tower point cloud data, and obtaining an iron tower PS pixel point scattering intensity value through weighted summation in combination with a weight adjustment parameter; constructing a complex value of the point cloud data of the iron tower according to the interferogram parameters and a radar scattering weight model, and obtaining a simulation interference elevation phase of the PS pixel point of the iron tower by using complex field summation; differential interferometric phases of iron tower PS pixel points in the interferogram are obtained through difference; and performing statistical analysis on the differential phase map of the iron tower PS pixel point set region in the time sequence interferogram to obtain a deformation result.

A method performed by a wireless device includes receiving ephemeris data from a network node associated with a first cell. The ephemeris data is associated with a satellite serving a second cell, and the ephemeris data is coarse ephemeris data comprising less than full accuracy than full ephemeris data. The wireless device uses the coarse ephemeris data to locate a beam associated with the second cell and synchronizes with the second cell.

The invention provides an earth surface deformation inversion method based on an overlapping region splicing InSAR deformation field. The earth surface deformation inversion method comprises the following steps: firstly, acquiring multiple frames of LOS images, image frame information and GNSS observation data in a research region range; classifying the LOS images to obtain a same-track data set, and screening overlapped regions to obtain the overlapped regions of the same track and corresponding GNSS station data; an observation equation is constructed to solve and eliminate errors of the same-rail data, and the same-rail data are spliced; performing spatial interpolation and projection on the GNSS three-dimensional deformation observation data to obtain GNSS three-dimensional deformation projection data; in combination with GNSS three-dimensional deformation projection data, weighted balancing of an error plane is carried out, and adjacent orbit data splicing is executed; earth surface three-dimensional deformation inversion is carried out, earth surface three-dimensional deformation data are obtained, strain rate data are calculated, and an evaluation result of earth surface three-dimensional deformation inversion is obtained. According to the method, the precision of three-dimensional deformation inversion of the earth surface is improved.

The invention relates to the technical field of side slope deformation space-time prediction and multi-source data fusion, in particular to a side slope deformation space-time prediction method and device based on multi-scale data fusion. The method comprises the following steps: acquiring slope synthetic aperture radar interference data and global navigation satellite data, determining a measuring point displacement time sequence and coordinates according to the data, and aligning a displacement direction with a time sequence window; extracting time sequence features of a displacement time sequence in the global navigation satellite data and generating a feature matrix, and fusing the feature matrix with the synthetic aperture radar interference data to obtain a feature fusion matrix; an intra-set adjacency matrix of the former measuring point is constructed according to the measuring point coordinates, an inter-set adjacency matrix of the two types of data measuring points is constructed through a Thiessen polygon graph, and the intra-set adjacency matrix and the inter-set adjacency matrix are fused to form a fused adjacency matrix. And based on the feature fusion matrix and the fusion adjacency matrix, aggregating space neighborhood information to generate a displacement time sequence prediction result of the slope global measurement point. Therefore, the problems of insufficient space-time correlation excavation, low slope global deformation prediction precision and the like in related technologies are solved.

The invention discloses a diffracted wave-based heterogeneous body prediction method, electronic equipment, a storage medium and a diffracted wave-based heterogeneous body prediction device. The method comprises the following steps: acquiring a diffracted wave three-dimensional seismic data volume, a horizon of a target interval and a comprehensive wavelet based on a full-wave-field three-dimensional seismic data volume; constructing a longitudinal wave velocity and density conversion relation based on the logging longitudinal wave velocity and the corresponding density, and further constructing a longitudinal wave velocity and longitudinal wave impedance conversion relation to convert the interval velocity of the full-wave-field three-dimensional seismic data volume into interval velocity low-frequency impedance; constructing a diffracted wave impedance inversion objective function based on the longitudinal wave diffraction coefficient and the comprehensive wavelet, and taking the interval velocity low-frequency impedance body as the background impedance of the inversion objective function; and solving the diffracted wave impedance inversion objective function based on an optimization algorithm, carrying out diffracted wave impedance inversion, and obtaining diffracted wave impedance to carry out heterogeneous body prediction. Quantitative prediction of the heterogeneous body is realized through the diffracted wave three-dimensional seismic data body, and the impedance inversion precision of the heterogeneous reservoir is improved.

The invention discloses a coherent radial velocity spectrum target detection method based on a backtracking filter. The method comprises the following steps: constructing a sea surface echo model and a multi-frame joint detection problem; the method comprises the following steps: based on a sea surface echo model and a multi-frame joint detection problem, performing coherent accumulation on a radar receiving pulse sequence, and calculating test statistics of intra-frame coherent accumulation under each frame, each distance unit and each radial speed; based on the inconsistency of the sea clutter structural texture between the slow time Doppler velocity and the frame time radial velocity and the velocity consistency of the target on two time scales, inter-frame accumulation is carried out on the test statistics of each frame, each distance unit and each radial velocity in a distance-frame time plane along the radial velocity direction; constructing a coherent radial velocity spectrum; and in an inter-frame accumulation process, replacing a linear search accumulation mode of the coherent radial velocity spectrum between scanning by using the backtracking filter bank to obtain a coherent radial velocity spectrum based on the backtracking filter bank, and carrying out target detection by using the coherent radial velocity spectrum.

The application provides a complex terrain area elevation inversion method based on small unmanned aerial vehicle borne interferometric SAR. N-track small unmanned aerial vehicle flight tracks are set, the X-axis and Y-axis coordinates of the starting point and ending point of each track small unmanned aerial vehicle are unchanged, and the value of the Z-axis is gradually increased along the vertical direction. N-track echo data are obtained by the small unmanned aerial vehicle borne SAR according to the set tracks, and BP algorithm is used for imaging processing to obtain N SAR images. The offset of the SAR images of adjacent short baseline small unmanned aerial vehicles is obtained. The offset of the SAR images of long baseline small unmanned aerial vehicles is calculated based on the offset of the SAR images of small unmanned aerial vehicles, and the SAR images of long baseline are image matched. The SAR images after image matching are processed to obtain corresponding interference phase images, and the interference phase images are filtered by using Goldstein algorithm. Multi-baseline phase unwrapping is performed based on the interference phase of the interference phase images. The complex terrain area elevation inversion method of SAR is realized.

The application discloses a method for suppressing ghost waves of a streamer based on frequency-slowness domain continuation and ghost wave parameter optimization. Considering that ghost wave parameters vary with time and space and affect the ghost wave removal effect, first, the marine seismic data is divided by using a sliding time window to reduce the influence of the time variation of the ghost wave parameters; the delay time and amplitude difference coefficient of the ghost wave and the primary wave are enumerated alternately; the data in the time window is continued by using a frequency-slowness domain wave field continuation operator, and then the original marine seismic data is added and subtracted by the continued record; the sum of the absolute values of the products of the new records obtained by operation is taken as the target, the optimal ghost wave delay time and amplitude difference coefficient are iteratively searched, and then the ghost wave is suppressed; after the data in each window is suppressed, all the window data are integrated into a complete record. The method can effectively suppress the ghost wave, widen the frequency band and improve the resolution of the seismic data.

The invention provides a method for obtaining a three-dimensional body through two-dimensional seismic reflection data gridding, and belongs to the field of seismic data processing. Through a technical means of constructing a virtual grid through rotation, a two-dimensional measuring line is converted into a three-dimensional data body which is uniformly distributed in space, and subsequent data processing is not limited to one direction any more. Specifically, latitude and longitude coordinates of each discrete reflection data are obtained and converted into UTM projection rectangular coordinate values, and coordinate values of a longitude minimum point are determined; calculating the optimal azimuth angle of the grid according to the plane distribution characteristics of the discrete points, and determining the grid spacing along the azimuth angle direction and the vertical azimuth angle direction; the grid is rotated to the three-dimensional grid direction, the coordinates of the point with the minimum longitude value after rotation are extended outwards by a half grid interval, and the coordinates of the point serve as base points of the grid; all the coordinate points are rotated in the same way, data falling into the same grid are superposed, then the coordinate number of each point in the grid is calculated, and a three-dimensional data volume is obtained accordingly.

The invention discloses a relative motion positioning method, system and device of a multi-seabed beacon network, a medium and a product, and relates to the technical field of ocean geodetic measurement data processing, and the method comprises the steps: determining a priori value of sound velocity deviation of the multi-seabed beacon network; respectively constructing a quadratic polynomial model of each seabed beacon and a grouping adjustment function model of each sampling observation moment; initializing a quadratic polynomial model of each seabed beacon and a grouping adjustment function model of all sampling observation moments; performing iterative solution on the quadratic polynomial models of all the seabed beacons and the grouping adjustment function models of all the sampling observation moments by using a least square method so as to obtain seabed beacon motion models of all the seabed beacons; and when the current seabed beacon is observed at the current moment, substituting the current moment into the seabed beacon motion model of the current seabed beacon to obtain a motion three-dimensional coordinate of the current seabed beacon at the current moment. According to the invention, the relative positioning precision of the seabed beacons is improved.

The invention discloses a depth domain non-zero phase seismic wavelet extraction method and device. The method may include: decomposing an amplitude spectrum of a depth domain seismic wavelet through a depth-wavenumber spectrum; estimating the phase spectrum of the seismic wavelet according to the phase shift characteristic of the depth domain seismic wavelet; and according to the phase spectrum, through amplitude phase matching fusion, depth domain non-zero phase seismic wavelets are extracted. The depth domain non-zero phase seismic wavelet extraction method can improve the accuracy of depth domain non-zero phase seismic wavelet extraction.

The invention discloses a long-distance and short-distance combined advanced geological forecasting method and system and storage equipment, belongs to the technical field of geological forecasting, and can solve the problem that long-distance accurate forecasting is difficult in the prior art. The method comprises the following steps: S1, constructing a three-dimensional geological model within a first length range from a tunnel face in a to-be-constructed area according to exposed geological information of a constructed area, and predicting a first heteroplasmon within the first length range based on the three-dimensional geological model; s2, predicting a second heteroplasmon within a first length range by using an electromagnetic prospecting method, and correcting the three-dimensional geologic model by using the second heteroplasmon; s3, determining a wave velocity parameter within a second length range from the tunnel face by using a seismic exploration method, and predicting a third heteroplasmon within the second length range based on the wave velocity parameter; and S4, correcting wave velocity parameters according to the first heteroplasmon, the second heteroplasmon and the third heteroplasmon, and carrying out geological forecast according to the corrected wave velocity parameters and the corrected three-dimensional geological model. The method is used for geological forecast.