30 results about "Kriging method" patented technology

The invention discloses a tunnel geological information modeling and dynamic updating method and system. The method comprises the following steps: firstly, collecting drilling data before construction and stratum information exposed by a tunnel face in a construction process, defining a lithology category set after standardization processing, and converting the lithology category set into an indication variable; empirical semi-variance is calculated based on the standardized data set, and variation parameters are obtained through variation function fitting; determining a neighborhood radius according to the variation parameter, and constructing an initial geologic model by adopting an indicating Kriging method; in the tunneling process, newly revealed stratum information is fused in an incremental mode, variation parameters are adjusted in a self-adaptive mode, a geologic model is updated, and meanwhile uncertainty quantification is conducted; and finally, determining the most probable lithology category of each spatial position based on the updated geologic model, and generating a lithology distribution diagram. According to the method, lithology category probability prediction and geologic model dynamic updating are realized, tunnel geologic spatial features can be accurately described, prediction uncertainty is quantified, and prediction precision and reliability are improved.

The invention relates to the technical field of offshore wind power equipment, and discloses an offshore floating type wind turbine mooring system optimization method based on a proxy model, and the method comprises the following steps: S1, determining key design parameters; s2, sampling to generate design points; s3, establishing a total index function of multi-target coupling; s4, generating a sample database through numerical simulation; s5, constructing an agent model by using a Kriging method; s6, calling the proxy model to search an optimal solution by adopting a differential evolution algorithm; s7, performing numerical simulation verification on the optimal solution and judging convergence; and S8, if not, updating the database and returning to S5 for iteration, and if so, outputting the optimal configuration. According to the method, the problems that a traditional optimization method is high in simulation calculation cost and too long in consumed time are solved, searching is carried out through the proxy model instead of high-precision simulation, the verification iteration step is combined, the optimization efficiency is greatly improved while the result accuracy is guaranteed, and the multi-target mooring scheme with the optimal comprehensive performance can be efficiently obtained.

The invention discloses a coal mine geological data three-dimensional modeling and transparent working face construction method and system. The method comprises the steps that multi-source geological and production data are collected and preprocessed; constructing a refined three-dimensional geologic body model based on multipoint geostatistics and a collaborative Kriging method; fusing the roadway model through Boolean operation; dynamically constructing a working face model by using real-time mining data; when a model error exceeds a threshold value, dynamic updating based on data assimilation is triggered; and finally, realizing model integration and transparent display in a three-dimensional visualization engine. According to the method, multi-point geostatistics and collaborative Kriging interpolation are fused, multi-source data such as geological drilling and three-dimensional earthquakes are integrated, the depicting precision of complex geological structures such as faults and collapse columns is greatly improved, the geometrical shape of the constructed three-dimensional geologic body model is more vivid, rock mechanics parameters are integrated, and the method is suitable for large-scale popularization and application. The problem that an existing three-dimensional geologic model is limited in precision is solved.

The invention relates to the technical field of geographic information systems and disaster monitoring, in particular to a multi-disaster historical data space-time interpolation method, which comprises a space-time collaborative interpolation module, a disaster-inducing factor identification module and a quantitative characterization module. Wherein the space-time collaborative interpolation module constructs a space-time variation function based on an optimized Kriging method and spline interpolation, and realizes collaborative interpolation of space and time dimensions; the disaster-inducing factor identification module is used for carrying out system identification and weight distribution on multiple disaster-inducing factors based on a variable fuzzy optimization theory and an entropy evaluation method; and the quantitative characterization module realizes dynamic quantization and hierarchical characterization of disaster-inducing factors based on a natural breakpoint grading method and an AHP weighting method. According to the method, high-precision space-time reconstruction and disaster-inducing factor dynamic quantification of historical disaster data in a multi-disaster complex terrain region are realized, the problems of insufficient precision and insufficient dynamic factor capture of a traditional interpolation method in a space-time heterogeneity region are effectively solved, and a reliable data basis is provided for multi-disaster risk assessment and disaster chain analysis.

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 belongs to the technical field of machine tool control, and particularly discloses a machine tool hydrostatic guideway feeding system motion error compensation method and a storage medium. The method comprises the following steps: measuring and acquiring Z-direction error data of the guide rail; processing the data by Fourier series and fitting a guide rail contour error; inputting load, oil pad geometric parameters and guide rail working parameters; a workbench-guide rail balance equation based on consideration of guide rail contour errors is established, and a Z-direction motion error model perpendicular to the feeding direction of a workbench is calculated to serve as a low-fidelity model; the Z-direction deviation ez'and the coordinate Y of the workbench are measured and obtained to serve as a high-fidelity model; constructing a multi-fidelity model by combining the low-fidelity model and the high-fidelity model based on a CO-Kriging method; and inputting the multi-fidelity model into a numerical control system, reading the position of the workbench in real time, substituting the position into the multi-fidelity model, and calculating and compensating the Z-direction motion error of the workbench. The method has the advantages of being easy to operate, high in compensation precision and good in reliability and flexibility.

The invention relates to a control method for a road construction wind shielding device. Collected wind condition information is subjected to filtering preprocessing; constructing wind field profile information continuously spatially distributed in the vertical direction through a common Kriging method, and outputting a wind field feature vector; calculating a vortex-induced resonance risk coefficient according to the wind field feature vector, judging whether the adjustment action is good or not through a reward function, finally training and optimizing by adopting a PPO-Clip algorithm, and outputting final action control parameters of the left wing plate, the right wing plate and the top plate; the movable wind shield is adopted, real-time changes of external airflow can be adapted, and a continuous and stable working environment is provided for a rear construction area; the system no longer responds to "arrived" wind, but predicts the morphology of "upcoming" wind. And the wing plate is adjusted to a pressure relief angle in advance, so that sudden impact is avoided.

The invention provides a joint optimization method for geometric parameters of an outer rotor brushless direct current motor. According to the joint optimization method, average torque maximization, torque pulsation minimization and cogging torque minimization are taken as optimization targets; the method comprises the following steps: firstly, constructing a motor performance meta-model based on AMOP sampling and a moving least square-isotropic Kriging method, and establishing comprehensive sensitivity through Pearson correlation to distinguish primary sensitive and secondary sensitive variables; and then, a q-ROFS-TOPSIS method is introduced, so that a unified closeness degree evaluation index can be obtained. The optimization process adopts a two-layer staged optimization strategy: in the first stage, sub-sensitive variables are dominated, main sensitive is allowed to perform fine tuning, the relevance between main and sub parameters is emphasized, and optimization is performed through a particle swarm optimization algorithm; and in the second stage, the main sensitive variable is adjusted on the basis of the sub-sensitive optimization, and optimization is carried out through an evolutionary algorithm. According to the method, the optimal motor geometric parameter combination with high average torque, low torque ripple and low cogging torque can be obtained under the conditions of incomplete samples and multiple working conditions.

The invention discloses a fuel gear pump behavior prediction method based on mechanism-data fusion. The method comprises the following steps: 1, establishing a fuel gear pump internal flow mechanism model; the data acquisition module is used for acquiring low-confidence data based on physical constraints; 2, acquiring multi-working-condition experimental data of the fuel gear pump, and using the experimental data as high-confidence data for training a fuel gear pump mechanism-data fusion prediction model; and step 3, establishing a fuel gear pump mixed behavior prediction model based on a collaborative Kriging method, and utilizing the auxiliary training data set obtained based on the mechanism model and the high-confidence data set obtained based on the multi-working-condition test bed to combine with the collaborative Kriging method to predict the mixed behavior of the fuel gear pump. And the physical boundary constraint of the mechanism model and the strong nonlinear fitting capability of data driving are fused, and a fuel gear pump mechanism-data fusion prediction model is trained, so that high-precision behavior prediction of the fuel gear pump is realized. According to the invention, the prediction precision and generalization ability are significantly improved.

This application discloses a lithofacies analysis system and method based on a shale geological-mechanical coupled lithofacies classification system, relating to the field of shale lithofacies analysis. The system construction module determines the geological-mechanical coupled lithofacies classification system based on information data from the data acquisition module; the prediction module determines single-well lithofacies prediction information; the feedback adjustment module adjusts the single-well lithofacies prediction information based on set rock mechanics and geostress parameter ranges to obtain adjusted information; the three-dimensional model construction module uses a genetic algorithm to perform seismic attribute embedding and integrated inversion based on seismic attribute data and adjusted information, and uses a co-kriging method to construct a three-dimensional model of well-seismic co-modeling of geological-mechanical coupled lithofacies; the parameter model determination module uses a layer-control-facies dual-control method to control constraints and determine the heterogeneity of the geological-mechanical characteristics of shale layers in three-dimensional space, providing a basis for establishing production areas for the exploration and development of artificial oil and gas reservoirs.

The invention discloses a method, a device and equipment for establishing a dynamic permeability model and a readable storage medium. The method comprises the following steps: selecting a target production well meeting a preset condition from each production well; calculating the dynamic permeability of each perforation section in the target production well; calculating a permeability dimensionless ratio based on the dynamic permeability of each perforation section; based on a preset logging curve and the permeability dimensionless ratio, obtaining a permeability ratio model through a Kriging method; and multiplying the static permeability model by the permeability ratio model to obtain a dynamic permeability model. According to the invention, the technical problem of low accuracy of the correction result of the permeability in the prior art is solved.

The application relates to the technical field of data processing, and discloses a bridge displacement remote intelligent monitoring method and system and computer equipment. The method comprises the following steps: setting multiple trigger thresholds to perform grouped polling sampling to obtain displacement data; adopting dynamic weight fusion to remove noise, identify abnormalities, and generate a digital twin model by using a Kriging method; extracting features to classify and distinguish environmental and structural response displacements; calculating a deviation index to estimate fatigue damage formation and generate a health report; calculating a risk entropy value to set an early warning and generate a treatment scheme; and evaluating an improvement rate to update and configure a closed-loop optimization. By fusing environmental data and historical monitoring data, the application establishes an adaptive sampling strategy and an intelligent analysis and processing mechanism, accurately distinguishes displacement changes caused by environmental factors and structural factors, and performs accurate risk quantification and automatic decision support.

The invention relates to a reliability analysis method for a missile wing unfolding mechanism, and belongs to the field of design and reliability analysis of missile structures. Comprising the following steps: step 1, establishing a main failure mode of the missile wing unfolding mechanism according to a failure condition and a failure mechanism of the missile wing unfolding mechanism; and 2, determining the distribution type and distribution parameters of the input variables, and establishing a reliability analysis model of the missile wing unfolding mechanism. And 3, solving the proxy model of the missile wing unfolding mechanism by adopting a self-adaptive Kriging method. 4, reliability analysis is conducted on the missile wing unfolding mechanism through the agent model and the Monte Carlo method. The method provided by the invention has strong operability, can control the uncertainty factors causing the reliability problem in the design link, and has a good application prospect.

The invention discloses a simulation prediction method based on reduced-order model machine learning, and relates to the technical field of crossing of computer-aided engineering simulation and machine learning, in particular to a simulation prediction method based on reduced-order model machine learning, which comprises the following steps: firstly, performing finite element simulation on automobile parts (such as an anti-collision beam) to obtain simulation data; secondly, selecting sample points by adopting a Latin hypercube sampling method, and verifying that a data set meets an order reduction requirement; then, carrying out dimension reduction processing on the data by applying an order reduction model (such as intrinsic orthogonal decomposition (POD) or singular value decomposition (SVD)); and finally, predicting the data after order reduction by using a Kriging method machine learning model, and evaluating the prediction precision by using a decision coefficient R2 and a root mean square error RMSE. According to the method, through combination of the reduced-order model and machine learning, the simulation calculation efficiency is effectively improved, meanwhile, the accuracy of a prediction result is ensured, and the method is suitable for rapid simulation optimization in automobile design.

The invention provides a river topography construction method and system fusing multiple nonlinear trends, and relates to the field of hydraulic engineering digitization. The curve and width change characteristics of a complex river are described through a dynamic coordinate conversion mode based on a normal vector field, a multiple nonlinear trend function parameter space is constructed, an optimal function combination is automatically optimized to generate a trend model, and spatial topology constraints are utilized to automatically adapt to the characteristics of different terrain sections of any river channel. The defects that the boundary of the traditional Kriging interpolation method is not limited, the terrain uniformity assumption is insufficient, and the system error is high when the section spacing is large are overcome. According to the method, normal vector field coordinate conversion is utilized to automatically adapt to river boundaries, multiple nonlinear trend function combinations describe comprehensive features of the river section, a Kriging method is combined to generate a high-precision interpolation result, and the method is suitable for sparse section-based river terrain reconstruction.

The application discloses a tunnel geological information modeling and dynamic updating method and system. The method first collects drilling data before construction and stratum information exposed in a working face during construction, defines a lithology category set and converts the lithology category set into an indicator variable after standardization processing; calculates an empirical semi-variance based on the standardized data set, obtains a variation parameter through variogram fitting; determines a neighborhood radius according to the variation parameter, and constructs an initial geological model by using an indicator Kriging method; in the tunneling process, newly exposed stratum information is incrementally integrated, the variation parameter is adaptively adjusted, the geological model is updated, and uncertainty quantification is simultaneously performed; finally, the most possible lithology category of each spatial position is determined based on the updated geological model, and a lithology distribution map is generated. The application realizes lithology category probability prediction and geological model dynamic updating, can accurately depict tunnel geological spatial characteristics, quantitatively predict uncertainty, and improves prediction accuracy and reliability.

The application provides a method and system for obtaining offshore wind farm site parameters, the method comprising: planning the distribution and number of overall wind farm fan positions, and selecting representative sampling points by using an optimization method; sampling static pressure at the representative sampling points to obtain low-disturbance in-situ samples; performing static sounding tests on all planned fan positions to obtain static sounding data of all positions; testing the low-disturbance in-situ samples to obtain soil strength interpretation parameters of the representative sampling points; and using the Kriging method to expand the interpretation parameters of a limited number of representative sampling points to all fan positions to interpret the static sounding data. The application can effectively improve the accuracy of site parameter prediction and evaluation, reduce the cost and workload of offshore survey, and reduce potential risks and uncertainties in the process of offshore wind power engineering design and construction.

A method of synthesizing reactor core power distribution is disclosed. According to the method, the power of all fuel assemblies in a reactor core is calculated from powers of fuel assemblies in which in-core instruments are located using the ordinary kriging methodology. Also, a hot-pin power distribution of each fuel assembly is synthesized from the power of all fuel assemblies calculated, whereby more accurate hot-pin axial power distribution, rather than pseudo hot-pin axial power distribution, may be synthesized.

The application discloses a kind of near-surface stratum Q value modeling methods for micro-logging data in loess plateau region, which includes: step one: using improved micro-logging first arrival spectrum ratio method to obtain the near-surface micro-logging Q value of logging point;Step two: using Kriging method to carry out spatial interpolation on micro-logging Q value data, obtain low-frequency near-surface Q value model;Step three: using machine learning to predict the rest of the spatial position Q value through known logging Q value data, obtain high-frequency near-surface Q value model;Step four: scale fusion is carried out to low-frequency, high-frequency near-surface Q value model, and a fine near-surface Q value model is established;Step five: using fine near-surface Q value model to carry out absorption and attenuation compensation processing on seismic data.The application solves the problem that the extracted logging Q value is not accurate due to the poor coupling relationship of shotpoint and receiver, provides reliable key Q value data for solving the absorption and attenuation compensation problem in complex near-surface area, and can be used in industrial production of seismic exploration.

The application provides a control surface structure flutter margin optimization method based on a Kriging surrogate model, comprising the following steps: establishing a control surface structure reference model; performing quality discrete processing on the reference model, updating the quality of each point of the control surface structure reference model according to the perturbation method, performing control surface structure flutter analysis, performing sensitivity analysis on each quality point of the control surface, obtaining the positions of a plurality of quality points with the greatest influence on the quality point discrete flutter speed, and marking the plurality of quality points as flutter significant points; updating the quality data of the plurality of flutter significant points in batches, establishing a flutter speed surrogate model based on the Kriging method; testing the flutter speed surrogate model by using a test data set; and performing rudder surface flutter optimization analysis based on the control surface flutter speed surrogate model and a genetic algorithm. The technical scheme of the application is used to solve the technical problem that the traditional design method usually depends on experience, experiments or a large amount of calculation, and the design cycle is long and the cost is high.

This invention discloses a high-precision 3D geological modeling method for chamber hydrogen storage using multi-source data collaborative calibration, belonging to the field of chamber hydrogen storage and 3D modeling technology. Addressing the shortcomings of existing methods, such as poor multi-source data fusion, lack of dynamic calibration mechanisms, and insufficient grid adaptability, this invention comprises the following steps: First, data from seismic, borehole, outcrop, and 3D laser scanning sources are collected and preprocessed to construct a standardized database. Then, a reliability evaluation system is established, and contradictory data are iteratively calibrated using an improved Kriging method. Next, key geological areas are identified, and an adaptive hybrid grid partitioning strategy is adopted. Subsequently, an initial model is constructed and iteratively optimized through cross-validation and borehole matching verification. Finally, the model is output, and a dynamic update mechanism is established. This invention effectively improves data reliability and modeling accuracy, reduces computational costs, and provides more accurate 3D geological modeling results for related fields of chamber hydrogen storage engineering.

The invention relates to the technical field of mineral resource quantity calculation, in particular to a polymetallic mineral resource quantity estimation method and system, and the method comprises the steps: carrying out the unified constraint of the spatial correlation of polymetallic grades in adjacent calculation units through a cooperative Kriging method, enabling the consistent relation of different metals in the trend and tendency change direction to be brought into a same judgment system, and carrying out the estimation of the polymetallic mineral resource quantity. The space splitting risk generated when multiple metals are calculated separately is reduced, contribution sorting is conducted on multiple metal grades in the same calculation unit through the Spearman grade correlation coefficient, the relative change relation between the metal grades is converted into a comparable grade sequence, it is avoided that low-contribution metals amplify unstable factors in numerical value conversion, and the calculation efficiency is improved. Mutational metal items are identified and eliminated through spatial correlation constraints, interference of local abnormal grades on an overall result is reduced, and the formed multi-metal resource quantity result is enhanced in the aspects of spatial consistency, numerical comparability and overall reliability.

The invention relates to the technical field of mining, in particular to a resource evaluation method and system for graphite mine mining, and the method comprises the steps: building a geological resource database of a graphite mine based on dynamically updated drilling data; the method comprises the following steps: determining the boundary grade of an ore body in DIMINE software by importing drilling data from the earth surface to an interval of-600m, connecting the ore body and modifying an ore body interpretation line based on an exploration line section, and forming a three-dimensional entity model of a graphite mine; building a block model; performing value estimation on the block model according to a Kriging value estimation method; determining the resource quantity loss and dilution rate and the mining grade of the graphite mine in an interval from the earth surface to-600m based on DIMINE software, and converting the resource quantity of the graphite mine meeting a preset conversion condition into reserves; and generating and outputting a table according to the resource quantity, the resource quantity loss and dilution rate, the mining grade and the reserves of the graphite mine. The production cost and mining dilution loss of the graphite mine can be reduced, and resource waste and environmental pollution are reduced.

The invention discloses a calculation method of an airplane under the action of a steady-state non-uniform temperature field. Effective temperature measuring point data in a flight test is screened according to aircraft structural components to serve as field data; an initial non-uniform temperature field is formed through a heat conduction correction Kriging method interpolation method according to the spatial position; an airplane finite element model is adopted, the initial non-uniform temperature field is used as a boundary condition, heat transfer analysis calculation of the airplane in the heat convection environment is carried out, and a steady-state non-uniform temperature field is obtained; in addition to the field data, temperature measuring points used for verifying the steady-state non-uniform temperature field are arranged, temperature data are obtained to serve as verification data, correction of the steady-state non-uniform temperature field is completed, and verification data measuring points are preferentially arranged on the typical tangent plane; and carrying out material / geometric nonlinear simulation calculation on the airplane finite element model loading temperature field and the severe working condition aerodynamic load. According to the method, the calculation precision of the high-stress risk position can be improved, the calculation reliability is improved, and the safety risk is reduced.

The invention discloses a soil thickness prediction method and system, and relates to the technical field of transmission tower construction. The method comprises the following steps: acquiring digital elevation model data, multispectral remote sensing image data and geological lithology data, and respectively calculating a terrain humidity index, a lithology weighting coefficient, a terrain roughness index, a slope index and a normalized vegetation index based on the data. According to the method, a digital elevation model, a multispectral remote sensing image and geological lithology data are systematically integrated, a multi-factor evaluation system comprising a terrain humidity index, a lithology weighting coefficient, a terrain roughness index, a slope index and a normalized vegetation index is constructed, and the weight of each factor is scientifically determined by adopting principal component analysis or multiple linear regression; the accumulation potential index comprehensively reflecting the soil is generated through weighted fusion, and the inherent limitation that a traditional Kriging method only depends on the spatial distance is overcome.

The application discloses a spatial fine interpolation method for sea ice density based on a collaborative Kriging method, and comprises the following steps: S1, dividing an interpolation area and obtaining sea ice density observation values and climate parameter observation values at observation points; S2, taking the sea ice density of the observation points as a main variable and taking the climate parameter with the highest correlation with the sea ice density as a covariant; S3, calculating the distance between all the observation points, the main variable variation function value and the covariant variation function value; calculating the covariant variation function value between the observation points; S4, obtaining the autocorrelation distance of the covariant variation function value based on a fitting curve; S5, calculating the distance between the observation points and the interpolation point within the autocorrelation distance range of each interpolation point and obtaining the covariant variation function value between the observation points and the interpolation point; S6, constructing a Lagrange equation and solving the Lagrange equation to obtain a weight coefficient; and S7, constructing a quality control equation of the collaborative Kriging method and calculating a sea ice density prediction value of the interpolation point.

The application discloses a kind of nuclear fuel rod wear depth efficient high-precision calculation method, device, equipment and storage medium, it is related to computer technical field, the method includes: obtaining initial sample data set containing vibration parameter and corresponding wear depth;Based on initial sample data set, initial Kriging method Kriging prediction model is constructed;The predicted response value of sample point is calculated using initial prediction model, and the learning function value of sample point is calculated according to predicted response value and preset failure threshold value;Judge whether to meet iteration termination condition;If not, determine the least uncertain point according to learning function value, supplement simulation wear depth data to update training set sample or introduce new experimental wear depth data to update training set sample;According to the training set sample after updating, a new Kriging prediction model is constructed;If yes, output current prediction model as target prediction model.The method can improve the accuracy of nuclear fuel rod positioning grid wear prediction.

The application relates to the technical field of intelligent agriculture and intelligent regulation and control of agricultural environment, and discloses an intelligent agriculture environment intelligent regulation and control method based on space-time big data modeling, which comprises the following steps: collecting multi-dimensional space-time data of a target planting area and performing pretreatment to obtain a training data set; based on a gradient-enhanced Taylor Kriging method, sample point observation values and environmental gradient information in the training data set are used to construct a proxy model; real-time collection of soil, meteorological and crop online data of the target planting area is performed, environmental disturbance events are identified from the data, and the analytical gradient output by the proxy model is used to quantize the estimated influence of environmental disturbance on crop attention indexes; based on the estimated influence and / or multi-target preferences preset by a user, yield or quality prediction values output by the proxy model are called to generate an optimal environmental regulation strategy vector; the optimal environmental regulation strategy vector is executed and the proxy model is updated incrementally online, so that precise regulation and control of agricultural environment parameters are realized.

Provided are a lithofacies analysis system and method based on a shale geological-mechanical coupled lithofacies classification system, which relate to the field of shale lithofacies analysis. The lithofacies analysis system includes a system construction module that determines a geological-mechanical coupled lithofacies classification system based on information data of a data acquisition module. A prediction module determines single-well lithofacies prediction information. A feedback adjustment module performs feedback adjustment on the single-well lithofacies prediction information based on a set range of a rock mechanics and geostress parameter and acquire adjustment information. A three-dimensional model construction module performs, by a genetic algorithm, seismic attribute embedding and integrated inversion based on seismic attribute data and the adjustment information, and constructs, by a collaborative Kriging method, a three-dimensional model for well-seismic collaborative modeling of geological-mechanical coupled lithofacies.

The present application relates to a kind of gravity background field construction method based on iterative optimization kriging, comprising the following steps: S1. obtain local gravity measurement data as known point set;S2. based on the known point set, preliminary interpolation is carried out to the grid point of grid background map using kriging method, and the kriging estimated value of each grid point is obtained;S3. by introducing local gravity anomaly covariance and kriging variation function, the random uncertainty of kriging estimated value in interpolation process is evaluated, and each grid point is divided into reference point and to be modified point based on random uncertainty;S4. reference point is combined with known point set to construct continuous distribution constraint surface, and the continuous distribution constraint surface is locally continuously modified based on the obtained reference point, and the continuous distribution constraint surface is iteratively extended based on the locally modified constraint surface until the continuous distribution constraint surface is continuously modified, to realize the construction of gravity background field.