96 results about "Gravity anomaly" patented technology

The invention discloses a submarine topography super-resolution reconstruction method based on window displacement multi-source fusion, and aims to solve the problem of precision attenuation caused by absence of multi-beam truth value evaluation and resolution improvement of an existing neural network submarine topography reconstruction method. The method comprises the following steps of: obtaining gravity anomaly, gravity vertical gradient, vertical line deviation component, submarine topography background and ship survey water depth multi-source data of a target sea area; constructing a residual U-Net attention neural network model, and performing training by taking multi-source data extracted by a 11 * 11 window and position codes as input; moving the trained model sampling window in a staggered manner according to a target resolution step length; extracting data points in the window during movement and injecting position codes; and outputting a high-resolution water depth predicted value based on the position code and the window data. The reconstruction precision is improved by fusing physical quantities, the multi-resolution robustness is guaranteed by a window displacement mechanism, and the medium-frequency feature reconstruction capability is enhanced by a residual U-Net attention model.

The invention relates to the technical field of earthquake monitoring, and discloses a non-constructional gravity effect correction system applied to flowing gravity observation, which comprises the following steps of: calculating a diffusion-non-diffusion deviation ratio representing dynamic properties by inverting a global optimal diffusivity adaptive to a current field; the diffusion non-diffusion deviation ratio quantifies the proportion of non-diffusion components which cannot be explained by the diffusion model in the gravity field change, so that whether the dominant mechanism of gravity anomaly is hydrological noise residue or real construction activity is quantitatively judged. According to the invention, self-inspection of the credibility of the correction result is realized, false signals are effectively eliminated, and the reliability of earthquake precursor monitoring is improved.

The invention relates to the technical field of underground medium density inversion, in particular to an underground medium density inversion method and system based on high-precision microgravity, and the method comprises the steps: obtaining the gravity data of each observation point in an underground medium density inversion ground surface measurement region; calculating a superposition abnormal value of each gravity subsequence of each observation point, extracting a superposition characteristic value of each modal component, obtaining the superposition interference confidence of each gravity subsequence of each observation point, obtaining the superposition influence persistence of each observation point, adjusting the wavelet threshold of the wavelet transform of the gravity data on each observation point, and calculating the superposition abnormal value of each gravity subsequence of each observation point. And correcting the gravity data of each observation point after denoising processing to extract Bouguer gravity anomaly of each observation point, and obtaining an inversion result of the underground medium density through a Parker-Oldenburg density interface inversion algorithm. According to the method, the underground medium density can be effectively and accurately inverted through the high-precision microgravity data.

The invention discloses a geologic body boundary extraction method based on gravity total horizontal derivative anomaly, and belongs to the technical field of geophysical data processing and interpretation. Calculating gravity total horizontal derivative anomaly based on actually measured gravity anomaly data; in the sliding window, detecting whether the gravity total horizontal derivative abnormity meets a preset extreme value structure condition in multiple directions, and screening out boundary points based on an amplitude threshold value and a slope threshold value; according to the invention, a geologic body boundary extraction method based on gravity total horizontal derivative anomaly is established by using directional extremum features of gravity total horizontal derivative anomaly; compared with a traditional boundary recognition method, the method does not depend on prior information, is high in anti-noise performance, can quickly and accurately recognize the boundary position of the geologic body under a complex background, and has higher boundary resolution.

The invention belongs to the technical field of geophysical inversion, and relates to a gravity anomaly three-dimensional inversion method based on multi-scale feature fusion, which comprises the following steps: randomly generating density anomaly bodies with different numbers, densities, sizes and shapes to obtain an underground medium three-dimensional density model; carrying out forward modeling on the three-dimensional density model of the underground medium to obtain two-dimensional gravity abnormal data of the earth surface; a multi-scale feature fusion model is trained with the underground medium three-dimensional density model as output and the two-dimensional gravity anomaly data as input, and the multi-scale feature fusion model comprises a two-dimensional encoder, a three-dimensional decoder, a dimension converter and a multi-scale feature fusion device, wherein the dimension converter and the multi-scale feature fusion device are located between the two-dimensional encoder and the three-dimensional decoder; and inputting actually measured two-dimensional gravity anomaly data into the trained multi-scale feature fusion model, and performing inversion to obtain a three-dimensional density body constrained by a physical range. According to the method, the expression capability of the multi-scale features of the underground structure is enhanced, and the multi-solution problem and the deep resolution reduction problem of inversion can be relieved.

