CN110927823B - Grid value correction method applied to shipborne gravimeter - Google Patents

Grid value correction method applied to shipborne gravimeter Download PDF

Info

Publication number
CN110927823B
CN110927823B CN201911227736.1A CN201911227736A CN110927823B CN 110927823 B CN110927823 B CN 110927823B CN 201911227736 A CN201911227736 A CN 201911227736A CN 110927823 B CN110927823 B CN 110927823B
Authority
CN
China
Prior art keywords
gravity
value
gravimeter
wharf
grid
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.)
Active
Application number
CN201911227736.1A
Other languages
Chinese (zh)
Other versions
CN110927823A (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.)
First Institute of Oceanography MNR
Original Assignee
First Institute of Oceanography MNR
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 First Institute of Oceanography MNR filed Critical First Institute of Oceanography MNR
Priority to CN201911227736.1A priority Critical patent/CN110927823B/en
Publication of CN110927823A publication Critical patent/CN110927823A/en
Application granted granted Critical
Publication of CN110927823B publication Critical patent/CN110927823B/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
    • G01V13/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00
    • 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
    • G01V7/16—Measuring gravitational fields or waves; Gravimetric prospecting or detecting specially adapted for use on moving platforms, e.g. ship, aircraft

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Manufacturing & Machinery (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Ocean & Marine Engineering (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention discloses a grid value correction method applied to a ship-borne gravimeter, which is characterized in that a ship-borne gravimeter is used for reciprocating departure from a harbor and stopping at a replenishing harbor in midway to serve as a basic mathematical model, and an error compensation correction model for observing a gravity value is established so as to carry out error compensation on the observed gravity value; calculating a grid value based on the corresponding relation between the compensated gravity observation value and the absolute gravity value, and compensating the data of the shipborne gravimeter; the method can eliminate the system error caused by the lattice value error, obtain high-precision voyage survey result data, and has important guiding value and practical significance for marine geophysical survey, especially marine gravity measurement of deep and distant voyages.

Description

Grid value correction method applied to shipborne gravimeter
Technical Field
The invention belongs to the field of application of difference calibration and data processing of an observation data system of a ship-borne gravimeter, and particularly relates to a grid value correction method applied to the ship-borne gravimeter.
Background
As an important component of the research of the earth gravimetry, the marine gravimetry provides gravity data for researching the earth spherical shape and the earth internal structure, guaranteeing the spaceflight and the remote weapon launching, exploring marine mineral resources and the like. The shipborne gravimeter is used as a key core device for marine geophysical survey, and besides various system and accidental errors caused by the influence of factors such as external environment, the level of field operators and the like, the grid value precision of the shipborne gravimeter can also be changed in the using process of the device.
Before the influence of factors such as shipbuilding technology and ship endurance, marine gravity survey is mainly focused on offshore areas, a port is unified to and fro or a nearby port, and the grid value of the gravimeter is calibrated by a gravimeter manufacturer at the beginning of purchase. In the actual survey process, the zero drift condition of the shipborne gravimeter can be easily judged by a mode of going to and fro a mother port, but whether the grid value is correct and how the accuracy of the grid value is often ignored subjectively. With the annual increase of global deep and far sea voyages, more and more voyages need to stop at a plurality of wharfs in the investigation process, and objectively have the condition and the requirement of long-baseline observation.
At present, domestic marine gravimeters are still widely used in domestic marine gravimeter surveys, the conventional method is to transport the gravimeter back to a production place for calibration or select two long base line points on land for calibration, however, as a fine instrument, a ship-borne gravimeter is very easy to be influenced by factors such as collision during transportation, and the secondary change of the grid value is caused. With the continuous development of scientific technology, deep and open sea gravity investigation is more and more common, and one or more wharfs are often parked in the middle of the deep and open sea gravity investigation.
Disclosure of Invention
The invention provides a grid value correction method applied to a shipborne gravimeter, which is characterized in that the grid value of the shipborne gravimeter in a single voyage is corrected and calculated based on the fact that the shipborne gravimeter comes and goes to and from a departure/return port and the gravity data of a midway stop replenishment port are taken as mathematical models, so that errors caused by the grid value errors are effectively eliminated, and high-precision voyage survey result data are obtained.
The invention is realized by adopting the following technical scheme: a grid value correction method applied to a ship-borne gravimeter comprises the following steps:
step S1, establishing an error compensation correction model of the observed gravity value by taking the gravity data of the ship-borne gravimeter which comes and goes to and fro the departure/return port and stops midway at the replenishment port as a basic mathematical model so as to carry out error compensation on the observed gravity value;
and step S2, obtaining a grid value based on the corresponding relation between the compensated observed gravity value and the absolute gravity value, and compensating the data of the shipborne gravimeter.
Further, the step S1 of performing error compensation on the observed gravity value specifically includes the following steps:
step S11, analyzing and respectively obtaining the gravity observation points of the departure wharf, the supply wharf and the return wharf and the changes of the normal gravity field values of the corresponding wharf base points caused by the latitude changes;
and step S12, calculating the gravity observation data of the starting wharf to the elevation surface of the gravity base point to obtain the calculated gravimeter data at the gravity base point of the starting wharf.
Further, the step S11 obtains the normal gravity field value variation due to the latitude variation specifically by the following means:
Figure BDA0002302711340000021
gL＝0.000814*sin 2L*dL
wherein: wherein: dLThe distance between the shipborne gravimeter and the base point in the north-south direction is the unit of meter (m); a is an included angle between a wharf (an sailing wharf, a supply wharf and a returning wharf) and the north direction, and rad; gLFor normal gravity field variations, 10-5m/s2(ii) a L is the gravity observation point latitude, rad.
Further, in step S12, the following formula is used for calculation:
Figure BDA0002302711340000022
Gbase opener＝GOpener+0.3086*h/K0+gL/K0
Wherein: h is the height from the center of gravity of the shipborne gravimeter to a base point, and the unit is meter (m); a is the distance from the ship port to the water surface, and the unit is meter (m); b is the distance from the starboard of the ship to the water surface, and the unit is meter (m); c is the distance from the center of gravity of the shipborne gravimeter to the deck, and the unit is meter (m); d is the distance from the gravity base point of the wharf to the water surface, and the unit is meter (m); gOpenerThe gravity observed value of the starting wharf is checked; gBase openerCalculating gravimeter data at the gravity base point; k0Original grid value of a shipborne gravimeter, 10-5m/s2And/grid.
Further, the step S2 specifically includes the following steps:
step S21, analyzing and obtaining gravimeter data G of the replenishment dock and the return dock to the gravity base point based on the step S1 respectivelyRadical supplementAnd GRadical reversionCalculating to obtain a zero drift rate by adopting gravity observation data of the starting wharf and the returning wharf, and compensating an observation gravity value G of the replenishment wharf;
step S22, calculating to obtain a shipborne gravimeter grid value K:
Figure BDA0002302711340000023
wherein: k is the on-board gravimeter grid value after this calculation is performed, 10-5m/s2A/grid; gInsulation opener、GJu BuThe absolute gravity values at the base points of the starting wharf and the replenishment wharf are directly measured by adopting an absolute gravimeter.
Further, the step S21 is specifically analyzed in the following manner:
Figure BDA0002302711340000024
wherein: t is t1、t2、t3Respectively as gravity observed values GOpener、GSupplement device、GReturn toObservation time of (1), GSupplement deviceGravity observation for supply terminals, GReturn toAnd the gravity observed value of the return wharf is obtained.
Further, the step S2 is followed by:
step S3: and based on the lattice value obtained in the step S2, obtaining an updated lattice value by iterating the steps S1 to S2 for multiple times until the absolute value of the difference between the lattice values K obtained in two adjacent iterations is smaller than an allowable error, and further obtaining a corrected lattice value.
Compared with the prior art, the invention has the advantages and positive effects that:
the scheme takes the gravity data of a ship-borne gravimeter which comes and goes to and fro the sailing/returning port and stops midway at a replenishing port as a basic mathematical model, and establishes an error compensation correction model for observing the gravity value so as to carry out error compensation on the observed gravity value; calculating a grid value based on the corresponding relation between the compensated gravity observation value and the absolute gravity value, and compensating the data of the shipborne gravimeter; the method can eliminate the system error caused by the lattice value error, obtain high-precision voyage survey result data, and has important guiding value and practical significance for marine geophysical survey, especially marine gravity measurement of deep and distant voyages.
Drawings
FIG. 1 is a schematic diagram illustrating a comparison between a ship's reciprocating motion and a base point according to an embodiment of the present invention;
FIG. 2 is a graph comparing pre-calibration ship data and satellite data in accordance with an embodiment of the present invention;
FIG. 3 is a graph comparing ship data and satellite data after correcting the lattice values in the embodiment of the present invention.
Detailed Description
In order to clearly understand the above objects and advantages of the present invention, the following detailed description of the embodiments of the present invention is made with reference to the accompanying drawings, and it should be noted that in the present embodiment, the departure port and the return port may be the same port, or two different ports with smaller absolute gravity difference; the departure wharf is a certain wharf in the departure wharf, and the return wharf and the supply wharf are in the same way.
The embodiment provides a grid value correction method applied to a ship-borne gravimeter, which comprises the following steps:
the method comprises the following steps that firstly, an error compensation correction model for observing a gravity value is established by taking gravity data of a shipborne gravimeter which comes and goes to and from a port of departure and stops midway at a port of replenishment as a basic mathematical model, so that the error compensation is carried out on the observed gravity value, and the method specifically comprises the following steps:
(1) analyzing and respectively obtaining the gravity observation points of the starting dock, the supply dock and the returning dock and the normal field value change of each corresponding dock base point due to the latitude change:
Figure BDA0002302711340000031
gL＝0.000814*sin 2L*dL
wherein: dLThe distance between the shipborne gravimeter and the base point in the north-south direction is the unit of meter (m); a is an included angle between a wharf (an sailing wharf, a supply wharf and a returning wharf) and the north direction, and rad; gLFor normal gravity field variations, 10-5m/s2(ii) a L is the gravity observation point latitude, rad.
(2) Wharf to be started (gravity observation value G)Opener) The gravity observation data is reduced to the elevation surface of the gravity base point, and the following formula is adopted for calculation:
Figure BDA0002302711340000041
Gbase opener＝GOpener+0.3086*h/K0+gL/K0
Wherein: h is the height from the gravity center of the gravity meter to a base point, and is meter (m); as shown in fig. 1, a is the ship port to surface distance, meters (m); b is the distance from the starboard of the ship to the water surface, meter (m); c is the distance from the gravity center of the gravity meter to the deck, meter (m); d is the distance from the gravity base point to the water surface, namely meter (m); gBase openerCalculating gravimeter data at the gravity base point; k0Original grid value of gravimeter, 10-5m/s2And/grid.
Step two, calculating a grid value based on the corresponding relation between the compensated observation gravity value and the absolute gravity value, and compensating the data of the shipborne gravimeter, wherein the method specifically comprises the following steps:
(3) respectively calculating a replenishment dock (gravity observation value G) according to the step (1) and the step (2)Supplement device) Return dock (gravity observation value G)Return to) Gravimeter data G reduced to the base point of gravityRadical supplementAnd GRadical reversionCalculating to obtain the zero drift rate by adopting the gravity observation data of the departure port and the return port, and compensating and calculating the observation gravity value G at the supply dock, wherein the method specifically comprises the following steps:
Figure BDA0002302711340000042
wherein: t is t1、t2、t3Respectively as gravity observed values GOpener、GSupplement device、GReturn toThe observation time of (2).
(4) Calculating the grid value of the shipborne gravimeter:
Figure BDA0002302711340000043
wherein: k is the on-board gravimeter grid value after this calculation is performed, 10-5m/s2A/grid; gInsulation opener、GJu BuThe absolute gravity values at the base points of the starting wharf and the replenishment wharf respectively can be directly measured by adopting an absolute gravimeter.
Step three, step (3) calculating the original gravimeter grid value K adopted when calculating the gravimeter data reduced to the gravity base point0And (4) updating the gravimeter grid value to K after calculation in the step (4). Based on the obtained latest lattice value K, adopting a mode of iterating the steps (1) to (4) for multiple times until the absolute value of the difference value of the lattice value K calculated twice before and after is less than 10-5×10-5m/s2Stopping iteration when/lattice.
The method can be applied to zero drift linearity of the gravimeter which meets the national standard survey requirements and stops two or more shipborne gravity surveys or grid value correction voyages. In practical application, when the data processing of a survey area is completed through zero drift control of a starting and returning wharf, the problem that an overlarge system difference exists between the data and the satellite gravity data is found, and the situation that a gravity instrument is inaccurate is deduced by the system difference. By the method, the corresponding system difference can be obviously removed, and the defect of the voyage data can be effectively made up.
According to the scheme, the systematic deviation of the gravimeter grid value possibly existing in the deep sea gravity investigation process is considered, the shipborne gravimeter grid value is corrected by utilizing the gravity data difference between the shipborne gravimeter which executes the field investigation voyage and the wharf, the rationality of the method is verified by utilizing the 1 '× 1' gridding data (version V27.1) of Sandwell and Smith, and the calibration of the shipborne data grid value of a certain voyage is taken as an example and is shown in the attached figures 2 and 3. Before the grid value is corrected, obvious system difference exists between ship measurement data and satellite data (figure 2), the ship measurement data and the satellite gravity data are well fitted after the grid value is corrected by the method (figure 3), and the ship measurement data carry a large amount of high-frequency information missing from the satellite data. The test shows that the method can obviously remedy systematic deviation caused by grid value errors, is easy to realize along with the annual increase of ocean global voyages, and can meet the requirements of shipborne gravity measurement users.
The above description is only a preferred embodiment of the present invention, and not intended to limit the present invention in other forms, and any person skilled in the art may apply the above modifications or changes to the equivalent embodiments with equivalent changes, without departing from the technical spirit of the present invention, and any simple modification, equivalent change and change made to the above embodiments according to the technical spirit of the present invention still belong to the protection scope of the technical spirit of the present invention.

Claims (4)

1. A lattice value correction method applied to a ship-borne gravimeter is characterized by comprising the following steps:
step S1, establishing an error compensation correction model of the observed gravity value by taking the gravity data of the ship-borne gravimeter which comes and goes to and fro the departure/return port and stops midway at the replenishment port as a basic mathematical model so as to carry out error compensation on the observed gravity value;
step S2, obtaining a grid value based on the corresponding relation between the compensated observation gravity value and the absolute gravity value, and compensating the data of the shipborne gravimeter;
the step S2 specifically includes the following steps:
step S21, analyzing and obtaining gravimeter data G of the replenishment dock and the return dock to the gravity base point based on the step S1 respectivelyRadical supplementAnd GRadical reversionCalculating to obtain a zero drift rate by adopting gravity observation data of the starting wharf and the returning wharf, and compensating an observation gravity value G of the replenishment wharf;
Figure FDA0002656681460000011
wherein: t is t1、t2、t3Respectively as gravity observed values GOpener、GSupplement device、GReturn toObservation time of (1), GSupplement deviceGravity observation for supply terminals, GReturn toThe gravity observed value of the return wharf is obtained;
step S22, calculating to obtain a shipborne gravimeter grid value K:
Figure FDA0002656681460000012
wherein G isInsulation opener、GJu BuAbsolute gravity values at the base points of the starting wharf and the replenishment wharf are respectively;
and step S3, based on the grid values obtained in step S2, calculating the updated grid values by adopting a mode of iterating steps S1 to S2 for multiple times until the absolute value of the difference value of the grid values K obtained in two adjacent iterations is smaller than an allowable error, and further obtaining the corrected grid values.
2. The grid value correction method applied to the ship-borne gravimeter according to claim 1, characterized in that: the error compensation of the observed gravity value in step S1 specifically includes the following steps:
step S11, analyzing and respectively obtaining the gravity observation points of the departure wharf, the supply wharf and the return wharf and the changes of the normal gravity field values of the corresponding wharf base points caused by the latitude changes;
and step S12, the gravity observation data of each wharf are reduced to the elevation surface of the gravity base point, and the gravity meter data at the gravity base point of each wharf after reduction are obtained.
3. The grid value correction method applied to the ship-borne gravimeter according to claim 2, characterized in that: the step S11 specifically obtains the normal gravity field value change by:
Figure FDA0002656681460000013
gL＝0.000814*sin2L*dL
wherein: dLThe distance between the shipborne gravimeter and the base point in the north-south direction is the distance between the shipborne gravimeter and the base point; a is an included angle between a wharf and the due north direction, and the wharf comprises a sailing wharf, a supply wharf and a return wharf; gLNormal gravity field changes; l is the gravity observation point latitude.
4. The grid value correction method applied to the ship-borne gravimeter according to claim 3, characterized in that: in step S12, the reduced gravimeter data at the gravity base point of the starting dock is calculated by using the following formula:
Figure FDA0002656681460000021
Gbase opener＝GOpener+0.3086*h/K0+gL/K0
Wherein: h is the height from the center of gravity of the shipborne gravimeter to a base point; a is the distance from the ship port to the water surface; b is the distance from the starboard of the ship to the water surface; c is the distance from the center of gravity of the shipborne gravimeter to the deck; d is the distance from the gravity base point of the wharf to the water surface; gOpenerThe gravity observed value of the sailing wharf is obtained; gBase openerCalculating gravimeter data at the gravity base point; k0Original grid value of the shipborne gravimeter.
CN201911227736.1A 2019-12-04 2019-12-04 Grid value correction method applied to shipborne gravimeter Active CN110927823B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911227736.1A CN110927823B (en) 2019-12-04 2019-12-04 Grid value correction method applied to shipborne gravimeter

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911227736.1A CN110927823B (en) 2019-12-04 2019-12-04 Grid value correction method applied to shipborne gravimeter

Publications (2)

Publication Number Publication Date
CN110927823A CN110927823A (en) 2020-03-27
CN110927823B true CN110927823B (en) 2021-03-16

Family

ID=69856644

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911227736.1A Active CN110927823B (en) 2019-12-04 2019-12-04 Grid value correction method applied to shipborne gravimeter

Country Status (1)

Country Link
CN (1) CN110927823B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112415634B (en) * 2020-10-27 2021-12-07 青岛海洋地质研究所 Dynamic gravimeter zero drift compensation method based on satellite gravity anomaly information
CN112487604B (en) * 2020-10-27 2022-08-16 青岛海洋地质研究所 Long-time nonlinear drift compensation method for output data of marine gravimeter
CN112731545B (en) * 2020-12-18 2022-07-22 自然资源部第一海洋研究所 Ocean irregular net measurement data adjustment method

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6615139B1 (en) * 2002-03-28 2003-09-02 Council Of Scientific & Industrial Research Digitally implemented method for automatic optimization of gravity fields obtained from three-dimensional density interfaces using depth dependent density
CN101201412A (en) * 2006-12-13 2008-06-18 吴海林 Instrument for standardizing gravimeter lattice value
CN103376471B (en) * 2012-04-13 2015-08-19 中国石油天然气集团公司 A kind of waters gravity prospecting acquisition method on ice
CN108387951B (en) * 2018-01-19 2019-11-12 中国人民解放军92859部队 A kind of new algorithm based on T1 Repeated Line Tl correction air-sea calibration of gravimeter

