CN102944897A - Correction method for sea well shock speed scissors difference based on standard reference layer - Google Patents
Correction method for sea well shock speed scissors difference based on standard reference layer Download PDFInfo
- Publication number
- CN102944897A CN102944897A CN2012104683075A CN201210468307A CN102944897A CN 102944897 A CN102944897 A CN 102944897A CN 2012104683075 A CN2012104683075 A CN 2012104683075A CN 201210468307 A CN201210468307 A CN 201210468307A CN 102944897 A CN102944897 A CN 102944897A
- Authority
- CN
- China
- Prior art keywords
- average velocity
- reference layer
- correction
- layer
- point set
- 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 invention provides a correction method for sea well shock speed scissors difference based on a standard reference layer. The method comprises the following steps of pre-collecting a layer speed control point set and a VSP (Vertical Seismic Profile) average speed data point set on a well; the method further comprises the following steps of converting the layer speed control point into average speed control point, performing water depth correction, extracting the VSP data of the well based on layer positions, extracting average speed of the control point according to the standard reference layer, correcting layer by layer and performing tunnel correction and anti-water depth correction according to the standard reference layer. By adopting the method provided by the invention, the correction subjective of the existing prior on the position away from the well point is overcome, and the lateral deviation sensitivity to stratum speed is small; as the water depth correction and anti-water depth correction are adopted for the measured sea blocks, the influence of the seawater layer to the measurement deviation is eliminated.
Description
Technical field
The invention belongs to geological exploration field, relate to a kind of sea well shake speed price scissors bearing calibration based on the canonical reference layer.
Background technology
Seismic velocity is related to accurately degree of depth structural map accurately, and seismic velocity is that fine structures is explained, the key of one-tenth figure and ask for accurately.
VSP (Vertical Seismic Profiling) is vertical seismic profile (VSP), is a kind of seismological observation method.The vertical seismic profiling (VSP) method is earthquake-wave-exciting near some points the earth's surface, observes at some multistage multicomponent geophone stations of arranging along the wellhole different depth.In vertical seismic profiling (VSP), because wave detector places the stratum inner by well, so can not only receive up compressional wave and the up transformed wave propagated from bottom to top, also can receive descending compressional wave and the down conversion ripple propagated from top to bottom, even can receive shear wave.This is that vertical seismic profiling (VSP) is compared most important characteristics with the surface seismic section.
Speed and moving school speed that well logging obtains derive from the different measuring method, and different physical significances and accuracy standard are arranged, so the price scissors phenomenon between their average velocity of converting to is ubiquitous.The price scissors of well shake speed is proofreaied and correct, and exactly VSP is logged well average velocity as target, realizes that earthquake average velocity data volume is to the correction of the log data with clear and definite geological Significance.
The wellblock is being arranged, can utilize the logging speed of a large amount of fixed wells to come interpolation to set up calibration model, proofreading and correct earthquake average velocity, thereby obtain accurately earthquake average velocity, and then obtain accurately degree of depth structural map.But the few well in ground or marine block carries out reference without the wellblock owing to lacking enough log datas, proofreaies and correct according to traditional method of interpolation, will be difficult to guarantee accuracy.Therefore must be for the special circumstances of few well without the wellblock, the bearing calibration of reliable applicable well shake velocity contrast is developed in further investigation.
Traditional bearing calibration implementation procedure is to utilize the t-v curve of well in the work area or VSP(speed as standard, take layer position as constraint, laterally carries out interpolation and extrapolation obtains the standard speed body; T-v curve with each road of velocity spectrum compares again, proofreaies and correct.This bearing calibration is simple, workable, but also has larger limitation, at first has larger subjectivity in the position away from the well point, is difficult to judge the order of accuarcy of proofreading and correct the result; Secondly the method is only applicable to the zone that stratigraphic dip is less and lateral speed change is comparatively mild, can there be larger error at higher formation clination or the large zone of speed horizontal change, the confidence level of correcting value is lower between well, has especially destroyed the horizontal relativeness of original velocity field.
So-called canonical reference layer is that sequence top, the end and the intermediate interface of selecting to satisfy the continuous trace ability of earthquake vertical resolution and space constructed a series of satisfy tautochronism between the well shake and the aspects of condition for consistence.The canonical reference layer has played the effect of bridge between time and the degree of depth, so selecting properly canonical reference layer, necessary condition and the basis of dark conversion when being the well shake.
The prior art scheme starting point is placed on the layering statistical law based on the canonical reference layer, thinks in certain specific zone, and the t-V of each layer relation has apparent in view regularity; Rule from each canonical reference layer of minority well data statistics, again based on the rule of each layer, seismic velocity spectrum is proofreaied and correct in layering, trend and the well rule of its each layer are coincide, at last in the reference layer data after to proofread and correct on each road as benchmark, vertically by all the other reference mark of channel correction, realize the correction of all reference mark data.The method has kept the horizontal relativeness of original velocity field when realizing correction, avoid well is shaken the difficulty that velocity contrast is classified and predicted, is a kind of succinct reliably disposal route.
But above-mentioned existing velocity correction method is equivalent in fact the match t-v curve of the minority well correction target as velocity spectrum t-v curve; Have larger subjectivity in the position away from the well point, be difficult to judge the order of accuarcy of proofreading and correct the result; Be only applicable to the comparatively mild zone of the less and lateral speed change of stratigraphic dip.
Summary of the invention
For overcome prior art to the minority well area to the subjectivity of well shake velocity correction be only applicable to the limitation in the less and comparatively mild zone of lateral speed change of inclination layer, the invention provides a kind of sea well shake speed price scissors bearing calibration based on the canonical reference layer.
Sea well shake speed price scissors bearing calibration based on the canonical reference layer of the present invention comprises in advance acquisition layer speed control point set A1 and aboveground VSP average velocity data point set B1, it is characterized in that, also comprises the steps:
Step 1. interval velocity reference mark transfers the average velocity reference mark to
To the some canonical reference layers of definition in the stratum in the work area, interval velocity is controlled among the point set A1 every each reference mark interval velocity together average speed by 1. formula and transform;
1. in
Be the interval velocity of each reference mark at i canonical reference layer,
For this canonical reference layer to two-way travel time on sea level, after each reference mark data of interval velocity control point set A1 are done above-mentioned average velocity and transform, obtain average velocity control point set A2;
Step 2. water deep correction
The average velocity control point set A2 that obtains in the aboveground VSP average velocity data point set B1 that gathers in advance and the step 1 is carried out water deep correction according to 2. formula
2. middle V ' and T ' are respectively average velocity and the two-way travel time after the correction, and V and T are respectively average velocity and the seismic event two-way time before proofreading and correct, Z
0And t
0Be respectively the degree of depth and the seismic event two-way time in seabed.
The average velocity control point set that obtains in aboveground VSP average velocity data point set and the step 1 is carried out obtaining VSP average velocity correction data point set B2 and average velocity correction control point set A3 behind the above-mentioned water deep correction;
Step 3. is by each aboveground VSP data of layer bit extraction
According to the time of canonical reference layer at well location, the draw standard reference layer is in the data point of each well among the VSP average velocity correction data point set B2 that obtains in step 2, utilize linear fit method, simulate the VSP average velocity tendency equation F1 of this canonical reference layer, each canonical reference layer is carried out aforesaid operations;
Step 4. is pressed the canonical reference layer and is extracted reference mark average velocity
According to the time of canonical reference layer on the road, the draw standard reference layer is in the data point in each road, reference mark among the average velocity correction control point set A3 that obtains in step 2, utilize linear fit method, simulate the reference mark average velocity tendency equation F2 of this canonical reference layer, each canonical reference layer is carried out aforesaid operations;
Step 5. is carried out layering by the canonical reference layer and is proofreaied and correct
On each canonical reference layer, the average velocity that obtains in the step 2 is proofreaied and correct each reference mark that belongs to this layer position among the control point set A3 adopts the integral shift method to move the VSP average velocity tendency equation F1 that obtains in the step 3, point after moving is utilized linear fit method, and reference mark average velocity tendency equation F3 is proofreaied and correct in the layering that obtains this layer; Each canonical reference layer is carried out aforesaid operations;
Extract the point of this road on each canonical reference layer by the road the point after above-mentioned moving, average velocity tendency equation F4 is proofreaied and correct in the road layering that utilizes linear fit method to obtain this road; To whenever finishing together aforesaid operations;
Step 6. is by channel correction
The reference mark employing integral shift method that does not belong among the point set A3 on this canonical reference layer is proofreaied and correct in the average velocity control that obtains in the step 2 move the road layering correction average velocity tendency equation F4 that obtains in the step 5;
Step 7. is turned one's coat and is deeply proofreaied and correct
To the average velocity reference mark after the velocity correction that obtains in the step 6, according to 3. formula dark correction of turning one's coat;
3. V ' and T ' are respectively average velocity and two-way travel time after the correction in the formula, and V and T are respectively average velocity and the seismic event two-way time before proofreading and correct, Z
0And t
0Be respectively the degree of depth and the seismic event two-way time in seabed.
Concrete, 2. formula and/or 3. in the formula
, V wherein
0Be work area seawater speed.
Concrete, described integral shift method comprises the steps:
Step 31. is got and is moved a little concentrated any two points, determines straight-line equation t=av+b, and wherein a, b are coefficient, and t, v are variable;
The quilt of obtaining in the step 32. pair step 31 is moved two definite point coordinate a little concentrating, utilizes speed to move equation:
v1’=v1-[(t1-b)/a-(t1-B)/A]?----------④
V1 wherein, t1 are above-mentioned two coordinates of determining any point in the point, and v1 ' is the speed coordinate of this point after moving, and 4. A, B are the coefficient of moving the target equation in the formula, and the described target equation of moving is t=Av+B, and wherein A, B are coefficient, and t, v are variable;
The v1 ' that utilizes 4. formula to obtain is the speed coordinate after moving;
In the described integral shift method, the time coordinate of being moved a little remains unchanged.
Concrete, described linear fit method is least square method.
Sea well shake speed price scissors bearing calibration based on the canonical reference layer of the present invention, seismic velocity time relationship to each well and outer each road of well is measured, employing is based on aboveground VSP average velocity data pick-up and the earthquake reference mark average velocity data pick-up of canonical reference layer, take aboveground data as the basis, match and correction are carried out in reference mark in conjunction with each road take average velocity as sample, overcome prior art to the correction subjectivity away from the position, well point, each reference mark, road is proofreaied and correct the edge layer integral shift that adopts based on the canonical reference layer and proofreaied and correct, and is little to the lateral deviation susceptibility of formation velocity.The present invention has adopted water deep correction and dark correction of turning one's coat to measuring the sea block, has eliminated the impact of sea water layer on measured deviation.
Description of drawings
Fig. 1 illustrates the position relationship schematic diagram of a kind of embodiment on canonical reference layer of the present invention, road, VSP well and sea level;
Fig. 2 illustrates a kind of schematic diagram of embodiment of the sample point mathematics match in well of the present invention and road;
Before Fig. 3 illustrates the present invention fit equation is proofreaied and correct and the contrast schematic diagram of a kind of embodiment after proofreading and correct.
Embodiment
Below in conjunction with accompanying drawing, the specific embodiment of the present invention is described in further detail.
Sea well shake speed price scissors bearing calibration based on the canonical reference layer of the present invention, at first want the canonical reference layer in objective definition work area, definition and structure about the canonical reference layer, on the books in the prior art, for example publish " dark switch technology research when the well in the reservoir description shakes " in " petroleum exploration " phase in January, 2011.Selected canonical reference layer had both satisfied earthquake information space and vertical explanation requirement, can satisfy again that well logging information and geologic objective explain to require unified etc. the time geology cross section, therefore the chosen distance formation at target locations is nearest, and the subsequence group top, the end or the intermediate interface that satisfy the continuous trace ability of earthquake vertical resolution and space are as the canonical reference layer.The canonical reference layer has played the effect of bridge between time and the degree of depth, so selecting properly canonical reference layer, necessary condition and the basis of dark conversion when being the well shake.
The present invention is after making up the canonical reference layer, again the sea-bed area in target work area carried out data acquisition and analysis, comprise the interval velocity control point set A1 and the aboveground VSP average velocity data point set B1 that obtain in this marine site, about above-mentioned data acquisition method and data handling procedure, in the prior art existing description, for example " seismic data analysis: seismic data processing, inverting and explanation " (U.S., the author: the Ilyushin horse hereby).Fig. 1 illustrates the position relationship schematic diagram on road, VSP well, canonical reference layer and sea level in the above-mentioned target work area.The plane of figure middle and upper part is the sea level, and each bar vertical line vertical with the sea level represent and well, wherein overstriking be well, each road of thinner expression represents the canonical reference layer of constructing with several oblique lines under the sea level.The intersection point of each road and canonical reference layer is the reference mark.Part road, well, canonical reference layer have only been drawn among Fig. 1.
Sea well shake speed price scissors bearing calibration based on the canonical reference layer of the present invention also comprises the steps:
Step 1. interval velocity reference mark transfers the average velocity reference mark to
To the some canonical reference layers of definition in the stratum in the work area, interval velocity is controlled among the point set A1 every each reference mark interval velocity together average speed by 1. formula and transform;
1. in
Be the interval velocity of each reference mark at i canonical reference layer,
For this canonical reference layer to two-way travel time on sea level, after each reference mark data of interval velocity control point set A1 are done above-mentioned average velocity and transform, obtain average velocity control point set A2;
1. the equal sign of formula the right minute subitem represents the round trip mileage of each layer from the sea to i canonical reference layer cumulative, i.e. total round trip mileage length of i layer, denominator term represents that the two-way travel time of each layer from the sea to i canonical reference layer is cumulative, i.e. the two-way time length of i layer.The two-way travel time here refers to seismic event two-way travel time, and by artificial manufacturing earthquake, measuring that the seismic event related data analyzes target work area rock mass information is technological means commonly used in this area.
Step 2. water deep correction
The average velocity control point set A2 that obtains in the aboveground VSP average velocity data point set B1 that gathers in advance and the step 1 is carried out water deep correction according to 2. formula
2. middle V ' and T ' are respectively average velocity and the two-way travel time after the correction, and V and T are respectively average velocity and the seismic event two-way time before proofreading and correct, Z
0And t
0Be respectively the degree of depth and the seismic event two-way time in seabed.
The average velocity control point set that obtains in aboveground VSP average velocity data point set and the step 1 is carried out obtaining VSP average velocity correction data point set B2 and average velocity correction control point set A3 behind the above-mentioned water deep correction;
Step 3. is by each aboveground VSP data of layer bit extraction
According to the time of canonical reference layer at well location, the draw standard reference layer is in the data point of each well among the VSP average velocity correction data point set B2 that obtains in step 2, utilize linear fit method, simulate the VSP average velocity tendency equation F1 of this canonical reference layer, each canonical reference layer is carried out aforesaid operations;
Step 4. is pressed the canonical reference layer and is extracted reference mark average velocity
According to the time of canonical reference layer on the road, the draw standard reference layer is in the data point in each road, reference mark among the average velocity correction control point set A3 that obtains in step 2, utilize linear fit method, simulate the reference mark average velocity tendency equation F2 of this canonical reference layer, each canonical reference layer is carried out aforesaid operations;
Each canonical reference layer time on well location and road in step 3 and the step 4 is the time coordinate of having set during self-defined canonical reference layer in advance.
Step 5. is carried out layering by the canonical reference layer and is proofreaied and correct
On each canonical reference layer, the average velocity that obtains in the step 2 is proofreaied and correct each reference mark that belongs to this layer position among the control point set A3 adopts the integral shift method to move the VSP average velocity tendency equation F1 that obtains in the step 3, point after moving is utilized linear fit method, and reference mark average velocity tendency equation F3 is proofreaied and correct in the layering that obtains this layer; Each canonical reference layer is carried out aforesaid operations;
Extract the point of this road on each canonical reference layer by the road the point after above-mentioned moving, average velocity tendency equation F4 is proofreaied and correct in the road layering that utilizes linear fit method to obtain this road; To whenever finishing together aforesaid operations; An embodiment who for example provides among Fig. 2, point set before the expression of Fig. 2 the right is moved, point set and tendency equation after the left side represents to move, concrete equation and variance parameter before and after moving have been provided among Fig. 2, horizontal ordinate represents the speed dimension, ordinate represents the time dimension, only represents inventor's a specific embodiment.
Step 6. is by channel correction
The reference mark employing integral shift method that does not belong among the point set A3 on this canonical reference layer is proofreaied and correct in the average velocity control that obtains in the step 2 move the road layering correction average velocity tendency equation F4 that obtains in the step 5; An embodiment who for example provides among Fig. 3, horizontal ordinate represents the speed dimension, ordinate represents the time dimension, the point set before Fig. 3 the right expression is moved, point set and tendency equation after the left side represents to move only represent inventor's a specific embodiment.
Step 7. is turned one's coat and is deeply proofreaied and correct
To the average velocity reference mark after the velocity correction that obtains in the step 6, according to 3. formula dark correction of turning one's coat;
3. V ' and T ' are respectively average velocity and two-way travel time after the correction in the formula, and V and T are respectively average velocity and the seismic event two-way time before proofreading and correct, Z
0And t
0Be respectively the degree of depth and the seismic event two-way time in seabed.
Sea well shake speed price scissors bearing calibration based on the canonical reference layer of the present invention, the following degree of depth in seabed that is used for Shao Jing work area in the analysis marine site is constructed.
In the 3. formula of the 2. formula of above-mentioned steps 2 and step 7, Z
0And t
0Be respectively the degree of depth in seabed and seabed to the seismic event two-way travel time on sea, Z
0And t
0Relation and seismic event travelling speed V in seawater
0Relevant, when measuring Z
0And t
0Among one the time, can utilize
Ask for another, here hypothetically seismic wave in the seabed in the whole sea water layer on sea velocity of propagation equate everywhere, can suppose V
0=1500 meter per seconds.
Concrete, comprise the steps: in the method for integral shift described in step 5 and the step 6
Step 31. is got and is moved a little concentrated any two points, determines straight-line equation t=av+b, and wherein a, b are coefficient, and t, v are variable;
The quilt of obtaining in the step 32. pair step 31 is moved two definite point coordinate a little concentrating, utilizes speed to move equation:
v1’=v1-[(t1-b)/a-(t1-B)/A]?----------④
V1 wherein, t1 are above-mentioned two coordinates of determining any point in the point, and v1 ' is the speed coordinate of this point after moving, 4. A, B are the coefficient of moving the target equation in the formula, the described target equation of moving is t=Av+B, and wherein A, B are coefficient, and time t, speed v are variable;
The v1 ' that utilizes 4. formula to obtain is the speed coordinate after moving;
In the described integral shift method, the time coordinate of being moved a little remains unchanged.
For example in the step 5, described VSP average velocity tendency equation F1 is for moving the target equation, average velocity is proofreaied and correct each reference mark that belongs to this layer position among the control point set A3 and is moved point set for quilt, in the step 6, the layering of described road is proofreaied and correct average velocity tendency equation F4 for moving the target equation, and average velocity control is proofreaied and correct the reference mark that does not belong among the point set A3 on this canonical reference layer and moved point set for quilt.
The linear fit method of utilizing point set to draw linear equation in each step of sea well shake speed price scissors bearing calibration based on the canonical reference layer of the present invention has several different methods, preferred least square method among the present invention in the prior art.
Sea well shake speed price scissors bearing calibration based on the canonical reference layer of the present invention, seismic velocity time relationship to each well and outer each road of well is measured, employing is based on aboveground VSP average velocity data pick-up and the earthquake reference mark average velocity data pick-up of canonical reference layer, take aboveground data as the basis, match and correction are carried out in reference mark in conjunction with each road take average velocity as sample, overcome prior art to the correction subjectivity away from the position, well point, each reference mark, road is proofreaied and correct the edge layer integral shift that adopts based on the canonical reference layer and proofreaied and correct, and is little to the lateral deviation susceptibility of formation velocity.The present invention has adopted water deep correction and dark correction of turning one's coat to measuring the sea block, has eliminated the impact of sea water layer on measured deviation.
Previously described is each preferred embodiment of the present invention; design parameter among described embodiment and the embodiment or Figure of description only are the invention proof procedures for clear statement inventor; be not to limit scope of patent protection of the present invention; scope of patent protection of the present invention still is as the criterion with its claims; the equivalent structure that every utilization instructions of the present invention and accompanying drawing content are done changes, and in like manner all should be included in protection scope of the present invention.
Claims (4)
1. based on the sea well shake speed price scissors bearing calibration of canonical reference layer, comprise in advance acquisition layer speed control point set A1 and aboveground VSP average velocity data point set B1, it is characterized in that, also comprise the steps:
Step 1. interval velocity reference mark transfers the average velocity reference mark to
To the some canonical reference layers of definition in the stratum in the work area, interval velocity is controlled among the point set A1 every each reference mark interval velocity together average speed by 1. formula and transform;
----------①
1. in
Be the interval velocity of each reference mark at i canonical reference layer,
For this canonical reference layer to two-way travel time on sea level, after each reference mark data of interval velocity control point set A1 are done above-mentioned average velocity and transform, obtain average velocity control point set A2;
Step 2. water deep correction
The average velocity control point set A2 that obtains in the aboveground VSP average velocity data point set B1 that gathers in advance and the step 1 is carried out water deep correction according to 2. formula
2. middle V ' and T ' are respectively average velocity and the two-way travel time after the correction, and V and T are respectively average velocity and the seismic event two-way time before proofreading and correct, Z
0And t
0Be respectively the degree of depth and the seismic event two-way time in seabed.
The average velocity control point set that obtains in aboveground VSP average velocity data point set and the step 1 is carried out obtaining VSP average velocity correction data point set B2 and average velocity correction control point set A3 behind the above-mentioned water deep correction;
Step 3. is by each aboveground VSP data of layer bit extraction
According to the time of canonical reference layer at well location, the draw standard reference layer is in the data point of each well among the VSP average velocity correction data point set B2 that obtains in step 2, utilize linear fit method, simulate the VSP average velocity tendency equation F1 of this canonical reference layer, each canonical reference layer is carried out aforesaid operations;
Step 4. is pressed the canonical reference layer and is extracted reference mark average velocity
According to the time of canonical reference layer on the road, the draw standard reference layer is in the data point in each road, reference mark among the average velocity correction control point set A3 that obtains in step 2, utilize linear fit method, simulate the reference mark average velocity tendency equation F2 of this canonical reference layer, each canonical reference layer is carried out aforesaid operations;
Step 5. is carried out layering by the canonical reference layer and is proofreaied and correct
On each canonical reference layer, the average velocity that obtains in the step 2 is proofreaied and correct each reference mark that belongs to this layer position among the control point set A3 adopts the integral shift method to move the VSP average velocity tendency equation F1 that obtains in the step 3, point after moving is utilized linear fit method, and reference mark average velocity tendency equation F3 is proofreaied and correct in the layering that obtains this layer; Each canonical reference layer is carried out aforesaid operations;
Extract the point of this road on each canonical reference layer by the road the point after above-mentioned moving, average velocity tendency equation F4 is proofreaied and correct in the road layering that utilizes linear fit method to obtain this road; To whenever finishing together aforesaid operations;
Step 6. is by channel correction
The reference mark employing integral shift method that does not belong among the point set A3 on this canonical reference layer is proofreaied and correct in the average velocity control that obtains in the step 2 move the road layering correction average velocity tendency equation F4 that obtains in the step 5;
Step 7. is turned one's coat and is deeply proofreaied and correct
To the average velocity reference mark after the velocity correction that obtains in the step 6, according to 3. formula dark correction of turning one's coat;
3. V ' and T ' are respectively average velocity and two-way travel time after the correction in the formula, and V and T are respectively average velocity and the seismic event two-way time before proofreading and correct, Z
0And t
0Be respectively the degree of depth and the seismic event two-way time in seabed.
2. the sea well shake speed price scissors bearing calibration based on the canonical reference layer as claimed in claim 1 is characterized in that: 2. formula and/or 3. in the formula
, V wherein
0Be work area seawater speed.
3. the sea well shake speed price scissors bearing calibration based on the canonical reference layer as claimed in claim 1, it is characterized in that: described integral shift method comprises the steps:
Step 31. is got and is moved a little concentrated any two points, determines straight-line equation t=av+b, and wherein a, b are coefficient, and t, v are variable;
The quilt of obtaining in the step 32. pair step 31 is moved two definite point coordinate a little concentrating, utilizes speed to move equation:
v1’=v1-[(t1-b)/a-(t1-B)/A]?----------④
V1 wherein, t1 are above-mentioned two coordinates of determining any point in the point, and v1 ' is the speed coordinate of this point after moving, and 4. A, B are the coefficient of moving the target equation in the formula, and the described target equation of moving is t=Av+B, and wherein A, B are coefficient, and t, v are variable;
The v1 ' that utilizes 4. formula to obtain is the speed coordinate after moving;
A little concentrated all in twos groupings of point to carry out aforesaid operations to moving;
In the described integral shift method, the time coordinate of being moved a little remains unchanged.
4. the sea well shake speed price scissors bearing calibration based on the canonical reference layer as claimed in claim 1, it is characterized in that: described linear fit method is least square method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210468307.5A CN102944897B (en) | 2012-11-20 | 2012-11-20 | Correction method for sea well shock speed scissors difference based on standard reference layer |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210468307.5A CN102944897B (en) | 2012-11-20 | 2012-11-20 | Correction method for sea well shock speed scissors difference based on standard reference layer |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102944897A true CN102944897A (en) | 2013-02-27 |
| CN102944897B CN102944897B (en) | 2015-07-01 |
Family
ID=47727858
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210468307.5A Active CN102944897B (en) | 2012-11-20 | 2012-11-20 | Correction method for sea well shock speed scissors difference based on standard reference layer |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102944897B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2018099880A1 (en) * | 2016-12-02 | 2018-06-07 | Spotlight | Method for improving seismic acquisitions utilising active ultralight seismic detection systems |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4884247A (en) * | 1987-03-09 | 1989-11-28 | Mobil Oil Company | Method of processing geophysical data to compensate for earth filter attenuation |
| CN102279416A (en) * | 2010-12-30 | 2011-12-14 | 北京诺克斯达石油科技有限公司 | Method for matching far-and-near offset-distance stacked profile phase frequencies |
-
2012
- 2012-11-20 CN CN201210468307.5A patent/CN102944897B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4884247A (en) * | 1987-03-09 | 1989-11-28 | Mobil Oil Company | Method of processing geophysical data to compensate for earth filter attenuation |
| CN102279416A (en) * | 2010-12-30 | 2011-12-14 | 北京诺克斯达石油科技有限公司 | Method for matching far-and-near offset-distance stacked profile phase frequencies |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2018099880A1 (en) * | 2016-12-02 | 2018-06-07 | Spotlight | Method for improving seismic acquisitions utilising active ultralight seismic detection systems |
| CN110050204A (en) * | 2016-12-02 | 2019-07-23 | 聚光灯公司 | A method of earthquake-capturing is improved using active ultralight seismic acquisition system |
| JP2019536050A (en) * | 2016-12-02 | 2019-12-12 | スポットライト | A method for improving seismic acquisition utilizing an active ultralight seismic detection system |
| RU2751573C2 (en) * | 2016-12-02 | 2021-07-15 | Спотлайт | Method for improving seismic data collection using ultra light active seismic control systems |
| CN110050204B (en) * | 2016-12-02 | 2021-12-03 | 聚光灯公司 | Method for improving seismic acquisition using active ultra-light seismic detection system |
| US11435492B2 (en) | 2016-12-02 | 2022-09-06 | Spotlight | Method for improving seismic acquisitions utilising active ultralight seismic detection systems |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102944897B (en) | 2015-07-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102759745B (en) | Method for forecasting carbonate reservoir based on forward modeling of digital geological outcrop model | |
| CN100487488C (en) | Zero hypocentral distance vertical seismic section compressional-shear wave data depth field corridor stacked section processing method | |
| CN106837297A (en) | A kind of method for recognizing inter well connectivity and profit dynamic prediction | |
| CN106875471A (en) | Coal measures contains or water barrier Visualization Modeling method | |
| CN103046868A (en) | Integrated horizontal well geosteering method | |
| CN103245971A (en) | Well breaking point guided earthquake minor fault interpretation method and device | |
| CN109870719B (en) | Well position layout method, device and system for carbonate rock tight thin reservoir | |
| CN103114840A (en) | Calculating method and device of organic carbon content of high-over high mature shale | |
| CN102272414A (en) | Method for optimizing well production in reservoirs having flow barriers | |
| CN105510993A (en) | Foreland basin deep buried and compressed type complex gypsum-salt rock identification and distribution prediction method | |
| CN106501854A (en) | Oil gas drilling target geological quantifying risk evaluation methodology | |
| CN105242307B (en) | Carbonate complex seismic reservoir porosity acquisition methods and device | |
| CN105911592B (en) | A kind of real three dimensional seismic data long wavelength's static correcting method for boring constraint | |
| Chen et al. | Study on the application of a comprehensive technique for geological prediction in tunneling | |
| CN104133250A (en) | Gypsiferous salt stratum geological layering method | |
| CN103821505A (en) | Sandstone oil and gas carrier bed geophysics-geology-geochemistry detection method and device | |
| CN103088803A (en) | Method of distinguishing collapse causing factors of karst collapse | |
| CN104899928A (en) | Three-directional geology modeling method based on sparse borehole points | |
| CN103792578B (en) | A kind of Fluid Identification Method of frequency dispersion AVO association attributes inverting | |
| CN105551083A (en) | Method for depicting relationship between horizontal well track and strata | |
| US20130327125A1 (en) | Method for geochemical grandient exploration | |
| Qi et al. | Hierarchy and subsurface correlation of muddy baffles in lacustrine delta fronts: a case study in the X Oilfield, Subei Basin, China | |
| Agoha et al. | Weathered layer determination and its effects on engineering structures: case study of parts of the Niger delta | |
| Liu et al. | The Control Theory and Application for Well Pattern Optimization of Heterogeneous Sandstone Reservoirs | |
| CN114910953B (en) | Method for explaining irregular section of complex fracture zone |
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 |