The invention relates to a gravity background field construction method based on iterative optimization Kriging. The method comprises the following steps: S1, obtaining in-situ gravity measurement data as a known point set; s2, based on the known point set, performing preliminary interpolation on grid points of the gridding background image by adopting a Kriging method to obtain a Kriging estimated value of each grid point; s3, by introducing a local gravity anomaly covariance and a Kriging variation function, evaluating the random uncertainty of a Kriging estimated value in an interpolation process, and dividing each grid point into a reference point and a point to be corrected based on the random uncertainty; and S4, combining the reference point and the known point set to construct a continuous distribution constraint curved surface, carrying out local continuity correction on the continuous distribution constraint curved surface based on the obtained reference point, and carrying out iterative expansion on the basis of the corrected local constraint curved surface until the continuous correction of the whole continuous distribution constraint curved surface is completed. And the construction of the gravity background field is realized.

The invention discloses a gravity anomaly upward continuation calculation method based on a central grid data block separation method, belongs to the technical field of geophysical exploration data processing, and solves the problems that when an existing method is used for processing an ultra-low altitude upward continuation problem, singularity occurs when an integral kernel function coincides with a calculation point and a data point; the method comprises the following steps: stripping a central data grid where a calculation point is located from an integral domain based on a block separation method, and carrying out practical modification processing on a global integral model based on the block separation method to obtain a modified model without integral singularity; according to the method, a practical modification scheme for converting a gravity anomaly upward continuation global integral model into a local numerical integral model is provided, numerical calculation verification is carried out on a kernel function singularity solution and the practical modification scheme of a calculation model by utilizing an ultrahigh-order bit model EGM2008, and model internal coincidence precision superior to 1mGal can be obtained; therefore, the high-precision earth external gravitational field assignment requirement is met.

The invention belongs to the field of gravity exploration, and particularly relates to a method for acquiring gravity anomaly blank area data based on potential field characteristics, which comprises the following steps: acquiring gravity and topographic data of an exploration area and a neighborhood, and generating standard grid data and identifying a blank area through preprocessing and gridding; multiple interpolation algorithms are combined for parallel calculation, cross validation quantitative evaluation is carried out, and an optimal model is selected to obtain a reliable Bouguer gravity abnormal value of the blank area; mathematical analysis and filtering smoothing processing are carried out on the interpolation result, noise is eliminated, and a Bouguer gravity abnormal value conforming to the physical continuity and smoothness of the gravity field is obtained; and superposing the processed Bouguer gravity anomaly with a terrain gravity effect to obtain a complete space gravity anomaly field. Through systematic processing of existing Bouguer gravity anomaly and topographic data, filling and balancing of gravity anomaly values in areas with complex terrains and sparse data are successfully realized by using technical means of interpolation, analysis, smoothing, reduction and the like.

A method for 3D inversion of Full Tensor Gradiometry (FTG) data comprising receiving observed FTG data, performing a kernel matrix calculation on a subset of the observed FTG data and the gravity anomaly data, performing one-time forward modeling in a wavenumber domain to produce predicted FTG data and predicted gravity anomaly data for the reference density model, performing a residual between the observed data and predicted data, performing a depth-weighting function, a model covariance matrix, and data error covariance matrix on the observed FTG data and observed gravity anomaly data, obtaining a model update based on the depth-weighting function, kernel matrix, the model covariance matrix and the data error covariance matrix, and the residual between the observed FTG data and the predicted FTG data in the wavenumber domain, and performing inversion by directly obtaining an inverted model based on the model update and reference model.

The present invention relates to the field of marine surveying and mapping technology, and more specifically, to an underwater navigation adaptability improvement method based on curvature direction feature adaptation area aggregation, which is beneficial for eliminating isolated adaptation areas and disconnected adaptation blocks and filling internal voids of adaptation blocks, thereby improving directional adaptability and regional connectivity. The method comprises the following steps: obtaining a gravity anomaly reference map: gridding gravity anomaly data and dividing it using a sliding window method to obtain a gravity anomaly reference map; designing curvature direction features: combining the directional characteristics of maximum and minimum curvatures and the characteristics of Gaussian curvature for evaluating regional gravity anomaly changes to construct curvature direction features; and pre-calibrating adaptation area clustering and implementing an adaptation area reclosing mechanism to obtain a final calibrated adaptation area. Compared with the prior art, the selected area not only has stronger connectivity but also has good directional adaptability, thereby ensuring the positioning accuracy of an underwater vehicle in multiple headings.

This invention discloses a high-precision density reversal method and system for gravity anomalies based on regular, stable, and high-speed solution. [Solution] Based on the above-mentioned field source regularization solution, we propose using data with false anomalies at the field source imaging depth as boundary value conditions. Again, by solving the upper space of the first class boundary value problem of Laplace's equation, we calculate gravity anomalies on the ground and perform least-squares optimization iterative solution with the actual ground anomalies to obtain high-precision anomaly imaging data at the field source depth that does not contain false anomalies due to the Gibbs effect. This depth is used as the top interface of the field source, and the top interface is divided into grid prisms to construct the field source morphology. Density inversion is performed on each prism top surface in the frequency domain to obtain the density distribution of the field source morphology. Since this algorithm is implemented in the frequency domain and the inversion is performed at a known field source depth, the present invention is a gravity density inversion method with high speed, high precision, and high resolution characteristics, and provides a new concept for realizing the solution of gravity field source density parameters.

The invention provides an island small landslide early warning method based on gravity sensing and a related device. The method comprises the following steps: respectively collecting gravity data and environment interference data of a target island in real time by using a gravity sensing sensor array and an environment sensing sensor; performing interference correction on the gravity data based on the environment interference data to obtain corrected gravity data; mapping the corrected gravity data into a three-dimensional coordinate system of the target island, and generating a gravity field space distribution matrix; based on the gravitational field spatial distribution matrix, marking an area with gravity anomaly in the target island as a potential anomaly area; and performing prediction based on gravity data corresponding to the potential abnormal area in the corrected gravity data by using the trained landslide early warning model, and performing small landslide early warning on the target island when determining that the predicted landslide occurrence probability is greater than a preset early warning threshold. According to the invention, high-sensitivity and extreme-weather-resistant real-time monitoring and accurate early warning of small landslides of islands and reefs can be realized.

This invention discloses a method, apparatus, and equipment for determining formation thickness based on combined gravity and seismic data. The method includes: extracting the vertical derivative of the Bouguer gravity anomaly in the target basin; regularizing the vertical derivative of the Bouguer gravity anomaly to obtain a regularized value; determining the formation thickness at several well points and in three-dimensional regions within the target basin based on seismic and well data; and establishing a formation thickness prediction function for the target basin based on the regularized value of the vertical derivative of the Bouguer gravity anomaly and the formation thicknesses at several well points and in three-dimensional regions, thereby determining the formation thickness of the entire target basin based on the formation thickness prediction function. This method not only reduces the uncertainty of difference caused by data sparsity when using only seismic data but also solves the drawback of existing gravity interface inversion techniques lacking seismic constraints.

The invention provides a microgravity identification method and system for an underground space. The method belongs to the technical field of gravity recognition, and comprises the following steps: carrying out forward modeling calculation on each group of geometric parameter combination by utilizing a gravity anomaly analysis formula, and generating a synthetic sample data set formed by a gravity anomaly sequence and corresponding geometric parameters; carrying out tidal correction and null drift correction on the actually measured microgravity original data to obtain a basic gravity abnormal value, carrying out terrain correction, height correction and interlayer correction on the basic gravity abnormal value, and calculating Bouguer gravity anomaly after correction; after determining the influence range of the superficial interference, stripping the superficial interference field to obtain target gravity abnormal data for inversion; training a deep learning inversion model by using the synthetic sample data set; and inputting target gravity anomaly data used for inversion into the trained deep learning inversion model, and outputting a prediction result of the geometric parameters and the spatial position of the underground target body. The microgravity recognition precision can be improved in the urban shallow complex environment.

The application belongs to the field of gravity anomaly detection, and in particular to a gravity anomaly target body boundary estimation method and system. Based on the comparison and calculation of boundary functions and forward data sets, the characteristic line segment length vector is obtained from the extreme points of the boundary functions, a reference model is set in the observation range, the characteristic line segment length of each reference model is calculated, the boundary length and the characteristic line segment vector forward data set are established, the target to be measured is matched and compared with the data set, and the boundary estimation of the target body to be measured in the medium field gravity environment is realized. The problems of data attenuation, divergence and low resolution in gravity detection in the medium field environment are solved.

The invention discloses a gravitational field adaptation degree calculation method based on course optimization. The gravitational field adaptation degree calculation method comprises the following steps: converting a gravity anomaly database into a gravity anomaly grayscale image; an improved Sobel operator template is adopted to carry out convolution on the gravity anomaly grayscale image; convolution factors in different directions in the improved Sobel operator template are determined according to the course; constructing a gray level co-occurrence matrix of each sample matching area according to a convolution result; calculating characteristic parameters according to the gray level co-occurrence matrix, and carrying out normalization processing on the characteristic parameters; and calculating the characteristic parameters of different sample matching areas through a principal component analysis method to obtain the fitness of each sample matching area. According to the method, the influence of the course is considered when the gravitational field adaptation degree is calculated, and the precision and accuracy of gravity matching navigation are improved.

This invention relates to the field of geophysical exploration, specifically disclosing a method for identifying geological body boundaries based on adaptive matching of multi-source gravity and magnetic data. The method includes: acquiring gravity and magnetic anomaly data; estimating the equivalent source depth; constructing a depth-response mismatch correction field based on the depth attenuation difference between gravity and magnetic potential fields; compensating for and registering the magnetic anomaly with remanent magnetization and demagnetization offset; constructing an adaptive coupling direction tensor based on the consistency of the local gradient structures of the two fields; performing anisotropic boundary enhancement on the two fields using the principal direction of the adaptive coupling direction tensor as a constraint and the coupling strength as a weight; and iteratively fusing the two fields through mutual verification using the adaptive coupling direction tensor as a constraint and the depth-response mismatch correction field as a weight, outputting the fine spatial location of the geological body boundary and an adaptive confidence evaluation result. This invention achieves adaptive registration and direction constraint enhancement in the depth domain of gravity and magnetic data, effectively improving the accuracy and reliability of geological body boundary identification under complex geological conditions.

The invention discloses a short critical strong earthquake occurrence time prediction method and device, a medium and equipment, and relates to the technical field of earthquake prediction. The method comprises the following steps: performing standardized preprocessing on an atmospheric tidal gravity signal in a preset time window, and after an atmospheric tidal gravity abnormal event is detected, inputting a plurality of target days into a pre-constructed prediction model respectively to obtain the probability of strong shock in each target day after the atmospheric tidal gravity abnormal event occurs; the prediction model is obtained by training according to statistical correlation characteristics between the occurrence time of the atmospheric tidal gravity abnormal event and the occurrence time of the strong earthquake; the strong earthquakes represent earthquakes of which the magnitude is greater than a preset threshold value; and determining a strong earthquake occurrence time prediction result according to the probability of occurrence of the strong earthquake in each target day after the abnormality occurs. According to the method, quantitative probability prediction of the strong earthquake occurrence time is realized, and the method has the advantages of interpretability, calibration and online updating.

The application belongs to the technical field of geophysical inversion, and relates to a gravity anomaly three-dimensional inversion method based on multi-scale feature fusion, comprising the following steps: randomly generating density anomaly bodies with different number, density, size and shape to obtain a three-dimensional density model of underground medium; forward modeling the three-dimensional density model of underground medium to obtain two-dimensional gravity anomaly data of the ground surface; training a multi-scale feature fusion model with the three-dimensional density model of underground medium as output and the two-dimensional gravity anomaly data as input, wherein the multi-scale feature fusion model comprises a two-dimensional encoder and a three-dimensional decoder, and a dimension converter and a multi-scale feature fusion device located between the two-dimensional encoder and the three-dimensional decoder; inputting the measured two-dimensional gravity anomaly data into the trained multi-scale feature fusion model to obtain a three-dimensional density body constrained by a physical range through inversion. The application enhances the expression capability of multi-scale features of underground structures, and can relieve the multi-solution problem of inversion and the problem of resolution decline in deep parts.

The invention discloses an edge detection filter based on a TANH function and an arcsine function, and the filter determines Bouguer gravity anomaly data through obtaining free air gravity anomaly data, terrain correction values and Bouguer correction values of a plurality of gravity points of a research region. Calculating Gaussian envelope gradient data based on the gravity gradient tensor component data and Gaussian envelope, and obtaining an edge detection filtering result according to the Gaussian envelope gradient data and TANH function filtering; noise interference can be effectively suppressed by adopting a double-constraint mechanism, a generated boundary result is clear and continuous, and the problems of boundary divergence, positioning fuzziness, form distortion and the like existing in a traditional method are effectively solved; anomalous bodies with different depths are synchronously and accurately identified, stable detection performance can be maintained under complex geological conditions, and generation of false boundaries under the condition of coexistence of positive and negative anomalies is reduced; the reliable boundary recognition capability can be kept under the strong noise background, manual intervention in the preprocessing link is reduced, and the practical application convenience and the result reliability are improved.

The invention relates to the technical field of data processing, in particular to a mining area geological information comprehensive analysis method and system based on data integration, and the method comprises the steps: obtaining and preprocessing multi-source geological data, carrying out the feature correlation analysis and normalization, and obtaining a standardized data set; classifying the magnetic anomaly and gravity anomaly data through a support vector machine to generate a fluid dynamic classification result; performing interpolation processing on the classification result to generate continuous fluid dynamic data, and analyzing the speed and direction change of the fluid to generate a permeability distribution diagram; extracting an abnormal region according to the dynamic change trend, matching the abnormal region with historical data, marking geological abnormal points and generating an abnormal distribution diagram; according to the geological flow monitoring report, the underground water body type or the change of the oil-gas mixing proportion is judged, a pressure fluctuation mode is recognized, and a refined geological abnormal point distribution diagram is generated. According to the invention, mining area resource management and potential risk early warning can be effectively supported.

The invention discloses a method for multi-source data joint identification and inversion of a middle and lower earth crust high conductor, and relates to the technical field of geophysical exploration. The method for multi-source data joint identification and inversion of the middle and lower earth crust high conductor mainly comprises the following steps: carrying out potential field separation on Bouguer gravity anomaly data to obtain residual gravity anomaly, regional gravity anomaly and field source depth, and further determining earth crust residual gravity anomaly and middle and lower earth crust density disturbance characteristics; and calculating residual gravity low anomalies reflecting the middle and lower earth crust high conductors in the research area, delineating plane distribution of the middle and lower earth crust high conductors, and performing random sub-domain three-dimensional automatic inversion for the residual gravity low anomalies to obtain a three-dimensional space distribution form of the middle and lower earth crust high conductors. By implementing the multi-source data combined identification and inversion method for the mid-lower earth crust high conductor provided by the invention, the mid-lower earth crust high conductor can be rapidly and effectively identified and inverted, and the three-dimensional space distribution of the mid-lower earth crust high conductor can be obtained.

The invention discloses a mining area stress-strain monitoring method based on gravity abnormal change, and the method comprises the specific steps: laying an observation network, and obtaining the gravity abnormal data of a measurement point; preprocessing the obtained gravity anomaly data of the measurement point to obtain processed gravity anomaly data; constructing a linear observation equation set, and introducing a regularization item to solve an underground disturbance density vector; obtaining a mapping relation between the underground disturbance density and the stress tensor through a relation between the stress tensor and the strain tensor, and solving the underground disturbance density through a gravity anomaly difference value between different time points so as to obtain the stress tensors at different moments; and superposing the obtained stress tensor variation result to a mining area mining map, a geological structure map and a historical deformation map to form a comprehensive analysis map, and if the stress tensor corresponding to a certain working face or fracture structure area is obviously enhanced, prompting that the area has a goaf dilatation or ground pressure display risk. According to the method, the measured space coverage degree is high, and maintenance and remeasurement are facilitated.

The application discloses a kind of estimation method and system based on high frequency information gravity-magnetic anomaly lateral resolution, comprising the following steps: obtaining gravity-magnetic anomaly data, griding gravity-magnetic anomaly data, obtaining griding gravity-magnetic anomaly data;Low-pass filtering is carried out to the griding gravity-magnetic anomaly data, and regional gravity-magnetic anomaly data are obtained;Based on the griding gravity-magnetic anomaly data and regional gravity-magnetic anomaly data, subtract to obtain residual gravity-magnetic anomaly data;Change geometric scale factor, obtain residual gravity-magnetic anomaly data of different scales;The residual gravity-magnetic anomaly data of different scales are compared and analyzed, and the preset scale that presents residual gravity-magnetic anomaly distribution characteristics is obtained, and the estimation of gravity-magnetic anomaly lateral resolution is realized.The lateral resolution of gravity anomaly or magnetic anomaly estimated by the method can evaluate the related geological problems that can be solved by existing gravity-magnetic data, and can improve the efficiency of gravity-magnetic data processing and interpretation.

A method for 3D inversion of Full Tensor Gradiometry (FTG) data comprising receiving observed FTG data, performing a kernel matrix calculation on a subset of the observed FTG data and the gravity anomaly data, performing one-time forward modeling in a wavenumber domain to produce predicted FTG data and predicted gravity anomaly data for the reference density model, performing a residual between the observed data and predicted data, performing a depth-weighting function, a model covariance matrix, and data error covariance matrix on the observed FTG data and observed gravity anomaly data, obtaining a model update based on the depth-weighting function, kernel matrix, the model covariance matrix and the data error covariance matrix, and the residual between the observed FTG data and the predicted FTG data in the wavenumber domain, and performing inversion by directly obtaining an inverted model based on the model update and reference model.

The invention relates to the technical field of shipborne heavy force measurement, in particular to a shipborne heavy force abnormity measurement method, a medium and equipment. The method comprises the following steps: performing combined navigation fusion on a shipborne inertial navigation system and a satellite navigation system through Kalman filtering to obtain corrected angular velocity and specific force; calculating a gravity anomaly measurement coarse value and a virtual vertical accelerometer specific force; enabling the gravity anomaly measurement coarse value, the virtual vertical accelerometer specific force and the actual vertical accelerometer specific force to respectively pass through an FIR low-pass filter to obtain a gravity anomaly measurement fine value, a virtual vertical accelerometer specific force filtering value and an actual vertical accelerometer specific force filtering value; and correcting the gravity anomaly measurement precise value to obtain a final gravity anomaly measurement result. According to the ship-borne gravity anomaly measurement method, gravity anomaly measurement errors generated by actual vertical accelerometer offset errors can be restrained, and the gravity anomaly measurement accuracy is improved; and meanwhile, the precision requirement on a gravity sensitive instrument such as an accelerometer is reduced.

The present disclosure provides a gravity anomaly extraction method, device, equipment and storage medium, and belongs to the technical field of gravity exploration. The method comprises: performing terrain correction and curving processing on the Bouguer gravity anomaly of a target region to obtain a first gravity anomaly; performing constant-density gravity stripping processing on the first gravity anomaly to obtain a second gravity anomaly; performing variable-density gravity stripping processing on the second gravity anomaly to obtain a third gravity anomaly, wherein the density of the multiple layers of strata between the stratum where the target geological body is located and the surface is variable in the along-layer lateral direction; determining a first gravity regional field of the target region; and extracting the gravity anomaly of the target geological body according to the third gravity anomaly and the first gravity regional field. The variable-density gravity stripping processing on the second gravity anomaly can eliminate the influence of the variable density of the strata in the along-layer lateral direction on the gravity anomaly of the target geological body, thereby improving the accuracy of the extracted target gravity anomaly.

The present invention discloses a method and system for high-precision density inversion of gravity anomalies based on regular stable fast solution. [Solution] Based on the above source regularization solution, we propose using data containing spurious anomalies at the source imaging depth as boundary value conditions. Then, through the upper half-space solution of the first-class boundary value problem of Laplace's equation, we calculate the gravity anomaly on the ground, and perform a least-squares optimization iteration with the actual ground anomaly to obtain high-precision anomaly imaging data at the source depth without spurious anomalies caused by the Gibbs effect. This depth is used as the top interface of the source, and the top interface is divided into grid prisms to construct a source configuration. Density inversion is performed on each prism top surface in the frequency domain to obtain the density distribution of the source configuration. Because this algorithm is implemented in the frequency domain and the inversion is performed at a known source depth, the present invention provides a gravity density inversion method with high speed, high precision, and high resolution, and provides a new concept for solving gravity source density parameters.

The invention discloses a method and device for calculating a full-crustal extension coefficient based on subcrustal flow and magma invasion recovery, and the method comprises the steps: obtaining geological comprehensive data, such as a seabed seismic survey line and deep reflection multichannel seismic section data, of a target region, carrying out the stratigraphic division according to the geological comprehensive data, constructing a time-depth conversion formula through the drilling and logging data, and carrying out the calculation of the full-crustal extension coefficient. And converting the time profile data into depth domain data. And establishing an earth crust structure model in combination with a stratigraphic division result and related data, calculating a gravity anomaly curve through forward modeling and correcting the model until a preset condition is met so as to obtain density structure distribution of each part of the profile, and extracting the thickness of the earth crust high-density body. And then, according to the abnormal settlement after cracking, the earth crust loss amount is extracted, the current earth crust thickness is recalculated by removing the earth crust high-density body thickness and supplementing the earth crust loss thickness, the corrected earth crust thickness is obtained, and finally the whole earth crust extension coefficient is calculated. According to the method, the influence of bringing-in and taking-out of crustal substances is eliminated, and the accuracy of evaluating the crustal extension change rule is greatly improved.

The present application relates to a strapdown dynamic gravity measurement method and device for mixed surface and underwater operations. The method comprises: determining the current state of a gravimeter through an autonomous detection system during navigation, and constructing a surface or underwater strapdown dynamic gravity measurement system based on the determination result; performing gravity measurement using the strapdown dynamic gravity measurement system, and calculating raw gravity anomaly values ​​based on high-precision velocity, position, and attitude information; filtering using a low-pass filter to obtain effective gravity anomaly results; constructing a celestial accelerometer drift model for the current voyage based on a segment of initial celestial relative force and a segment of final celestial relative force; and calibrating the effective gravity anomaly results using the initial celestial relative force, the celestial accelerometer drift model, and the gravity value of the gravimeter's onshore reference point to obtain the final gravity measurement result. This method can realize strapdown dynamic gravity measurement for mixed surface / underwater operations.