Also Published As

Publication number Publication date
CN110927823A (en) 2020-03-27

Similar Documents

Publication Publication Date Title
CN110927823B (en) Grid value correction method applied to shipborne gravimeter
US6543146B2 (en) Electronic compass and compensation of large magnetic errors for operation over all orientations
US10215564B2 (en) Automatic compass calibration systems and methods
Hare Depth and position error budgets for mulitbeam echosounding
CN100535683C (en) Calibration method for ultra-short baseline acoustic positioning system
CN109085655B (en) Underwater platform gravity measurement scheme and verification method
CN110220491B (en) Method for estimating installation error angle of optical pod of unmanned aerial vehicle
CN108802737A (en) A kind of ultra short base line and calibration method for shallow water area Underwater Navigation
CN103017755A (en) Measuring method for underwater navigation attitudes
CN110132308B (en) Attitude determination-based USBL installation error angle calibration method
CN103925904B (en) A kind of ultra-short baseline setting angle deviation unbiased esti-mator method based on symmetrical survey line
Mohammadloo et al. Correcting multibeam echosounder bathymetric measurements for errors induced by inaccurate water column sound speeds
Roggenbuck et al. Sea surface heights retrieval from ship-based measurements assisted by GNSS signal reflections
EP0907898A2 (en) Method of locating hydrophones
CN110109167B (en) Offshore precision positioning method based on elevation constraint
CN104965102A (en) Method of effectively improving flow direction measurement precision of acoustic doppler current profiler
CN104567802A (en) Survey line land-sea elevation transfer method employing integrated shipborne gravity and GNSS
CN113758469A (en) Sea level monitoring method and system based on multimode multi-frequency GNSS receiver
Bjorke Computation of calibration parameters for multibeam echo sounders using the least squares method
CN113819882B (en) Method for calculating gravity potential difference between cross-sea elevation points
CN211977915U (en) Real-time tide level determining device
CN114994649A (en) Multi-beam probe correction method and multi-beam sounding data correction method
Zhang et al. Experimental study on underwater moving gravity measurement by using strapdown gravimeter based on AUV platform
Ardalan et al. Improved vessel squat modeling for hydrographic and navigation applications using kinematic GNSS positioning
Riley et al. An accuracy statement for the buoy heading component of NDBC directional wave measurements

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