CN109143375A - A kind of method and device in determining drilling well orientation - Google Patents
A kind of method and device in determining drilling well orientation Download PDFInfo
- Publication number
- CN109143375A CN109143375A CN201810946533.7A CN201810946533A CN109143375A CN 109143375 A CN109143375 A CN 109143375A CN 201810946533 A CN201810946533 A CN 201810946533A CN 109143375 A CN109143375 A CN 109143375A
- Authority
- CN
- China
- Prior art keywords
- target
- subterranean formation
- velocity field
- well
- formation zone
- 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
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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2200/00—Details of seismic or acoustic prospecting or detecting in general
- G01V2200/10—Miscellaneous details
- G01V2200/16—Measure-while-drilling or logging-while-drilling
Abstract
The embodiment of the present application discloses a kind of method and device in determining drilling well orientation.The described method includes: being based on vertical seismic profiling (VSP) velocity information, processing is corrected to the initial vertical misalignment velocity field of the first subterranean formation zone in initial depth migration velocity field, obtains the target vertical misalignment velocity field of the first subterranean formation zone；The well individual-layer data that place is set according to standard well location is corrected processing to the initial lateral migration velocity field of the first subterranean formation zone, obtains the target lateral migration velocity field of the first subterranean formation zone；According to the target longitudinal wave migration velocity field and target lateral migration velocity field of the first subterranean formation zone, and second subterranean formation zone initial vertical misalignment velocity field and initial lateral migration velocity field, pre-stack depth migration processing is carried out to primary earthquake data, obtains the target depth migration imaging result in purpose work area；Based on target depth migration imaging as a result, determining drilling well orientation of the target well on target zone.The accuracy in identified drilling well orientation can be improved.
Description
Technical field
This application involves Seismic Data Processing Technique field, in particular to a kind of method and device in determining drilling well orientation.
Background technique
Reservoir space in carbonates is mainly made of secondary hole, hole, seam, referred to as fracture hole type reservoir.This kind high-quality
Reservoir is usually expressed as " beading strong reflection " on seismic profile, is the main probing mesh of current carbonate exploration exploitation
Mark.In order to be accurately positioned the position of fracture hole type reservoir, and for probing drilling well orientation is provided, needs to fold seismic data
Preceding depth migration processing.Studies have shown that carbonate rock fractured cave body is very sensitive to migration velocity, the speed of general 30m~50m is missed
Difference will influence the imaging precision of a beading.In the exploration and development of early stage, target mainly based on large-scale cave, miss by small speed
Difference can be tolerated, however in the exploration of the new period, carbonate reservoir probing feature shows themselves in that 1, ultra-deep, often beyond
7000 meters；2, target spot is small, and area is in 20 meters or so of circumference, in this way, wanting when determining drilling well orientation to the precision of seismic data
Asking becomes higher, is drilled so as to cause using traditional depth migration result, and risk can become higher.Therefore, at present urgently
A kind of method for needing effective determining drilling well orientation, to improve the accuracy in identified drilling well orientation, to reduce drilling risk.
Summary of the invention
The purpose of the embodiment of the present application is to provide a kind of method and device in determining drilling well orientation, to improve identified bore
The accuracy in well orientation, to reduce drilling risk.
In order to solve the above technical problems, the method and device that the embodiment of the present application provides a kind of determining drilling well orientation is in this way
It realizes:
A kind of method in determining drilling well orientation, provides the primary earthquake data and initial depth migration velocity in purpose work area
Target well location sets the vertical seismic profiling (VSP) velocity information at place in field and the purpose work area；Wherein, the target well is drilled into institute
State the target zone coboundary in purpose work area；The described method includes:
Based on the vertical seismic profiling (VSP) velocity information, to the first subterranean formation zone in the initial depth migration velocity field
Initial vertical misalignment velocity field is corrected processing, obtains the target vertical misalignment velocity field of first subterranean formation zone；Wherein,
First subterranean formation zone indicates to be located at the subterranean formation zone above the target zone in the purpose work area；
For the standard well for being located at finishing drilling around the target well in the purpose work area, obtains the standard well location and set
The well individual-layer data at place, and according to the well individual-layer data, to the first subterranean formation zone in the initial depth migration velocity field
Initial transverse direction migration velocity field is corrected processing, obtains the target lateral migration velocity field of first subterranean formation zone；
According to the target longitudinal wave migration velocity field of first subterranean formation zone and target lateral migration velocity field and described
The initial vertical misalignment velocity field of second subterranean formation zone and initial lateral migration velocity field in initial depth migration velocity field, to institute
It states primary earthquake data and carries out pre-stack depth migration processing, obtain the target depth migration imaging result in the purpose work area；Its
In, second subterranean formation zone indicates target zone described in the purpose work area and the stratigraphic region below the target zone
Domain；
Based on the target depth migration imaging as a result, determining drilling well orientation of the target well on the target zone.
In preferred embodiment, to the initial vertical misalignment velocity field of the first subterranean formation zone in the initial depth migration velocity field
It is corrected processing, comprising:
By the way of Event tracking, picked up from the primary earthquake data first subterranean formation zone multiplely
The seismic profile layer position of layer；
From the initial vertical misalignment velocity field of first subterranean formation zone, mesh described in seismic profile layer position is obtained
It marks well and bores the initial vertical misalignment speed met at position, and the seismic profile described in the vertical seismic profiling (VSP) velocity information
Target well described in layer position bores the vertical seismic profiling (VSP) speed met at position, and based on target described in seismic profile layer position
Well bores the initial vertical misalignment speed and vertical seismic profiling (VSP) speed met at position, determines that seismic profile layer position is corresponding vertical
To correction coefficient；
By the corresponding longitudinal correction coefficient in seismic profile layer position respectively with each sampling on seismic profile layer position
The product of initial vertical misalignment speed at point position, as the target at each sampling point position on seismic profile layer position
Vertical misalignment speed；
According to the target vertical misalignment speed at the sampling point position on determining each seismic profile layer position, determine
The target vertical misalignment velocity field of first subterranean formation zone.
In preferred embodiment, the corresponding longitudinal correction coefficient in seismic profile layer position is determined, comprising:
By target well described in seismic profile layer position bore meet position at initial vertical misalignment speed average value with
The ratio of the average value of vertical seismic profiling (VSP) speed, as the corresponding longitudinal correction coefficient in seismic profile layer position.
In preferred embodiment, to the initial lateral migration velocity field of the first subterranean formation zone in the initial depth migration velocity field
It is corrected processing, comprising:
For in multiple stratum of the first subterranean formation zone specify stratum, obtained from the well individual-layer data it is described specifiedly
Standard well described in layer bores the well layer position met at position, and the earthquake from multiple stratum of first subterranean formation zone of pickup
Standard well described in the specified stratum is obtained in section layer position bores the seismic profile layer position met at position；Wherein, described specified
The absolute value of difference in the log data on stratum between the log data value of two neighboring logging point and the specified stratum
The ratio of the average value of the log data value of each logging point is less than or equal to prediction proportion threshold value；
According to standard well described in the specified stratum bore meet position at well layer position and seismic profile layer position, determine described in
The corresponding grid scale factor of each sampled point in specified stratum；
According to the corresponding grid scale factor of sampled point each in the specified stratum, determine in first subterranean formation zone
The corresponding grid scale factor of each sampled point；
Based on the corresponding grid scale factor of sampled point each in first subterranean formation zone, to first subterranean formation zone
Initial lateral migration velocity field be corrected processing, obtain the target lateral migration velocity field of first subterranean formation zone.
In preferred embodiment, the corresponding grid scale factor of each sampled point in the specified stratum is determined, comprising:
The standard well described in the specified stratum bores the depth met between well layer position and seismic profile layer position at position
When the absolute value of difference is less than or equal to predetermined depth threshold value, determines that standard well described in the specified stratum bores and meet at position
The grid scale factor is 1；
The standard well described in the specified stratum bores the depth met between well layer position and seismic profile layer position at position
When the absolute value of difference is greater than predetermined depth threshold value, by standard described in the absolute value of the depth difference and the specified stratum
Well bore meet position at well layer position to earth's surface distance ratio, as standard well described in the specified stratum bore meet position at
The grid scale factor；
The grid scale factor met at position is bored based on standard well described in the specified stratum, it is true by the way of interpolation
The corresponding grid scale factor of each sampled point in the fixed specified stratum.
In preferred embodiment, processing is corrected to the initial lateral migration velocity field of first subterranean formation zone, comprising:
When the corresponding grid scale factor of the specified sampled point in first subterranean formation zone is 1, described specify is adopted
The corresponding grid scale factor of sampling point is as the corresponding cross-level coefficient of the specified sampled point；
It, will be described specified when the corresponding grid scale factor of the specified sampled point in first subterranean formation zone is not 1
The corresponding grid scale factor of sampled point adds 1 to be used as the corresponding cross-level coefficient of the specified sampled point；
The corresponding cross-level coefficient of the specified sampled point and the initial transverse direction at the specified sampling point position is inclined
The product for moving speed, as the target lateral migration velocity at the specified sampling point position；
According to the target lateral migration velocity at the sampling point position in determining first subterranean formation zone, determine described in
The target lateral migration velocity field of first subterranean formation zone.
In preferred embodiment, before carrying out pre-stack depth migration processing to the primary earthquake data, the method is also wrapped
It includes:
It is less than the first specified region of specified degree of covering for degree of covering in the purpose work area, is based on described first
The degree of covering in specified region, repairs the primary earthquake subdata in the first specified region described in the primary earthquake data
Positive processing, obtains the seismic target earthquakes subdata in the described first specified region；
According to the second specified region in the seismic target earthquakes subdata in the described first specified region and the primary earthquake data
Primary earthquake subdata, determine the seismic target earthquakes data in the purpose work area；Wherein, described in the described second specified region expression
Degree of covering is greater than or equal to the region of specified degree of covering in purpose work area；
Correspondingly, pre-stack depth migration processing is carried out to the primary earthquake data, comprising: to the seismic target earthquakes data
Carry out pre-stack depth migration processing.
In preferred embodiment, the primary earthquake subdata in the first specified region described in the primary earthquake data is repaired
Positive processing, comprising:
The initial amplitude value of sampled point will be specified divided by the finger in the primary earthquake subdata in the described first specified region
The square root for determining the degree of covering at sampling point position obtains the target amplitude value of the specified sampled point；
According to the target amplitude value of the sampled point in the described first determining specified region, the described first specified region is determined
Seismic target earthquakes subdata.
A kind of device in determining drilling well orientation, the primary earthquake data and initial depth in described device offer purpose work area are inclined
Move the vertical seismic profiling (VSP) velocity information that target well location in velocity field and the purpose work area sets place；Wherein, the target well
It is drilled into the target zone coboundary in the purpose work area；Described device includes: target vertical misalignment velocity field determining module, target cross
To migration velocity field determining module, migration imaging result determining module and drilling well orientation determining module；Wherein,
The target vertical misalignment velocity field determining module, for being based on the vertical seismic profiling (VSP) velocity information, to institute
The initial vertical misalignment velocity field for stating the first subterranean formation zone in initial depth migration velocity field is corrected processing, obtains described
The target vertical misalignment velocity field of one subterranean formation zone；Wherein, first subterranean formation zone indicates to be located at institute in the purpose work area
State the subterranean formation zone above target zone；
The target lateral migration velocity field determining module, for all for the target well is located in the purpose work area
The standard well for enclosing finishing drilling obtains the well individual-layer data that the standard well location sets place, and according to the well individual-layer data, to described
The initial lateral migration velocity field of the first subterranean formation zone is corrected processing in initial depth migration velocity field, obtains described first
The target lateral migration velocity field of subterranean formation zone；
The migration imaging result determining module, for the target longitudinal wave migration velocity field according to first subterranean formation zone
With the initial vertical misalignment speed of the second subterranean formation zone in target lateral migration velocity field and the initial depth migration velocity field
Field and initial lateral migration velocity field are spent, pre-stack depth migration processing is carried out to the primary earthquake data, obtains the purpose
The target depth migration imaging result in work area；Wherein, second subterranean formation zone indicates target zone described in the purpose work area
With the subterranean formation zone being located at below the target zone；
Drilling well orientation determining module is used for based on the target depth migration imaging as a result, determining the target well
Drilling well orientation on the target zone.
A kind of device in determining drilling well orientation, including memory, processor and the calculating being stored on the memory
Machine program, the memory storage have the primary earthquake data and initial depth migration velocity field and the mesh in purpose work area
Work area in target well location set the vertical seismic profiling (VSP) velocity information at place；Wherein, the target well is drilled into the purpose work area
Target zone coboundary；The computer program executes following steps when being run by the processor:
Based on the vertical seismic profiling (VSP) velocity information, to the first subterranean formation zone in the initial depth migration velocity field
Initial vertical misalignment velocity field is corrected processing, obtains the target vertical misalignment velocity field of first subterranean formation zone；Wherein,
First subterranean formation zone indicates to be located at the subterranean formation zone above the target zone in the purpose work area；
For the standard well for being located at finishing drilling around the target well in the purpose work area, obtains the standard well location and set
The well individual-layer data at place, and according to the well individual-layer data, to the first subterranean formation zone in the initial depth migration velocity field
Initial transverse direction migration velocity field is corrected processing, obtains the target lateral migration velocity field of first subterranean formation zone；
According to the target longitudinal wave migration velocity field of first subterranean formation zone and target lateral migration velocity field and described
The initial vertical misalignment velocity field of second subterranean formation zone and initial lateral migration velocity field in initial depth migration velocity field, to institute
It states primary earthquake data and carries out pre-stack depth migration processing, obtain the target depth migration imaging result in the purpose work area；Its
In, second subterranean formation zone indicates target zone described in the purpose work area and the stratigraphic region below the target zone
Domain；
Based on the target depth migration imaging as a result, determining drilling well orientation of the target well on the target zone.
As can be seen from the technical scheme provided by the above embodiments of the present application, determining drilling well orientation provided by the embodiments of the present application
Method and device can be based on the vertical seismic profiling (VSP) velocity information, in the initial depth migration velocity field the first
The initial vertical misalignment velocity field of one subterranean formation zone is corrected processing, obtains the target vertical misalignment of first subterranean formation zone
Velocity field；Wherein, first subterranean formation zone indicates to be located at the subterranean formation zone above the target zone in the purpose work area.It connects
, for the standard well for being located at finishing drilling around the target well in the purpose work area, the available standard well location is set
The well individual-layer data at place, and can be according to the well individual-layer data, to the first stratigraphic region in the initial depth migration velocity field
The initial lateral migration velocity field in domain is corrected processing, obtains the target lateral migration velocity field of first subterranean formation zone.
It is then possible to according to the target longitudinal wave migration velocity field and target lateral migration velocity field of first subterranean formation zone, Yi Jisuo
The initial vertical misalignment velocity field of the second subterranean formation zone and initial lateral migration velocity field in initial depth migration velocity field are stated, it is right
The primary earthquake data carry out pre-stack depth migration processing, obtain the target depth migration imaging result in the purpose work area；
Wherein, second subterranean formation zone indicates target zone described in the purpose work area and the stratigraphic region below the target zone
Domain.Finally, can be based on the target depth migration imaging as a result, determining drilling well of the target well on the target zone
Orientation.In this way, by the determining higher target depth migration imaging of precision as a result, to instruct the target well in the purpose
The accuracy in identified drilling well orientation can be improved, to reduce drilling risk in drilling well orientation on layer.
Detailed description of the invention
In order to illustrate the technical solutions in the embodiments of the present application or in the prior art more clearly, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
The some embodiments recorded in application, for those of ordinary skill in the art, in the premise of not making the creative labor property
Under, it is also possible to obtain other drawings based on these drawings.
Fig. 1 is a kind of flow chart of the embodiment of the method in determining drilling well orientation of the application；
Fig. 2 is each seismic profile layer position longitudinal direction correction coefficient in the first subterranean formation zone determined in the embodiment of the present application
Schematic diagram；
Fig. 3 is showing of showing on seismic profile of the well individual-layer data of the standard well in the embodiment of the present application around target well
It is intended to；
Fig. 4 is a kind of composed structure schematic diagram of the embodiment for the device that the application determines drilling well orientation；
Fig. 5 is the composed structure schematic diagram of another embodiment for the device that the application determines drilling well orientation.
Specific embodiment
The embodiment of the present application provides a kind of method and device in determining drilling well orientation.
In order to make those skilled in the art better understand the technical solutions in the application, below in conjunction with the application reality
The attached drawing in example is applied, the technical scheme in the embodiment of the application is clearly and completely described, it is clear that described implementation
Example is merely a part but not all of the embodiments of the present application.Based on the embodiment in the application, this field is common
The application protection all should belong in technical staff's every other embodiment obtained without creative efforts
Range.
This application provides a kind of methods in determining drilling well orientation.The method in the determining drilling well orientation can be provided with mesh
Work area primary earthquake data and initial depth migration velocity field and the purpose work area in target well location set the vertical of place
Seismic profile velocity information；Wherein, the target well is drilled into the target zone coboundary in the purpose work area.
In the present embodiment, the purpose work area can be work area still leaved for development or being developed, such as in tower
Tub ground TZ area etc..Target well location in the purpose work area, which is set, can refer to by opening up in the purpose work area
Drilling well, be able to detect that the position of vertical seismic profiling (VSP) (Vertical Seismic Profiling, VSP) velocity information.
Wherein, the target well is drilled into the coboundary of the target zone in the purpose work area from earth's surface, does not bore also and meets the target zone.
In the present embodiment, the target zone can refer to the reservoir for not yet determining drilling well orientation.The reservoir for example may be used
To be fracture hole type reservoir.The depth of the target zone is usually relatively deep, for example, the depth of the target zone can be more than 7000
Rice.
In the present embodiment, can by way of seismic prospecting and data acquisition mode, obtain the purpose work
The primary earthquake data and initial depth migration velocity field in area.It may include initial longitudinal direction in the initial depth migration velocity field
Migration velocity field and initial lateral migration velocity field.Wherein, horizontal and vertical is respectively horizontal direction and vertical with horizontal direction
Direction.
Fig. 1 is a kind of flow chart of the embodiment of the method in determining drilling well orientation of the application.As shown in Figure 1, the determining brill
The method in well orientation, includes the following steps.
Step S101: it is based on the vertical seismic profiling (VSP) velocity information, in the initial depth migration velocity field first
The initial vertical misalignment velocity field of subterranean formation zone is corrected processing, obtains the target vertical misalignment speed of first subterranean formation zone
Spend field；Wherein, first subterranean formation zone indicates to be located at the subterranean formation zone above the target zone in the purpose work area.
In the present embodiment, it is based on the vertical seismic profiling (VSP) velocity information, speed can be deviated to the initial depth
The initial vertical misalignment velocity field of the first subterranean formation zone is corrected processing in degree field, obtains the target of first subterranean formation zone
Vertical misalignment velocity field.Wherein, first subterranean formation zone indicates to be located at the ground above the target zone in the purpose work area
Layer region.Specifically, first ground can be picked up from the primary earthquake data by the way of Event tracking first
The seismic profile layer position on multiple stratum of layer region.It then, can be from the initial vertical misalignment speed of first subterranean formation zone
In, obtains target well described in seismic profile layer position and bore the initial vertical misalignment speed met at position, and from described
Target well described in seismic profile layer position described in vertical seismic profiling (VSP) velocity information bores the vertical seismic profiling (VSP) speed met at position
Degree, and the initial vertical misalignment speed and vertically met at position can be bored based on target well described in seismic profile layer position
Section speed is shaken, determines the corresponding longitudinal correction coefficient in seismic profile layer position.It is then possible to by seismic profile layer position
It is corresponding longitudinal direction correction coefficient respectively with the initial vertical misalignment speed at each sampling point position on seismic profile layer position
Product, as the target vertical misalignment speed at each sampling point position on seismic profile layer position.Finally, can basis
The target vertical misalignment speed at sampling point position on determining each seismic profile layer position, determines first stratum
The target vertical misalignment velocity field in region.
In the present embodiment, by the way of Event tracking, described first is picked up from the primary earthquake data
The seismic profile layer position on multiple stratum of subterranean formation zone, can specifically include, and can compare adjacent in the primary earthquake data
Two seismic channel datas are referring to corresponding wave character in timing window, and according to the similarity of wave character, pick up described first
The seismic profile layer position on multiple stratum of subterranean formation zone.Seismic profile layer position can be used for characterizing first subterranean formation zone
Depth of stratum of the middle Different Strata on the seismic profile of the primary earthquake data characterization.
In the present embodiment, it is inclined that the initial longitudinal direction met at position is bored based on target well described in seismic profile layer position
Speed and vertical seismic profiling (VSP) speed are moved, the corresponding longitudinal correction coefficient in seismic profile layer position is determined, can specifically include,
Can by target well described in seismic profile layer position bore meet position at initial vertical misalignment speed average value with it is vertical
The ratio of the average value of seismic profile speed, as the corresponding longitudinal correction coefficient in seismic profile layer position.Wherein, describedly
The corresponding longitudinal correction coefficient of each sampled point at shake section layer position can be identical.The sampled point can refer to described
Sampled point in primary earthquake data.For example, Fig. 2 is that each earthquake is cutd open in the first subterranean formation zone determined in the embodiment of the present application
The schematic diagram of surface layer position longitudinal direction correction coefficient.Abscissa and ordinate in Fig. 2 are respectively earthquake Taoist monastic name and correction coefficient.
Step S102: for the standard well for being located at finishing drilling around the target well in the purpose work area, described in acquisition
Standard well location sets the well individual-layer data at place, and according to the well individual-layer data, in the initial depth migration velocity field first
The initial lateral migration velocity field of subterranean formation zone is corrected processing, obtains the target lateral offset speed of first subterranean formation zone
Spend field.
It in the present embodiment, can for the standard well for being located at finishing drilling around the target well in the purpose work area
To obtain the well individual-layer data that the standard well location sets place.
It in the present embodiment, can be in the initial depth migration velocity field first according to the well individual-layer data
The initial lateral migration velocity field of subterranean formation zone is corrected processing, obtains the target lateral offset speed of first subterranean formation zone
Field is spent, can specifically include, firstly, can be layered from the well for stratum is specified in multiple stratum of the first subterranean formation zone
Standard well described in the specified stratum is obtained in data bores the well layer position met at position.For example, Fig. 3 is in the embodiment of the present application
The schematic diagram that the well individual-layer data of standard well around target well is shown on seismic profile.Abscissa and ordinate in Fig. 3
Respectively earthquake Taoist monastic name and sampling time, wherein the unit in sampling time is millisecond.As shown in figure 3, around the target well
Standard well for example can be the standard well of ZG7-5, ZG702, ZG7, TZ and ZG512, and in these standard wells, stratum TC is corresponding
Well layer position is respectively TC-TC3069, TC-TC3053, TC-TC3035, TC-TC2966 and TC-TC3025 and TO3s pairs of stratum
The well layer position answered be respectively TO3s-TO3s4261, TO3s-TO3s4168, TO3s-TO3s4149, TO3s-TO3s3776 and
TO3s-TO3s3715.Described in being obtained from the seismic profile layer position on multiple stratum of first subterranean formation zone of pickup
Standard well described in specified stratum bores the seismic profile layer position met at position.Wherein, phase in the log data on the specified stratum
The well logging number of each logging point of the absolute value and specified stratum of difference between the log data value of adjacent two logging points
It is less than or equal to prediction proportion threshold value according to the ratio of the average value of value.Wherein, the log data specifically can be resistivity survey
Well curve data, density log curve data etc..Correspondingly, the log data value can be resistivity value, density etc..
The value range of the prediction proportion threshold value can be 0~0.2, specifically, can be set according to practical logging situation, herein
It is not construed as limiting.Then, can the standard well according to the specified stratum bore meet position at well layer position and seismic profile layer
Position determines the corresponding grid scale factor of each sampled point in the specified stratum.It is then possible to according in the specified stratum
The corresponding grid scale factor of each sampled point, determine in first subterranean formation zone corresponding grid scale of each sampled point because
Son.Finally, can be based on the corresponding grid scale factor of sampled point each in first subterranean formation zone, to first stratum
The initial lateral migration velocity field in region is corrected processing, obtains the target lateral migration velocity of first subterranean formation zone
?.
In the present embodiment, the well layer position at the brill of the standard well according to specified stratum chance position and earthquake are cutd open
Surface layer position determines the corresponding grid scale factor of each sampled point in the specified stratum, can specifically include, when described specified
Standard well described in stratum bore meet position at well layer position and seismic profile layer position between depth difference absolute value be less than or
When equal to predetermined depth threshold value, it can determine that standard well described in the specified stratum bores the grid scale factor met at position and is
1.The standard well described in the specified stratum, which bores, meets the depth difference between well layer position and seismic profile layer position at position
It, can be by standard well described in the absolute value of the depth difference and the specified stratum when absolute value is greater than predetermined depth threshold value
Bore meet position at well layer position to earth's surface distance ratio, as standard well described in the specified stratum bore meet position at
The grid scale factor.The grid scale factor met at position can be bored based on standard well described in the specified stratum, using slotting
The mode of value determines the corresponding grid scale factor of each sampled point in the specified stratum.For example, can be based on described specified
Standard well described in stratum bores the grid scale factor met at position and determines the finger using the interpolation method of inverse distance-weighting
Determine in stratum except the standard well bores the corresponding grid scale factor of other sampled points in addition to meeting the sampled point at position.
In the present embodiment, according to the corresponding grid scale factor of sampled point each in the specified stratum, institute is determined
The corresponding grid scale factor of each sampled point in the first subterranean formation zone is stated, can specifically include, it can be by the specified stratum
In the corresponding grid scale factor of a sampled point, respectively as in other stratum of first subterranean formation zone with the sampled point
The corresponding grid scale factor of the identical each sampled point in horizontal position.Wherein, other described stratum can indicate described
Stratum in one subterranean formation zone in addition to the specified stratum.In this way, each in first subterranean formation zone adopt can be determined
The corresponding grid scale factor of sampling point.
In the present embodiment, right based on the corresponding grid scale factor of sampled point each in first subterranean formation zone
The initial lateral migration velocity field of first subterranean formation zone is corrected processing, and the target for obtaining first subterranean formation zone is horizontal
It to migration velocity field, can specifically include, when the corresponding grid scale factor of specified sampled point in first subterranean formation zone
It, can be using the corresponding grid scale factor of the specified sampled point as the corresponding cross-level of the specified sampled point when being 1
Coefficient.When the corresponding grid scale factor of the specified sampled point in first subterranean formation zone is not 1, described can will specify
The corresponding grid scale factor of sampled point adds 1 to be used as the corresponding cross-level coefficient of the specified sampled point.It then, can be by institute
The product of the initial lateral shift speed at the corresponding cross-level coefficient of specified sampled point and the specified sampling point position is stated,
As the target lateral migration velocity at the specified sampling point position.Finally, can be according to determining first stratigraphic region
The target lateral migration velocity at each sampling point position in domain determines the target lateral offset speed of first subterranean formation zone
Spend field.
Step S103: according to the target longitudinal wave migration velocity field and target lateral migration velocity of first subterranean formation zone
The initial vertical misalignment velocity field and initial lateral shift of second subterranean formation zone in field and the initial depth migration velocity field
Velocity field carries out pre-stack depth migration processing to the primary earthquake data, obtains the target depth offset in the purpose work area
Imaging results；Wherein, second subterranean formation zone indicates target zone described in the purpose work area and is located under the target zone
The subterranean formation zone of side.
In the present embodiment, it is deviated according to the target longitudinal wave migration velocity field of first subterranean formation zone and target lateral
The initial vertical misalignment velocity field and target lateral of second subterranean formation zone in velocity field and the initial depth migration velocity field
Migration velocity field can carry out pre-stack depth migration processing to the primary earthquake data, obtain the target in the purpose work area
Depth migration imaging result.Wherein, second subterranean formation zone indicates target zone described in the purpose work area and is located at described
Subterranean formation zone below target zone.Specifically, it is possible, firstly, to according to the target longitudinal wave migration velocity field of first subterranean formation zone
With the initial vertical misalignment speed of the second subterranean formation zone in target lateral migration velocity field and the initial depth migration velocity field
Field and initial lateral migration velocity field are spent, determines the target offset velocity field in the purpose work area.It is then possible to according to the mesh
Migration velocity field is marked, pre-stack depth migration processing, the mesh in the available purpose work area are carried out to the primary earthquake data
Mark depth migration imaging result.Wherein it is possible to identify the seam in the target zone from the target depth migration imaging result
The position of hole body.
Step S104: based on the target depth migration imaging as a result, determining the target well on the target zone
Drilling well orientation.
In the present embodiment, based on the target depth migration imaging as a result, determining the target well in the purpose
Drilling well orientation on layer, can specifically include, and can identify in the target zone from the target depth migration imaging result
Fracture hole body position.For example, the beading that can be showed in the target depth migration imaging result according to the fracture hole body
Shape strong reflection feature is identified.It, can be according to the target after the position for identifying the fracture hole body in the target zone
Well is drilled into the position of the fracture hole body in the position and the target zone of the target zone coboundary, determines the target well described
Drilling well orientation on target zone.
In one embodiment of the application, in order to avoid because the target well location is in the discontented covering in the purpose work area
Region, caused by fracture hole body position drift condition, wherein the discontented overlay area can indicate to cover in the purpose work area
Lid number is less than the region of specified degree of covering, can handle it carrying out pre-stack depth migration to the primary earthquake data
Before, the method can also include that be less than the specified degree of covering for degree of covering in the purpose work area first refers to
Region is determined, based on the degree of covering in the described first specified region, to the first specified region described in the primary earthquake data
Primary earthquake subdata is modified processing, obtains the seismic target earthquakes subdata in the described first specified region.It then, can basis
Primary earthquake in the second specified region in the seismic target earthquakes subdata and the primary earthquake data in the first specified region
Data determine the seismic target earthquakes data in the purpose work area；Wherein, the described second specified region indicates to cover in the purpose work area
Lid number is greater than or equal to the region of specified degree of covering.Correspondingly, pre-stack depth migration is carried out to the primary earthquake data
Processing, can specifically include, and carry out pre-stack depth migration processing to the seismic target earthquakes data.In the present embodiment, described
Specified degree of covering specifically can be the degree of covering that observation system designs during seismic prospecting.
In the present embodiment, the degree of covering based on the described first specified region, to institute in the primary earthquake data
The primary earthquake subdata for stating the first specified region is modified processing, obtains the seismic target earthquakes subnumber in the described first specified region
According to, can specifically include, can will in the primary earthquake subdata in the described first specified region specify sampled point initial amplitude
Value obtains the target amplitude value of the specified sampled point divided by the square root of the degree of covering at the specified sampling point position.
Wherein, the specified sampled point can indicate any one sampled point in the primary earthquake subdata.It can be according to determining
The target amplitude value of sampled point in the first specified region, determines the seismic target earthquakes subdata in the described first specified region.
The embodiment of the method in the determining drilling well orientation can be based on the vertical seismic profiling (VSP) velocity information first, right
The initial vertical misalignment velocity field of the first subterranean formation zone is corrected processing in the initial depth migration velocity field, obtains described
The target vertical misalignment velocity field of first subterranean formation zone；Wherein, first subterranean formation zone indicates to be located in the purpose work area
Subterranean formation zone above the target zone.Then, for the mark for being located at finishing drilling around the target well in the purpose work area
Quasi- well, the available standard well location set the well individual-layer data at place, and can be according to the well individual-layer data, to described initial
The initial lateral migration velocity field of the first subterranean formation zone is corrected processing in depth migration velocity field, obtains first stratum
The target lateral migration velocity field in region.It is then possible to according to the target longitudinal wave migration velocity field of first subterranean formation zone and
The initial vertical misalignment speed of second subterranean formation zone in target lateral migration velocity field and the initial depth migration velocity field
The primary earthquake data are carried out pre-stack depth migration processing, obtain the purpose work by field and initial lateral migration velocity field
The target depth migration imaging result in area；Wherein, second subterranean formation zone indicate target zone described in the purpose work area and
Subterranean formation zone below the target zone.Finally, can be based on the target depth migration imaging as a result, described in determining
Drilling well orientation of the target well on the target zone.In this way, by the determining higher target depth migration imaging of precision as a result,
It instructs drilling well orientation of the target well on the target zone, the accuracy in identified drilling well orientation can be improved, from
And reduce drilling risk.
Fig. 4 is a kind of composed structure schematic diagram of the embodiment for the device that the application determines drilling well orientation.The determining brill
The device in well orientation can provide the primary earthquake data and initial depth migration velocity field and the purpose work in purpose work area
Target well location sets the vertical seismic profiling (VSP) velocity information at place in area；Wherein, the target well is drilled into the purpose in the purpose work area
Layer coboundary.As shown in figure 4, determine drilling well orientation device may include: target vertical misalignment velocity field determining module 100,
Target lateral migration velocity field determining module 200, migration imaging result determining module 300 and drilling well orientation determining module 400.
The target vertical misalignment velocity field determining module 100 can be used for believing based on the vertical seismic profiling (VSP) speed
Breath, is corrected processing to the initial vertical misalignment velocity field of the first subterranean formation zone in the initial depth migration velocity field, obtains
To the target vertical misalignment velocity field of first subterranean formation zone；Wherein, first subterranean formation zone indicates the purpose work area
In be located at the target zone above subterranean formation zone.
The target lateral migration velocity field determining module 200 can be used for described for being located in the purpose work area
The standard well of finishing drilling around target well obtains the well individual-layer data that the standard well location sets place, and according to the well hierarchy number
According to being corrected processing to the initial lateral migration velocity field of the first subterranean formation zone in the initial depth migration velocity field, obtain
To the target lateral migration velocity field of first subterranean formation zone.
The migration imaging result determining module 300 can be used for inclined according to the target longitudinal wave of first subterranean formation zone
The initial of the second subterranean formation zone in velocity field and target lateral migration velocity field and the initial depth migration velocity field is moved to indulge
To migration velocity field and initial lateral migration velocity field, pre-stack depth migration processing is carried out to the primary earthquake data, is obtained
The target depth migration imaging result in the purpose work area；Wherein, second subterranean formation zone indicates institute in the purpose work area
State target zone and the subterranean formation zone below the target zone.
Drilling well orientation determining module 400 can be used for based on the target depth migration imaging as a result, described in determining
Drilling well orientation of the target well on the target zone.
Fig. 5 is the composed structure schematic diagram of another embodiment for the device that the application determines drilling well orientation.Such as Fig. 5 institute
Show, the device in the determining drilling well orientation may include memory, processor and the calculating being stored on the memory
Machine program, the memory storage have the primary earthquake data and initial depth migration velocity field and the mesh in purpose work area
Work area in target well location set the vertical seismic profiling (VSP) velocity information at place；Wherein, the target well is drilled into the purpose work area
Target zone coboundary；The computer program executes following steps when being run by the processor:
Step S101: it is based on the vertical seismic profiling (VSP) velocity information, in the initial depth migration velocity field first
The initial vertical misalignment velocity field of subterranean formation zone is corrected processing, obtains the target vertical misalignment speed of first subterranean formation zone
Spend field；Wherein, first subterranean formation zone indicates to be located at the subterranean formation zone above the target zone in the purpose work area；
Step S102: for the standard well for being located at finishing drilling around the target well in the purpose work area, described in acquisition
Standard well location sets the well individual-layer data at place, and according to the well individual-layer data, in the initial depth migration velocity field first
The initial lateral migration velocity field of subterranean formation zone is corrected processing, obtains the target lateral offset speed of first subterranean formation zone
Spend field；
Step S103: according to the target longitudinal wave migration velocity field and target lateral migration velocity of first subterranean formation zone
The initial vertical misalignment velocity field and initial lateral shift of second subterranean formation zone in field and the initial depth migration velocity field
Velocity field carries out pre-stack depth migration processing to the primary earthquake data, obtains the target depth offset in the purpose work area
Imaging results；Wherein, second subterranean formation zone indicates target zone described in the purpose work area and is located under the target zone
The subterranean formation zone of side；
Step S104: based on the target depth migration imaging as a result, determining the target well on the target zone
Drilling well orientation.
The Installation practice in the determining drilling well orientation is corresponding with the embodiment of the method in the determining drilling well orientation, can be with
Realize the technical solution for determining the embodiment of the method in drilling well orientation, and the technical effect of adquisitiones embodiment.
In the 1990s, the improvement of a technology can be distinguished clearly be on hardware improvement (for example,
Improvement to circuit structures such as diode, transistor, switches) or software on improvement (improvement for method flow).So
And with the development of technology, the improvement of current many method flows can be considered as directly improving for hardware circuit.
Designer nearly all obtains corresponding hardware circuit by the way that improved method flow to be programmed into hardware circuit.Cause
This, it cannot be said that the improvement of a method flow cannot be realized with hardware entities module.For example, programmable logic device
(Programmable Logic Device, PLD) (such as field programmable gate array (Field Programmable Gate
Array, FPGA)) it is exactly such a integrated circuit, logic function determines device programming by user.By designer
Voluntarily programming comes a digital display circuit " integrated " on a piece of PLD, designs and makes without asking chip maker
Dedicated IC chip.Moreover, nowadays, substitution manually makes IC chip, this programming is also used instead mostly " is patrolled
Volume compiler (logic compiler) " software realizes that software compiler used is similar when it writes with program development,
And the source code before compiling also write by handy specific programming language, this is referred to as hardware description language
(Hardware Description Language, HDL), and HDL is also not only a kind of, but there are many kind, such as ABEL
(Advanced Boolean Expression Language)、AHDL(Altera Hardware Description
Language)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL
(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(Ruby
Hardware Description Language) etc., VHDL (Very-High-Speed is most generally used at present
Integrated Circuit Hardware Description Language) and Verilog.Those skilled in the art also answer
This understands, it is only necessary to method flow slightly programming in logic and is programmed into integrated circuit with above-mentioned several hardware description languages,
The hardware circuit for realizing the logical method process can be readily available.
It is also known in the art that other than realizing controller in a manner of pure computer readable program code, it is complete
Entirely can by by method and step carry out programming in logic come so that controller with logic gate, switch, specific integrated circuit, programmable
Logic controller realizes identical function with the form for being embedded in microcontroller etc..Therefore this controller is considered one kind
Hardware component, and the structure that the device for realizing various functions for including in it can also be considered as in hardware component.Or
Even, can will be considered as realizing the device of various functions either the software module of implementation method can be Hardware Subdivision again
Structure in part.
Device that above-described embodiment illustrates, module can specifically realize by computer chip or entity, or by having certain
The product of function is planted to realize.
For convenience of description, it is divided into various modules when description apparatus above with function to describe respectively.Certainly, implementing this
The function of each module can be realized in the same or multiple software and or hardware when application.
As seen through the above description of the embodiments, those skilled in the art can be understood that the application can
It realizes by means of software and necessary general hardware platform.Based on this understanding, the technical solution essence of the application
On in other words the part that contributes to existing technology can be embodied in the form of software products, in a typical configuration
In, calculating equipment includes one or more processors (CPU), input/output interface, network interface and memory.The computer is soft
Part product may include that some instructions are used so that a computer equipment (can be personal computer, server or network
Equipment etc.) execute method described in certain parts of each embodiment of the application or embodiment.The computer software product can
To be stored in memory, memory may include the non-volatile memory in computer-readable medium, random access memory
(RAM) and/or the forms such as Nonvolatile memory, such as read-only memory (ROM) or flash memory (flash RAM).Memory is computer
The example of readable medium.Computer-readable medium includes that permanent and non-permanent, removable and non-removable media can be by
Any method or technique come realize information store.Information can be computer readable instructions, data structure, the module of program or its
His data.The example of the storage medium of computer includes, but are not limited to phase change memory (PRAM), static random access memory
(SRAM), dynamic random access memory (DRAM), other kinds of random access memory (RAM), read-only memory
(ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory techniques, CD-ROM are read-only
Memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassettes, tape magnetic disk storage or
Other magnetic storage devices or any other non-transmission medium, can be used for storage can be accessed by a computing device information.According to
Herein defines, and computer-readable medium does not include of short duration computer readable media (transitory media), such as modulation
Data-signal and carrier wave.
All the embodiments in this specification are described in a progressive manner, same and similar portion between each embodiment
Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.Especially for device reality
For applying example, since it is substantially similar to the method embodiment, so being described relatively simple, related place is referring to embodiment of the method
Part explanation.
The application can be used in numerous general or special purpose computing system environments or configuration.Such as: personal computer, clothes
Business device computer, handheld device or portable device, laptop device, multicomputer system, microprocessor-based system, set
Top box, programmable consumer-elcetronics devices, network PC, minicomputer, mainframe computer including any of the above system or equipment
Distributed computing environment etc..
The application can describe in the general context of computer-executable instructions executed by a computer, such as program
Module.Generally, program module includes routines performing specific tasks or implementing specific abstract data types, programs, objects, group
Part, data structure etc..The application can also be practiced in a distributed computing environment, in these distributed computing environments, by
Task is executed by the connected remote processing devices of communication network.In a distributed computing environment, program module can be with
In the local and remote computer storage media including storage equipment.
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application there are many deformation and
Variation is without departing from spirit herein, it is desirable to which the attached claims include these deformations and change without departing from the application's
Spirit.
Claims (10)
1. a kind of method in determining drilling well orientation, which is characterized in that provide the primary earthquake data and initial depth in purpose work area
Target well location sets the vertical seismic profiling (VSP) velocity information at place in degree migration velocity field and the purpose work area；Wherein, the mesh
Mark well is drilled into the target zone coboundary in the purpose work area；The described method includes:
Based on the vertical seismic profiling (VSP) velocity information, in the initial depth migration velocity field the first subterranean formation zone it is initial
Vertical misalignment velocity field is corrected processing, obtains the target vertical misalignment velocity field of first subterranean formation zone；Wherein, described
First subterranean formation zone indicates to be located at the subterranean formation zone above the target zone in the purpose work area；
For the standard well for being located at finishing drilling around the target well in the purpose work area, obtains the standard well location and set place
Well individual-layer data, and according to the well individual-layer data, in the initial depth migration velocity field the first subterranean formation zone it is initial
Lateral shift speed field is corrected processing, obtains the target lateral migration velocity field of first subterranean formation zone；
According to the target longitudinal wave migration velocity field of first subterranean formation zone and target lateral migration velocity field and described initial
The initial vertical misalignment velocity field of second subterranean formation zone and initial lateral migration velocity field in depth migration velocity field, to described first
Beginning seismic data carries out pre-stack depth migration processing, obtains the target depth migration imaging result in the purpose work area；Wherein, institute
Stating the second subterranean formation zone indicates target zone described in the purpose work area and the subterranean formation zone below the target zone；
Based on the target depth migration imaging as a result, determining drilling well orientation of the target well on the target zone.
2. the method according to claim 1, wherein to the first stratigraphic region in the initial depth migration velocity field
The initial vertical misalignment velocity field in domain is corrected processing, comprising:
By the way of Event tracking, multiple stratum of first subterranean formation zone are picked up from the primary earthquake data
Seismic profile layer position；
From the initial vertical misalignment velocity field of first subterranean formation zone, target well described in seismic profile layer position is obtained
Bore the initial vertical misalignment speed met at position, and the seismic profile layer position described in the vertical seismic profiling (VSP) velocity information
Described in target well bore the vertical seismic profiling (VSP) speed met at position, and bored based on target well described in seismic profile layer position
The initial vertical misalignment speed and vertical seismic profiling (VSP) speed at position are met, determines the corresponding longitudinal school in seismic profile layer position
Positive coefficient；
By the corresponding longitudinal correction coefficient in seismic profile layer position respectively with each sample on seismic profile layer position
The product for setting the initial vertical misalignment speed at place, it is longitudinal as the target at each sampling point position on seismic profile layer position
Migration velocity；
According to the target vertical misalignment speed at the sampling point position on determining each seismic profile layer position, determine described in
The target vertical misalignment velocity field of first subterranean formation zone.
3. according to the method described in claim 2, it is characterized in that, determining the corresponding longitudinal correction system in seismic profile layer position
Number, comprising:
By target well described in seismic profile layer position bore meet position at initial vertical misalignment speed average value with it is vertical
The ratio of the average value of seismic profile speed, as the corresponding longitudinal correction coefficient in seismic profile layer position.
4. according to the method described in claim 2, it is characterized in that, to the first stratigraphic region in the initial depth migration velocity field
The initial lateral migration velocity field in domain is corrected processing, comprising:
For stratum is specified in multiple stratum of the first subterranean formation zone, obtained in the specified stratum from the well individual-layer data
The standard well bores the well layer position met at position, and the seismic profile from multiple stratum of first subterranean formation zone of pickup
Standard well described in the specified stratum, which is obtained, in layer position bores the seismic profile layer position met at position；Wherein, the specified stratum
Log data in two neighboring logging point log data value between difference absolute value and the specified stratum it is each
The ratio of the average value of the log data value of logging point is less than or equal to prediction proportion threshold value；
The well layer position and seismic profile layer position met at position are bored according to standard well described in the specified stratum, is determined described specified
The corresponding grid scale factor of each sampled point in stratum；
According to the corresponding grid scale factor of sampled point each in the specified stratum, determine each in first subterranean formation zone
The corresponding grid scale factor of sampled point；
Based on the corresponding grid scale factor of sampled point each in first subterranean formation zone, to the first of first subterranean formation zone
Beginning lateral shift speed field is corrected processing, obtains the target lateral migration velocity field of first subterranean formation zone.
5. according to the method described in claim 4, it is characterized in that, determining the corresponding cross of each sampled point in the specified stratum
To scale factor, comprising:
The standard well described in the specified stratum bores the depth difference met between well layer position and seismic profile layer position at position
Absolute value be less than or equal to predetermined depth threshold value when, determine standard well described in the specified stratum bore meet position at transverse direction
Scale factor is 1；
The standard well described in the specified stratum bores the depth difference met between well layer position and seismic profile layer position at position
Absolute value when being greater than predetermined depth threshold value, standard well described in the absolute value of the depth difference and the specified stratum is bored
Meet position at well layer position to earth's surface distance ratio, as standard well described in the specified stratum bore meet position at cross
To scale factor；
The grid scale factor met at position is bored based on standard well described in the specified stratum, and institute is determined by the way of interpolation
State the corresponding grid scale factor of each sampled point in specified stratum.
6. according to the method described in claim 4, it is characterized in that, initial lateral shift speed to first subterranean formation zone
Field is corrected processing, comprising:
When the corresponding grid scale factor of the specified sampled point in first subterranean formation zone is 1, by the specified sampled point
The corresponding grid scale factor is as the corresponding cross-level coefficient of the specified sampled point；
When the corresponding grid scale factor of the specified sampled point in first subterranean formation zone is not 1, by the specified sampling
The corresponding grid scale factor of point adds 1 to be used as the corresponding cross-level coefficient of the specified sampled point；
By the initial lateral shift speed at the specified corresponding cross-level coefficient of sampled point and the specified sampling point position
The product of degree, as the target lateral migration velocity at the specified sampling point position；
According to the target lateral migration velocity at the sampling point position in determining first subterranean formation zone, described first is determined
The target lateral migration velocity field of subterranean formation zone.
7. the method according to claim 1, wherein carrying out pre-stack depth migration to the primary earthquake data
Before processing, the method also includes:
It is less than the first specified region of specified degree of covering for degree of covering in the purpose work area, it is specified based on described first
The degree of covering in region is modified place to the primary earthquake subdata in the first specified region described in the primary earthquake data
Reason, obtains the seismic target earthquakes subdata in the described first specified region；
According in the seismic target earthquakes subdata in the described first specified region and the primary earthquake data the second specified region just
Beginning earthquake subdata determines the seismic target earthquakes data in the purpose work area；Wherein, the described second specified region indicates the purpose
Degree of covering is greater than or equal to the region of specified degree of covering in work area；
Correspondingly, pre-stack depth migration processing is carried out to the primary earthquake data, comprising: carry out to the seismic target earthquakes data
Pre-stack depth migration processing.
8. the method according to the description of claim 7 is characterized in that the first specified region described in the primary earthquake data
Primary earthquake subdata be modified processing, comprising:
The initial amplitude value of sampled point will be specified to adopt in the primary earthquake subdata in the described first specified region divided by described specify
The square root of degree of covering at sampling point position obtains the target amplitude value of the specified sampled point；
According to the target amplitude value of the sampled point in the described first determining specified region, the mesh in the described first specified region is determined
Mark earthquake subdata.
9. a kind of device in determining drilling well orientation, which is characterized in that described device provide purpose work area primary earthquake data and
Target well location sets the vertical seismic profiling (VSP) velocity information at place in initial depth migration velocity field and the purpose work area；Wherein,
The target well is drilled into the target zone coboundary in the purpose work area；Described device includes: that target vertical misalignment velocity field determines
Module, target lateral migration velocity field determining module, migration imaging result determining module and drilling well orientation determining module；Wherein,
The target vertical misalignment velocity field determining module, for being based on the vertical seismic profiling (VSP) velocity information, to described first
The initial vertical misalignment velocity field of the first subterranean formation zone is corrected processing in beginning depth migration velocity field, obtains first ground
The target vertical misalignment velocity field of layer region；Wherein, first subterranean formation zone indicates to be located at the mesh in the purpose work area
Layer above subterranean formation zone；
The target lateral migration velocity field determining module, for being located at around the target well in the purpose work area
The standard well of finishing drilling obtains the well individual-layer data that the standard well location sets place, and according to the well individual-layer data, to described initial
The initial lateral migration velocity field of the first subterranean formation zone is corrected processing in depth migration velocity field, obtains first stratum
The target lateral migration velocity field in region；
The migration imaging result determining module, for the target longitudinal wave migration velocity field and mesh according to first subterranean formation zone
Mark the initial vertical misalignment velocity field of the second subterranean formation zone in lateral shift speed field and the initial depth migration velocity field
With initial lateral migration velocity field, pre-stack depth migration processing is carried out to the primary earthquake data, obtains the purpose work area
Target depth migration imaging result；Wherein, second subterranean formation zone indicates target zone and position described in the purpose work area
Subterranean formation zone below the target zone；
Drilling well orientation determining module is used for based on the target depth migration imaging as a result, determining the target well in institute
State the drilling well orientation on target zone.
10. a kind of device in determining drilling well orientation, including memory, processor and the calculating being stored on the memory
Machine program, which is characterized in that the memory storage has the primary earthquake data and initial depth migration velocity field in purpose work area,
And target well location sets the vertical seismic profiling (VSP) velocity information at place in the purpose work area；Wherein, the target well is drilled into described
The target zone coboundary in purpose work area；The computer program executes following steps when being run by the processor:
Based on the vertical seismic profiling (VSP) velocity information, in the initial depth migration velocity field the first subterranean formation zone it is initial
Vertical misalignment velocity field is corrected processing, obtains the target vertical misalignment velocity field of first subterranean formation zone；Wherein, described
First subterranean formation zone indicates to be located at the subterranean formation zone above the target zone in the purpose work area；
For the standard well for being located at finishing drilling around the target well in the purpose work area, obtains the standard well location and set place
Well individual-layer data, and according to the well individual-layer data, in the initial depth migration velocity field the first subterranean formation zone it is initial
Lateral shift speed field is corrected processing, obtains the target lateral migration velocity field of first subterranean formation zone；
According to the target longitudinal wave migration velocity field of first subterranean formation zone and target lateral migration velocity field and described initial
The initial vertical misalignment velocity field of second subterranean formation zone and initial lateral migration velocity field in depth migration velocity field, to described first
Beginning seismic data carries out pre-stack depth migration processing, obtains the target depth migration imaging result in the purpose work area；Wherein, institute
Stating the second subterranean formation zone indicates target zone described in the purpose work area and the subterranean formation zone below the target zone；
Based on the target depth migration imaging as a result, determining drilling well orientation of the target well on the target zone.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810946533.7A CN109143375B (en) | 2018-08-20 | 2018-08-20 | Method and device for determining drilling direction |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810946533.7A CN109143375B (en) | 2018-08-20 | 2018-08-20 | Method and device for determining drilling direction |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109143375A true CN109143375A (en) | 2019-01-04 |
| CN109143375B CN109143375B (en) | 2020-06-09 |
Family
ID=64790203
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810946533.7A Active CN109143375B (en) | 2018-08-20 | 2018-08-20 | Method and device for determining drilling direction |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109143375B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111624662A (en) * | 2019-02-28 | 2020-09-04 | 中国石油天然气集团有限公司 | Speed field correction method, device and equipment |
| CN112083484A (en) * | 2019-06-12 | 2020-12-15 | 中国石油天然气集团有限公司 | Method and device for detecting true thickness of stratum |
| CN112433247A (en) * | 2020-11-17 | 2021-03-02 | 中国石油化工股份有限公司 | Method and device for adjusting position of stratum to be drilled while drilling |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105842736A (en) * | 2016-05-27 | 2016-08-10 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for building stratum velocity model |
| US20160356909A1 (en) * | 2014-12-23 | 2016-12-08 | Hall Iburton Energy Service, Inc. | Determining p-wave azimuthal anisotropy from walkaround vsp with offset dependent slowness corrections |
| WO2017007738A1 (en) * | 2015-07-05 | 2017-01-12 | Schlumberger Technology Corporation | Borehole trajectory via multi-component borehole seismic receiver |
| WO2017044103A1 (en) * | 2015-09-10 | 2017-03-16 | Halliburton Energy Services, Inc. | Passive ranging using acoustic energy originating from a target wellbore |
| CN107576987A (en) * | 2017-08-07 | 2018-01-12 | 中国石油天然气股份有限公司 | A kind of design method of carbonate beading reservoir-level well track |
| CN108008469A (en) * | 2016-10-28 | 2018-05-08 | 中石化石油工程技术服务有限公司 | The horizontal well geosteering modeling method that well shake combines |
| CN108397130A (en) * | 2018-02-02 | 2018-08-14 | 中国石油天然气集团有限公司 | Boring method and device |
-
2018
- 2018-08-20 CN CN201810946533.7A patent/CN109143375B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20160356909A1 (en) * | 2014-12-23 | 2016-12-08 | Hall Iburton Energy Service, Inc. | Determining p-wave azimuthal anisotropy from walkaround vsp with offset dependent slowness corrections |
| WO2017007738A1 (en) * | 2015-07-05 | 2017-01-12 | Schlumberger Technology Corporation | Borehole trajectory via multi-component borehole seismic receiver |
| WO2017044103A1 (en) * | 2015-09-10 | 2017-03-16 | Halliburton Energy Services, Inc. | Passive ranging using acoustic energy originating from a target wellbore |
| CN105842736A (en) * | 2016-05-27 | 2016-08-10 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for building stratum velocity model |
| CN108008469A (en) * | 2016-10-28 | 2018-05-08 | 中石化石油工程技术服务有限公司 | The horizontal well geosteering modeling method that well shake combines |
| CN107576987A (en) * | 2017-08-07 | 2018-01-12 | 中国石油天然气股份有限公司 | A kind of design method of carbonate beading reservoir-level well track |
| CN108397130A (en) * | 2018-02-02 | 2018-08-14 | 中国石油天然气集团有限公司 | Boring method and device |
Non-Patent Citations (1)
| Title |
|---|
| 史鸿祥 等: "基于随钻地震测井的地震导向钻井技术", 《石油勘探与开发》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111624662A (en) * | 2019-02-28 | 2020-09-04 | 中国石油天然气集团有限公司 | Speed field correction method, device and equipment |
| CN112083484A (en) * | 2019-06-12 | 2020-12-15 | 中国石油天然气集团有限公司 | Method and device for detecting true thickness of stratum |
| CN112433247A (en) * | 2020-11-17 | 2021-03-02 | 中国石油化工股份有限公司 | Method and device for adjusting position of stratum to be drilled while drilling |
| CN112433247B (en) * | 2020-11-17 | 2022-09-02 | 中国石油化工股份有限公司 | While-drilling adjusting method and device for position of stratum to be drilled |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109143375B (en) | 2020-06-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106324676B (en) | A kind of method and device of determining fault sealing property | |
| CN104818982B (en) | A kind of fractured-vuggy reservoir connectedness quantitative analysis method and device | |
| CN108252710B (en) | Method and device for determining water saturation of oil reservoir | |
| CN109143375A (en) | A kind of method and device in determining drilling well orientation | |
| CN104755960A (en) | Improving velocity models for processing seismic data based on basin modeling | |
| US10928533B2 (en) | Identifying potential hydrocarbon traps in a subterranean region using recursive anisotropic erosion of seismic data | |
| CN106761672B (en) | Every interlayer inter-wells predicting method and device | |
| CN106873030B (en) | The determination method and apparatus of Fault Seal | |
| CN104391319B (en) | A kind of determination method and device of seismic data acquisition system | |
| CN104375178B (en) | A kind of carbonate rock fractured cave method for predicting reservoir and device | |
| CN106970423B (en) | A kind of method and device of determining shale gas stratigraphic structure model | |
| CN107179546B (en) | A kind of method and device of determining velocity field | |
| US8593462B2 (en) | Method and system of correlating a measured log to a predicted log | |
| CN107144878A (en) | A kind of fault recognizing method and device | |
| CN108035708A (en) | A kind of method and device for removing bed boundary back wave | |
| CN104614762A (en) | Loose sandstone gas reservoir boundary determination method and device | |
| CN107656309B (en) | A kind of method and device of determining fault throw | |
| CN106772599A (en) | A kind of method and device for calculating formation shear speed | |
| CN106772596A (en) | A kind of method and device for determining pre-stack time migration velocity field | |
| CN107765318B (en) | A kind of method and device of determining geologic horizon | |
| CN107368688A (en) | The determination method and apparatus of marine facies single sand body | |
| CN106908842A (en) | A kind of method and device for determining individual well sedimentary micro curve | |
| BR112019006417B1 (en) | METHOD TO AVOID GEOLOGICAL FORMATION LIMITS, AND INSTRUCTION STORAGE DEVICE | |
| CN109143338A (en) | A kind of method and device that seismic data merges | |
| EP3685190B1 (en) | Processing a 4d seismic signal based on noise model |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |