CN106405643A - Well-to-seismic integration velocity field construction method and device - Google Patents
Well-to-seismic integration velocity field construction method and device Download PDFInfo
- Publication number
- CN106405643A CN106405643A CN201610767702.1A CN201610767702A CN106405643A CN 106405643 A CN106405643 A CN 106405643A CN 201610767702 A CN201610767702 A CN 201610767702A CN 106405643 A CN106405643 A CN 106405643A
- Authority
- CN
- China
- Prior art keywords
- point
- time
- well
- velocity
- contact
- 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
- 238000010276 construction Methods 0.000 title claims abstract description 27
- 230000010354 integration Effects 0.000 title abstract 4
- 238000012937 correction Methods 0.000 claims abstract description 51
- 230000015572 biosynthetic process Effects 0.000 claims description 79
- 238000000034 method Methods 0.000 claims description 32
- 238000006243 chemical reaction Methods 0.000 claims description 12
- 238000005553 drilling Methods 0.000 claims description 11
- 238000005516 engineering process Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000013500 data storage Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000011514 reflex Effects 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Abstract
The invention discloses a well-to-seismic integration velocity field construction method and device. The well-to-seismic integration velocity field construction method comprises the following steps that (A) seismic stack velocity of each time point on the common depth point (CDP) gather is converted into first stratum velocity of each time point through a generalized Dix formula algorithm; (B) the correction coefficient of each time point is acquired by using logging information and the first stratum velocity of each time point; and (C) the correction coefficient of each time point is multiplied with the first stratum velocity of each time point so that well-to-seismic integration velocity of each time point is acquired.
Description
Technical field
The present invention relates to seismic data process field, more particularly, be related to a kind of by seismic velocity and well point speed phase
In conjunction with well shake joint velocity field construction method.
Background technology
The complex area velocity field construction method commonly used at present has two kinds：One is using drilling depth conversion average speed,
Then velocity field structure is carried out with the variation tendency of depth according to speed；Two is to be carried out using the speed being converted to of normal-moveout spectrum
Velocity field builds.Above two method respectively has shortcoming, and the former can not effectively utilize many well informations and seismic velocity information, from
The precision of constraint well location speed farther out is very low；The latter only using earthquake average speed desin speed field, due to seismic velocity
Influence factor a lot, and precision is relatively low, is difficult to obtain accurate structural map.
Content of the invention
Therefore, it is an object of the invention to provide one kind can effectively utilize many well informations and seismic velocity information thus
Build the well shake joint velocity field construction method of accurate velocity field.
According to an aspect of the present invention, provide a kind of well shake joint velocity field construction method, comprise the following steps：(A) lead to
Cross broad sense Dix formula algorithm by the earthquake overlap rate conversion of the Each point in time in common depth point (CDP) road collection be each when
Between point the first formation velocity；(B) well-log information and the first formation velocity of Each point in time is utilized to obtain Each point in time
Correction coefficient；(C) pass through the correction coefficient of Each point in time is multiplied with the first formation velocity of Each point in time, obtain each
The well shake joint speed of time point.
Preferably, step (A) includes：(a1) Each point in time in CDP road collection and corresponding earthquake overlap speed are read
Degree；(a2) for first time point in CDP road collection, will earthquake overlap speed corresponding with described time point as described when
Between point the first formation velocity；(a3) for any time point in the other times point in CDP road collection, with described any time
The time t of the previous time point of pointoldAnd its corresponding earthquake overlap speed VoldTime t with described any time pointnewAnd its
Corresponding earthquake overlap speed VnewSet up the first Hyperbolic EquationWith the second Hyperbolic Equation(a4) choose specified point on the second Hyperbolic Equation as the second point of contact, and ask for the second point of contact
Tangent slope K2；(a5) translate described tangent line tangent with the first Hyperbolic Equation, so that it is determined that the first tangent point of contact；(a6)
Calculated the slope K of the straight line at the first point of contact and the second point of contact1；(a7) pass through equationWhen calculating described arbitrary
Between point the first formation velocity.
Preferably, step (a4) includes：Equally spaced choose multiple points on the second Hyperbolic Equation as multiple second
Point of contact, and asked for a plurality of tangent slope K at the plurality of second point of contact2i, wherein, i=1 ..., N, N are whole more than 2
Number；Wherein, step (a5) includes：Translate described a plurality of tangent line tangent with the first Hyperbolic Equation, so that it is determined that tangent is multiple
First point of contact；Wherein, step (a6) includes：Calculated each first point of contact in the plurality of first point of contact and the plurality of
The slope K of the straight line at corresponding second point of contact in two point of contacts1i；Wherein, step (a7) includes：By equationCalculate multiple speed, and by calculating the meansigma methodss of the plurality of speed as described any time point
First formation velocity.
Preferably, step (B) includes：(b1) the second formation velocity of each well point is asked for using drilling data；(b2) obtain
Take the first formation velocity of time point corresponding with each well point in the first formation velocity of conversion in step (A), as each
First formation velocity of well point；(b3) pass through the first formation velocity phase of the second formation velocity of each well point and each well point
Remove, obtain the correction coefficient of each well point；(b4) by row interpolation is entered to the correction coefficient of each well point, obtain Each point in time
Correction coefficient.
Preferably, in step (b4), using Kriging regression method, row interpolation is entered to the correction coefficient of each well point.
According to a further aspect in the invention, a kind of well shake joint velocity field construction device is provided, including：Stack velocity is changed
Unit, is configured to broad sense Dix formula algorithm the earthquake overlap of the Each point in time in common depth point (CDP) road collection is fast
Degree is converted to the first formation velocity of Each point in time；Correction coefficient acquiring unit, be configured to, with well-log information and each
First formation velocity of time point obtains the correction coefficient of Each point in time；Well shakes joint speed acquiring unit, is configured to lead to
Cross and the correction coefficient of Each point in time is multiplied with the first formation velocity of Each point in time, obtain the well shake connection of Each point in time
Sum velocity.
Preferably, stack velocity converting unit is configured to：Read the Each point in time and accordingly in CDP road collection
Shake stack velocity；For first time point in CDP road collection, will earthquake overlap speed corresponding with described time point as institute
State the first formation velocity of time point；For any time point in the other times point in CDP road collection, with described any time
The time t of the previous time point of pointoldAnd its corresponding earthquake overlap speed VoldTime t with described any time pointnewAnd its
Corresponding earthquake overlap speed VnewSet up the first Hyperbolic EquationWith the second Hyperbolic EquationChoose specified point on the second Hyperbolic Equation as the second point of contact, and ask for the tangent line at the second point of contact
Slope K2；Translate described tangent line tangent with the first Hyperbolic Equation, so that it is determined that the first tangent point of contact；Calculated first to cut
The slope K of the straight line at point and the second point of contact1；By equationCalculate the first stratum speed of described any time point
Degree.
Preferably, stack velocity converting unit is additionally configured to：Equally spaced choose multiple on the second Hyperbolic Equation
Point is as multiple second point of contacts, and asks for a plurality of tangent slope K at the plurality of second point of contact2i, wherein, i=1 ...,
N, N are the integer more than 2；Translate that described a plurality of tangent line is tangent with the first Hyperbolic Equation, so that it is determined that tangent multiple first
Point of contact；Corresponding second calculating in each first point of contact and the plurality of second point of contact in the plurality of first point of contact cuts
The slope K of the straight line of point1i；By equationCalculate multiple speed, and by calculating the plurality of speed
Meansigma methodss are as the first formation velocity of described any time point.
Preferably, correction coefficient acquiring unit is configured to：Ask for the second stratum speed of each well point using drilling data
Degree；Obtain the first stratum speed of time point corresponding with each well point in the first formation velocity of stack velocity converting unit conversion
Degree, as the first formation velocity of each well point；By by the first ground of the second formation velocity of each well point and each well point
Interval velocity is divided by, and obtains the correction coefficient of each well point；By row interpolation is entered to the correction coefficient of each well point, when obtaining each
Between point correction coefficient.
Preferably, correction coefficient acquiring unit is configured with Kriging regression method the correction coefficient of each well point is entered
Row interpolation.
Brief description
By the description carrying out below in conjunction with accompanying drawing, the above and other aspect of the exemplary embodiment of the present invention, feature
Will be more readily apparent from advantage, in the accompanying drawings：
Fig. 1 is to illustrate the flow chart that well shakes joint velocity field construction method according to an embodiment of the invention；
Fig. 2 is to illustrate using broad sense Dix formula algorithm, stack velocity to be converted to stratum according to an embodiment of the invention
The flow chart of the method for speed；
Fig. 3 is the method illustrating to ask for correction coefficient using well-log information formation velocity according to an embodiment of the invention
Flow chart；
Fig. 4 is the block diagram illustrating well shake joint velocity field construction device according to an embodiment of the invention.
In the accompanying drawings, identical label will be understood to refer to identical element, feature and structure.
Specific embodiment
The present invention describing to help comprehensive understanding to be limited by claim and its equivalent referring to the drawings is provided
Exemplary embodiment.Description referring to the drawings includes various specific detail to help understand, but described specific detail
It only is seen as exemplary.Therefore, it will be appreciated by those of ordinary skill in the art that without departing from scope and spirit of the present invention
In the case of, the embodiments described herein can be made various changes and modifications.Additionally, in order to clear and brief, public affairs can be omitted
Know the description of function and structure.
Used in following described and claimed, term and word are not limited to literal meaning, but are only used by inventor
Enable to clear and as one man understand the present invention.Therefore, it should be appreciated by the person skilled in the art that providing the present invention's
The following description of exemplary embodiment is only descriptive purpose, rather than limits this that limited by claim and its equivalent
Bright purpose.
Hereinafter, describe first and ask for by using Dix formula (hereinafter referred to as narrow sense Dix formula) commonly used in the art
The method of formation velocity.Shown in narrow sense Dix formula such as below equation (1).
Wherein, vjFor the formation velocity of jth layer, T0,jFor the double-pass reflection time of jth layer, vrms,jRoot-mean-square for jth layer
Speed.The supposed premise of narrow sense Dix formula is medium is horizontal layer, and seismic wave vertical interface is incident, and (i.e. angle of incidence is equal to
Zero, ray parameter P=0).In actual applications, root mean sequare velocity v in equation (1)rmsStack velocity v can only be usedmcsReplace
Interval velocity v can be calculatedj.However, under the premise of being assumed as above, vrmsWith vmcsBetween error larger.Further, medium is water
The practical situation that flat bed shape is constructed with complicated geology has greater difference, and seismic ray angle of incidence less (geophone offset is less) and height
Precision stack velocity calculates the cannon offset (seismic ray angle of incidence is larger) requiring and just constitutes a pair implacable contradiction.
As can be seen here, asking for formation velocity by using narrow sense Dix formula has great limitation.
According to embodiments of the invention, for a horizontal interface, be can determine as follows by narrow sense Dix formula
Reflection time curve.
Wherein, t0It it is the round trip vertical reflection time of this horizontal interface, x is receiving point and shot point on this horizontal interface
The distance between, t is the time receiving echo at x, vrmsRoot mean sequare velocity for this horizontal interface.
So, in t2-x2In coordinate system, reflection time curve as implied above is expressed as straight line, and its slope can pass through
Below equation determines.
Above equation shows t2-x2The slope of straight line be root mean sequare velocity square inverse, thus inspire and can define all
Root speed is t2-x2The inverse of straight slope.Therefore, the rewritable equation (2) being as follows of above equation.
If in equation (2) formula, t is certain in the layer reflex time, the level that x propagates for seismic wave in this layer away from
From then the root mean sequare velocity of at this moment this layer can be changed into the formation velocity as represented by below equation (3).
Wherein, Δ xj=xj+1-xj,ΔTj=Tj+1-Tj.
Dx/dt is to be 1 second upper and lower two reflection line-ups the time in coordinate points (xj,Tj) and (xj+1,Tj+1) go up and reflect together
The tangent slope of phase axle, Δ xj/ΔTjIt is the straight slope connecting two point of contacts.According to push over process above it was determined that according to
The ultimate principle of the Slope Method for calculating formation velocity of embodiments of the invention is that top bottom reflection line-ups two are parallel cuts
Line slope and the evolution of the straight slope product being connected two point of contacts, are the formation velocity of this layer.According to embodiments of the invention,
Due to obtaining above-mentioned Slope Method by improving to narrow sense Dix formula, therefore Slope Method can be referred to as broad sense Dix formula and calculate
Method.Below with reference to accompanying drawing, this is explained in more detail.
Fig. 1 is to illustrate the flow chart that well shakes joint velocity field construction method according to an embodiment of the invention.
With reference to Fig. 1, in step S101, during by broad sense Dix formula algorithm by each in common depth point (CDP) road collection
Between point earthquake overlap rate conversion be Each point in time formation velocity.
Fig. 2 is to illustrate using broad sense Dix formula algorithm, stack velocity to be converted to stratum according to an embodiment of the invention
The flow chart of the method for speed.
With reference to Fig. 2, in step s 201, read the time point in CDP road collection and corresponding earthquake overlap speed.According to
Embodiments of the invention, last time point inverted order from CDP road collection can read Each point in time and accordingly
Shake stack velocity is it is also possible to first time dot sequency from CDP road collection reads Each point in time and corresponding earthquake is folded
Acceleration.However, the invention is not restricted to this.For example, it is possible to randomly read time point and the corresponding earthquake in CDP road collection
Stack velocity.The time of the Each point in time in CDP road collection and corresponding earthquake overlap speed can be passed through of the prior art
Various methods obtain, and repeat no more here.In step S202, determine that whether the time point reading is the in this CDP road collection
One time point.If the time point reading is first time point in this CDP road collection, in step S203, will be with institute
State time point corresponding earthquake overlap speed as the formation velocity of described time point.Subsequently, return to step S201, reads another
Individual time point and corresponding earthquake overlap speed.If however, when the time point reading is not first in this CDP road collection
Between point, then in step S204, with the time t of the previous time point of described time pointoldAnd its corresponding earthquake overlap speed
VoldTime t with described time pointnewAnd its corresponding earthquake overlap speed VnewSet up following two Hyperbolic Equations (4) and
(5).
Next, in step S205, a bit (X, T) in selection the second Hyperbolic Equation (5) is as point of contact, and asks for
Cross the tangent slope K at this point of contact2.In step S206, the described tangent line of translation is tangent with the first Hyperbolic Equation (4), thus
Determine tangent point of contact (X', T').In step S207, calculated the slope of the straight line of point of contact (X, T) and point of contact (X', T')
K1.In step S208, by equationCalculate the formation velocity of described time point.In step S209, really
Determine whether all time points in CDP road collection to be processed.If clicked through to all times in CDP road collection
Go process, then the method terminating for stack velocity to be converted to formation velocity.Without to all time points in CDP road collection
Processed, then return to step S201 is to read next time point.
According to embodiments of the invention, for not being the time point of first time point in CDP road collection, in step
In S205, the multiple point of contacts of the upper selection of the second Hyperbolic Equation (5) can be equally spaced on, and asked for the many of the plurality of point of contact
Bar tangent slope K2i.Here, i=1 ..., N, N are the integer more than 2.Then, in step S206, translation is described a plurality of
Tangent line is tangent with the first Hyperbolic Equation, so that it is determined that the multiple point of contacts on the first Hyperbolic Equation (4), and in step S207
In, calculated multiple on each point of contact and the second Hyperbolic Equation (5) in the multiple point of contacts on the first Hyperbolic Equation (4)
The slope K of the straight line at corresponding point of contact in point of contact1i.Finally, in step S207, by equationCalculate
Multiple speed, ask for the meansigma methodss of the plurality of speed, and using the meansigma methodss asked for as described time point formation velocity.
Referring back to Fig. 1, in step s 102, when obtaining each using the formation velocity of well-log information and Each point in time
Between point correction coefficient.The method being more fully described the correction coefficient obtaining Each point in time hereinafter with reference to Fig. 3.
Fig. 3 is the method illustrating to ask for correction coefficient using well-log information formation velocity according to an embodiment of the invention
Flow chart.
With reference to Fig. 3, in step S301, ask for the formation velocity of each well point using drilling data.Here, can be by existing
The various methods in technology are had to ask for the formation velocity of each well point using drilling data.For example, can be by drilling well be obtained
It is divided by obtain the formation velocity of well point reflex time in the stratum obtaining in formation thickness and geological data.In step s 302,
Obtain the formation velocity of time point corresponding with each well point in the formation velocity of conversion in step S101, as each well point
Formation velocity.In other words, in step s 302, the time in CDP road collection corresponding with each well point can be determined first
Point, then using the formation velocity of the corresponding time point being converted in step S101 as each well point formation velocity.
Next, in step S303, by by the formation velocity of each well point asked in step S301 with step s 302
The formation velocity of each well point obtaining is divided by, and to obtain the correction coefficient of each well point.Finally, in step s 304, pass through
Row interpolation is entered to the correction coefficient of each well point, obtains the correction coefficient of Each point in time.According to embodiments of the invention, can make
With Kriging regression method, row interpolation is entered to the correction coefficient of each well point.However, the invention is not restricted to this, it is possible to use existing skill
Various interpolation methods in art enter row interpolation to the correction coefficient of each well point.
Referring back to Fig. 1, in step s 103, by by the stratum of the correction coefficient of Each point in time and Each point in time
Speed is multiplied, and obtains the well shake joint speed of Each point in time.
Fig. 4 is the block diagram illustrating well shake joint velocity field construction device according to an embodiment of the invention.
With reference to Fig. 4, joint velocity field construction device 400 may include stack velocity converting unit 401, correction coefficient obtains for well shake
Take unit 402 and well shake joint speed acquiring unit 403.According to embodiments of the invention, well shakes joint velocity field construction device
400 can be realized by various computing devices (for example, computer, server, work station etc.).Specifically, stack velocity turns
Changing unit 401 can be by broad sense Dix formula algorithm by the earthquake overlap speed of the Each point in time in common depth point (CDP) road collection
Be converted to the formation velocity of Each point in time.Correction coefficient acquiring unit 402 can be using the ground of well-log information and Each point in time
Interval velocity obtains the correction coefficient of Each point in time.Well shake joint speed acquiring unit 403 can pass through the school of Each point in time
Positive coefficient is multiplied with the formation velocity of Each point in time, obtains the well shake joint speed of Each point in time.
As described above, stack velocity converting unit 401 can read Each point in time and corresponding earthquake in CDP road collection
Stack velocity.For first time point in CDP road collection, stack velocity converting unit 401 can will be corresponding to described time point
Earthquake overlap speed as described time point formation velocity.For any time in the other times point in CDP road collection
Point, stack velocity converting unit 401 can be with the time t of the previous time point of described any time pointoldAnd its earthquake accordingly
Stack velocity VoldTime t with described any time pointnewAnd its corresponding earthquake overlap speed VnewSet up the first hyperbola side
JourneyWith the second Hyperbolic EquationAfter setting up two Hyperbolic Equations, stack velocity
Converting unit 401 can choose a bit (X, T) on the second Hyperbolic Equation as point of contact, and asked for this point of contact tangent line oblique
Rate K2.Then, the translatable described tangent line of stack velocity converting unit 401 is tangent with the first Hyperbolic Equation, so that it is determined that tangent
Point of contact (X', T'), and then calculated the slope K of the straight line of point of contact (X, T) and point of contact (X', T')1.Finally, stack velocity turns
Change unit 401 and can pass through equationCalculate the formation velocity of described any time point.
Selectively, for any time point not being first time point in CDP road collection, stack velocity converting unit
401 can be equally spaced on the second Hyperbolic Equation and choose multiple point of contacts, and ask for a plurality of tangent line at the plurality of point of contact
Slope K2i.Here, i=1 ..., N, N are the integer more than 2.Then, stack velocity converting unit 401 is translatable described many
Bar tangent line is tangent with the first Hyperbolic Equation, so that it is determined that the multiple point of contacts on the first Hyperbolic Equation, and calculated first
The corresponding point of contact in the multiple point of contacts on each point of contact and the second Hyperbolic Equation in multiple point of contacts on Hyperbolic Equation
Straight line slope K1i.Finally, stack velocity converting unit 401 can pass through equationCalculate multiple speed,
Ask for the meansigma methodss of the plurality of speed, and using the meansigma methodss asked for as described any time point formation velocity.
Additionally, correction coefficient acquiring unit 402 can ask for the formation velocity of each well point first with drilling data, then
Obtain the formation velocity of time point corresponding with each well point in the formation velocity of stack velocity converting unit conversion.Then, school
Positive coefficient acquiring unit 402 can be by the formation velocity of each well point that will be asked for using drilling data and the stack velocity obtaining
In the formation velocity of converting unit 401 conversion, the formation velocity of time point corresponding with each well point is divided by, and to obtain each well
The correction coefficient of point.Finally, (for example, correction coefficient acquiring unit 402 can be by entering row interpolation to the correction coefficient of each well point
By using Kriging regression method, row interpolation is entered to the correction coefficient of each well point), obtain the correction coefficient of Each point in time.
Well shake joint velocity field construction method can effectively utilize many well informations and ground according to an embodiment of the invention
Shake velocity information, is applied to multiplexed address zone velocity field and builds, thus the mapping precision of structural map and becoming of drilling effectively
Power.
Well shake joint velocity field construction method can achieve as computer readable recording medium storing program for performing according to an embodiment of the invention
On computer-readable code, or can be sent by transmission medium.Computer readable recording medium storing program for performing is can to store hereafter may be used
The arbitrary data storage device of the data being read by computer system.The example of computer readable recording medium storing program for performing includes read-only storage
Device (ROM), random access memory (RAM), CD (CD)-ROM, digital versatile disc (DVD), tape, floppy disk, optical data
Storage device, but not limited to this.Transmission medium may include the carrier wave sending by network or various types of communication port.Calculate
Machine readable medium recording program performing also can be distributed in the computer system connecting network, thus computer-readable code is deposited in a distributed fashion
Storage and execution.
Although the certain exemplary embodiments with reference to the present invention have been shown and described the present invention, art technology
Personnel will be understood that, in the case of without departing from the spirit and scope of the present invention being limited by claim and its equivalent, can enter
Various changes on row various forms and details.
Claims (10)
1. a kind of well shakes joint velocity field construction method, comprises the following steps：
(A) pass through broad sense Dix formula algorithm by the earthquake overlap rate conversion of the Each point in time in common depth point (CDP) road collection
The first formation velocity for Each point in time；
(B) utilize the correction coefficient of the first formation velocity acquisition Each point in time of well-log information and Each point in time；
(C) pass through the correction coefficient of Each point in time is multiplied with the first formation velocity of Each point in time, obtain each time
The well shake joint speed of point.
2. well as claimed in claim 1 shakes joint velocity field construction method, and wherein, step (A) includes：
(a1) read Each point in time and the corresponding earthquake overlap speed in CDP road collection；
(a2) for first time point in CDP road collection, will earthquake overlap speed corresponding with described time point as described
First formation velocity of time point；
(a3) for any time point in the other times point in CDP road collection, with the previous time point of described any time point
Time toldAnd its corresponding earthquake overlap speed VoldTime t with described any time pointnewAnd its corresponding earthquake overlap
Speed VnewSet up the first Hyperbolic EquationWith the second Hyperbolic Equation
(a4) choose specified point on the second Hyperbolic Equation as the second point of contact, and ask for the tangent slope at the second point of contact
K2；
(a5) translate described tangent line tangent with the first Hyperbolic Equation, so that it is determined that the first tangent point of contact；
(a6) calculated the slope K of the straight line at the first point of contact and the second point of contact1；
(a7) pass through equationCalculate the first formation velocity of described any time point.
3. well as claimed in claim 2 shakes joint velocity field construction method, and wherein, step (a4) includes：Equally spaced choose
Multiple points on second Hyperbolic Equation are as multiple second point of contacts, and ask for a plurality of tangent line at the plurality of second point of contact
Slope K2i, wherein, i=1 ..., N, N are the integer more than 2；
Wherein, step (a5) includes：Translate described a plurality of tangent line tangent with the first Hyperbolic Equation, so that it is determined that tangent is multiple
First point of contact；
Wherein, step (a6) includes：Calculated each first point of contact in the plurality of first point of contact and the plurality of second to cut
The slope K of the straight line at corresponding second point of contact in point1i；
Wherein, step (a7) includes：By equationCalculate multiple speed, and by calculating the plurality of speed
The meansigma methodss of degree are as the first formation velocity of described any time point.
4. well as claimed in claim 1 shakes joint velocity field construction method, and wherein, step (B) includes：
(b1) the second formation velocity of each well point is asked for using drilling data；
(b2) in the first formation velocity of conversion in obtaining step (A), the first stratum of time point corresponding with each well point is fast
Degree, as the first formation velocity of each well point；
(b3) by the first formation velocity of the second formation velocity of each well point and each well point is divided by, obtain each well point
Correction coefficient；
(b4) by the correction coefficient of each well point is entered with row interpolation, obtain the correction coefficient of Each point in time.
5. well as claimed in claim 4 shakes joint velocity field construction method, wherein, in step (b4), is inserted using Ke Lijin
Value method enters row interpolation to the correction coefficient of each well point.
6. a kind of well shakes joint velocity field construction device, including：
Stack velocity converting unit, when being configured to broad sense Dix formula algorithm by each in common depth point (CDP) road collection
Between point earthquake overlap rate conversion be Each point in time the first formation velocity；
Correction coefficient acquiring unit, when being configured to, with the first formation velocity of well-log information and Each point in time and obtaining each
Between point correction coefficient；
Well shakes joint speed acquiring unit, is configured to the correction coefficient of Each point in time and the first of Each point in time
Formation velocity is multiplied, and obtains the well shake joint speed of Each point in time.
7. well as claimed in claim 6 shakes joint velocity field construction device, and wherein, stack velocity converting unit is configured to：
Read the Each point in time in CDP road collection and corresponding earthquake overlap speed；For first time point in CDP road collection,
Using earthquake overlap speed corresponding with described time point as described time point the first formation velocity；For in CDP road collection
Any time point in other times point, with the time t of the previous time point of described any time pointoldAnd its earthquake accordingly
Stack velocity VoldTime t with described any time pointnewAnd its corresponding earthquake overlap speed VnewSet up the first hyperbola side
JourneyWith the second Hyperbolic EquationChoose the specified point conduct on the second Hyperbolic Equation
Second point of contact, and asked for the tangent slope K at the second point of contact2；Translate described tangent line tangent with the first Hyperbolic Equation, thus
Determine the first tangent point of contact；Calculated the slope K of the straight line at the first point of contact and the second point of contact1；By equationCalculate the first formation velocity of described any time point.
8. well as claimed in claim 7 shakes joint velocity field construction device, and wherein, stack velocity converting unit is also configured
For：Equally spaced choose multiple points on the second Hyperbolic Equation as multiple second point of contacts, and ask for the plurality of second
The a plurality of tangent slope K at point of contact2i, wherein, i=1 ..., N, N are the integer more than 2；Translate described a plurality of tangent line and first
Hyperbolic Equation is tangent, so that it is determined that tangent multiple first point of contacts；Calculated in the plurality of first point of contact each first
The slope K of the straight line at corresponding second point of contact in point of contact and the plurality of second point of contact1i；By equation
Calculate multiple speed, and by calculating first formation velocity as described any time point for the meansigma methodss of the plurality of speed.
9. well as claimed in claim 6 shakes joint velocity field construction device, and wherein, correction coefficient acquiring unit is configured to：
Ask for the second formation velocity of each well point using drilling data；Obtain the first formation velocity of stack velocity converting unit conversion
In time point corresponding with each well point the first formation velocity, as the first formation velocity of each well point；By by each
Second formation velocity of well point is divided by with the first formation velocity of each well point, obtains the correction coefficient of each well point；By right
The correction coefficient of each well point enters row interpolation, obtains the correction coefficient of Each point in time.
10. well as claimed in claim 9 shakes joint velocity field construction device, and wherein, correction coefficient acquiring unit is configured to
Using Kriging regression method, row interpolation is entered to the correction coefficient of each well point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610767702.1A CN106405643B (en) | 2016-08-30 | 2016-08-30 | Well shakes joint velocity field construction method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610767702.1A CN106405643B (en) | 2016-08-30 | 2016-08-30 | Well shakes joint velocity field construction method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106405643A true CN106405643A (en) | 2017-02-15 |
| CN106405643B CN106405643B (en) | 2018-09-25 |
Family
ID=58003648
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610767702.1A Active CN106405643B (en) | 2016-08-30 | 2016-08-30 | Well shakes joint velocity field construction method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106405643B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107315194A (en) * | 2017-05-24 | 2017-11-03 | 中国石油天然气股份有限公司 | The determination method and apparatus of seismic pre-stack time migration speed |
| CN109188522A (en) * | 2018-10-09 | 2019-01-11 | 中国石油天然气股份有限公司 | Velocity field construction method and device |
| CN109459790A (en) * | 2017-09-06 | 2019-03-12 | 中国石油化工股份有限公司 | For coal measure strata seismic velocity field method for building up and system |
| CN111060969A (en) * | 2019-12-25 | 2020-04-24 | 恒泰艾普（北京）能源科技研究院有限公司 | Well control Q compensation method |
| CN111208558A (en) * | 2018-11-22 | 2020-05-29 | 中国石油天然气股份有限公司 | Method and device for establishing ultra-deep low-amplitude three-dimensional geological structure |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5663928A (en) * | 1994-10-19 | 1997-09-02 | Elf Aquitaine Production | Method for analysing and processing seismic reflection data for the determination of a high resolution spatial velocity field for hyperbolicity correction |
| CN102129084A (en) * | 2010-12-17 | 2011-07-20 | 中国石油天然气股份有限公司 | Method and device for acquiring seismic velocity in thin reservoir layer through well control |
| CN102967882A (en) * | 2012-11-16 | 2013-03-13 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for building layer velocity model of stratum |
| CN104122582A (en) * | 2014-06-24 | 2014-10-29 | 夏正元 | Method for calculating high-accuracy seismic wave velocity by utilizing stacking velocity |
| CN105842736A (en) * | 2016-05-27 | 2016-08-10 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for building stratum velocity model |
-
2016
- 2016-08-30 CN CN201610767702.1A patent/CN106405643B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5663928A (en) * | 1994-10-19 | 1997-09-02 | Elf Aquitaine Production | Method for analysing and processing seismic reflection data for the determination of a high resolution spatial velocity field for hyperbolicity correction |
| CN102129084A (en) * | 2010-12-17 | 2011-07-20 | 中国石油天然气股份有限公司 | Method and device for acquiring seismic velocity in thin reservoir layer through well control |
| CN102967882A (en) * | 2012-11-16 | 2013-03-13 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for building layer velocity model of stratum |
| CN104122582A (en) * | 2014-06-24 | 2014-10-29 | 夏正元 | Method for calculating high-accuracy seismic wave velocity by utilizing stacking velocity |
| CN105842736A (en) * | 2016-05-27 | 2016-08-10 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for building stratum velocity model |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107315194A (en) * | 2017-05-24 | 2017-11-03 | 中国石油天然气股份有限公司 | The determination method and apparatus of seismic pre-stack time migration speed |
| CN107315194B (en) * | 2017-05-24 | 2019-04-12 | 中国石油天然气股份有限公司 | The determination method and apparatus of seismic pre-stack time migration speed |
| CN109459790A (en) * | 2017-09-06 | 2019-03-12 | 中国石油化工股份有限公司 | For coal measure strata seismic velocity field method for building up and system |
| CN109188522A (en) * | 2018-10-09 | 2019-01-11 | 中国石油天然气股份有限公司 | Velocity field construction method and device |
| CN109188522B (en) * | 2018-10-09 | 2020-02-14 | 中国石油天然气股份有限公司 | Velocity field construction method and device |
| CN111208558A (en) * | 2018-11-22 | 2020-05-29 | 中国石油天然气股份有限公司 | Method and device for establishing ultra-deep low-amplitude three-dimensional geological structure |
| CN111060969A (en) * | 2019-12-25 | 2020-04-24 | 恒泰艾普（北京）能源科技研究院有限公司 | Well control Q compensation method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106405643B (en) | 2018-09-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106405643A (en) | Well-to-seismic integration velocity field construction method and device | |
| WO2017024702A1 (en) | Inversion system for ray elastic parameter | |
| CN104375188A (en) | Seismic wave transmission attenuation compensation method and device | |
| CN103365916A (en) | Earthquake event parameter estimation and acquirement method and system and earthquake event search engine | |
| CN110389382B (en) | Oil and gas reservoir characterization method based on convolutional neural network | |
| CN105093319A (en) | Ground micro-seismic static correction method based on three-dimensional seismic data | |
| CN105242315B (en) | A kind of method and device of geological structure imaging | |
| CN104865597A (en) | Modeling method of depth domain interval velocity initial model | |
| CN115373024A (en) | Method and device for inverting passive land edge crustal structure based on stratum recording settlement | |
| CN104570122A (en) | Seismic data static correction method and device based on reflected waves | |
| NO20101734A1 (en) | Speed models for a single well and for a set of wells | |
| CN102866422A (en) | Depth domain geological entity model generation method | |
| CN108594300B (en) | Landform imaging method and device and computer storage medium | |
| CN105717538B (en) | Relief surface earthquake data offset reference plane conversion method and device | |
| CN103777242A (en) | Speed discrimination method with combination of depth focusing and gather event flattening | |
| CN102053275B (en) | Method for calculating relative statics correction value of combination in single-point earthquake chamber | |
| CN106443829B (en) | A kind of near-surface model construction method and device | |
| CN111208558B (en) | Method and device for establishing ultra-deep low-amplitude three-dimensional geological structure | |
| CN105891887B (en) | Speed based on superposition of data vertically and horizontally high density assay method | |
| CN104181610A (en) | Ray path elasticity inversion method and system | |
| CN103901475A (en) | Method and device for drawing attribute contour map | |
| CN114841586A (en) | Oil gas scenic spot delineation method and system | |
| CN107544092A (en) | Data processing method and device in vertical seismic profiling (VSP) collection | |
| Biasi et al. | Disaggregating UCERF3 for Site-Specific Applications | |
| CN104199107B (en) | Depth prediction approach and system before brill based on the many wave datum of vertical seismic |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: No. 216, No. 216, Huayang Avenue, Tianfu New District, Sichuan, Sichuan Applicant after: China National Petroleum Corporation Chuanqing Drilling Engineering Geophysical Exploration Company Ltd. Address before: Shuangliu County Huayang Huayang Road in Chengdu city of Sichuan Province in 610213 section of No. 216, Igawa geophysical exploration company of the Ministry of science and technology Applicant before: China National Petroleum Corporation Chuanqing Drilling Engineering Geophysical Exploration Company Ltd. |
|
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20180514 Address after: 100007 Dongzhimen North Street, Dongcheng District, Dongcheng District, Beijing Applicant after: China Petroleum and Natural Gas Group Co., Ltd. Applicant after: Dongfang Geophysical Exploration Co., Ltd., China Petrochemical Corp. Address before: 610213 No. 216, Huayang Road, Tianfu New District, Chengdu, Sichuan Applicant before: China National Petroleum Corporation Chuanqing Drilling Engineering Geophysical Exploration Company Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |