阅读说明：本技术 一种海底重力测量中区—远一区地形改正计算方法 (Submarine gravity measurement middle-far one-region terrain correction calculation method ) 是由 刘洪波 孔祥明 张一� 崔灿 齐永亮 于嘉宾 于 2019-10-12 设计创作，主要内容包括：本发明涉及一种海底重力测量中区—远一区地形改正计算方法,属于重力测量领域,方法包括建立海底地形节点网、选取海底重力测点在海底地形节点网中最近的四个节点、利用改正计算公式计算每个节点的地形改正值,再通过X、Y方向双线性插值法插值到重力测点位置上作为重力测点地形改正值,本发明针对海底重力测量的特点,充分考虑海水引力对重力测量的影响,提出了修正海水引力影响的海底重力地形改正计算方法,推导了可适用于从海岸线到深海全海域海底重力测量的中区—远一区地形改正计算公式,以提高海底重力测量地形改正的精度,进而提高了海底重力测量异常精度。(The invention relates to a middle area-far area terrain correction calculation method in submarine gravity measurement, belonging to the gravity measurement field, the method comprises establishing a submarine terrain node network, selecting four nearest nodes of a submarine gravity measuring point in the submarine terrain node network, calculating the terrain correction value of each node by using a correction calculation formula, and interpolating to the gravity measuring point position by a bilinear interpolation method in direction X, Y to be used as the terrain correction value of the gravity measuring point, aiming at the characteristics of submarine gravity measurement, fully considering the influence of seawater gravity on gravity measurement, providing a submarine gravity terrain correction calculation method for correcting the influence of seawater gravity, deducing a middle area-far area terrain correction calculation formula which can be suitable for the gravity measurement from a coastline to the submarine all-sea area, the method improves the accuracy of correcting the topography of the submarine gravity measurement, and further improves the abnormal accuracy of the submarine gravity measurement.)

1. A submarine gravity measurement middle-far area terrain correction calculation method is characterized by comprising the following steps:

(1) dividing the terrain around the gravity measuring point as the center into small blocks by using the collected submarine topography map, and establishing a submarine topography node network;

(2) determining the position of a gravity measuring point in a submarine topography node network according to the submarine gravity measuring point coordinates;

(3) selecting four nearest nodes of the seabed gravity measuring point in a seabed terrain node network, wherein the four nodes are named as A, B, C, D;

(4) calculating the terrain correction value delta g of each node according to the submarine terrain node network by using a formula (2)_A、△g_B、△g_C、△g_D；

Correcting a calculation formula:

the values of the parameters in formula (2) are as follows:

ρ₀the density of the seawater is 1.03 x 10³g/cm³；

h₀-height values of gravity measurement points;

S_ijunderwater topographic coefficient when (h)_ij+h₀) When the value is negative, the value is 1, (h)_ij+h₀) Taking 0 when the value is not negative;

T_ijthe water topographic coefficient when (h)_ij+h₀) When the value is negative, 0 is taken, (h)_ij+h₀) Taking 1 when the value is not negative;

ρ₁-average density of crust, 2.67 of10³g/cm³；

l is integration grid distance;

C_ij-integrating constants, selecting trapezoidal coefficients;

r_ij-the distance between the integration node (i, j) and the calculation point;

h_ij-the elevation difference between the integration node (i, j) and the calculation point;

g-universal constant of gravity (6.67X 10)^-11m³/(kg·s²))；

The left delta g of the formula 2 is a terrain correction value, and the elevation value h of the gravity measurement point is used when the terrain correction value of the node is calculated₀Replacing the elevation values of the four nodes;

(5) interpolating to the position of the gravity measuring point by an X, Y direction bilinear interpolation method according to the topographic correction value of the peripheral node of the gravity measuring point to be used as the topographic correction value of the gravity measuring point;

(6) and (5) comprehensively correcting the influence of the rocks and the seawater on the gravity measurement through the gravity measuring point terrain correction value calculated in the step (5), wherein the obtained gravity measuring point terrain correction value is used for accurately calculating the grid abnormal value.

2. The correction and calculation method for the topography of the middle-distant area in the submarine gravity measurement according to claim 1, wherein in the step (1), the topography within 20km around the gravity measurement point is divided into small blocks to establish the submarine topography node network.

3. The method for correcting and calculating the terrain in the middle-far area in the seafloor gravimetry according to claim 1, wherein in the step (1), the terrain-divided small blocks are square small blocks.

4. The method for correcting and calculating the topography of the middle-to-far area in the submarine gravity measurement according to claim 1, wherein in the step (1), more than 1: the mesh degree of the established node network is more than 5m multiplied by 5 m.

5. The method for correcting and calculating the topography of the middle-to-far area in the submarine gravity measurement according to claim 1, wherein in the step (5), the bilinear interpolation formula of the bilinear interpolation method is as follows:

in formula (3): delta g_p-a topographical correction value for the gravity measurement point P; delta g_A、△g_B、△g_C、△g_D-a terrain correction value for node A, B, C, D; delta x, delta y-node grid distance; (X, Y) -coordinates of a point P of the gravity measuring point; (X)_i,Y_j) -the coordinates of point a; (X)_i,Y_j+1) -coordinates of point B; (X)_i+1,Y_j) -coordinates of point C; (X)_i+1,Y_j+1) -coordinates of the D-point.

Technical Field

The invention relates to a method for correcting and calculating the topography of a middle area and a far area in submarine gravity measurement, in particular to the method for correcting and calculating the topography within the range of 100 m-20000 m in submarine high-precision gravity measurement, and belongs to the field of gravity measurement.

Background

At present, in the various modifications of the submarine gravity measurement, the terrain correction value is almost the same order of magnitude as the gravity abnormal value generated by an abnormal body, and the accuracy of the terrain correction plays a key role in the accuracy of gravity exploration and interpretation. The terrain correction calculation formula of the middle area-the far area specified by the national gravity survey standard is as follows:

the values of the parameters in the formula are as follows:

ρ₁-average density of crust of 2.67X 10³g/cm³；

l is integration grid distance;

C_ij-integrating constants, selecting trapezoidal coefficients;

r_ij-the distance between the integration node (i, j) and the calculation point;

h_ij-the difference in elevation between the integration node (i, j) and the calculation point

When the submarine gravity measurement work is carried out, the influence of seawater is not considered when the terrain correction is carried out by adopting the formula (1), and the water body correction is not carried out when the terrain is corrected. Along with the expansion of the gravity measurement work to the sea area, when the submarine gravity measurement is carried out, especially when the submarine high-precision gravity measurement is carried out, the correction method of the topography of the middle area to the far area directly influences the accuracy of the gravity anomaly calculation and influences the precision of the gravity anomaly.

Disclosure of Invention

Aiming at the defects of the prior art, the invention provides a correction calculation method for the topography of a middle area-a far area in the submarine gravity measurement.

The technical scheme of the invention is as follows:

a terrain correction calculation method for a middle area and a far area in submarine gravity measurement comprises the following steps:

(1) dividing the terrain around the gravity measuring point as the center into a plurality of small blocks by using the collected submarine topography map, and establishing a submarine topography node network;

(2) determining the position of a gravity measuring point in a submarine topography node network according to the submarine gravity measuring point coordinates;

(3) selecting four nearest nodes of the seabed gravity measuring point in a seabed terrain node network, wherein the four nodes are named as A, B, C, D;

(4) calculating the terrain correction value delta g of each node according to the submarine terrain node network by using a formula (2)_A、△g_B、△g_C、△g_D。

Correcting a calculation formula:

the values of the parameters in formula (2) are as follows:

ρ₀the density of the seawater is 1.03 x 10³g/cm³；

h₀-height values of gravity measurement points;

S_ijunderwater topographic coefficient when (h)_ij+h₀) When the value is negative, the value is 1, (h)_ij+h₀) Taking 0 when the value is not negative;

T_ijthe water topographic coefficient when (h)_ij+h₀) When the value is negative, 0 is taken, (h)_ij+h₀) Taking 1 when the value is not negative;

ρ₁-average density of crust of 2.67X 10³g/cm³；

l is integration grid distance;

C_ij-integrating constants, selecting trapezoidal coefficients;

r_ij-the distance between the integration node (i, j) and the calculation point;

h_ij-the elevation difference between the integration node (i, j) and the calculation point;

g-universal constant of gravity (6.67X 10)^-11m³/(kg·s²))；

The left delta g of the formula 2 is a terrain correction value, and the elevation value h of the gravity measurement point is used when the terrain correction value of the node is calculated₀Instead of elevation values for four nodes, i.e. h in the formula when calculating the terrain correction value for each node₀All the gravity measuring points have the same elevation value;

(5) interpolating to the position of the gravity measuring point by an X, Y direction bilinear interpolation method according to the topographic correction value of the peripheral node of the gravity measuring point to be used as the topographic correction value of the gravity measuring point;

(6) and (5) comprehensively correcting the influence of the rocks and the seawater on the gravity measurement through the gravity measuring point terrain correction value calculated in the step (5), wherein the obtained gravity measuring point terrain correction value is used for accurately calculating the grid abnormal value.

Preferably, in the step (1), the terrain within the range of 20km around the gravity measuring point as the center is divided into a plurality of small blocks, and the submarine terrain node network is established.

Preferably, in step (1), the terrain-segmented small blocks are square small blocks.

Preferably, in step (1), more than 1: the mesh degree of the established node network is more than 5m multiplied by 5 m.

Preferably, in step (5), the bilinear interpolation formula is:

in formula (3): delta g_p-a topographical correction value for the gravity measurement point P; delta g_A、△g_B、△g_C、△g_D-a terrain correction value for node A, B, C, D; delta x, delta y-node grid distance; (X, Y) -coordinates of a point P of the gravity measuring point; (X)_i,Y_j) -the coordinates of point a; (X)_i,Y_j+1) -coordinates of point B; (X)_i+1,Y_j) -coordinates of point C; (X)_i+1,Y_j+1) -coordinates of the D-point.

The invention has the beneficial effects that:

the invention provides a submarine gravity terrain correction calculation method for correcting the influence of seawater gravity on gravity measurement by fully considering the influence of the seawater gravity on the gravity measurement aiming at the characteristics of submarine gravity measurement, and deduces a middle-far-one-region terrain correction calculation formula applicable to submarine gravity measurement from a coastline to a deep sea whole sea area so as to improve the precision of submarine gravity measurement terrain correction and further improve the abnormal precision of submarine gravity measurement.

Drawings

FIG. 1 is a schematic diagram of a terrain partitioning domain according to the present invention;

fig. 2 is a schematic view of the terrain effects.

Detailed Description

The present invention will be further described by way of examples, but not limited thereto, with reference to the accompanying drawings.