CN107589457B - The method and apparatus of well-log information scaleup - Google Patents
The method and apparatus of well-log information scaleup Download PDFInfo
- Publication number
- CN107589457B CN107589457B CN201710598821.3A CN201710598821A CN107589457B CN 107589457 B CN107589457 B CN 107589457B CN 201710598821 A CN201710598821 A CN 201710598821A CN 107589457 B CN107589457 B CN 107589457B
- Authority
- CN
- China
- Prior art keywords
- well
- scaleup
- log information
- lithology
- lift height
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention provides a kind of method and apparatus of well-log information scaleup, wherein this method comprises: the well-log information based on reservoir to be measured, obtains at least one lithology lift height of reservoir to be measured；Using lithology lift height as the time window length of scaleup, Bacchus average computation is carried out to the well-log information within the scope of lithology lift height, obtains the well-log information within the scope of lithology lift height after scaleup.In embodiments of the present invention, it can preferably match with the earthquake scale feature at given depth position, the numberical range of elastic parameter is more smoothly concentrated after scaleup.In conjunction with the seismic data of well-log information and the reservoir to be measured after scaleup, well shake calibration can be carried out to the reservoir.
Description
Technical field
The present invention relates to technical field of geological exploration, in particular to a kind of method and apparatus of well-log information scaleup.
Background technique
Well-log information is the capsule information that can not be lacked in geologic information analytic process, and the well shake based on well-log information is fine
Calibration, reservoir modeling and Rock physical analysis are prestack Geological deduction, reservoir prediction is laid a good foundation.However, since reservoir is non-
The influence of homogenieity, working frequency and observed pattern, such as: the working frequency of Sonic Logging Data is about 2-20KHz,
Wavelet wavelength is very short, can detect than relatively thin stratum；The dominant frequency of surface-seismic data is generally 10-100Hz, earthquake wavelength
Very long, range can only detect very thick stratum about between several meters to tens of rice, and the formation thickness of seismic wave detection is far longer than
Log well the formation thickness of detection, therefore, well logging scale and earthquake scale both different measurement scales can be to rock physics
The analysis of feature has large effect.Since reflection shooting is the effective way for grasping subsurface three-dimensional spatial information, and survey
Well data can accurately describe the petrophysics property at wellbore locations, therefore, can by the data measured by seismic prospecting with
And data measured by well-log information combine and are analyzed.Therefore, well-log information can be subjected to scaleup, then will
Well-log information and seismic exploration data after scaleup are combined, the matching to well shake relationship is carried out under earthquake scale
Analysis, i.e. well shake calibration, this is of great significance to development reservoir prediction under earthquake scale.
Summary of the invention
The present invention provides a kind of method and apparatus of well-log information scaleup, and well is carried out under earthquake scale to reach
The matching analysis of shake relationship and the purpose of reservoir prediction.
It may include: based on reservoir to be measured the embodiment of the invention provides a kind of method of well-log information scaleup
Well-log information, obtains the result of log interpretation of the reservoir to be measured, and the result of log interpretation includes: the layering of at least one lithology
Thickness；Using the lithology lift height as the time window length of scaleup, to the well logging within the scope of the lithology lift height
Data carries out Bacchus average computation, obtains the well-log information within the scope of the lithology lift height after scaleup.
In one embodiment, after obtaining the well-log information within the scope of the lithology lift height after scaleup,
It can also include: the seismic data in conjunction with well-log information and the reservoir to be measured after the scaleup, to the storage to be measured
Layer carry out well shake calibration.
In one embodiment, the well-log information within the scope of the lithology lift height after scaleup is obtained, comprising:
Determine the sampling number of the lithology lift height；
Well-log information after each sampled point scaleup is calculated in the following way: using current sampling point as institute
The central point for stating time window length is adopted using the sampling number as the time window length according to each in the time window length
Elastic parameter corresponding to sampling point carries out Bacchus average computation to the well-log information of the current sampling point, obtains described current
Well-log information after sampled point scaleup.
In one embodiment, the sampling number of the lithology lift height is determined, comprising:
According to the lithology lift height and the depth sampling interval of well-log information, it is thick that the lithology layering is calculated
The sampling number of degree.
In one embodiment, Bacchus is carried out according to well-log information of the following formula to the current sampling point averagely to count
It calculates, the well-log information after obtaining the current sampling point scaleup:
ρB=<ρ>
In above formula,<>indicates to elastic parameter according to percent by volume
It is weighted and averaged；
In above formula, ρBAverag density after indicating the current sampling point scaleup, VSBIndicate the current sampling point
After scaleup it is vertical on shear wave velocity, VPBVelocity of longitudinal wave on indicating vertical after the current sampling point scaleup, ρ table
Averag density before showing the current sampling point scaleup, VSCross on indicating vertical before the current sampling point scaleup
Wave velocity, VPVelocity of longitudinal wave on indicating vertical before the current sampling point scaleup.
The embodiment of the invention also provides a kind of devices of well-log information scaleup, may include:
Explanation results determining module obtains the well logging solution of the reservoir to be measured for the well-log information based on reservoir to be measured
It releases as a result, the result of log interpretation includes: at least one lithology lift height；
Well-log information determining module, for the time window length using the lithology lift height as scaleup, to described
Well-log information within the scope of lithology lift height carries out Bacchus average computation, obtains scale within the scope of the lithology lift height
Well-log information after roughening.
In one embodiment, described device can also include:
Well shake demarcating module, for obtain the well-log information within the scope of the lithology lift height after scaleup it
Afterwards, in conjunction with the seismic data of well-log information and the reservoir to be measured after the scaleup, well is carried out to the reservoir to be measured
Shake calibration.
In one embodiment, the well-log information determining module includes: sampling number determination unit, described in determining
The sampling number of lithology lift height；Average calculation unit, it is thick for each sampling point scale to be calculated in the following way
Well-log information after change: using current sampling point as the central point of the time window length, using the sampling number as it is described when
Window length provides the well logging of the current sampling point according to elastic parameter corresponding to each sampled point in the time window length
Material carries out Bacchus average computation, the well-log information after obtaining the current sampling point scaleup.
In one embodiment, the sampling number determination unit is used to be provided according to the lithology lift height and well logging
The sampling number of the lithology lift height is calculated in the depth sampling interval of material.
In one embodiment, the well-log information determining module is specifically used for according to following formula to the present sample
The well-log information of point carries out Bacchus average computation, the well-log information after obtaining the current sampling point scaleup:
ρB=<ρ>
In above formula,<>indicates to elastic parameter according to percent by volume
It is weighted and averaged；
In above formula, ρBAverag density after indicating the current sampling point scaleup, VSBIndicate the current sampling point
After scaleup it is vertical on shear wave velocity, VPBVelocity of longitudinal wave on indicating vertical after the current sampling point scaleup, ρ table
Averag density before showing the current sampling point scaleup, VSCross on indicating vertical before the current sampling point scaleup
Wave velocity, VPVelocity of longitudinal wave on indicating vertical before the current sampling point scaleup.
In embodiments of the present invention, due to the seimic wave velocity in stratum from shallowly to it is deep be variation, in certain depth
In range, the variation of seismic wave is mainly as caused by variation of lithological.Therefore, in this application, it is layered according to lithology to well logging
Curve carries out roughening treatment, specifically, the well-log information based on reservoir to be measured, obtains at least one lithology layering of reservoir to be measured
Thickness；Using lithology lift height as the time window length of scaleup, the well-log information within the scope of lithology lift height is carried out
Bacchus average computation obtains the well-log information within the scope of lithology lift height after scaleup.Using aforesaid way, Neng Gougeng
It matches with the earthquake scale feature at given depth position well, the numberical range of elastic parameter is more smooth after scaleup
It concentrates.In conjunction with the seismic data of well-log information and the reservoir to be measured after scaleup, well shake can be carried out to the reservoir
Calibration.
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 the method flow diagram of well-log information scaleup provided by the present application；
Fig. 2 is the log and its lithology layering before and after scaleup provided by the present application, wherein first row indicates deep
Schematic diagram is spent, secondary series indicates that the density schematic diagram before scaleup, third column indicate the density schematic diagram after scaleup, the
Four column indicate scaleups before velocity of longitudinal wave schematic diagram, the 5th column indicate scaleup after velocity of longitudinal wave schematic diagram, the 6th
Column indicate that the shear wave velocity schematic diagram before scaleup, the 7th column indicate the shear wave velocity schematic diagram after scaleup, the 8th column
Indicate lithology hierarchical diagram, the 9th column indicate reservoir lithology schematic diagram to be measured；
Fig. 3 is a kind of structural block diagram of the device of well-log information scaleup provided by the present application.
Specific embodiment
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.
Log is carried out in general, can use Bacchus (can hereinafter be referred to as Backus) average method
Scaleup.Scaleup refers to that carrying out Backus to the data in certain thickness is averaged, thus by the thin layer of smaller scale
Medium is equivalent to the elastic fluid of large scale.When selecting thickness parameter, traditionally there are two types of ways.First method is to adopt
It is handled with fixed window, that is, identical thickness is all made of to whole log；Second method be using sliding window scan come
Processing, that is, seismic wave wavelength is determined according to the velocity amplitude at log current depth, determines that scale is thick according to the earthquake wavelength
The time window length of change, that is, be weighted and averaged the length of window.
Since seimic wave velocity to deep from shallowly gradually increasing, the wavelength of seismic wave can also become at different depth
Change.Thus, in the first method, ruler all is carried out using identical thickness to whole log under the conditions of different depth
Degree roughening is unreasonable.In the second approach, the speed of a seismic data reference frequency and current point is needed to calculate
Seismic wave wavelength since frequency of seismic wave is variation, thus selects a fixed reference frequency only to determine time window length
It is unreasonable；In addition, adjacent spots value changes greatly, and the when window of adjacent spots is long since log is influenced by noise
Degree differs greatly, and therefore, second method is also unreasonable.
Defect when in view of using above two method choice thickness parameter, inventors herein proposes based on well log interpretation knot
Lithology lift height in fruit carries out Backus average computation to well-log information, carries out scaleup to well-log information to reach
Purpose.Based on this, a kind of method of well-log information scaleup is proposed, as shown in Figure 1, may comprise steps of:
S101: the well-log information based on reservoir to be measured obtains the result of log interpretation of the reservoir to be measured, the well logging solution
Releasing result may include: at least one lithology lift height.
Seimic wave velocity in stratum from shallowly to it is deep be variation, in certain depth bounds, the change of seimic wave velocity
Change mainly as caused by variation of lithological.It, can be according to the layering pair of the lithology of well log interpretation in one embodiment of the application
The division of log progress different lithology thickness range.
Specifically, lithology to sentence knowledge standard as follows:
(1) coal seam: the interval transit time in coal seam is bigger, generally in 300us/m or more；
(2) limestone: the interval transit time of limestone is smaller, in 190us/m hereinafter, and photoelectric coefficient > 3.5Ba/e；
(3) dolomite: dolomitic interval transit time is close with limestone, in 190us/m hereinafter, but photoelectric coefficient < 3.5Ba/
e；
(4) dry sand rock: the interval transit time of dry sand rock is between 235~300.Natural gamma value is generally less than 85API.Shale
Content < 25%；
(5) gas sand: the interval transit time of gas sand is more slightly lower than dry sand rock, generally between 240~300.Natural gal
Horse value is generally less than 85API, shale content < 25%.Gas saturation < 50%；
(6) gas sandstone: the interval transit time of gas sandstone is more slightly lower than gas sand, generally between 250~300.Natural gamma
Value is generally less than 85API, shale content < 25%.Gas saturation < 50%；
(7) argillaceous sandstone: the interval transit time and dry sand rock of argillaceous sandstone are close, generally between 235~300.Natural gal
Horse value is between 85~120API, shale content 25~50%；
(8) Sandy Silt: the interval transit time and dry sand rock of Sandy Silt are close, generally between 235~300.Natural gal
Horse value is in 120API or more, shale content 50~75%；
(9) mud stone: the interval transit time and sandstone of mud stone have greater overlap, generally between 250~300.Natural gamma value
In 120API or more, shale content is greater than 75%；
Therefore, in log, with interval transit time curve can Direct Recognition coal seam (time difference be greater than 300us/m), but
Sand, mud stone are (time difference is between 230~300us/m) stacked on top of each other, and limestone and dolomite are also stacked on top of each other, and (time difference is respectively less than
230us/m), therefore single all lithology cannot be distinguished with the curve.
Under the premise of determining that lithology is sand shale with interval transit time curve, knowledge can be directly sentenced with shale content curve
Sand shale；It, can be straight with photoelectric coefficient curve under the premise of determining that lithology is limestone or is dolomitic with interval transit time curve
It connects and distinguishes both curves.
In this application, it is layered according to lithology and roughening treatment is carried out to log, it is flat can preferably to meet Backus
Calculate layered thickness choose primary condition, can preferably with the earthquake scale feature kissing at given depth position
It closes, the numberical range of elastic parameter is more smoothly concentrated after scaleup.Compared with the conventional method, this method does not need earthquake number
It is a more reasonably selection according to reference frequency.
S102: using the lithology lift height as the time window length of scaleup, to the lithology lift height range
Interior well-log information carries out Backus average computation, obtains the well logging money within the scope of the lithology lift height after scaleup
Material.
It, can be using the lithology lift height as scaleup after obtaining at least one described lithology lift height
Time window length, within the scope of the lithology lift height well-log information carry out Backus average computation, obtain the lithology
Well-log information within the scope of lift height after scaleup.Backus averaging method is: assuming that stratum is by a series of thin layer
A series of thin layer can be averagely equivalent to a uniform layer by composition, Backus, it can the thin layer of smaller scale is situated between
Matter is equivalent to the elastic fluid of large scale.Backus averagely obtains one with a series of elastic parameter of thin layers and density
Elastic constant even, after weighted average, recycles the elastic constant to calculate equivalent velocity of longitudinal wave, shear wave velocity and arithmetic
Density after average.Specifically, may comprise steps of:
S2-1: the sampling number of the lithology lift height is determined.
The lithology point can be calculated according to the lithology lift height and the depth sampling interval of well-log information
The sampling number of thickness degree.
Some lithology is layered, using the depth sampling interval d of its thickness h and log, calculates the lithology point
Backus sampling number L in layer, specifically, formula is as follows:
In this application, above-mentioned depth sampling interval d can be 0.125m.
S2-2: the well-log information after each sampled point scaleup is calculated in the following way: with current sampling point
As the central point of the time window length, using the sampling number as the time window length, according in the time window length
Elastic parameter corresponding to each sampled point carries out Backus average computation to the well-log information of the current sampling point, obtains institute
Well-log information after stating current sampling point scaleup.
In the present embodiment, it is adopted according to the point centered on current sampling point, described in the current sampling point is corresponding up and down
It is average to carry out Backus to the well-log information of the current sampling point for elastic parameter corresponding to the sampled point of the half of number of samples
It calculates, the well-log information after obtaining the current sampling point scaleup.
Specifically, Backus average computation can be carried out according to well-log information of the following formula to the current sampling point,
Well-log information after obtaining the current sampling point scaleup:
In above formula,<>indicates to elastic parameter according to percent by volume
It is weighted and averaged；
In above formula, ρBAverag density after indicating the current sampling point scaleup, VSBIndicate the current sampling point
After scaleup it is vertical on shear wave velocity, VPBVelocity of longitudinal wave on indicating vertical after the current sampling point scaleup, ρ table
Averag density before showing the current sampling point scaleup, VSCross on indicating vertical before the current sampling point scaleup
Wave velocity, VPVelocity of longitudinal wave on indicating vertical before the current sampling point scaleup.
In the present embodiment,<>can indicate to the point centered on current sampling point, above and below the current sampling point
Value corresponding to the sampled point of the half of the corresponding sampling number is weighted and averaged, such as: when the sampling number is 50
When, centered on current sampling point, upward 25 sampled points of current sampling point, corresponding to downward 25 sampled points of current sampling point
Elastic parameter, Backus average computation is carried out to the well-log information of the current sampling point, the survey after obtaining scaleup
Well profile.
It further, can be with after obtaining the well-log information within the scope of the lithology lift height after scaleup
In conjunction with the seismic data of well-log information and the reservoir to be measured after the scaleup, well shake mark is carried out to the reservoir to be measured
It is fixed.
It is specifically described below with reference to method of the specific embodiment to above-mentioned well-log information scaleup, however
It is worth noting that, the specific embodiment merely to the present invention is better described, does not constitute improper limitations of the present invention.
Step 1: the log of Different Strata attribute is obtained in field use logger measurement, according to conventional well logging
Explain that process obtains well log interpretation achievement.Well log interpretation achievement include lithology layering and its different lithology layering initial depth and
Lithology lift height.Log and its lithology layering being illustrated in figure 2 before and after scaleup.The 1st column are depth in Fig. 2,
2nd, 4,6 column are density before scaleup, velocity of longitudinal wave, Shear Wave Velocity Well Logging curve respectively, and the 8th, 9 column are lithology respectively
Layering and its corresponding lithology.
Step 2: the initial depth and thickness of different lithology layering are obtained according to the 8th, 9 column, according to the layering of current lithology
The depth of thickness and log calculates the Backus average length parameter in current lithology layering using interval；
Step 3: Backus average computation being carried out to log using formula (1), in current lithology depth of seam division range
It is interior, using the average length parameter in this depth bounds, obtain log after the scaleup in this depth bounds.
Step 4: Backus average computation is carried out to the log in the layering of all lithology, it is entire bent after being roughened
Line.As the 3rd, 5,7 column are density after Backus average, velocity of longitudinal wave, Shear Wave Velocity Well Logging curve respectively in Fig. 2.
Based on the same inventive concept, a kind of device of well-log information scaleup is additionally provided in the embodiment of the present invention, such as
Described in the following examples.The principle and well-log information scaleup solved the problems, such as due to the device of well-log information scaleup
Method is similar, therefore the implementation of the device of well-log information scaleup may refer to the reality of the method for well-log information scaleup
It applies, overlaps will not be repeated.Used below, the software of predetermined function may be implemented in term " unit " or " module "
And/or the combination of hardware.Although device described in following embodiment is preferably realized with software, hardware or soft
The realization of the combination of part and hardware is also that may and be contemplated.Fig. 3 is the well-log information scaleup of the embodiment of the present invention
A kind of structural block diagram of device, as shown in figure 3, may include: explanation results determining module 301, well-log information determining module
302, the structure is illustrated below.
Explanation results determining module 301 can be used for the well-log information based on reservoir to be measured, obtain the reservoir to be measured
Result of log interpretation, the result of log interpretation include: at least one lithology lift height；
Well-log information determining module 302 can be used for the time window length using the lithology lift height as scaleup,
Backus average computation is carried out to the well-log information within the scope of the lithology lift height, obtains the lithology lift height range
Well-log information after interior scaleup.
In one embodiment, described device can also include: well shake demarcating module, can be used for obtaining the lithology
After well-log information within the scope of lift height after scaleup, in conjunction with well-log information after the scaleup and described to be measured
The seismic data of reservoir carries out well shake calibration to the reservoir to be measured.
In one embodiment, the well-log information determining module may include: sampling number determination unit, can be used for
Determine the sampling number of the lithology lift height；Average calculation unit can be used for being calculated in the following way each
Well-log information after sampled point scaleup: using current sampling point as the central point of the time window length, with the sampled point
Number is used as the time window length, is currently adopted according to elastic parameter corresponding to each sampled point in the time window length to described
The well-log information of sampling point carries out Backus average computation, the well-log information after obtaining the current sampling point scaleup.
In one embodiment, the sampling number determination unit can be used for according to the lithology lift height and survey
The sampling number of the lithology lift height is calculated in the depth sampling interval of well data.
In one embodiment, the well-log information determining module specifically can be used for according to following formula to described current
The well-log information of sampled point carries out Backus average computation, the well-log information after obtaining the current sampling point scaleup:
ρB=<ρ>
In above formula,<>indicates to elastic parameter according to percent by volume
It is weighted and averaged；
In above formula, ρBAverag density after indicating the current sampling point scaleup, VSBIndicate the current sampling point
After scaleup it is vertical on shear wave velocity, VPBVelocity of longitudinal wave on indicating vertical after the current sampling point scaleup, ρ table
Averag density before showing the current sampling point scaleup, VSCross on indicating vertical before the current sampling point scaleup
Wave velocity, VPVelocity of longitudinal wave on indicating vertical before the current sampling point scaleup.
It, can be with automatic implementation using the embodiment of the device of well-log information scaleup provided by the various embodiments described above
The processing method of the well-log information scaleup carries out well shake calibration to reservoir to be measured, it may not be necessary to implement the tool of personnel
Body participates in, and directly can shake calibration result by output well, simple and quick, effectively increases user experience.
In the device of the well-log information scaleup, the survey within the scope of the lithology lift height after scaleup is obtained
The extension of well data, the sampling number for determining the lithology lift height, Backus average computation mode is referred to aforementioned side
The associated description of method.
Although mentioning the mode of well-log information scaleup, the method for determination of sampling number, Backus in teachings herein
The description such as average computation mode, still, the application is not limited to be situation described in the embodiment of the present application.Certain rows
Embodiment modified slightly can also be real in industry standard or the practice processes described using customized mode or embodiment
Existing above-described embodiment is identical, after equivalent or close or deformation it is anticipated that implementation result.Using these modifications or deformed number
According to the embodiment of the acquisitions such as roughening/determination/calculating, still may belong within the scope of the optional embodiment of the application.
Although this application provides the method operating procedure as described in embodiment or flow chart, based on conventional or noninvasive
The means for the property made may include more or less operating procedure.The step of enumerating in embodiment sequence is only numerous steps
One of execution sequence mode, does not represent and unique executes sequence.It, can be with when device in practice or end product execute
It is executed according to embodiment or method shown in the drawings sequence or parallel executes (such as parallel processor or multiple threads
Environment, even distributed data processing environment).The terms "include", "comprise" or its any other variant are intended to
Non-exclusive inclusion, so that process, method, product or equipment including a series of elements are not only wanted including those
Element, but also including other elements that are not explicitly listed, or further include for this process, method, product or equipment
Intrinsic element.In the absence of more restrictions, be not precluded include the process, method of the element, product or
There is also other identical or equivalent elements in person's equipment.
Unit, device or module that above-described embodiment illustrates etc. can specifically realize by computer chip or entity, or
It is realized by the product with certain function.For convenience of description, various modules point are divided into function when describing apparatus above
It does not describe.It certainly, when implementing the application can the function of each module is real in the same or multiple software and or hardware
It is existing, the module for realizing same function can also be realized by the combination of multiple submodule or subelement etc..Dress described above
Set that embodiment is only schematical, for example, the division of the unit, only a kind of logical function partition, in actual implementation
There may be another division manner, such as multiple units or components can be combined or can be integrated into another system or one
A little features can be ignored, or not execute.Another point, shown or discussed mutual coupling or direct-coupling or communication link
Connecing can be through some interfaces, the indirect coupling or communication connection of device or unit, can be electrical property, mechanical or other shapes
Formula.
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 that its inside includes can also be considered as in hardware component.Or
Person even, can will be considered as realizing the device of various functions either the software module of implementation method can be hardware again
Structure in component.
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, class etc..The application can also be practiced in a distributed computing environment, in these distributed computing environments,
By executing task by the connected remote processing devices of communication network.In a distributed computing environment, program module can
To be located in the local and remote computer storage media including storage equipment.
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, the computer software product
It can store in storage medium, such as ROM/RAM, magnetic disk, CD, including some instructions are used so that a computer equipment
(can be personal computer, mobile terminal, server or the network equipment etc.) executes each embodiment of the application or implementation
Method described in certain parts of example.
Each embodiment in this specification is described in a progressive manner, the same or similar portion between each embodiment
Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.The application can be used for crowd
In mostly general or special purpose computing system environments or configuration.Such as: personal computer, server computer, handheld device or
Portable device, laptop device, multicomputer system, microprocessor-based system, set top box, programmable electronics set
Standby, network PC, minicomputer, mainframe computer, distributed computing environment including any of the above system or equipment etc..
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 of well-log information scaleup characterized by comprising
Based on the well-log information of reservoir to be measured, the result of log interpretation of the reservoir to be measured, the result of log interpretation packet are obtained
It includes: at least one lithology lift height；
Using the lithology lift height as the time window length of scaleup, the well logging within the scope of the lithology lift height is provided
Material carries out Bacchus average computation, obtains the well-log information within the scope of the lithology lift height after scaleup.
2. the method as described in claim 1, which is characterized in that obtaining within the scope of the lithology lift height after scaleup
Well-log information after, the method also includes:
In conjunction with the seismic data of well-log information and the reservoir to be measured after the scaleup, well is carried out to the reservoir to be measured
Shake calibration.
3. the method as described in claim 1, which is characterized in that obtain within the scope of the lithology lift height after scaleup
Well-log information, comprising:
Determine the sampling number of the lithology lift height；
Well-log information after each sampled point scaleup is calculated in the following way: using current sampling point as it is described when
The central point of window length, using the sampling number as the time window length, according to each sampled point in the time window length
Corresponding elastic parameter carries out Bacchus average computation to the well-log information of the current sampling point, obtains the present sample
Well-log information after point scale roughening.
4. method as claimed in claim 3, which is characterized in that determine the sampling number of the lithology lift height, comprising:
According to the lithology lift height and the depth sampling interval of well-log information, the lithology lift height is calculated
Sampling number.
5. method as claimed in claim 3, which is characterized in that according to following formula to the well-log information of the current sampling point
Bacchus average computation is carried out, the well-log information after obtaining the current sampling point scaleup:
ρB=<ρ>
In above formula,<>expression carries out elastic parameter according to percent by volume
Weighted average；
In above formula, ρBAverag density after indicating the current sampling point scaleup, VSBIndicate the present sample point scale
After roughening it is vertical on shear wave velocity, VPBVelocity of longitudinal wave on indicating vertical after the current sampling point scaleup, ρ indicate institute
Averag density before stating current sampling point scaleup, VSShear wave speed on indicating vertical before the current sampling point scaleup
Degree, VPVelocity of longitudinal wave on indicating vertical before the current sampling point scaleup.
6. a kind of device of well-log information scaleup characterized by comprising
Explanation results determining module obtains the well log interpretation knot of the reservoir to be measured for the well-log information based on reservoir to be measured
Fruit, the result of log interpretation include: at least one lithology lift height；
Well-log information determining module, for the time window length using the lithology lift height as scaleup, to the lithology
Well-log information within the scope of lift height carries out Bacchus average computation, obtains scaleup within the scope of the lithology lift height
Well-log information afterwards.
7. device as claimed in claim 6, which is characterized in that further include:
Well shakes demarcating module, for tying after obtaining the well-log information within the scope of the lithology lift height after scaleup
The seismic data of well-log information and the reservoir to be measured after closing the scaleup carries out well shake mark to the reservoir to be measured
It is fixed.
8. device as claimed in claim 6, which is characterized in that the well-log information determining module includes:
Sampling number determination unit, for determining the sampling number of the lithology lift height；
Average calculation unit, for the well-log information after each sampled point scaleup to be calculated in the following way: to work as
Central point of the preceding sampled point as the time window length, using the sampling number as the time window length, according to it is described when window
Elastic parameter corresponding to each sampled point in length carries out Bacchus to the well-log information of the current sampling point and averagely counts
It calculates, the well-log information after obtaining the current sampling point scaleup.
9. device as claimed in claim 8, which is characterized in that the sampling number determination unit is used for according to the lithology point
The sampling number of the lithology lift height is calculated in the depth sampling interval of thickness degree and well-log information.
10. device as claimed in claim 8, which is characterized in that the well-log information determining module is specifically used for according to following
Formula carries out Bacchus average computation to the well-log information of the current sampling point, after obtaining the current sampling point scaleup
Well-log information:
ρB=<ρ>
In above formula,<>expression carries out elastic parameter according to percent by volume
Weighted average；
In above formula, ρBAverag density after indicating the current sampling point scaleup, VSBIndicate the present sample point scale
After roughening it is vertical on shear wave velocity, VPBVelocity of longitudinal wave on indicating vertical after the current sampling point scaleup, ρ indicate institute
Averag density before stating current sampling point scaleup, VSShear wave speed on indicating vertical before the current sampling point scaleup
Degree, VPVelocity of longitudinal wave on indicating vertical before the current sampling point scaleup.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710598821.3A CN107589457B (en) | 2017-07-21 | 2017-07-21 | The method and apparatus of well-log information scaleup |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710598821.3A CN107589457B (en) | 2017-07-21 | 2017-07-21 | The method and apparatus of well-log information scaleup |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107589457A CN107589457A (en) | 2018-01-16 |
| CN107589457B true CN107589457B (en) | 2019-08-06 |
Family
ID=61041693
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710598821.3A Active CN107589457B (en) | 2017-07-21 | 2017-07-21 | The method and apparatus of well-log information scaleup |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107589457B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110646864B (en) * | 2019-10-16 | 2021-06-25 | 中国电建集团成都勘测设计研究院有限公司 | Complex engineering geological rock group classification method for inserting irregular invaded rocks and veins |
-
2017
- 2017-07-21 CN CN201710598821.3A patent/CN107589457B/en active Active
Non-Patent Citations (1)
| Title |
|---|
| "基于Backus等效平均的测井资料尺度粗化方法研究";曹丹平;《石油物探》;20150131;第54卷(第1期);第105-109页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107589457A (en) | 2018-01-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107817535B (en) | The determination method and apparatus of short lap | |
| CN103258091B (en) | Unconventionaloil pool hides the method and device that net horizontal section three-dimensional mechanical models for rock mass is set up | |
| CN112703429A (en) | Machine learning based seismic attribute analysis | |
| Close et al. | Integrated workflows for shale gas and case study results for the Horn River Basin, British Columbia, Canada | |
| US10359529B2 (en) | Singularity spectrum analysis of microseismic data | |
| CN106597544B (en) | The brittle prediction technique of fine and close oil-gas reservoir and device | |
| EA014144B1 (en) | Method and apparatus for predicting the hydrocarbon production of a well location | |
| US20140052377A1 (en) | System and method for performing reservoir stimulation operations | |
| CN103256046A (en) | Unconventional oil and gas reservoir horizontal well section full-fracture-length fracturing parameter analog method and device | |
| US11879316B2 (en) | Geostatistical analysis of microseismic data in fracture modeling | |
| US20220291418A1 (en) | An integrated geomechanics model for predicting hydrocarbon and migration pathways | |
| CN108397130A (en) | Boring method and device | |
| CN107450108A (en) | The determination method and apparatus in dessert area | |
| Askari et al. | A fully integrated method for dynamic rock type characterization development in one of Iranian off-shore oil reservoir | |
| Schneider et al. | Interpretation of fractured zones using seismic attributes—Case study from Teapot Dome, Wyoming, USA | |
| CN107817519A (en) | Basement lithology recognition methods and system based on weight magnetic shake Comprehensive interpretation | |
| CN112114380A (en) | Prediction method for conglomerate oil reservoir | |
| Yousefzadeh et al. | Microseismic 101: monitoring and evaluating hydraulic fracturing to improve the efficiency of oil and gas recovery from unconventional reservoirs | |
| EP3301485B1 (en) | Crosswell microseismic system | |
| CN107589457B (en) | The method and apparatus of well-log information scaleup | |
| CN106842316B (en) | Crack determines method and apparatus | |
| CN108008464A (en) | Crack anisotropism quantitatively characterizing method and its system | |
| CN107807409B (en) | The determination method and apparatus of density of earth formations and resistivity relation | |
| CN107576984B (en) | The determination method and apparatus of window when during seismic amplitude attributes extraction | |
| CN108072901A (en) | A kind of method and system for obtaining accurate hydrostatic pressure and overlying formation pressure |
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 |