CN113960690A - A method and device for calculating the influence of sea surface gravity data measurement accuracy on the results of seabed topography inversion - Google Patents

A method and device for calculating the influence of sea surface gravity data measurement accuracy on the results of seabed topography inversion Download PDF

Info

Publication number
CN113960690A
CN113960690A CN202111033017.3A CN202111033017A CN113960690A CN 113960690 A CN113960690 A CN 113960690A CN 202111033017 A CN202111033017 A CN 202111033017A CN 113960690 A CN113960690 A CN 113960690A
Authority
CN
China
Prior art keywords
sea
gravity
depth
field element
data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202111033017.3A
Other languages
Chinese (zh)
Other versions
CN113960690B (en
Inventor
范雕
李姗姗
赵东明
李新星
张金辉
范昊鹏
冯进凯
单建晨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Information Engineering University Of Chinese People's Liberation Army Cyberspace Force
Original Assignee
PLA Information Engineering University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by PLA Information Engineering University filed Critical PLA Information Engineering University
Priority to CN202111033017.3A priority Critical patent/CN113960690B/en
Publication of CN113960690A publication Critical patent/CN113960690A/en
Application granted granted Critical
Publication of CN113960690B publication Critical patent/CN113960690B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention discloses a method and a device for calculating the influence of sea surface gravity data measurement precision on a submarine topography inversion result, wherein the method comprises the following steps: taking the sea depth as input data, resolving gravity field element information generated by the sea depth by using a strict prism volume division method, and obtaining an effective integral radius; adding sea depth Gaussian white noises corresponding to sea depths with different relative accuracies on the basis of the error-free sea depth, and further generating sea depth data with different relative accuracies; based on the obtained effective integral radius, respectively taking the error-free sea depth and the obtained sea depth data with different relative precisions as input, and forward calculating to obtain gravity field element results corresponding to the error-free sea depth and the sea depth data with different relative precisions; the forward gravity field element result of the sea depth data with different relative precisions is subtracted from the error-free forward gravity field element result of the sea depth to obtain a forward gravity field element difference value; and taking the difference value as input, and counting the difference result index. The invention discloses an internal correlation between gravity data and seafloor topography data.

Description

Method and device for calculating influence of sea surface gravity data measurement precision on submarine topography inversion result
Technical Field
The invention belongs to the technical field of submarine topography inversion, and particularly relates to a method and a device for calculating influence of sea surface gravity data measurement precision on a submarine topography inversion result.
Background
The measurement of the submarine topography (sea depth) is used as a basic means and a system engineering for observing the sea and recognizing the sea, and can play an irreplaceable role in the aspects of sea resource development, sea ecological environment protection, sea technological innovation, sea equity maintenance and the like. Currently, the submarine topography measuring technology mainly includes a ship-based submarine topography measuring technology, a submarine topography measuring technology, an Airborne laser radar depth measurement (ALB) technology, a satellite-based submarine topography measuring technology, and the like. The space grid size of the firstly issued grid is a 15 'global submarine topography model (GEBCO _2019), and the space grid size recovered by organizations such as SIO, NOAA, NGA and the like based on the SRTM15+ V1.0 and other terrain models is a 15' global topography model (SRTM15+ V2.0), wherein most submarine topography data of the ocean area are recovered by means of satellite altimetry gravity data inversion.
Generally, the accuracy of inverting a submarine topography model by using a satellite height measurement technology is to be improved. The inversion of the nature of the submarine topography by the multi-source sea surface gravity data (satellite altimetry gravity data, ship survey gravity data, aviation gravity data and the like) is the inverse solution of the problem of the forward disturbance field element of the submarine topography, so that the accuracy of constructing a submarine topography model by the gravity data is mainly influenced by two factors, namely the calculation model error of the submarine topography by gravity inversion; secondly, gravity input data error. The sea surface gravity input data quality will directly affect the sea surface gravity data based seafloor terrain modeling effect.
Disclosure of Invention
The invention provides a method and a device for calculating the influence of sea surface gravity data measurement precision on a sea surface topography inversion result, aiming at the problem that the sea surface gravity data inversion precision is influenced by sea surface input gravity field metadata quality by utilizing satellite height measurement gravity data.
In order to achieve the purpose, the invention adopts the following technical scheme:
the invention provides a method for calculating the influence of sea surface gravity data measurement precision on a submarine topography inversion result, which comprises the following steps:
step 1: taking the sea depth as input data, resolving gravity field element information generated by the sea depth by using a strict prism volume division method, and obtaining an effective integral radius;
step 2: adding sea depth Gaussian white noises corresponding to sea depths with different relative accuracies on the basis of the error-free sea depth, and further generating sea depth data with different relative accuracies;
and step 3: taking the effective integral radius obtained in the step 1 as a basis, respectively taking the error-free sea depth and the sea depth data with different relative precisions obtained in the step 2 as input, and forward-calculating to obtain gravity field element results corresponding to the error-free sea depth and the sea depth data with different relative precisions;
and 4, step 4: 3, the gravity field element result obtained by forward modeling of the sea depth data with different relative precisions in the step 3 is subtracted from the gravity field element result obtained by forward modeling of the sea depth without errors to obtain a forward gravity field element difference value;
and 5: and 4, taking the difference obtained in the step 4 as input, and counting the difference result indexes.
Further, the step 1 comprises:
resolving gravity field element information generated by sea depth by using a strict prism volume division method, wherein the gravity field elements comprise sea surface gravity anomaly and gravity anomaly vertical gradient:
Figure BDA0003245886070000021
in the formula (I), the compound is shown in the specification,
Figure BDA0003245886070000022
wherein G is the earth gravitational potential constant and takes the value of 6.672 multiplied by 10-8cm3/(g·s2)；ρcAnd ρwRespectively representing the density of the crust of the earth and the density of the seawater; (i)p,jp) To study the site location; (i, j) is the flow point position; h (i, j) is the corresponding sea depth of the (i, j) point; l is the distance between the point of interest and the flow point; (x, y, z) are flow point coordinates; sLAnd SBRespectively representing the points of the integral radius of the central area in the longitude and latitude directions; Δ x and Δ y represent grid sizes in the longitude and latitude directions, respectively; t iszFor disturbance gravity, the disturbance gravity is used to represent sea surface gravity anomaly; t iszzIs a gravity abnormal vertical gradient.
Further, the step 2 comprises:
sea depth data of different relative precisions are generated by formula (3):
herr＝h0+hno (3)
in the formula, herrDifferent relative precision depths; h is0The sea depth is error-free; h isnoThe sea depth Gaussian white noise standard deviation is equal to the product of the error-free sea depth average value and the relative precision.
Further, in step 4, the forward gravitational field component difference value is obtained according to the formula (4):
Figure BDA0003245886070000031
in the formula, gammaiGravity representing the ith forward performanceA field element difference value;
Figure BDA0003245886070000032
and
Figure BDA0003245886070000033
respectively representing the ith gravity field element values obtained by forward modeling of different relative precision sea depths and error-free sea depths; and N is the number of the statistical gravity difference values.
Further, in step 5, the difference result indicator includes:
Figure BDA0003245886070000034
in the formula, gammaiRepresenting the difference value of the ith forward acting gravity field element; delta gammamax、△γmin、△γmeanAnd Δ γstdThe maximum, minimum, mean and standard deviation of the difference are indicated, respectively.
The invention also provides a device for calculating the influence of the sea surface gravity data measurement precision on the inversion result of the submarine topography, which comprises:
the resolving module is used for resolving gravity field element information generated by the sea depth by using a strict prism volume division method by taking the sea depth as input data to obtain an effective integral radius;
the noise adding module is used for adding sea depth Gaussian white noise corresponding to sea depths with different relative accuracies on the basis of the error-free sea depth, and further generating sea depth data with different relative accuracies;
the forward modeling module is used for forward modeling to obtain gravity field element results corresponding to the error-free sea depth and the sea depth data with different relative precisions, which are obtained by the noise adding module, by taking the effective integral radius obtained by the resolving module as a basis and respectively taking the error-free sea depth and the sea depth data with different relative precisions obtained by the noise adding module as input;
the error making module is used for making an error between a gravity field element result obtained by forward modeling of the sea depth data with different relative precisions of the forward modeling module and a gravity field element result obtained by forward modeling of the error-free sea depth to obtain a forward gravity field element difference value;
and the statistical module is used for taking the difference value of the difference making module as input and counting the difference result index.
Further, the gravity field element comprises a sea surface gravity anomaly and a gravity anomaly vertical gradient.
Further, the difference result indicator includes a maximum value, a minimum value, an average value, and a standard deviation of the difference.
Compared with the prior art, the invention has the following beneficial effects:
the invention provides a method for calculating theoretical precision requirements of sea surface gravity anomaly and gravity anomaly vertical gradient measurement data under different inversion precision index requirements of submarine topography. The invention removes the entanglement influence of algorithm errors and data errors, completely realizes the requirement of calculating the inversion result quality of the submarine topography on the accuracy of the gravity input data from the input and output data level, and reveals the internal association between the gravity data (input data) and the submarine topography data (output data).
The method has beneficial reference values for researching the measurement precision requirements of different construction precision requirements of the submarine topography on the sea surface gravity data and utilizing the result to design a height measurement satellite constellation scheme for inverting the submarine topography by utilizing the satellite height measurement gravity data and recovering, analyzing and processing the marine gravity field data based on the satellite height measurement technology.
Drawings
FIG. 1 is a basic flow chart of a method for calculating the influence of sea surface gravity data measurement accuracy on a submarine topography inversion result according to an embodiment of the present invention;
FIG. 2 is an exemplary diagram of a noise-free sea depth model;
FIG. 3 is an exemplary diagram of the calculation results of different gravitational field components at an exemplary point;
FIG. 4 is an exemplary diagram of a model of different relative precisions of the sea depth of the target sea area;
fig. 5 is a schematic structural diagram of a device for calculating an influence of sea surface gravity data measurement accuracy on a submarine topography inversion result according to an embodiment of the present invention.
Detailed Description
The invention is further illustrated by the following examples in conjunction with the accompanying drawings:
as shown in fig. 1, a method for calculating the influence of the measurement accuracy of the sea surface gravity data on the inversion result of the submarine topography includes:
step S101: taking the sea depth as input data, resolving gravity field element information generated by the sea depth by using a strict prism volume division method, and obtaining an effective integral radius;
step S102: adding sea depth Gaussian white noises corresponding to sea depths with different relative accuracies on the basis of the error-free sea depth, and further generating sea depth data with different relative accuracies;
step S103: taking the effective integral radius obtained in the step S101 as a basis, respectively taking the error-free sea depth and the sea depth data with different relative precisions obtained in the step S102 as input, and forward-calculating to obtain gravity field element results corresponding to the error-free sea depth and the sea depth data with different relative precisions;
step S104: the gravity field element result obtained by forward modeling of the sea depth data with different relative precisions in the step S103 is subtracted from the gravity field element result obtained by forward modeling of the sea depth without error to obtain a forward gravity field element difference value; specifically, the forward gravity field element difference includes a forward gravity abnormal difference and a forward gravity abnormal vertical gradient difference;
step S105: and taking the difference obtained in the step S104 as an input, and counting the difference result index.
Further, the step S101 includes:
resolving gravity field element information generated by sea depth by using a strict prism volume division method, wherein the gravity field elements comprise sea surface gravity anomaly and gravity anomaly vertical gradient:
Figure BDA0003245886070000051
in the formula (I), the compound is shown in the specification,
Figure BDA0003245886070000052
wherein G is the earth gravitational potential constant, and takes the valueIs 6.672 x 10-8cm3/(g·s2)；ρcAnd ρwRespectively representing the density of the crust of the earth and the density of the seawater; (i)p,jp) To study the site location; (i, j) is the flow point position; h (i, j) is the corresponding sea depth of the (i, j) point; l is the distance between the point of interest and the flow point; (x, y, z) are flow point coordinates; sLAnd SBRespectively representing the points of the integral radius of the central area in the longitude and latitude directions; Δ x and Δ y represent grid sizes in the longitude and latitude directions, respectively; t iszFor disturbance gravity, the disturbance gravity is used to represent sea surface gravity anomaly; t iszzIs a gravity abnormal vertical gradient.
Further, the step S102 includes:
sea depth data of different relative precisions are generated by formula (3):
herr＝ho+hno (3)
in the formula, herrDifferent relative precision depths; h is0The sea depth is error-free; h isnoThe sea depth Gaussian white noise standard deviation is equal to the product of the error-free sea depth average value and the relative precision.
Further, in step S104, a forward gravitational field component difference value is obtained according to the formula (4):
Figure BDA0003245886070000061
in the formula, gammaiRepresenting the difference value of the ith forward acting gravity field element;
Figure BDA0003245886070000062
and
Figure BDA0003245886070000063
respectively representing the ith gravity field element values obtained by forward modeling of different relative precision sea depths and error-free sea depths; and N is the number of the statistical gravity difference values.
Further, in step S105, the difference result index includes:
Figure BDA0003245886070000064
in the formula, gammaiRepresenting the difference value of the ith forward acting gravity field element; delta gammamax、△γmin、△γmeanAnd Δ γstdThe maximum, minimum, mean and standard deviation of the difference are indicated, respectively.
Specifically, taking a 3 ° × 3 ° (111 ° E-114 ° E, 12 ° N-15 ° N) sea area in south china as an example, an S & S V19.1.1 error-free sea depth model is used as an error-free sea depth reference, and corresponding sea surface gravity anomaly and gravity anomaly vertical gradient accuracy requirements are respectively calculated under the conditions that the relative accuracies of the sea floor terrain are respectively 6%, 8% and 10%.
The method comprises the following steps: the calculation uses the central point (13 degrees 29 '30' N, 112 degrees 29 '30' E) of the south sea test sea area as an analysis sample point, and adopts the formula (1) to calculate the variation trend results of Tz and Tzz along with the integral radius under the noise-free south sea depth model (figure 2) (see figure 3). The effective integration radius is set to be about 61km, taking the figure 3 trend into account.
Step two: based on the S & S V19.1.1 error-free sea depth model, sea depth models with relative accuracies of 6%, 8%, and 10% were obtained using equation (3), and the results are shown in fig. 4.
Step three: and (3) calculating gravity field elements obtained by forward modeling of the sea depths with relative accuracies of 6%, 8% and 10% and gravity field element results obtained by forward modeling of the error-free sea depths respectively by adopting a formula (1) based on the effective integration radius of 61km obtained in the step one.
Step four: the gravity field elements obtained by forward performance of the sea depths with the relative precision of 6%, 8% and 10% obtained in the third step are subjected to error correction with the gravity field element result obtained by forward performance of the sea depths without errors through a formula (4), so that a forward gravity field element difference value is obtained;
step five: taking the difference obtained in the step four as input, and counting difference result indexes according to the formula (5); the difference statistics are shown in table 1.
TABLE 1 statistics of sea surface disturbance field element difference results corresponding to different relative precision sea depths
Figure BDA0003245886070000071
According to table 1, it can be known that the higher the sea depth relative accuracy is, the greater the corresponding sea surface gravity data accuracy requirement is, for example, when the sea depth relative accuracy in the south sea area is 6%, the corresponding sea surface Tz theoretical accuracy is 1.99 mGal; and when the relative precision of the sea depth of the sea area of the south sea is reduced to 8%, the theoretical precision of the corresponding sea surface Tz is reduced to 2.47 mGal. When the relative accuracy requirement of the sea depth of the south sea is within 10%, the corresponding sea surface Tz and Tzz accuracy at least reaches 3.23mGal and 18.25E theoretically.
On the basis of the foregoing embodiment, as shown in fig. 5, another aspect of the present invention provides an apparatus for calculating an influence of sea surface gravity data measurement accuracy on a submarine topography inversion result, including:
the calculating module 201 is used for calculating gravity field element information generated by the sea depth by using a strict prism volume division method by taking the sea depth as input data to obtain an effective integral radius;
the noise adding module 202 is used for adding sea depth Gaussian white noise corresponding to sea depths with different relative accuracies on the basis of the error-free sea depth, and further generating sea depth data with different relative accuracies;
the forward modeling module 203 is used for forward modeling to obtain gravity field element results corresponding to the error-free sea depth and the sea depth data with different relative accuracies, which are obtained by the noise adding module, by taking the effective integral radius obtained by the resolving module as a basis and respectively taking the error-free sea depth and the sea depth data with different relative accuracies obtained by the noise adding module as input;
the difference making module 204 is configured to make a difference between a gravity field element result obtained by forward modeling of the sea depth data with different relative accuracies of the forward modeling module and a gravity field element result obtained by forward modeling of the error-free sea depth, so as to obtain a forward-modeled gravity field element difference value;
and the statistical module 205 is configured to take the difference value of the difference module as an input, and count a difference result index.
Further, the resolving module 201 includes:
resolving gravity field element information generated by sea depth by using a strict prism volume division method, wherein the gravity field elements comprise sea surface gravity anomaly and gravity anomaly vertical gradient:
Figure BDA0003245886070000081
in the formula (I), the compound is shown in the specification,
Figure BDA0003245886070000082
wherein G is the earth gravitational potential constant and takes the value of 6.672 multiplied by 10-8cm3/(g·s2)；ρcAnd ρwRespectively representing the density of the crust of the earth and the density of the seawater; (i)p,jp) To study the site location; (i, j) is the flow point position; h (i, j) is the corresponding sea depth of the (i, j) point; l is the distance between the point of interest and the flow point; (x, y, z) is the flow point coordinate SLAnd SBRespectively representing the points of the integral radius of the central area in the longitude and latitude directions; Δ x and Δ y represent grid sizes in the longitude and latitude directions, respectively; t iszFor disturbance gravity, the disturbance gravity is used to represent sea surface gravity anomaly; t iszzIs a gravity abnormal vertical gradient.
Further, the noise adding module 202 includes:
sea depth data of different relative precisions are generated by formula (3):
herr＝ho+hno (3)
in the formula, herrDifferent relative precision depths; h is0The sea depth is error-free; h isnoThe sea depth Gaussian white noise standard deviation is equal to the product of the error-free sea depth average value and the relative precision.
Further, in the difference module 204, the forward gravitational field component difference value is obtained according to the formula (4):
Figure BDA0003245886070000083
in the formula, gammaiRepresenting the difference value of the ith forward acting gravity field element;
Figure BDA0003245886070000084
and
Figure BDA0003245886070000085
respectively representing the ith gravity field element values obtained by forward modeling of different relative precision sea depths and error-free sea depths; and N is the number of the statistical gravity difference values.
Further, in the statistical module 205, the difference result indicator includes:
Figure BDA0003245886070000091
in the formula, gammaiRepresenting the difference value of the ith forward acting gravity field element; delta gammamax、△γmin、△γmeanAnd Δ γstdThe maximum, minimum, mean and standard deviation of the difference are indicated, respectively.
In conclusion, the invention provides a calculation method for theoretical precision requirements of sea surface gravity anomaly and gravity anomaly vertical gradient measurement data under different inversion precision index requirements of submarine topography. The invention removes the entanglement influence of algorithm errors and data errors, completely realizes the requirement of calculating the inversion result quality of the submarine topography on the accuracy of the gravity input data from the input and output data level, and reveals the internal association between the gravity data (input data) and the submarine topography data (output data).
The method has beneficial values for exploring the measurement precision requirements of different construction precision requirements of the submarine topography on the sea surface gravity data, reflecting the underwater topography and topography as truly, completely, globally and accurately as possible according to the result, improving the global marine information acquisition, analysis and forecast capability, enhancing the marine cognition, marine forecast and marine information service capability, and realizing the forecast capability of hundreds of meters in the global marine mesoscale and key sea areas.
The above shows only the preferred embodiments of the present invention, and it should be noted that it is obvious to those skilled in the art that various modifications and improvements can be made without departing from the principle of the present invention, and these modifications and improvements should also be considered as the protection scope of the present invention.

Claims (8)

1.一种海面重力数据测量精度对海底地形反演结果影响计算方法，其特征在于，包括：1. a method for calculating the influence of sea surface gravity data measurement accuracy on the seabed topography inversion result, is characterized in that, comprises: 步骤1：将海深作为输入数据，利用严格棱柱体积分方法解算海深产生的重力场元信息，获取有效积分半径；Step 1: Take the sea depth as the input data, use the strict prism volume integration method to solve the gravity field element information generated by the sea depth, and obtain the effective integral radius; 步骤2：在无误差海深基础上添加不同相对精度海深对应的海深高斯白噪声，进而生成不同相对精度的海深数据；Step 2: Add sea depth Gaussian white noise corresponding to sea depths with different relative precisions on the basis of error-free sea depths, and then generate sea depth data with different relative precisions; 步骤3：以步骤1获得的有效积分半径为依据，分别以无误差海深、步骤2获得的不同相对精度的海深数据为输入，正演得到无误差海深及不同相对精度的海深数据对应的重力场元结果；Step 3: Based on the effective integral radius obtained in step 1, the error-free sea depth and the sea depth data with different relative accuracy obtained in step 2 are used as input, and the error-free sea depth and sea depth data with different relative accuracy are obtained through forward modeling. The corresponding gravity field element result; 步骤4：将步骤3不同相对精度的海深数据正演得到的重力场元结果与无误差海深正演得到的重力场元结果作差，得到正演的重力场元差值；Step 4: Make the difference between the gravity field element result obtained by forward modeling of ocean depth data with different relative precisions in step 3 and the gravity field element result obtained by the error-free ocean depth forward modeling to obtain the difference value of the forward modeling gravity field element; 步骤5：以步骤4得到的差值为输入，统计差值结果指标。Step 5: Take the difference obtained in step 4 as input, and count the difference result index. 2.根据权利要求1所述的一种海面重力数据测量精度对海底地形反演结果影响计算方法，其特征在于，所述步骤1包括：2. a kind of sea surface gravity data measurement accuracy according to claim 1 influences calculation method on bottom topography inversion result, is characterized in that, described step 1 comprises: 使用严格棱柱体积分方法解算海深产生的重力场元信息，所述重力场元包括海面重力异常和重力异常垂直梯度：Use the strict prismatic volume integral method to solve the gravity field element information generated by the ocean depth, the gravity field element includes the sea surface gravity anomaly and the gravity anomaly vertical gradient:
Figure FDA0003245886060000011
Figure FDA0003245886060000011
 式中，In the formula,
Figure FDA0003245886060000012
Figure FDA0003245886060000012
 其中，G为地球引力位常数，取值为6.672×10-8cm3/(g·s2)；ρc和ρw分别表示地壳密度和海水密度；(ip,jp)为研究点位置；(i,j)为流动点位置；h(i,j)为(i,j)点对应海深；l是研究点与流动点距离；(x,y,z)为流动点坐标；SL和SB分别表示经度和纬度方向中央区积分半径点数；△x和△y分别表示沿经度和纬度方向格网大小；Tz为扰动重力，此处用扰动重力表示海面重力异常；Tzz为重力异常垂直梯度。Among them, G is the gravitational potential constant of the earth, and the value is 6.672×10 -8 cm 3 /(g·s 2 ); ρ c and ρ w represent the density of the crust and seawater, respectively; ( ip , j p ) are the research points position; (i, j) is the position of the flow point; h(i, j) is the sea depth corresponding to the point (i, j); l is the distance between the research point and the flow point; (x, y, z) is the coordinate of the flow point; S L and S B represent the integral radius points in the central area in the longitude and latitude directions, respectively; △x and △y represent the grid size along the longitude and latitude directions, respectively; T z is the disturbance gravity, where the disturbance gravity is used to represent the sea surface gravity anomaly; T zz is the vertical gradient of gravity anomaly. 3.根据权利要求1所述的一种海面重力数据测量精度对海底地形反演结果影响计算方法，其特征在于，所述步骤2包括：3. a kind of sea surface gravity data measurement accuracy according to claim 1 influences calculation method on bottom topography inversion result, it is characterized in that, described step 2 comprises: 通过公式(3)生成不同相对精度的海深数据：The depth data of different relative precisions are generated by formula (3): herr＝h0+hno (3)h err = h 0 +h no (3) 式中，herr为不同相对精度海深；h0为无误差海深；hno是不同相对精度海深对应的海深高斯白噪声，海深高斯白噪声标准差等于无误差海深平均值与相对精度的乘积。In the formula, h err is the sea depth with different relative accuracy; h 0 is the error-free sea depth; h no is the sea depth Gaussian white noise corresponding to the sea depth with different relative accuracy, and the standard deviation of the sea depth Gaussian white noise is equal to the mean value of the sea depth without error. Product with relative precision. 4.根据权利要求1所述的一种海面重力数据测量精度对海底地形反演结果影响计算方法，其特征在于，所述步骤4中，按照(4)式得到正演的重力场元差值：4. The method for calculating the influence of sea surface gravity data measurement accuracy on seabed topography inversion results according to claim 1, wherein in the step 4, the forward modeled gravity field element difference is obtained according to formula (4) :
Figure FDA0003245886060000021
Figure FDA0003245886060000021
 式中，γi表示第i个正演的重力场元差值；
Figure FDA0003245886060000022
和
Figure FDA0003245886060000023
分别表示不同相对精度海深和无误差海深正演得到的第i个重力场元值；N为统计的重力差值个数。In the formula, γ i represents the gravitational field element difference of the i-th forward model;
Figure FDA0003245886060000022
and
Figure FDA0003245886060000023
Represents the i-th gravity field element value obtained by forward modeling of different relative precision ocean depths and error-free ocean depths respectively; N is the number of statistical gravity differences. 5.根据权利要求4所述的一种海面重力数据测量精度对海底地形反演结果影响计算方法，其特征在于，所述步骤5中，所述差值结果指标包括：5. The method for calculating the influence of sea surface gravity data measurement accuracy on seabed topography inversion results according to claim 4, wherein in the step 5, the difference result index comprises:
Figure FDA0003245886060000024
Figure FDA0003245886060000024
 式中，γi表示第i个正演的重力场元差值；△γmax、△γmin、△γmean和△γstd分别表示差值的最大值、最小值、平均值和标准差。In the formula, γ i represents the gravity field element difference of the i-th forward modeling; △γ max , △γ min , △γ mean and △γ std represent the maximum, minimum, average and standard deviation of the difference, respectively. 6.一种海面重力数据测量精度对海底地形反演结果影响计算装置，其特征在于，包括：6. A calculation device for the influence of sea surface gravity data measurement accuracy on seabed topography inversion results, characterized in that it comprises: 解算模块，用于将海深作为输入数据，利用严格棱柱体积分方法解算海深产生的重力场元信息，获取有效积分半径；The solution module is used to take the sea depth as input data, and use the strict prism volume integration method to solve the gravity field element information generated by the sea depth to obtain the effective integral radius; 加噪模块，用于在无误差海深基础上添加不同相对精度海深对应的海深高斯白噪声，进而生成不同相对精度的海深数据；The noise-adding module is used to add ocean depth Gaussian white noise corresponding to ocean depths with different relative precisions on the basis of error-free ocean depths, thereby generating ocean depth data with different relative precisions; 正演模块，用于以解算模块获得的有效积分半径为依据，分别以无误差海深、加噪模块获得的不同相对精度的海深数据为输入，正演得到无误差海深及不同相对精度的海深数据对应的重力场元结果；The forward modeling module is used to obtain the error-free sea depth and different relative accuracy based on the effective integral radius obtained by the solving module, and the error-free sea depth and the sea depth data of different relative precision obtained by the noise-added module as input respectively. Gravity field element results corresponding to accurate ocean depth data; 作差模块，用于将正演模块不同相对精度的海深数据正演得到的重力场元结果与无误差海深正演得到的重力重力场元结果作差，得到正演的重力场元差值；The difference module is used to make the difference between the gravity field element results obtained by forward modeling of ocean depth data with different relative precisions of the forward modeling module and the gravity field element results obtained by the error-free ocean depth forward modeling to obtain the forward modeling gravity field element difference. value; 统计模块，用于以作差模块的差值为输入，统计差值结果指标。The statistic module is used to take the difference value of the difference module as input, and count the difference value result indicators. 7.根据权利要求6所述的一种海面重力数据测量精度对海底地形反演结果影响计算装置，其特征在于，所述重力场元包括海面重力异常和重力异常垂直梯度。7 . The device for calculating the influence of sea surface gravity data measurement accuracy on seabed topography inversion results according to claim 6 , wherein the gravity field elements include sea surface gravity anomalies and gravity anomalies vertical gradients. 8 . 8.根据权利要求6所述的一种海面重力数据测量精度对海底地形反演结果影响计算装置，其特征在于，所述差值结果指标包括差值的最大值、最小值、平均值和标准差。8. The device for calculating the influence of sea surface gravity data measurement accuracy on seabed topography inversion results according to claim 6, wherein the difference result index comprises a maximum value, a minimum value, an average value and a standard value of the difference value. Difference.
CN202111033017.3A 2021-09-03 2021-09-03 Method and device for calculating influence of sea surface gravity data measurement accuracy on submarine topography inversion result Active CN113960690B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111033017.3A CN113960690B (en) 2021-09-03 2021-09-03 Method and device for calculating influence of sea surface gravity data measurement accuracy on submarine topography inversion result

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111033017.3A CN113960690B (en) 2021-09-03 2021-09-03 Method and device for calculating influence of sea surface gravity data measurement accuracy on submarine topography inversion result

Publications (2)

Publication Number Publication Date
CN113960690A true CN113960690A (en) 2022-01-21
CN113960690B CN113960690B (en) 2023-05-05

Family

ID=79460875

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111033017.3A Active CN113960690B (en) 2021-09-03 2021-09-03 Method and device for calculating influence of sea surface gravity data measurement accuracy on submarine topography inversion result

Country Status (1)

Country Link
CN (1) CN113960690B (en)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101416430B1 (en) * 2014-02-19 2014-07-10 한국지질자원연구원 Method for processing marine gravity data and Processing apparatus therefor
CN104035138A (en) * 2014-04-25 2014-09-10 西安测绘研究所 Accurate and rapid calculation method for global and local ocean disturbing gravity
US20150120196A1 (en) * 2013-10-29 2015-04-30 Weichang Li Method for Estimating Subsurface Properties from Geophysical Survey Data Using Physics-Based Inversion
CN106157368A (en) * 2016-06-25 2016-11-23 中国人民解放军信息工程大学 A kind of extensive area gravitational field Accurate Model and reconstructing method
CN108490496A (en) * 2018-03-26 2018-09-04 中国石油化工股份有限公司 Gravitational field inversion of Density method based on pseudo-radial basis function neural network
US20190004206A1 (en) * 2015-08-13 2019-01-03 Cgg Services Sas System and method for gravity and/or gravity gradient terrain corrections
CN110989021A (en) * 2019-12-03 2020-04-10 国家海洋局东海海洋环境调查勘察中心 Water depth inversion method and device and computer readable storage medium
CN112229403A (en) * 2020-08-31 2021-01-15 中国空间技术研究院 Method for improving marine gravity reconstruction precision based on geodetic level three-dimensional correction principle
CN112762927A (en) * 2020-12-18 2021-05-07 中国人民解放军战略支援部队信息工程大学 Semi-physical simulation method and system for underwater dynamic gravity data acquisition
CN112965124A (en) * 2021-02-08 2021-06-15 中国人民解放军92859部队 Method for calculating abnormal vertical gradient of external gravity by considering local guarantee conditions
CN113109883A (en) * 2021-05-26 2021-07-13 湖南科技大学 Method for forward modeling satellite gravity field based on isoparametric transformation global discrete grid spherical coordinates
CN113267822A (en) * 2021-04-30 2021-08-17 中国空间技术研究院 Method for improving inversion precision of submarine topography based on terrain constraint factor weight optimization

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150120196A1 (en) * 2013-10-29 2015-04-30 Weichang Li Method for Estimating Subsurface Properties from Geophysical Survey Data Using Physics-Based Inversion
KR101416430B1 (en) * 2014-02-19 2014-07-10 한국지질자원연구원 Method for processing marine gravity data and Processing apparatus therefor
CN104035138A (en) * 2014-04-25 2014-09-10 西安测绘研究所 Accurate and rapid calculation method for global and local ocean disturbing gravity
US20190004206A1 (en) * 2015-08-13 2019-01-03 Cgg Services Sas System and method for gravity and/or gravity gradient terrain corrections
CN106157368A (en) * 2016-06-25 2016-11-23 中国人民解放军信息工程大学 A kind of extensive area gravitational field Accurate Model and reconstructing method
CN108490496A (en) * 2018-03-26 2018-09-04 中国石油化工股份有限公司 Gravitational field inversion of Density method based on pseudo-radial basis function neural network
CN110989021A (en) * 2019-12-03 2020-04-10 国家海洋局东海海洋环境调查勘察中心 Water depth inversion method and device and computer readable storage medium
CN112229403A (en) * 2020-08-31 2021-01-15 中国空间技术研究院 Method for improving marine gravity reconstruction precision based on geodetic level three-dimensional correction principle
CN112762927A (en) * 2020-12-18 2021-05-07 中国人民解放军战略支援部队信息工程大学 Semi-physical simulation method and system for underwater dynamic gravity data acquisition
CN112965124A (en) * 2021-02-08 2021-06-15 中国人民解放军92859部队 Method for calculating abnormal vertical gradient of external gravity by considering local guarantee conditions
CN113267822A (en) * 2021-04-30 2021-08-17 中国空间技术研究院 Method for improving inversion precision of submarine topography based on terrain constraint factor weight optimization
CN113109883A (en) * 2021-05-26 2021-07-13 湖南科技大学 Method for forward modeling satellite gravity field based on isoparametric transformation global discrete grid spherical coordinates

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
MASOOD AZIZI: ""Integrating Gravity Data With Remotely Sensed Data for Structural Investigation of the Aynak-Logar Valley, Eastern Afghanistan, and the Surrounding Area"", 《IEEE JOURNAL OF SELECTED TOPICS IN APPLIED EARTH OBSERVATIONS AND REMOTE SENSING》 *
欧阳明达: ""利用海洋重力数据反演海底地形的理论与方法"", 《中国优秀硕士学位论文全文数据库》 *
范雕: ""海底地形高次项对海面重力信息影响分析"", 《武汉大学学报· 信息科学版》 *

Also Published As

Publication number Publication date
CN113960690B (en) 2023-05-05

Similar Documents

Publication Publication Date Title
Yang et al. Refraction correction of airborne LiDAR bathymetry based on sea surface profile and ray tracing
Hare et al. Modelling bathymetric uncertainty
EP3084350B1 (en) Method and system for generating a geoid via three computation spaces and airborne-acquired gravity data
CN114137624A (en) Method and system for inverting submarine topography based on satellite altimeter
CN109085655B (en) Underwater platform gravity measurement scheme and verification method
Zhang et al. Refraction and coordinate correction with the JONSWAP model for ICESat-2 bathymetry
CN113341476A (en) Method for improving ocean gravity spatial resolution based on submarine topography-gravity combination
CN101957193B (en) An optimization method for elevation transfer of islands and reefs
Yildiz et al. Determination and validation of the Turkish Geoid Model-2020 (TG-20)
Guo et al. Validation of just-released SWOT L2 KaRIn beta prevalidated data based on restore the marine gravity field and its application
Abdelaal et al. Assessing coastal heritage sustainability: Crustal deformation and sea-level trends at the Qaitbay Citadel in Alexandria, Egypt
Partama et al. A simple and empirical refraction correction method for UAV-based shallow-water photogrammetry
Sikder et al. Assessment of orthometric height derived from levelling, GNSS and EGM2008 geoid model in Bangladesh
Chu Two types of absolute dynamic ocean topography
Jiang et al. Enhanced gravity-geologic method to predict bathymetry by considering non-linear effects of surrounding seafloor topography
CN113960690A (en) A method and device for calculating the influence of sea surface gravity data measurement accuracy on the results of seabed topography inversion
Timilsina et al. Regional geoid for Nepal using least-squares collocation
CN108398126A (en) A kind of high-precision air-sea gravity measurement platform inclination correction model
Lee et al. Transformation of vertical datum surface in the coastal area using hybrid geoid models
Huang et al. Research and evaluation on key technological indicators for airborne and shipborne gravimetry
CN109297530A (en) A deep-sea full-water column multi-environmental element information fusion method and processing terminal
Wang et al. Estimating the true spatial resolution of satellite altimeter-derived gravity field models with shipborne data in the South China Sea
Khedr Sea level modeling and realization of lowest astronomical tide (LAT) for Alexandria Harbor, Egypt
CN113868838B (en) A method and system for evaluating the advantages of inverting seafloor terrain models using disturbance field elements
Ziebart et al. VORF-The UK Vertical Offshore Reference Frame: Enabling Real-time Hydrographic Surveying

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: 450000 Science Avenue 62, Zhengzhou High-tech Zone, Henan Province

Patentee after: Information Engineering University of the Chinese People's Liberation Army Cyberspace Force

Country or region after: China

Address before: No. 62 Science Avenue, High tech Zone, Zhengzhou City, Henan Province

Patentee before: Information Engineering University of Strategic Support Force,PLA

Country or region before: China