CN101086534A - Demodulator probe secondary localization method - Google Patents
Demodulator probe secondary localization method Download PDFInfo
- Publication number
- CN101086534A CN101086534A CN 200610083878 CN200610083878A CN101086534A CN 101086534 A CN101086534 A CN 101086534A CN 200610083878 CN200610083878 CN 200610083878 CN 200610083878 A CN200610083878 A CN 200610083878A CN 101086534 A CN101086534 A CN 101086534A
- Authority
- CN
- China
- Prior art keywords
- geophone station
- arrival
- grid
- geophone
- sdt
- 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
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The demodulator probe secondary positioning method of offshore oil earthquake starts with the looking for a demodulator probe location. Calibration of the demodulator probe based on the data of the location's data, if it is pulled in the best position, the coordinate of the point is the actual coordinate of the demodulator probe. It improves the existing speed over 7 times, with positioning precision significantly improved. It can reduce working load of relevant personnel, saving a lot of working time.
Description
Technical field
The present invention relates to the petroleum geology exploration operating technique, is a kind of demodulator probe secondary localization method in the petroleum seismic of neritic area, beach.
Background technology
In neritic area, beach oil geology seismic prospecting construction, need the wave detector sinking be arrived the seabed with ship, generally by the terrestrial coordinate of satnav (GPS) equipment records geophone station, this process is the one-time positioning of geophone station usually during sinking.If the position of geophone station is accurate, after the road collection data of geophone station were done linear normal moveout correction in the information of seismological observation collection, the first arrival ski-jump should be evened up (as Fig. 1).When seismological observation is gathered accurately the geophone station position be the necessary condition that data is prospected in high-quality earthquake, be the condition that obtains high-quality stacked section and end result section.
In actual observation, because the enforcement state of ship, sinking operation and stormy waves, morning and evening tides, influences such as trench can cause the wave detector offset, make that the first arrival ski-jump can not be evened up (as Fig. 2) after the linear normal moveout correction, like this, the reflected signal that can not guarantee same reflection horizon when road collection data are done overlap-add procedure reduces the resolution of stacked section with superimposed, has a strong impact on the quality of final processing achievement section.So need carry out second positioning to the geophone station position.
Demodulator probe secondary is positioned with two kinds of methods at present: radiocoustic position finding method and primary wave localization method, " petroleum exploration " 2003 the 4th phases disclose " the round crossing localization method " of employing, principle such as Fig. 3.R is the physical location of the geophone station that will ask for, and sonar is tied up on wave detector during radiocoustic position finding; S1 during the primary wave location, S2, S3 are the sp locations; During radiocoustic position finding, S1, S2, S3 are the positions of radiocoustic position finding ship.
With S1, S2 is the center of circle, with speed (V1 respectively, V2) multiply by the time (T1 T2) is that radius draws circle, when V1*T1+V2*T2 greater than the absolute value of S1S2 and V1*T1-V2*T2 during less than S1S2, two circles have two intersection points, and the point that has in two intersection points is exactly the physical location of geophone station; Equally, respectively with S1, S3 is the center of circle, (V1 multiply by the time V3) that (T1 T3) is that radius draws circle with speed, when V1*T1+V3*T3 greater than the absolute value of S1S3 and V1*T1-V3*T3 during less than S1S3, two circles have two intersection points, the point that has in two intersection points is exactly the physical location of geophone station; Work as S1, S2, the position of S3 is accurate, and V1, V2, V3, T1, T2, T3 be accurately the time, and S1 and S2 intersect, and S1 and S3 are crossing to have a public intersection point, and this public intersection point is exactly the physical location of geophone station.
In addition, when the absolute value that equals S1S2 or V1*T1-V2*T2 as V1*T1+V2*T2 equals S1S2, with S1, S2 is the center of circle, and (V1 V2) multiply by the time (T1 with speed, T2) draw circle for radius, two circles will have an intersection point, and this intersection point is exactly the physical location of geophone station.Wherein: when doing the primary wave location, speed V1, V2, V3 are the comprehensive speeds of seismic wave, and time T 1, T2, T3 are the first break times of seismic wave; When doing the sonar location, speed V1, V2, V3 are the velocity of propagation of sound wave in water, and time T 1, T2, T3 are the single way time of sound wave at water transmission.
This method in actual applications, S1, S2, the position of S3 and V1, V2, V3, T1, T2, T3 is not an absolutely accurate, so, three circles generally do not have public intersection point when intersecting in twos, can not be by the position of intersecting the second positioning point of determining geophone station of 3 circles, which also to finish: 1) all definite circles in the sp location relevant with this geophone station (radiocoustic position finding ship position) are carried out computing, will judge that simultaneously which circle has intersection point, the no intersection point of circle by following two steps, and will consider that two circles intersect the time, the relative position of two circles by the computing of intersecting in twos, obtains a plurality of intersection points (point looses) to influence of the error of calculation or the like; 2) a plurality of intersection points (point looses) that obtain being done the computing of asking " center of gravity ", should " center of gravity " position be exactly the position of the second positioning point of this geophone station.
As seen, the operand of existing " round crossing localization method " is big, and the man power and material who takies is a lot, and simultaneously, when the computing of location, the relative position of two circles can exert an influence to bearing accuracy, makes it to be difficult to obtain higher second positioning.
Summary of the invention
The object of the invention provides a kind of demodulator probe secondary localization method that can improve locating speed and precision effectively.
The invention provides following technical scheme, adopt following steps:
1) adopts conventional method that acquiring seismic data is observed in the construction area, carry out the first break time of shot gather data and pick up, obtain the first break time T of each shot point to geophone station;
2) adopt conventional method to set up the velocity field of construction area, obtain the seismic wave propagation speed V of shot point to geophone station;
3) the one-time positioning point position with geophone station is the center, and area dividing is on every side become grid, the length of side 2*DS of net region, and the grid number on every limit is: m=2*int (DS/d), d are the interval of grid, DS is the deviation range of geophone station;
4) with a grid node position as the geophone station position, the road collection data of this geophone station are done the line school, and calculate first arrival and even up effect;
With a grid node position XX[j], YY[k] as the geophone station position, calculate the linear dynamic correction value in each road and the difference of this road first break time
, after this difference was exactly linear normal moveout correction, the first arrival ski-jump in this road departed from the deviation of 0 datum line;
In the formula: j=1,2 ... m+1; K=1,2 ... m+1 represents (m+1) * (m+1) individual grid node is calculated;
T[i] to be seismic wave propagate into used time of this geophone station physical location, i.e. first break time from i shot point (X[i], Y[i]);
Effect SDT[j is evened up in first arrival] [k] be the absolute value sum that the first arrival ski-jump departs from the deviation of 0 datum line after the linear normal moveout correction in each road,
Wherein: i=1,2 ... n represents all data tracks of this geophone station are calculated;
5) obtain first arrival and even up effect S DT[j] [k] (j=1,2 ... m+1; K=1,2 ... m+1) a S DT[j0 of minimum in] [k0], its XX[j0], YY[k0] be the second positioning point position coordinates of geophone station.
The invention provides following technical scheme, in the described step 3), the interval d bearing accuracy as requested of grid is determined; The deviation range DS of geophone station determines according to the bearing accuracy of one-time positioning point.
Second positioning repeating step 3 to all geophone stations of construction area) to 5) step.
More traditional " the round crossing localization method " of speed that the present invention does the second positioning computing improved more than 7 times, bearing accuracy also is significantly improved simultaneously. for example: at the about big gun more than 3300 of a bunch, have in the construction of 960 geophone stations, carry out the second positioning computing with this method and need 60 minutes, need 420 minutes and utilize traditional " round crossing localization method " to do the location computing, alleviate related personnel's working strength, saved a large amount of working times.
Description of drawings
Fig. 1 for the geophone station position when accurate behind the line school first arrival even up design sketch;
Fig. 2 for the geophone station position when inaccurate behind the line school first arrival even up design sketch;
Fig. 3 is round crossing localization method schematic diagram.
Embodiment
The present invention is based on the demodulator probe secondary localization method that the first arrival of line school is evened up, find the position of a geophone station earlier, if after with the data of this location point the road collection data of this geophone station being done linear normal moveout correction, it is best that effect is evened up in first arrival, and then the coordinate of this point is exactly the actual position coordinate of geophone station.
Concrete steps are as follows:
1) adopts conventional method that acquiring seismic data is observed in the construction area, carry out the first break time of shot gather data and pick up, obtain the first break time T of each shot point to geophone station;
2) adopt conventional method to set up the velocity field of construction area, obtain the seismic wave propagation speed V of shot point to geophone station;
3) the one-time positioning point position with geophone station is the center, and area dividing is on every side become grid, the length of side 2*DS of net region, and the grid number on every limit is: m=2*int (DS/d), d are the interval of grid, DS is the deviation range of geophone station;
4). with a grid node hypothesis on location is the geophone station position, and the road collection data of this geophone station are done the line school, and calculates first arrival and even up effect;
With each grid node position XX[j], YY[k] be assumed to be the geophone station position, calculate the linear dynamic correction value in each road and the difference of this road first break time
After this difference was exactly linear normal moveout correction, the first arrival ski-jump in this road departed from the deviation of 0 datum line, and the absolute value sum of the deviate in all roads can represent that first arrival evens up effect SDT[j] [k].
SDT[j] [k] be the absolute value sum that the first arrival ski-jump departs from the deviation of 0 datum line after the linear normal moveout correction in each road.This value is more little, and more near 0, it is good more to illustrate that effect is evened up in first arrival.
In the expression formula (1):
SDT[j] [k] middle j=1,2 ... m+1; K=1,2 ... m+1 represents (m+1) * (m+1) individual grid node is calculated;
Be i shot point (X[i], Y[i]) to the distance of this geophone station assumed position (XX[j], YY[k]);
Being seismic wave propagates into used time of this geophone station assumed position (XX[j], YY[k]), promptly linear dynamic correction value from i shot point (X[i], Y[i]);
T[i] to be seismic wave propagate into used time of this geophone station physical location, i.e. first break time from i shot point (X[i], Y[i]);
T[i] 5) first arrival of alternative line school evens up effect SDT[j] [k] (j=1,2 ... m+1; K=1,2 ... m+1), find minimum value wherein, the position coordinates of the grid node of this minimum value correspondence is exactly the demodulator probe secondary elements of a fix.
Find SDT[j] [k] (j=1,2 ... m+1; K=1,2 ... m+1) SDTmin of minimum in;
SDTmin=SDT[1][1];
SDTmin=min(SDTmin,SDT[j][k])…………(2)
In the expression formula (2): j=1,2 ... m+1; K=1,2 ... m+1
Suppose: work as j=j0, during k=k0, SDT[j] [k] be worth minimum, so, XX[j0], YY[k0] be exactly the second positioning point position coordinates of geophone station.
6) all geophone stations are repeated 3,4,5 steps, finish the second positioning work of all geophone stations.
In certain shallow sea prospection project, a bunch has 3300 big guns, has 960 geophone stations, utilizes method of the present invention as follows to the process that geophone station carries out second positioning:
1. 3300 big guns are carried out first break time and pick up, obtained the first break time of each shot point to geophone station;
2. set up the velocity field of construction area, obtain the seismic wave propagation speed of shot point to geophone station;
3. the position area dividing on every side of once putting with a geophone station becomes grid;
As: this geophone station has received the data of n shot point, and n the data track of acceptance point altogether promptly arranged;
The coordinate of each shot point is: X[i], Y[i] (i=1,2 ... n);
Each shot point to the first break time of geophone station is: T[i] (i=1,2 ... n);
Each shot point to the seismic wave propagation speed of geophone station is: V[i] (i=1,2 ... n);
The coordinate of the one-time positioning point of this geophone station is: (X0, Y0);
The positioning error that requires is 1m;
The deviation range of geophone station is 50m, and promptly physical location and deviation once are in 50m, and (X0 is that grid is divided at the center Y0) with the one-time positioning point of this geophone station; The scope of net region is the Probability Area at the physical location place of this geophone station, so, the length of side 100m of net region, grid be spaced apart 1, so the grid number on every limit is: 100; Obtain the individual grid node of (101) * (101), XX[j], YY[k] (j=1,2 ... 101; K=1,2 ... 101);
4. be the geophone station position with a grid node hypothesis on location, the road collection data of this geophone station are done the line school and calculated first arrival and even up effect;
With each grid node position XX[j], YY[k] be assumed to be the geophone station position, road collection data are done the line school, obtain first arrival and even up effect SDT[j] [k], SDT[j] [k] be the absolute value sum that the first arrival ski-jump departs from the deviation of 0 datum line after the linear normal moveout correction in each road.This value is more little, and more near 0, it is good more to illustrate that effect is evened up in first arrival.
In the expression formula: j=1,2 ... 101; K=1,2 ... 101 expressions are calculated 101 * 101 grid nodes;
5. effect SDT[j is evened up in the first arrival of alternative line school] [k] (j=1,2 ... 101; K=1,2 ... 101), find minimum value wherein, the position coordinates of the grid node of this minimum value correspondence is exactly the demodulator probe secondary elements of a fix.
Find SDT[j] [k] (j=1,2 ... 101; K=1,2 ... 101) SDTmin of minimum in;
SDTmin=SDT[1][1];
SDTmin=min(SDTmin,SDT[j][k])
In the expression formula: j=1,2 ... 101; K=1,2 ... 101
Suppose: work as j=j0, during k=k0, SDT[j] [k] be worth minimum, so, XX[j0], YY[k0] be exactly the second positioning point position coordinates of geophone station.
6. all geophone stations are repeated 3,4,5 steps, finish the second positioning work of all geophone stations.
Claims (3)
1, a kind of demodulator probe secondary localization method is characterized in that adopting following steps:
1) adopts conventional method that acquiring seismic data is observed in the construction area, carry out the first break time of shot gather data and pick up, obtain the first break time T of each shot point to geophone station;
2) adopt conventional method to set up the velocity field of construction area, obtain the seismic wave propagation speed V of shot point to geophone station;
3) the one-time positioning point position with geophone station is the center, and area dividing is on every side become grid, the length of side 2*DS of net region, and the grid number on every limit is: m=2*int (DS/d), d are the interval of grid, DS is the deviation range of geophone station;
4) with a grid node position as the geophone station position, the road collection data of this geophone station are done the line school, and calculate first arrival and even up effect;
With a grid node position XX[j], YY[k] as the geophone station position, calculate the linear dynamic correction value in each road and the difference of this road first break time
After this difference was exactly linear normal moveout correction, the first arrival ski-jump in this road departed from the deviation of 0 datum line;
In the formula: j=1,2 ... m+1; K=1,2 ... m+1 represents (m+1) * (m+1) individual grid node is calculated;
T[i] to be seismic wave propagate into used time of this geophone station physical location, i.e. first break time from i shot point (X[i], Y[i]);
Effect SDT[j is evened up in first arrival] [k] be the absolute value sum that the first arrival ski-jump departs from the deviation of 0 datum line after the linear normal moveout correction in each road,
Wherein: i=1,2 ... n represents all data tracks of this geophone station are calculated;
5) obtain first arrival and even up effect SDT[j] [k] (j=1,2 ... m+1; K=1,2 ... m+1) SDT[j0 of minimum in] [k0], its XX[j0], YY[k0] be the second positioning point position coordinates of geophone station.
2, demodulator probe secondary localization method according to claim 1 is characterized in that: in the step 3), the interval d of grid determines according to the bearing accuracy that known exploration requires; The deviation range DS of geophone station determines according to the bearing accuracy of known one-time positioning point.
3, demodulator probe secondary localization method according to claim 1 is characterized in that: repeating step 3) to 5) step finishes the second positioning of all geophone stations of construction area.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006100838781A CN100549731C (en) | 2006-06-07 | 2006-06-07 | Demodulator probe secondary localization method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006100838781A CN100549731C (en) | 2006-06-07 | 2006-06-07 | Demodulator probe secondary localization method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101086534A true CN101086534A (en) | 2007-12-12 |
CN100549731C CN100549731C (en) | 2009-10-14 |
Family
ID=38937595
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2006100838781A Expired - Fee Related CN100549731C (en) | 2006-06-07 | 2006-06-07 | Demodulator probe secondary localization method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100549731C (en) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101551464B (en) * | 2008-03-31 | 2011-04-20 | 中国石油天然气集团公司 | Complex surface three-dimensional exploration first-motion wave time determining method |
CN102147395A (en) * | 2011-01-11 | 2011-08-10 | 长安大学 | Location method for massive cement in nondestructive testing |
CN102338887A (en) * | 2010-07-14 | 2012-02-01 | 中国石油天然气集团公司 | Irregular-size space-variant grid tomography imaging statics correction method |
CN101930080B (en) * | 2009-06-26 | 2012-08-15 | 中国石油集团东方地球物理勘探有限责任公司 | Secondary positioning method for submarine cable based on surface fitting |
CN103163544A (en) * | 2013-03-20 | 2013-06-19 | 中国海洋石油总公司 | Method and system for positioning shot point |
CN103941290A (en) * | 2013-01-22 | 2014-07-23 | 中国石油集团东方地球物理勘探有限责任公司 | Submarine cable movement locus simulation method and system |
CN105866839A (en) * | 2016-06-28 | 2016-08-17 | 中国矿业大学(北京) | Static correction method and static correction device on basis of common-attitude trace gathers |
CN106443776A (en) * | 2016-06-07 | 2017-02-22 | 中国地质大学(北京) | Relocation method of ocean bottom seismometer based on time slice method |
CN107843928A (en) * | 2016-09-21 | 2018-03-27 | 中国石油化工股份有限公司 | A kind of shot point method for correcting error based on grid-search algorithms |
CN107944064A (en) * | 2016-10-12 | 2018-04-20 | 中国石油化工股份有限公司 | The parallelization computational methods and system of big gun attribute fusion |
CN108121018A (en) * | 2017-11-30 | 2018-06-05 | 中国石油天然气集团公司 | Geophone station positional accuracy evaluation method and device |
CN109212598A (en) * | 2018-10-31 | 2019-01-15 | 中国石油化工股份有限公司 | Three-dimensional space secondary localization method based on direct wave inverting |
CN112764092A (en) * | 2019-10-21 | 2021-05-07 | 中国石油天然气集团有限公司 | Detection point position positioning method and system |
-
2006
- 2006-06-07 CN CNB2006100838781A patent/CN100549731C/en not_active Expired - Fee Related
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101551464B (en) * | 2008-03-31 | 2011-04-20 | 中国石油天然气集团公司 | Complex surface three-dimensional exploration first-motion wave time determining method |
CN101930080B (en) * | 2009-06-26 | 2012-08-15 | 中国石油集团东方地球物理勘探有限责任公司 | Secondary positioning method for submarine cable based on surface fitting |
CN102338887A (en) * | 2010-07-14 | 2012-02-01 | 中国石油天然气集团公司 | Irregular-size space-variant grid tomography imaging statics correction method |
CN102338887B (en) * | 2010-07-14 | 2013-09-04 | 中国石油天然气集团公司 | Irregular-size space-variant grid tomography imaging statics correction method |
CN102147395A (en) * | 2011-01-11 | 2011-08-10 | 长安大学 | Location method for massive cement in nondestructive testing |
CN102147395B (en) * | 2011-01-11 | 2012-08-15 | 长安大学 | Location method for massive cement in nondestructive testing |
CN103941290B (en) * | 2013-01-22 | 2017-04-05 | 中国石油集团东方地球物理勘探有限责任公司 | A kind of submarine cable movement locus analogy method and system |
CN103941290A (en) * | 2013-01-22 | 2014-07-23 | 中国石油集团东方地球物理勘探有限责任公司 | Submarine cable movement locus simulation method and system |
CN103163544A (en) * | 2013-03-20 | 2013-06-19 | 中国海洋石油总公司 | Method and system for positioning shot point |
CN106443776A (en) * | 2016-06-07 | 2017-02-22 | 中国地质大学(北京) | Relocation method of ocean bottom seismometer based on time slice method |
CN105866839A (en) * | 2016-06-28 | 2016-08-17 | 中国矿业大学(北京) | Static correction method and static correction device on basis of common-attitude trace gathers |
CN107843928A (en) * | 2016-09-21 | 2018-03-27 | 中国石油化工股份有限公司 | A kind of shot point method for correcting error based on grid-search algorithms |
CN107944064A (en) * | 2016-10-12 | 2018-04-20 | 中国石油化工股份有限公司 | The parallelization computational methods and system of big gun attribute fusion |
CN107944064B (en) * | 2016-10-12 | 2020-09-29 | 中国石油化工股份有限公司 | Parallel computing method and system for shot attribute fusion |
CN108121018A (en) * | 2017-11-30 | 2018-06-05 | 中国石油天然气集团公司 | Geophone station positional accuracy evaluation method and device |
WO2019105173A1 (en) * | 2017-11-30 | 2019-06-06 | 中国石油天然气集团有限公司 | Geophone locating accuracy evaluation method and device |
US11493595B2 (en) | 2017-11-30 | 2022-11-08 | China National Petroleum Corporation | Method and apparatus for evaluating accuracy in positioning a receiver point |
CN109212598A (en) * | 2018-10-31 | 2019-01-15 | 中国石油化工股份有限公司 | Three-dimensional space secondary localization method based on direct wave inverting |
CN109212598B (en) * | 2018-10-31 | 2021-05-04 | 中国石油化工股份有限公司 | Three-dimensional space secondary positioning method based on direct wave inversion |
CN112764092A (en) * | 2019-10-21 | 2021-05-07 | 中国石油天然气集团有限公司 | Detection point position positioning method and system |
CN112764092B (en) * | 2019-10-21 | 2023-11-28 | 中国石油天然气集团有限公司 | Method and system for positioning position of wave detection point |
Also Published As
Publication number | Publication date |
---|---|
CN100549731C (en) | 2009-10-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100549731C (en) | Demodulator probe secondary localization method | |
CN103345759B (en) | Accurate detection method for submarine large complex sandwave landforms | |
CN101852867B (en) | Vectorial resultant demodulator probe secondary positioning method | |
CN101980054B (en) | Method for establishing near-surface velocity model in high-density seismic static correction processing | |
CN101930080B (en) | Secondary positioning method for submarine cable based on surface fitting | |
Kent et al. | Distribution of magma beneath the East Pacific Rise near the 9 03′ N overlapping spreading center from forward modeling of common depth point data | |
CN104730579B (en) | A kind of joint static correcting method of ripple in length and breadth based on calculation of near surface shear velocity inverting | |
Shor Jr et al. | Mohole site selection studies north of Maui | |
Searle et al. | Gloria observations of the propagating rift at 95.50 W on the Cocos‐Nazca Spreading Center | |
Fricker et al. | Digital elevation models for the Lambert Glacier–Amery Ice Shelf system, East Antarctica, from ERS-1 satellite radar altimetry | |
CN102213769A (en) | Method for determining anisotropic parameters by utilizing data of three-dimensional VSP (Vertical Seismic Profile) | |
CN103389077B (en) | Seabed sand wave geomorphology movement detection method based on MBES (multi-beam echo sounding) | |
CN102590860B (en) | Seismic wave primary arrival information-based reflected wave modeling method | |
CN105259570A (en) | Seaborne time shifting earthquake receiving point displacement correction method | |
CN102879820A (en) | Three-dimensional surface layer model construction method based on triangular grid | |
CN102269811A (en) | Edge beam water depth data sound velocity correction method based on vertical beam water depth | |
CN104569988A (en) | Echo sounding-based correction method for great sounding | |
CN105093319A (en) | Ground micro-seismic static correction method based on three-dimensional seismic data | |
CN103954996B (en) | Method and device for determining fracture strike of stratum fracture based on travel time method | |
CN101551464B (en) | Complex surface three-dimensional exploration first-motion wave time determining method | |
CN113960532B (en) | Microseism positioning method based on secondary positioning calculation of virtual source | |
CN103744117B (en) | Non-contact underwater detection method of river channel underwater foundation stone | |
CN105487046A (en) | Large-incidence-angle sound ray tracking and positioning method | |
CN106959152A (en) | A kind of method of SVEL in measurement seawater | |
CN110967751B (en) | Positioning method of micro-seismic event based on ground shallow well monitoring and storage medium |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
EE01 | Entry into force of recordation of patent licensing contract |
Assignee: TIANJIN HAILONG PETROLEUM GEOPHYSICAL EXPLORATION Co.,Ltd. Assignor: BGP INC., CHINA NATIONAL PETROLEUM Corp. Contract record no.: 2010990000820 Denomination of invention: Demodulator probe secondary localization method Granted publication date: 20091014 License type: Exclusive License Open date: 20071212 Record date: 20101015 |
|
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20091014 |