[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN105388512B - A kind of calibration method and device of earthquake data acquisition exploration time - Google Patents

A kind of calibration method and device of earthquake data acquisition exploration time Download PDF

Info

Publication number
CN105388512B
CN105388512B CN201510705641.1A CN201510705641A CN105388512B CN 105388512 B CN105388512 B CN 105388512B CN 201510705641 A CN201510705641 A CN 201510705641A CN 105388512 B CN105388512 B CN 105388512B
Authority
CN
China
Prior art keywords
time
current
crystal oscillator
local crystal
count value
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
CN201510705641.1A
Other languages
Chinese (zh)
Other versions
CN105388512A (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.)
China National Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
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 China National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN201510705641.1A priority Critical patent/CN105388512B/en
Publication of CN105388512A publication Critical patent/CN105388512A/en
Application granted granted Critical
Publication of CN105388512B publication Critical patent/CN105388512B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/003Seismic data acquisition in general, e.g. survey design

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Electric Clocks (AREA)

Abstract

The application provides a kind of calibration method and device of earthquake data acquisition exploration time.This method includes:The first current time of satellite time service system is obtained as the local clock very first time, initializes local crystal oscillator count value;When reaching default prover time, judge whether satellite-signal is effective;When satellite-signal is invalid, local crystal oscillator counts, and is acquired time mark to geological data using count value;In default prover time and effective satellite-signal, the current time of time dissemination system second is obtained, local crystal oscillator current count value is recorded, determines the current local crystal oscillator time of current count value;The single counting error time is calculated according to the first current time, the very first time, the second current time, current local crystal oscillator time and current count value;Using the single counting error time is carried out to the local crystal oscillator time of geological data to be calibrated time calibration.The technical scheme provided using the embodiment of the present application can accurately determine the time of earthquake data acquisition exploration.

Description

A kind of calibration method and device of earthquake data acquisition exploration time
Technical field
The present invention relates to technical field of geophysical exploration, more particularly to the calibration of earthquake data acquisition exploration time a kind of Method and apparatus.
Background technology
In seismic exploration, it is desirable to must keep strict time consistency at epicenter excitation moment and data acquisition moment Property, typically require that error must not exceed a sampling interval.In recent years, Service of Timing by Satellite is progressively dissolved into earthquake data acquisition In exploration process, satellite time service system can provide fiducial time for earthquake data acquisition equipment.
When the earthquake data acquisition equipment with satellite time transfer function is mostly using the standard from satellite in the prior art Between the working method directly calibrated to local clock of information, reach the purpose of synchronous earthquake data acquisition exploration time.But Ground observation construction environment is complicated, and satellite-signal can not be received or satellite-signal is unstable by occurring in some construction areas Phenomenon.This, which may result in earthquake data acquisition equipment, can not carry out the calibration of local clock in time, to earthquake data markers Whether standard time accurately will depend entirely on local crystal oscillator precision, and high-precision local crystal oscillator cost is high, and works as long-time nothing After method receives the standard time information from satellite time service system, earthquake-capturing equipment will be unable to normal work.
Therefore, need a kind of calibration method of earthquake data acquisition exploration time badly in the prior art, can solve at low cost The problem of certainly earthquake data acquisition equipment is high to satellite-signal dependence.
The content of the invention
The purpose of the application is to provide a kind of calibration method and device of earthquake data acquisition exploration time, can be accurate Determination earthquake data acquisition exploration time.
What the calibration method and device for the earthquake data acquisition exploration time that the application provides were realized in:
A kind of calibration method of earthquake data acquisition exploration time, methods described include:
Obtain satellite time service system the first current time, using first current time as local clock first when Between, initialize the count value of local crystal oscillator;
When reaching default prover time, judge whether satellite-signal is effective;
When the result judged is invalid as satellite-signal, the local crystal oscillator is counted, and utilizes the local crystalline substance The count value shaken is acquired time mark to geological data;
When reaching the default prover time and effective satellite-signal, obtain the time dissemination system second is current Time, the current count value of the local crystal oscillator is recorded, determine the current local crystal oscillator time corresponding to the current count value;
During according to first current time, the very first time, second current time, the current local crystal oscillator Between and the current count value single counting error time is calculated;
The local crystal oscillator time of geological data to be calibrated is determined, utilizes single counting error time and described to be calibrated Count value corresponding to geological data carries out time calibration to the local crystal oscillator time of the geological data to be calibrated, obtains described treat Calibrate the collection exploration time after geological data calibration.
In a preferred embodiment, it is described after judging whether satellite-signal is effective when reaching default prover time Method also includes:
When the result judged is effective as satellite-signal, the current time of satellite time service system is obtained, described will be worked as First current time of the preceding time as the satellite time service system, first using the current time as the local clock Time, initialize the count value of the local crystal oscillator.
In a preferred embodiment, it is described to determine current local crystal oscillator time bag corresponding to the current count value Include:
During according to the first of the current count value of the local crystal oscillator, the frequency of the local crystal oscillator and the local clock Between the current local crystal oscillator time is calculated.
In a preferred embodiment, it is described to be worked as according to first current time, the very first time, described second The calculation formula of single counting error time is calculated in preceding time, the current local crystal oscillator time and the current count value It is as follows:
In above formula, t0Represent the single counting error time;ts2Represent second current time;ts1Represent described One current time;t'2Represent the current local crystal oscillator time;t1Represent the very first time;Count1Represent the current meter Numerical value.
In a preferred embodiment, the local crystal oscillator time for determining geological data to be calibrated includes:
According to count value corresponding to the geological data to be calibrated, the frequency of the local crystal oscillator and the local clock The local crystal oscillator time of geological data to be calibrated is calculated in the very first time.
In a preferred embodiment, it is described to utilize the single counting error time and the geological data to be calibrated Corresponding count value carries out time calibration to the local crystal oscillator time of the geological data to be calibrated, obtains the earthquake to be calibrated The calculation formula of collection exploration time after data calibration is as follows:
T=t0*Countd+td
In above formula, t represents the collection exploration time of the geological data to be calibrated;t0When representing the single counting error Between;Count1Represent that the geological data to be calibrated corresponds to count value, tdWhen representing the local crystal oscillator of the geological data to be calibrated Between.
A kind of calibrating installation of earthquake data acquisition exploration time, described device include:
First data processing module, for obtaining the first current time of satellite time service system, by described first it is current when Between very first time as local clock, initialize the count value of local crystal oscillator;
Judge module, during for reaching default prover time, judge whether satellite-signal is effective;
Second data processing module, for when the result judged is invalid as satellite-signal, the local crystal oscillator to enter Row is counted, and time mark is acquired to geological data using the count value of the local crystal oscillator;
3rd data processing module, for when reaching the default prover time and effective satellite-signal, obtaining Second current time of the time dissemination system, the current count value of the local crystal oscillator is recorded, determines the current count value pair The current local crystal oscillator time answered;
Computing module, for according to the first current time, very first time, second current time, described The single counting error time is calculated in current local crystal oscillator time and the current count value;
4th data processing module, for determining the local crystal oscillator time of geological data to be calibrated, utilize the single meter Local crystal oscillator time of the count value to the geological data to be calibrated corresponding to number error time and the geological data to be calibrated Time calibration is carried out, obtains the collection exploration time after the geological data calibration to be calibrated.
In a preferred embodiment, after judge module judges whether satellite-signal is effective, described device is also wrapped Include:
5th data processing module, for when the result judged is effective as satellite-signal, obtaining satellite time transfer system The current time of system, the first current time using the current time as the satellite time service system, by the current time As the very first time of the local clock, the count value of the initialization local crystal oscillator.
In a preferred embodiment, the 3rd data processing module includes:
First computing unit, for the current count value according to the local crystal oscillator, the frequency of the local crystal oscillator and institute The current local crystal oscillator time is calculated in the very first time for stating local clock.
In a preferred embodiment, the calculation formula of the computing module is as follows:
In above formula, t0Represent the single counting error time;ts2Represent second current time;ts1Represent described One current time;t'2Represent the current local crystal oscillator time;t1Represent the very first time;Count1Represent the current meter Numerical value.
In a preferred embodiment, the 4th data processing module includes:
Second computing unit, for according to corresponding to the geological data to be calibrated count value, the local crystal oscillator frequency The local crystal oscillator time of geological data to be calibrated is calculated in the very first time of rate and the local clock.
In a preferred embodiment, the calculation formula of the 4th data processing module is as follows:
T=t0*Countd+td
In above formula, t represents the collection exploration time of the geological data to be calibrated;t0When representing the single counting error Between;Count1Represent that the geological data to be calibrated corresponds to count value, tdWhen representing the local crystal oscillator of the geological data to be calibrated Between.
The technical scheme provided from above the embodiment of the present application, the embodiment of the present application connect in earthquake data acquisition equipment During continuous acquisition and recording, when satellite-signal can be received, can use the temporal information of satellite time service system to it is local when When clock is calibrated., can be in default prover time, using local crystal oscillator count value over the ground when satellite-signal can not be received Shake data and carry out interim acquisition time mark.Then in default prover time and when can receive satellite-signal, can use , can be according to a current and preceding satellite time transfer system while temporal information of satellite time service system is calibrated to local clock The single counting error time is calculated in the temporal information of system and corresponding local crystal oscillator count value etc..Finally, can utilize The geological data to be calibrated gathered when the single counting error time and invalid satellite-signal corresponds to count value to be calibrated The local crystal oscillator time for shaking data is carried out time calibration, the collection exploration time of the geological data to be calibrated after being calibrated.With Prior art is compared, can the inexpensive time accuracy for improving earthquake data acquisition exploration.
Brief description of the drawings
, below will be to embodiment or existing in order to illustrate more clearly of the embodiment of the present application or technical scheme of the prior art There is the required accompanying drawing used in technology description to be briefly described, it should be apparent that, drawings in the following description are only this Some embodiments described in application, for those of ordinary skill in the art, do not paying the premise of creative labor Under, other accompanying drawings can also be obtained according to these accompanying drawings.
Fig. 1 is a kind of flow chart of embodiment of the calibration method for the earthquake data acquisition exploration time that the application provides
Fig. 2 is the flow of another embodiment of the calibration method for the earthquake data acquisition exploration time that the application provides Figure;
Fig. 3 is a kind of schematic diagram of the calibrating installation for the earthquake data acquisition exploration time that the application provides;
Fig. 4 is another schematic diagram of the calibrating installation for the earthquake data acquisition exploration time that the application provides.
Embodiment
In order that those skilled in the art more fully understand the technical scheme in the application, it is real below in conjunction with the application The accompanying drawing in example is applied, the technical scheme in the embodiment of the present application is clearly and completely described, it is clear that described implementation Example only some embodiments of the present application, rather than whole embodiments.It is common based on the embodiment in the application, this area The every other embodiment that technical staff is obtained under the premise of creative work is not made, it should all belong to the application protection Scope.
Describe the specific implementation of the embodiment of the present application in detail with several specific examples below.
Introduce a kind of a kind of embodiment of the calibration method of earthquake data acquisition exploration time of the application first below.Fig. 1 It is a kind of flow chart of embodiment of the calibration method for the earthquake data acquisition exploration time that the application provides, this application provides Method operating procedure as described in embodiment or flow chart, but based on it is conventional or without performing creative labour can include it is more or The less operating procedure of person.The step of being enumerated in embodiment order is only a kind of mode in numerous step execution sequences, no Represent unique execution sequence., can be according to embodiment or accompanying drawing institute when device or client production in practice performs The method order shown performs either performs (such as environment of parallel processor or multiple threads) parallel.Specific such as Fig. 1 Shown, methods described can include:
S110:The first current time of satellite time service system is obtained, using first current time as local clock The very first time, initialize the count value of local crystal oscillator.
In this application, the first current time of satellite time service system can be obtained, using first current time as The very first time of local clock, initialize the count value of local crystal oscillator.
In actual applications, the satellite time service system can include but is not limited to satellite time service system GPS or the Big Dipper is defended Star time dissemination system.
In actual applications, the local crystal oscillator can include but is not limited to VCXO.
S120:When reaching default prover time, judge whether satellite-signal is effective.
In this application, when can reach default prover time, judge whether satellite-signal is effective.Specifically, in reality Can be previously according to the parameter setting calibration cycles such as the precision of the local crystal oscillator of selection, the default prover time in It can be determined according to the calibration cycle pre-set.For example the calibration cycle can be 10S, accordingly, reached every 10S During the default prover time, it can be determined that whether satellite-signal is effective.
S130:When the result judged is invalid as satellite-signal, the local crystal oscillator is counted, and utilizes described The count value of ground crystal oscillator is acquired time mark to geological data.
In this application, when the result that step S120 judges is invalid as satellite-signal, the local crystal oscillator can start Counted, and time mark is acquired to geological data using the count value of the local crystal oscillator.Specifically, the local The count value of crystal oscillator can be represented by the way of number and current count value is overflowed.
S140:When reaching the default prover time and effective satellite-signal, the of the time dissemination system is obtained Two current times, the current count value of the local crystal oscillator is recorded, determine current local crystal oscillator corresponding to the current count value Time.
In this application, after step s 130, can have in the arrival default prover time and the satellite-signal During effect, the second current time of the time dissemination system is obtained, records the current count value of the local crystal oscillator, is determined described current The current local crystal oscillator time corresponding to count value.
Specifically, described determine that the current local crystal oscillator time can include corresponding to the current count value:
During according to the first of the current count value of the local crystal oscillator, the frequency of the local crystal oscillator and the local clock Between the current local crystal oscillator time is calculated.
It is specifically, described according to the current count value of the local crystal oscillator, the frequency of the local crystal oscillator and the local The very first time of clock be calculated the current local crystal oscillator time according to calculate formula can be with as follows:
In above formula, t'2Represent the current local crystal oscillator time;ts1Represent first current time;Count1Represent The current count value;F represents the frequency of the local crystal oscillator.
Further, the frequency of the local crystal oscillator can use the temperature parameter of the local crystal oscillator to be adjusted.
S150:According to first current time, the very first time, second current time, the current local The single counting error time is calculated in crystal oscillator time and the current count value.
In this application, after step s 140, can be according to the first current time, very first time, described The single counting error time is calculated in second current time, the current local crystal oscillator time and the current count value.
Specifically, it is described according to first current time, the very first time, second current time, it is described work as The calculation formula that the single counting error time is calculated in preceding local crystal oscillator time and the current count value can be as follows:
In above formula, t0Represent the single counting error time;ts2Represent second current time;ts1Represent described One current time;t'2Represent the current local crystal oscillator time;t1Represent the very first time;Count1Represent the current meter Numerical value.
S160:The local crystal oscillator time of geological data to be calibrated is determined, utilizes single counting error time and described Count value corresponding to geological data to be calibrated carries out time calibration to the local crystal oscillator time of the geological data to be calibrated, obtains The collection exploration time after the geological data calibration to be calibrated.
In this application, after the step s 150, it may be determined that the local crystal oscillator time of geological data to be calibrated, utilize institute Count value corresponding to single counting error time and the geological data to be calibrated is stated to the local of the geological data to be calibrated The crystal oscillator time carries out time calibration, obtains the collection exploration time after the geological data calibration to be calibrated.
Specifically, the local crystal oscillator time for determining geological data to be calibrated can include:
According to count value corresponding to the geological data to be calibrated, the frequency of the local crystal oscillator and the local clock The local crystal oscillator time of geological data to be calibrated is calculated in the very first time.
Specifically, the count value according to corresponding to the geological data to be calibrated, the frequency of the local crystal oscillator and the institute The calculation formula for the local crystal oscillator time that geological data to be calibrated is calculated in the very first time for stating local clock can be as follows:
In above formula, tdRepresent the local crystal oscillator time of geological data to be calibrated;ts1Represent first current time; CountdRepresent count value corresponding to the geological data to be calibrated;F represents the frequency of the local crystal oscillator.
Further, the frequency of the local crystal oscillator can use the temperature parameter of the local crystal oscillator to be adjusted.
Specifically, described utilize count value pair corresponding to the single counting error time and the geological data to be calibrated The local crystal oscillator time of the geological data to be calibrated carries out time calibration, obtains adopting after the geological data calibration to be calibrated The calculation formula of collection exploration time can be as follows:
T=t0*Countd+td
In above formula, t represents the collection exploration time of the geological data to be calibrated;t0When representing the single counting error Between;Count1Represent that the geological data to be calibrated corresponds to count value, tdWhen representing the local crystal oscillator of the geological data to be calibrated Between.
Fig. 2 is the flow of another embodiment of the calibration method for the earthquake data acquisition exploration time that the application provides Figure, in certain embodiments, after step S120, as shown in Fig. 2 methods described also includes:
S170:When the result judged is effective as satellite-signal, the current time of satellite time service system is obtained, by institute First current time of the current time as the satellite time service system is stated, using the current time as the local clock The very first time, initialize the count value of the local crystal oscillator.
As can be seen here, the technical side that a kind of embodiment of the calibration method of earthquake data acquisition exploration time of the application provides Case, when that can receive satellite-signal, can use satellite time transfer in earthquake data acquisition equipment continuous acquisition recording process When the temporal information of system is calibrated to local clock;When satellite-signal can not be received, can in default prover time, Interim acquisition time mark is carried out to geological data using local crystal oscillator count value, then in default prover time and can be received , can basis while the temporal information of satellite time service system being used to calibrate local clock during to satellite-signal Single counting is calculated in current and a preceding satellite time service system temporal information and corresponding local crystal oscillator count value etc. Error time, finally, the earthquake number to be calibrated gathered when the single counting error time and invalid satellite-signal can be utilized Time calibration, the earthquake to be calibrated after being calibrated are carried out to the local crystal oscillator time of geological data to be calibrated according to corresponding count value The collection exploration time of data.Compared with prior art, can the inexpensive time accuracy for improving earthquake data acquisition exploration.
Fig. 3 is a kind of schematic diagram of the calibrating installation for the earthquake data acquisition exploration time that the application provides, and the application is another On the one hand a kind of calibrating installation of earthquake data acquisition exploration time is also provided, as shown in figure 3, described device 300 can include:
First data processing module 310, it can be used for the first current time for obtaining satellite time service system, by described first The very first time of the current time as local clock, initialize the count value of local crystal oscillator;
Judge module 320, when can be used for reaching default prover time, judge whether satellite-signal is effective;
Second data processing module 330, it can be used for when the result judged is invalid as satellite-signal, the local Crystal oscillator is counted, and time mark is acquired to geological data using the count value of the local crystal oscillator;
3rd data processing module 340, can be used for reach the default prover time and the satellite-signal it is effective When, the second current time of the time dissemination system is obtained, the current count value of the local crystal oscillator is recorded, determines the current meter The current local crystal oscillator time corresponding to numerical value;
Computing module 350, can be used for according to first current time, the very first time, it is described second it is current when Between, the current local crystal oscillator time and the current count value single counting error time is calculated;
4th data processing module 360, the local crystal oscillator time of geological data to be calibrated is determined for, using described Count value corresponding to single counting error time and the geological data to be calibrated is to the local brilliant of the geological data to be calibrated Shake progress time calibration time, obtain the collection exploration time after the geological data calibration to be calibrated.
Fig. 4 is another schematic diagram of the calibrating installation for the earthquake data acquisition exploration time that the application provides, at one In preferred embodiment, as shown in figure 4, judge module judge satellite-signal whether effectively after, described device 300 can be with Including:
It 5th data processing module 370, can be used for when the result judged is effective as satellite-signal, obtain satellite The current time of time dissemination system, the first current time using the current time as the satellite time service system, described it will work as The very first time of the preceding time as the local clock, initialize the count value of the local crystal oscillator.
In a preferred embodiment, the 3rd data processing module 340 can include:
First computing unit, it can be used for the current count value according to the local crystal oscillator, the frequency of the local crystal oscillator The current local crystal oscillator time is calculated with the very first time of the local clock.
In a preferred embodiment, the calculation formula of the computing module 350 can be as follows:
In above formula, t0Represent the single counting error time;ts2Represent second current time;ts1Represent described One current time;t'2Represent the current local crystal oscillator time;t1Represent the very first time;Count1Represent the current meter Numerical value.
In a preferred embodiment, the 4th data processing module 360 can include:
Second computing unit, for according to corresponding to the geological data to be calibrated count value, the local crystal oscillator frequency The local crystal oscillator time of geological data to be calibrated is calculated in the very first time of rate and the local clock.
In a preferred embodiment, the calculation formula of the 4th data processing module 360 can be as follows:
T=t0*Countd+td
In above formula, t represents the collection exploration time of the geological data to be calibrated;t0When representing the single counting error Between;Count1Represent that the geological data to be calibrated corresponds to count value, tdWhen representing the local crystal oscillator of the geological data to be calibrated Between.
As can be seen here, the calibration method of the application earthquake data acquisition exploration time a kind of and the embodiment of device provide Technical scheme is in earthquake data acquisition equipment continuous acquisition recording process, when satellite-signal can be received, can use defend When the temporal information of star time dissemination system is calibrated to local clock;, can be in default school when satellite-signal can not be received Between punctual, interim acquisition time mark is carried out to geological data using local crystal oscillator count value, then in default prover time and When can receive satellite-signal, while the temporal information of satellite time service system being used to calibrate local clock, It can be calculated according to the temporal information of a current and preceding satellite time service system and corresponding local crystal oscillator count value etc. The single counting error time, finally, it can utilize what is gathered when the single counting error time and invalid satellite-signal to treat school Quasi- geological data corresponds to count value and carries out time calibration to the local crystal oscillator time of geological data to be calibrated, is treated after being calibrated Calibrate the collection exploration time of geological data.Compared with prior art, can low cost improve earthquake data acquisition exploration when Between accuracy.
Each embodiment in this specification is described by the way of progressive, what each embodiment stressed be with The difference of other embodiment, between each embodiment identical similar part mutually referring to.It is real especially for system For applying example, because it is substantially similar to embodiment of the method, so description is fairly simple, related part is referring to embodiment of the method Part explanation.
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application have it is many deformation and Change is without departing from spirit herein, it is desirable to which appended claim includes these deformations and changed without departing from the application's Spirit.

Claims (10)

1. a kind of calibration method of earthquake data acquisition exploration time, it is characterised in that methods described includes:
First current time of acquisition satellite time service system, the very first time using first current time as local clock, Initialize the count value of local crystal oscillator;
When reaching default prover time, judge whether satellite-signal is effective;
When the result judged is invalid as satellite-signal, the local crystal oscillator is counted, and utilizes the local crystal oscillator Count value is acquired time mark to geological data;
When reaching the default prover time and effective satellite-signal, obtain the time dissemination system second it is current when Between, the current count value of the local crystal oscillator is recorded, determines the current local crystal oscillator time corresponding to the current count value;
According to first current time, the very first time, second current time, the current local crystal oscillator time and The single counting error time is calculated in the current count value;
The local crystal oscillator time of geological data to be calibrated is determined, utilizes the single counting error time and the earthquake to be calibrated Count value corresponding to data carries out time calibration to the local crystal oscillator time of the geological data to be calibrated, obtains described to be calibrated The collection exploration time after geological data calibration;
Wherein, it is described according to first current time, the very first time, second current time, the current local The calculation formula that the single counting error time is calculated in crystal oscillator time and the current count value is as follows:
<mrow> <msub> <mi>t</mi> <mn>0</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mrow> <mi>s</mi> <mn>2</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>t</mi> <mrow> <mi>s</mi> <mn>1</mn> </mrow> </msub> <mo>)</mo> <mo>-</mo> <mo>(</mo> <msubsup> <mi>t</mi> <mn>2</mn> <mo>&amp;prime;</mo> </msubsup> <mo>-</mo> <msub> <mi>t</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mrow> <msub> <mi>Count</mi> <mn>1</mn> </msub> </mrow> </mfrac> </mrow>
In above formula, t0Represent the single counting error time;ts2Represent second current time;ts1Represent that described first works as The preceding time;t'2Represent the current local crystal oscillator time;t1Represent the very first time;Count1Represent the current count Value.
2. according to the method for claim 1, it is characterised in that when reaching default prover time, whether judge satellite-signal After effectively, methods described also includes:
When the result judged is effective as satellite-signal, the current time of satellite time service system is obtained, when will be described current Between the first current time as the satellite time service system, using the current time as the local clock first when Between, the count value of the initialization local crystal oscillator.
3. method according to claim 1 or 2, it is characterised in that described to determine corresponding to the current count value currently The local crystal oscillator time includes:
Counted according to the very first time of the current count value of the local crystal oscillator, the frequency of the local crystal oscillator and the local clock Calculation obtains the current local crystal oscillator time.
4. method according to claim 1 or 2, it is characterised in that the local crystal oscillator for determining geological data to be calibrated Time includes:
According to the first of count value corresponding to the geological data to be calibrated, the frequency of the local crystal oscillator and the local clock The local crystal oscillator time of geological data to be calibrated is calculated in time.
5. method according to claim 1 or 2, it is characterised in that described to utilize the single counting error time and institute State count value corresponding to geological data to be calibrated and time calibration is carried out to the local crystal oscillator time of the geological data to be calibrated, obtain The calculation formula of collection exploration time after to the geological data calibration to be calibrated is as follows:
T=t0*Countd+td
In above formula, t represents the collection exploration time of the geological data to be calibrated;t0Represent the single counting error time; CountdRepresent that the geological data to be calibrated corresponds to count value, tdRepresent the local crystal oscillator time of the geological data to be calibrated.
6. a kind of calibrating installation of earthquake data acquisition exploration time, it is characterised in that described device includes:
First data processing module, for obtaining the first current time of satellite time service system, first current time is made For the very first time of local clock, the count value of local crystal oscillator is initialized;
Judge module, during for reaching default prover time, judge whether satellite-signal is effective;
Second data processing module, for when the result judged is invalid as satellite-signal, the local crystal oscillator to be counted Number, time mark is acquired to geological data using the count value of the local crystal oscillator;
3rd data processing module, for when reaching the default prover time and effective satellite-signal, described in acquisition Second current time of time dissemination system, the current count value of the local crystal oscillator is recorded, is determined corresponding to the current count value The current local crystal oscillator time;
Computing module, for according to the first current time, very first time, second current time, described current The single counting error time is calculated in local crystal oscillator time and the current count value;
4th data processing module, for determining the local crystal oscillator time of geological data to be calibrated, counted and missed using the single Count value corresponding to poor time and the geological data to be calibrated is carried out to the local crystal oscillator time of the geological data to be calibrated Time calibration, obtain the collection exploration time after the geological data calibration to be calibrated;
Wherein, the computing module be calculated the single counting error time use calculation formula it is as follows:
<mrow> <msub> <mi>t</mi> <mn>0</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mrow> <mi>s</mi> <mn>2</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>t</mi> <mrow> <mi>s</mi> <mn>1</mn> </mrow> </msub> <mo>)</mo> <mo>-</mo> <mo>(</mo> <msubsup> <mi>t</mi> <mn>2</mn> <mo>&amp;prime;</mo> </msubsup> <mo>-</mo> <msub> <mi>t</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mrow> <msub> <mi>Count</mi> <mn>1</mn> </msub> </mrow> </mfrac> </mrow>
In above formula, t0Represent the single counting error time;ts2Represent second current time;ts1Represent that described first works as The preceding time;t'2Represent the current local crystal oscillator time;t1Represent the very first time;Count1Represent the current count Value.
7. device according to claim 6, it is characterised in that after judge module judges whether satellite-signal is effective, Described device also includes:
5th data processing module, for when the result judged is effective as satellite-signal, obtaining satellite time service system Current time, the first current time using the current time as the satellite time service system, using the current time as The very first time of the local clock, initialize the count value of the local crystal oscillator.
8. the device according to claim 6 or 7, it is characterised in that the 3rd data processing module includes:
First computing unit, for the current count value according to the local crystal oscillator, the frequency of the local crystal oscillator and described The current local crystal oscillator time is calculated in the very first time of ground clock.
9. the device according to claim 6 or 7, it is characterised in that the 4th data processing module includes:
Second computing unit, for according to corresponding to the geological data to be calibrated count value, the frequency of the local crystal oscillator and The local crystal oscillator time of geological data to be calibrated is calculated in the very first time of the local clock.
10. the device according to claim 6 or 7, it is characterised in that the calculation formula of the 4th data processing module is such as Under:
T=t0*Countd+td
In above formula, t represents the collection exploration time of the geological data to be calibrated;t0Represent the single counting error time; CountdRepresent that the geological data to be calibrated corresponds to count value, tdRepresent the local crystal oscillator time of the geological data to be calibrated.
CN201510705641.1A 2015-10-27 2015-10-27 A kind of calibration method and device of earthquake data acquisition exploration time Active CN105388512B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510705641.1A CN105388512B (en) 2015-10-27 2015-10-27 A kind of calibration method and device of earthquake data acquisition exploration time

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510705641.1A CN105388512B (en) 2015-10-27 2015-10-27 A kind of calibration method and device of earthquake data acquisition exploration time

Publications (2)

Publication Number Publication Date
CN105388512A CN105388512A (en) 2016-03-09
CN105388512B true CN105388512B (en) 2018-01-05

Family

ID=55420993

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510705641.1A Active CN105388512B (en) 2015-10-27 2015-10-27 A kind of calibration method and device of earthquake data acquisition exploration time

Country Status (1)

Country Link
CN (1) CN105388512B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3580586B1 (en) 2017-02-09 2024-10-23 Services Pétroliers Schlumberger Geophysical deep learning
CN112564693B (en) * 2020-12-18 2024-01-05 北京自动化控制设备研究所 Self-adaptive time keeping time service method
CN115314142B (en) * 2022-06-29 2023-04-28 阿波罗智能技术(北京)有限公司 Time calibration method, device, electronic equipment and storage medium
CN118131308A (en) * 2022-12-01 2024-06-04 中国石油化工股份有限公司 Sampling data management method and device suitable for seismograph, seismograph and storage medium

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1294471A (en) * 1999-10-20 2001-05-09 索尼国际(欧洲)股份有限公司 Mobile terminal with accuracy real time generator for radio telecommunicating system
CN101655547A (en) * 2008-08-08 2010-02-24 凹凸电子(武汉)有限公司 Satellite navigation receiver and method for locating by utilizing satellite navigation receiver
CN102938220A (en) * 2012-11-28 2013-02-20 天津七一二通信广播有限公司 Method for correcting clock of automatic ship identifying terminal equipment
CN103532652A (en) * 2013-10-18 2014-01-22 杭州华三通信技术有限公司 Time synchronizing device and method
CN104460311A (en) * 2014-12-30 2015-03-25 四川九洲电器集团有限责任公司 Time calibration method and device

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050047275A1 (en) * 2003-09-01 2005-03-03 Geo-X Systems, Ltd. Synchronization and positioning of seismic data acquisition systems

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1294471A (en) * 1999-10-20 2001-05-09 索尼国际(欧洲)股份有限公司 Mobile terminal with accuracy real time generator for radio telecommunicating system
CN101655547A (en) * 2008-08-08 2010-02-24 凹凸电子(武汉)有限公司 Satellite navigation receiver and method for locating by utilizing satellite navigation receiver
CN102938220A (en) * 2012-11-28 2013-02-20 天津七一二通信广播有限公司 Method for correcting clock of automatic ship identifying terminal equipment
CN103532652A (en) * 2013-10-18 2014-01-22 杭州华三通信技术有限公司 Time synchronizing device and method
CN104460311A (en) * 2014-12-30 2015-03-25 四川九洲电器集团有限责任公司 Time calibration method and device

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
GPS 时钟同步法地震勘探仪器的设计与实现;易碧金;《石油仪器》;20050228;第19卷(第1期);第5-8页 *
地震数据采集无线同步技术研究;任家富等;《中国测试技术》;20081130;第34卷(第6期);第1-3页 *
时间统一技术研究及应用;刘凯;《中国优秀硕士学位论文全文数据库 信息科技辑》;20101115(第11期);第21-29页、53-54页 *

Also Published As

Publication number Publication date
CN105388512A (en) 2016-03-09

Similar Documents

Publication Publication Date Title
CN105388512B (en) A kind of calibration method and device of earthquake data acquisition exploration time
CN105444678B (en) Handset size measuring method and system
Samushia et al. Interpreting large-scale redshift-space distortion measurements
Scognamiglio et al. Fast determination of moment tensors and rupture history: What has been learned from the 6 April 2009 L'Aquila earthquake sequence
CN105204065B (en) A kind of method and apparatus of first break picking ripple
CN106097361A (en) A kind of defective area detection method and device
CN106896337A (en) For the method for Magnetic Sensor calibration
CN105549083B (en) A kind of method and apparatus for determining seismic observation system repeatability and integrally measuring
CN106324665A (en) Method and system of inverting fracture density
CN103499834B (en) The method and apparatus recovering analog detector geological data low frequency signal
Xu et al. Systemic comparison of seismometer horizontal orientations based on teleseismic earthquakes and ambient‐noise data
Peng et al. Application of a Threshold‐Based Earthquake Early Warning Method to the M w 6.6 Lushan Earthquake, Sichuan, China
CN106324671A (en) Method and device for checking shot point offset
CN106842302B (en) Method and device for editing first arrivals in batches
CN105547237B (en) A kind of shadow location technology based on least square method
CN108267782B (en) Earth&#39;s surface absolute displacement test macro and method based on machine vision and numeric value analysis
Telesca et al. Visibility graph analysis of seismicity around Enguri high arch dam, Caucasus
CN104316960B (en) A kind of Forecast Means of Reservoir Fractures and system based on VSP
CN106375055A (en) Method and device for measuring clock error of network device
CN107942391A (en) A kind of seabed geophone station localization method and device
Garcia et al. A nonlinear method to estimate source parameters, amplitude, and travel times of teleseismic body waves
CN111060416A (en) Method and system for dividing surface outcrop and rock core mechanical layer
CN109188515A (en) Micro-seismic monitoring crack focal point position calculating method and system
Materna et al. Slip deficit rates on southern Cascadia faults resolved with viscoelastic earthquake cycle modeling of geodetic deformation
Kelevitz et al. Performance of High‐Rate GPS Waveforms at Long Periods: Moment Tensor Inversion of the 2003 M w 8.3 Tokachi‐Oki Earthquake

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant