CN103869369B - Method for correcting interval velocity by using residue curvature - Google Patents

Method for correcting interval velocity by using residue curvature Download PDF

Info

Publication number
CN103869369B
CN103869369B CN201410092817.6A CN201410092817A CN103869369B CN 103869369 B CN103869369 B CN 103869369B CN 201410092817 A CN201410092817 A CN 201410092817A CN 103869369 B CN103869369 B CN 103869369B
Authority
CN
China
Prior art keywords
layer
interval velocity
velocity
updating
interval
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
CN201410092817.6A
Other languages
Chinese (zh)
Other versions
CN103869369A (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
Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd
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 Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd filed Critical Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd
Priority to CN201410092817.6A priority Critical patent/CN103869369B/en
Publication of CN103869369A publication Critical patent/CN103869369A/en
Application granted granted Critical
Publication of CN103869369B publication Critical patent/CN103869369B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Lining And Supports For Tunnels (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention provides a method for correcting an interval velocity by using a residue curvature. The method comprises the steps of extracting root-mean-square velocities of m intervals from a vertical speed function, and obtaining the interval velocities of the m intervals before updating according to the extracted root-mean-square velocities of the m intervals; solving the average interval velocity of the m intervals before updating; establishing a general function about the interval velocities of the m intervals and calculating the corresponding interval velocities of the m intervals when the general function is minimal as the interval velocities of the m intervals after updating; solving the average interval velocity of the m intervals after updating; calculating the residue curvature of the m intervals after updating according to the average interval velocity of the m intervals before updating and the average interval velocity of the m intervals after updating; correcting the interval velocities of the m layers before updating by using the residue curvature of the m intervals after updating. The stable interval velocity can be obtained by adopting the method and high-frequency oscillation due to residue curvature updating is prevented.

Description

Method interval velocity being corrected using remaining curvature
Technical field
The invention belongs to oil seismic exploration field, it is mainly used in seismic prospecting data and processes and explain, more specifically Say, be related to a kind of realize stable method interval velocity being corrected using remaining curvature.
Background technology
In process of seismic data processing, interval velocity serves very important effect, the direct shadow of precision of rate pattern Ring the quality of imaging results.
At present, the method for existing correction interval velocity is to extract stack velocity from velocity function spectrum, and will be superimposed speed Degree is converted to root-mean-square(RMS)Speed, then with Dix formula, root mean sequare velocity is converted into interval velocity.In theory, the party Method is only applicable to horizontal layered earth and the situation of little geophone offset, poor anti jamming capability, stack velocity or root mean sequare velocity micro- Little change all will make remaining curvature update and produce the higher-order of oscillation.
Content of the invention
It is an object of the invention to proposing a kind of method interval velocity being corrected using remaining curvature it is achieved that stablizing Reliable instantaneous velocity calculates, it is to avoid the HFO causing because interval velocity is unstable.
An aspect of of the present present invention provides a kind of method interval velocity being corrected using remaining curvature, methods described bag Include:(a)The root mean sequare velocity of m layer is extracted from vertical speed functionWherein,Mean square for the n-th layer in m layer Root speed, n ∈ [1, m], according to the root mean sequare velocity of the m layer extractingThe interval velocity of m layer before being updatedIts In,Interval velocity for the n-th layer in m layer before updating;(b)Interval velocity to m layer before updatingIt is averaging, obtain The average interval velocity of the m layer to before updateWherein,Average layer speed for the n-th layer in m layer before updating Degree;(c)Set up interval velocity V with regard to m layernFunctional, wherein, VnFor the interval velocity of the n-th layer in m layer, calculate Functional is interval velocity V of corresponding m layer during minimumn, by interval velocity V of calculated m layernAs m after updating The interval velocity of layerWherein,Interval velocity for the n-th layer in m layer after updating;(d)To m layer after updating Interval velocityIt is averaging, the average interval velocity of m layer after being updatedWherein,For in m layer after updating N-th layer average interval velocity;(e)Average interval velocity according to m layer before updatingWith update after m layer average Interval velocityCalculate remaining curvature γ of m layer after updatingn, wherein, γnFor the n-th layer in m layer after updating Remaining curvature;(f)Remaining curvature γ using m layer after updatingnInterval velocity to m layer before updatingIt is corrected.
Alternatively, step(a)The interval velocity of the n-th layer in the middle m layer calculating before updatingFormula as follows,
Wherein,For the root mean sequare velocity of (n-1)th layer in m layer, tnFor the time of the n-th layer in m layer, tn-1For The time of (n-1)th layer in m layer.
Alternatively, step(b)The average interval velocity of the n-th layer in the middle m layer calculating before updatingFormula as follows,
Wherein,For the interval velocity of i-th layer in n layer before updating, znLayer locating depth for the n-th layer in m layer Degree, zn-1Layer bit depth for (n-1)th layer in m layer.
Alternatively, described functional is that data item fitting expression, trend term fitting expression and smooth item matching are expressed The sum of formula.
Alternatively, the interval velocity according to m layer before updatingInterval velocity V with m layern, obtain data item matching table Reach formula.
Alternatively, data item fitting expression is,
Wherein, B represents data item fitting expression, Δ tnFor the time interval of the n-th layer in m layer and (n-1)th layer, Δ tn=tn-tn-1, tnFor the time of the n-th layer in m layer, tn-1For the time of (n-1)th layer in m layer,For in m layer The weight coefficient of the data item of n-layer, Vn-1Interval velocity for (n-1)th layer in m layer.
Alternatively, interval velocity V according to m layernAnd trend function, obtain trend term fitting expression.
Alternatively, the calculation expression of trend function is,
Wherein, VtrendT () represents trend function, VaFor the instantaneous velocity on ground top layer, VInstantaneous speed for infinite deep layer Degree, Δ V is the excursion of instantaneous velocity, kaFor previously given control parameter, t is the time.
Alternatively, wherein, trend term fitting expression is,
Wherein, C represents trend term fitting expression, Δ tnFor the time interval of the n-th layer in m layer and (n-1)th layer, Δ tn=tn-tn-1, tnFor the time of the n-th layer in m layer, tn-1For the time of (n-1)th layer in m layer,For in m layer The weight coefficient of the trend term of n-th layer,Interval velocity V for (n-1)th layer in m layern-11- (τ/Δ tn) power,For in m layer n-th layer interval velocity Vnτ/Δ tnPower, Vtrend(tn-1+ τ) it is time t=tn-1Becoming during+τ Potential function.
Alternatively, choose interval velocity V to m layer for the smoothing factornIt is smoothed, obtain smooth item matching expression Formula.
Alternatively, smooth item fitting expression is,
Wherein, D represents smooth item fitting expression,For unit characteristic number,Flat for the n-th layer in m layer The sliding factor, Δ tnFor the time interval of the n-th layer in m layer and (n-1)th layer, Δ tn=tn-tn-1, Δ tn+1For n-th in m layer + 1 layer of time interval with n-th layer, Δ tn+1=tn+1-tn, tnFor the time of the n-th layer in m layer, tn-1For n-th in m layer- 1 layer of time, tn+1For the time of (n+1)th layer in m layer, VnFor the interval velocity of the n-th layer in m layer, Vn+1For in m layer (n+1)th layer of interval velocity.
Alternatively, calculate remaining curvature γ of the n-th layer in m layer after updatingnFormula be,
The interval velocity of m layer is carried out using the method that remaining curvature is corrected to interval velocity using of the present invention Correction, makes the interval velocity of m layer after correction more accurate, is conducive to the judgement to subsurface geologic structures and pre- to oil-gas Layer Survey.
Brief description
By the detailed description carrying out below in conjunction with the accompanying drawings, above and other objects of the present invention, feature and advantage will Become more fully apparent, wherein:
Fig. 1 is the flow process of the method illustrating using remaining curvature, interval velocity to be corrected according to an embodiment of the invention Figure.
Specific embodiment
Describe embodiments of the invention next, with reference to accompanying drawing in detail.
Offer following description referring to the drawings is to help the enforcement to the present invention being limited by claim and its equivalent The comprehensive understanding of example.Including various specific detail to help understand, but these details are considered only as being exemplary.Therefore, originally Field skilled artisan will realize that without departing from the scope and spirit of the present invention, can be to reality described herein Apply example to make various changes and modifications.Additionally, for clarity and brevity, omit the description to known function and structure.
Fig. 1 is the flow process of the method illustrating using remaining curvature, interval velocity to be corrected according to an embodiment of the invention Figure.
With reference to Fig. 1, from vertical speed function, in a step 101, extract the root mean sequare velocity of m layerWherein, For the root mean sequare velocity of the n-th layer in m layer, n ∈ [1, m], according to the root mean sequare velocity of the m layer extractingUpdated The interval velocity of m front layerWherein,Interval velocity for the n-th layer in m layer before updating.
For example, various existing methods are can be utilized to extract the root mean sequare velocity of m layer from vertical speed functionWith m The time t of layern, tnFor the time of the n-th layer in m layer, wherein, as n=1, in m layer (n-1)th layer of root mean sequare velocityFor the root mean sequare velocity on ground top layer, the time t of (n-1)th in m layer layer0For the time on ground top layer, in m layer Layer bit depth z of n-1 layer0For the layer bit depth on ground top layer, interval velocity V of (n-1)th in m layer layer0Layer for ground top layer Speed.Here, the root mean sequare velocity on ground top layerTime t with ground top layer0Can extract from vertical speed function and obtain.
For example, calculate the interval velocity of the n-th layer in m layer before updatingFormula as follows,
Formula(1)In,For the root mean sequare velocity of (n-1)th layer in m layer, tn-1For (n-1)th layer in m layer when Between.
In a step 102, the interval velocity to m layer before updatingBe averaging, m layer before being updated average Interval velocityWherein,Average interval velocity for the n-th layer in m layer before updating.
For example, calculate the average interval velocity of the n-th layer in m layer before updatingFormula as follows
Formula(2)In,For the interval velocity of i-th layer in n layer before updating, znLayer for the n-th layer in m layer Bit depth, zn-1Layer bit depth for (n-1)th layer in m layer.
In step 103, set up interval velocity V with regard to m layernFunctional, wherein, VnFor n-th in m layer The interval velocity of layer, calculates interval velocity V that functional is corresponding m layer during minimumn, by interval velocity V of calculated m layern Interval velocity as m layer after updatingWherein,Interval velocity for the n-th layer in m layer after updating.Here, The minimum of functional can be calculated using various existing methods, for example, calculate the minimum of functional using quasi-Newton method.
For example, functional can be data item fitting expression, trend term fitting expression and smooth item fitting expression With.
For example, can be according to the interval velocity of m layer before updatingInterval velocity V with m layern, obtain data item matching table Reach formula.For example, data item fitting expression can be,
Formula(3)In, B represents data item fitting expression, Δ tnBetween the time for the n-th layer in m layer and (n-1)th layer Every Δ tn=tn-tn-1,For the weight coefficient of the data item of the n-th layer in m layer, Vn-1Layer speed for (n-1)th layer in m layer Degree.
For example, can be according to interval velocity V of m layernAnd trend function, obtain trend term fitting expression.For example, trend letter Number calculation expression can be,
Formula(4)In, VtrendT () represents trend function, VaFor the instantaneous velocity on ground top layer, VWink for infinite deep layer Shi Sudu, Δ V are the excursion of instantaneous velocity, kaFor previously given control parameter, t is the time.
For example, trend term fitting expression can be,
Formula(5)In, C represents trend term fitting expression,For the weight coefficient of the trend term of the n-th layer in m layer,Interval velocity V for (n-1)th layer in m layern-11- (τ/Δ tn) power,For in m layer n-th layer Interval velocity Vnτ/Δ tnPower, Vtrend(tn-1+ τ) it is time t=tn-1Trend function during+τ.
For example, choose smoothing factor to m layer interval velocity VnIt is smoothed, obtain smooth item matching expression Formula.For example, smooth item fitting expression is,
Formula(6)In, D represents smooth item fitting expression,For unit characteristic number,For the n-th layer in m layer Smoothing factor, Δ tn+1For (n+1)th layer of time interval with n-th layer in m layer, Δ tn+1=tn+1-tn, tn+1For in m layer (n+1)th layer of time, Vn+1Interval velocity for (n+1)th layer in m layer.
The weight coefficient of the data item of the n-th layer in m layerThe weight coefficient of the trend term of the n-th layer in m layerAnd m The smoothing factor of the n-th layer in individual layerValue can artificially be set(For example, ), only need to meet?.
At step 104, the interval velocity to m layer after updatingBe averaging, m layer after being updated average Interval velocityWherein,Average interval velocity for the n-th layer in m layer after updating.
For example, calculate the average interval velocity of the n-th layer in m layer after updatingFormula as follows,
Formula(7)In,For update after m layer in n-th layer average interval velocity,For m after updating The interval velocity of i-th layer in layer.
In step 105, the average interval velocity according to m layer before updatingAverage layer with m layer after renewal SpeedCalculate remaining curvature γ of m layer after updatingn, wherein, γnSurplus for the n-th layer in m layer after updating Remaining curvature.
For example, calculate remaining curvature γ of the n-th layer in m layer after updatingnFormula be,
In step 106, using remaining curvature γ of m layer after updatingnInterval velocity to m layer before updating It is corrected.
Specifically, remaining curvature γ of m layer after being updated according to method of the present inventionn, typically recognize It is the interval velocity of m layer after the renewal being obtained using method of the present inventionFor the interval velocity of accurate m layer, count Remaining curvature γ of m layer after the renewal obtainingnAverage interval velocity for m layer before updatingWith the m after renewal The average interval velocity of individual layerRatio, therefore can be utilized this ratio to update before m layer interval velocityCarry out school Just, specifically, by the interval velocity of m layer before updatingIt is multiplied by the remaining curvature of m layer after respective renewal respectively γnInverse, the interval velocity V ' of m layer after being correctedn.For example, calculate the interval velocity of the n-th layer in m layer after correction V′nFormula as follows:
Through formula(9)The interval velocity V ' of the m layer obtaining after correctionnMore stable reliability, it is to avoid oscillatory occurences.For example, As γ 5=0.5, after correction the 5th layer of interval velocity
According to the interval velocity of m layer after the renewal that method of the present invention obtains, can be with m layer before updating Interval velocity keeps less deviation, the value of functional can be made less, such interval velocity is more reliable and more stable again.
After the interval velocity of m layer before updating being corrected using remaining curvature using method of the present invention, right Answer the layer bit depth on time point also can be adjusted, for example, due to the time of interval velocity and this layer that layer bit depth is this layer Product, therefore, after interval velocity changes, corresponding layer bit depth also can change, can according to correction after m layer interval velocity and The product of the time of m layer corrected after m layer layer bit depth, m layer after remaining curvature correction layer speed The layer bit depth of degree and m layer is more accurate, is conducive to the judgement to subsurface geologic structures and the prediction to oil-gas Layer.
Method of the present invention obtains reliable and stable instantaneous velocity based on the root mean sequare velocity of pickup, by belt restraining Optimal way realize stable remaining curvature speed and update.
Although be particularly shown and described the present invention, those skilled in the art with reference to its exemplary embodiment It should be understood that in the case of the spirit and scope of the present invention being limited without departing from claim, can be so that form be carried out to it With the various changes in details.

Claims (11)

1. a kind of method interval velocity being corrected using remaining curvature, methods described includes:
A () extracts the root mean sequare velocity of m layer from vertical speed functionWherein,Mean square for the n-th layer in m layer Root speed, n ∈ [1, m],
Root mean sequare velocity according to the m layer extractingThe interval velocity of m layer before being updatedWherein,For updating The interval velocity of the n-th layer in m front layer;
(b) interval velocity to m layer before updatingIt is averaging, the average interval velocity of m layer before being updatedIts In,Average interval velocity for the n-th layer in m layer before updating;
C () sets up interval velocity V with regard to m layernFunctional, wherein, VnFor the interval velocity of the n-th layer in m layer, count Calculate interval velocity V that functional is corresponding m layer during minimumn, by interval velocity V of calculated m layernAs the m after updating The interval velocity of individual layerWherein,Interval velocity for the n-th layer in m layer after updating;
(d) interval velocity to m layer after updatingIt is averaging, the average interval velocity of m layer after being updatedIts In,Average interval velocity for the n-th layer in m layer after updating;
E () is according to the average interval velocity of m layer before updatingAverage interval velocity with m layer after renewalCalculate more Remaining curvature γ of m layer after newn, wherein, γnRemaining curvature for the n-th layer in m layer after updating;
F () is using remaining curvature γ of m layer after updatingnInterval velocity to m layer before updatingIt is corrected,
Wherein, calculate the average interval velocity of the n-th layer in m layer before updating in step (b)Formula as follows,
V n o l d ‾ = z n Σ i = 1 n z n - z n - 1 V i o l d ,
Wherein, Vi oldFor the interval velocity of i-th layer in n layer before updating, znFor the layer bit depth of the n-th layer in m layer, zn-1 Layer bit depth for (n-1)th layer in m layer.
2. method according to claim 1, wherein, calculates the layer speed of the n-th layer in m layer before updating in step (a) DegreeFormula as follows,
V n o l d = V ~ n 2 t n - V ~ n - 1 2 t n - 1 t n - t n - 1 ,
Wherein,For the root mean sequare velocity of (n-1)th layer in m layer, tnFor the time of the n-th layer in m layer, tn-1For m layer In (n-1)th layer of time.
3. method according to claim 1, wherein, described functional is data item fitting expression, trend term matching table Reach formula and the sum of smooth item fitting expression.
4. method according to claim 3, wherein, according to the interval velocity of m layer before updatingLayer speed with m layer Degree Vn, obtain data item fitting expression.
5. method according to claim 4, wherein, data item fitting expression is,
B = 1 2 Σ n = 1 m Δt n · w n B [ V n 2 - V n - 1 2 2 l n ( V n / V n - 1 ) - V n o l d ] 2 ,
Wherein, B represents data item fitting expression, Δ tnFor the time interval of the n-th layer in m layer and (n-1)th layer, Δ tn= tn-tn-1, tnFor the time of the n-th layer in m layer, tn-1For the time of (n-1)th layer in m layer,For n-th in m layer The weight coefficient of the data item of layer, Vn-1Interval velocity for (n-1)th layer in m layer.
6. method according to claim 3, wherein, according to interval velocity V of m layernAnd trend function, obtain trend term and intend Close expression formula.
7. method according to claim 6, wherein, the calculation expression of trend function is,
V t r e n d ( t ) = V a · V ∞ V a + Δ V · exp ( - k a t · V ∞ / Δ V ) ,
Wherein, VtrendT () represents trend function, VaFor the instantaneous velocity on ground top layer, VFor the instantaneous velocity of infinite deep layer, Δ V is the excursion of instantaneous velocity, kaFor previously given control parameter, t is the time.
8. method according to claim 7, wherein, trend term fitting expression is,
C = 1 2 Σ n = 1 m ∫ 0 Δt n w n C [ V n - 1 1 - ( τ / Δt n ) · V n τ / Δt n - V t r e n d ( t n - 1 + τ ) ] 2 d τ ,
Wherein, C represents trend term fitting expression, Δ tnFor the time interval of the n-th layer in m layer and (n-1)th layer, Δ tn= tn-tn-1, tnFor the time of the n-th layer in m layer, tn-1For the time of (n-1)th layer in m layer,For n-th in m layer The weight coefficient of the trend term of layer,Interval velocity V for (n-1)th layer in m layern-11- (τ/Δ tn) power, For in m layer n-th layer interval velocity Vnτ/Δ tnPower, Vtrend(tn-1+ τ) it is time t=tn-1Trend letter during+τ Number.
9. method according to claim 3, wherein, chooses interval velocity V to m layer for the smoothing factornIt is smoothed, Obtain smooth item fitting expression.
10. method according to claim 9, wherein, smooth item fitting expression is,
D = V 2 Δ t ‾ 2 Σ n = 1 m w n D · Δt n · Δt n + 1 [ l n ( V n + 1 / V n ) Δt n + 1 · l n ( V n / V n - 1 ) Δt n ] 2 ,
Wherein, D represents smooth item fitting expression,For unit characteristic number,For the smoothing factor of the n-th layer in m layer, ΔtnFor the time interval of the n-th layer in m layer and (n-1)th layer, Δ tn=tn-tn-1, Δ tn+1For (n+1)th layer in m layer with The time interval of n-th layer, Δ tn+1=tn+1-tn, tnFor the time of the n-th layer in m layer, tn-1For (n-1)th layer in m layer Time, tn+1For the time of (n+1)th layer in m layer, VnFor the interval velocity of the n-th layer in m layer, Vn+1For n-th in m layer + 1 layer of interval velocity, Vn-1Interval velocity for (n-1)th layer in m layer.
11. methods according to claim 1, wherein, calculate remaining curvature γ of the n-th layer in m layer after updatingn's Formula is,
γ n = V n o l d ‾ V n n e w ‾ .
CN201410092817.6A 2014-03-13 2014-03-13 Method for correcting interval velocity by using residue curvature Active CN103869369B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410092817.6A CN103869369B (en) 2014-03-13 2014-03-13 Method for correcting interval velocity by using residue curvature

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410092817.6A CN103869369B (en) 2014-03-13 2014-03-13 Method for correcting interval velocity by using residue curvature

Publications (2)

Publication Number Publication Date
CN103869369A CN103869369A (en) 2014-06-18
CN103869369B true CN103869369B (en) 2017-02-15

Family

ID=50908102

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410092817.6A Active CN103869369B (en) 2014-03-13 2014-03-13 Method for correcting interval velocity by using residue curvature

Country Status (1)

Country Link
CN (1) CN103869369B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101315427A (en) * 2007-05-29 2008-12-03 中国石油天然气集团公司 Method and system for processing seismic exploration data of complex area
US8363509B2 (en) * 2006-09-04 2013-01-29 Daniele Colombo Method for building velocity models for pre-stack depth migration via the simultaneous joint inversion of seismic, gravity and magnetotelluric data
CN102901985A (en) * 2011-07-25 2013-01-30 中国石油化工股份有限公司 Depth domain layer speed correcting method suitable for undulating surface

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8363509B2 (en) * 2006-09-04 2013-01-29 Daniele Colombo Method for building velocity models for pre-stack depth migration via the simultaneous joint inversion of seismic, gravity and magnetotelluric data
EP2062071B1 (en) * 2006-09-04 2014-10-22 Geosystem S.r.l. Method for building velocity models for pre-stack depth migration via the simultaneous joint inversion of seismic, gravity and magnetotelluric data
CN101315427A (en) * 2007-05-29 2008-12-03 中国石油天然气集团公司 Method and system for processing seismic exploration data of complex area
CN102901985A (en) * 2011-07-25 2013-01-30 中国石油化工股份有限公司 Depth domain layer speed correcting method suitable for undulating surface

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Velocity analysis by iterative profile migration;Kamal AI-Yahya;《Geophysics》;19890630;第54卷(第6期);第718-729页 *
南海北部地质结构复杂区域地层速度求取研究;周洪生;《中国优秀硕士学位论文全文数据库·基础科学辑》;20131215(第12期);A011-223 *
时间域剩余曲率偏移速度分析技术在iCluster软件中的实现;张兵 等;《勘探地球物理进展》;20100630;第33卷(第3期);第168-173页 *
角度域共成像点道集剩余曲率法偏移速度分析;张凯 等;《应用地球物理》;20100331;第7卷(第1期);第49-56页 *

Also Published As

Publication number Publication date
CN103869369A (en) 2014-06-18

Similar Documents

Publication Publication Date Title
CN104133245B (en) The static correcting method and system of a kind of seismic data
CN104914473B (en) A kind of inversion method and device of transient electromagnetic resistivity
Marto et al. A correlation of shear wave velocity and standard penetration resistance
CN105386756B (en) A method of brittle formation porosity is calculated using dependent variable
CN103837895B (en) Matching preliminary wave obtains the method for short-period static corrections
CN108089227B (en) Novel stratum pore pressure prediction method based on three-dimensional seismic data
CN103412333B (en) A kind of static correction base level defining method
CN102901985B (en) A kind of Depth Domain interval velocity modification method being applicable to relief surface
CN104483705A (en) Three-dimensional residual static correction method
CN105607146A (en) Quantitative characterization method of meandering river sand body scale
CN109375253A (en) Ground motion parameter evaluation method based on whole seismic structure maximum credible earthquakes
CN104375181B (en) A kind of method quickly generating three-dimensional common receiver first arrival road collection
CN103913770B (en) Method geological data processed based on VSP data
CN110197318A (en) Surface mine Grade Model update method
CN105093318B (en) A kind of adaptive wave equation wave field extrapolation static correcting method
CN103984016B (en) Converted wave anisotropy Amplitudeversusangle road collection abstracting method
CN106569259A (en) Regularized tomographic velocity inversion method and device based on structural inclination angle
CN104914469B (en) Static correcting method and device in a kind of converted shear wave
CN105093280B (en) Surface-level model is to the low frequency of earthquake data influence and the decomposition method of radio-frequency component
CN105242315A (en) Geologic structure imagining method and device
CN103869369B (en) Method for correcting interval velocity by using residue curvature
WO2015171215A1 (en) Efficient line search methods for multi-parameter full wavefield inversion
CN104459800A (en) Sand body pinch-out prediction method and device
CN105116466B (en) Method and device for determining oil field reservoir physical property characteristics
US20190106974A1 (en) Systems and methods for horizontal well geosteering

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20180202

Address after: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189

Patentee after: BGP INC., CHINA NATIONAL PETROLEUM Corp.

Address before: No. 216, No. 216, Huayang Avenue, Huayang Town, Shuangliu County, Shuangliu County, Sichuan

Patentee before: CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20200916

Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No.

Co-patentee after: BGP Inc., China National Petroleum Corp.

Patentee after: CHINA NATIONAL PETROLEUM Corp.

Address before: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189

Patentee before: BGP Inc., China National Petroleum Corp.