CN104793263B - A kind of AUTOMATIC SORTING IN LAYERS BY WELL-LOGGING CURVES value and evaluation method - Google Patents
A kind of AUTOMATIC SORTING IN LAYERS BY WELL-LOGGING CURVES value and evaluation method Download PDFInfo
- Publication number
- CN104793263B CN104793263B CN201510222975.3A CN201510222975A CN104793263B CN 104793263 B CN104793263 B CN 104793263B CN 201510222975 A CN201510222975 A CN 201510222975A CN 104793263 B CN104793263 B CN 104793263B
- Authority
- CN
- China
- Prior art keywords
- curve
- value
- automatic
- evaluation
- well
- 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
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a kind of AUTOMATIC SORTING IN LAYERS BY WELL-LOGGING CURVES value and evaluation method, carries out AUTOMATIC ZONING, automatic value and automatic Evaluation for well logging primitive curve and achievement curve, the method achieve 3 objectives:1. a plurality of log variation characteristic of synthesis and boundary condition carry out AUTOMATIC ZONING;2. according to tracing pattern, the representative value of each bar curve is successively obtained;3. according to value result and discrimination standard, automatically formation properties are made with evaluation.Present invention is primarily aimed at initial results are provided for follow-up well logging interactive interpretation and evaluation, to improve explanation operating efficiency.The present invention is log is squared, resistivity logging 2D invertings initialization etc. more can be applied directly.
Description
Technical field
The present invention relates to geophysical well logging technology field, specifically a kind of AUTOMATIC SORTING IN LAYERS BY WELL-LOGGING CURVES value and comment
Valency method.
Background technology
With the fast development of human civilization, the demand of the energy is increasingly increased, and modern society's consumer demand
Also become more diverse, demand of industrial production product can fast transformation, this all increases the demand to oil and natural gas.For
Effective detection oil-gas Layer, during oilfield prospecting developing, well-log information digital processing turns into vast with integrated interpretation
One of main task of well log interpretation personnel, in the process, explain personnel by the response according to original log, Yi Jijing
The achievement curve crossed after well-log information digital processing, with reference to other data such as drilling well, geology, well logging, rock core, carry out stratum
Evaluation, it is intended to determine reservoir property(Oil reservoir, gas-bearing formation, water layer or dried layer etc.), and obtain reservoir lithology, physical property(Porosity, ooze
Saturating rate)And oiliness(Oil saturation)Deng reservoir parameter.
Past is generally divided into three due to the reasons such as the limitation of computer software technology, traditional logging Reservoir Evaluation method
Step, 1. carry out digital processing to log data using computer and export drafting pattern;2. in a manner of desktop drawing operation
Reservoir division, curve readings and reservoir property evaluation are carried out to log;3. manual entry the above results form well logging
Interpretation results form simultaneously prints out.The logging evaluation work of a bite well is completed, generally requires and spends the 2 people time of more than 3 days,
Work period length be present, poor in timeliness, precision is not high, the problems such as easily make mistakes.Nowadays logging evaluation work is at data
Completion can be interacted under computer graphical interface completely by managing integrated interpretation, though but the interactive interpretation under graphical interfaces
Say and be extremely improved in mode of operation, but be still difficult to avoid that most of problem present in explanation by hand.
The content of the invention
Present applicant proposes a kind of AUTOMATIC SORTING IN LAYERS BY WELL-LOGGING CURVES value and evaluation method, it is intended to is well-log information digital processing
A set of algorithm logged well and explained automatically is provided with integrated interpretation software platform, a set of initial automatic explanation knot is provided for interactive interpretation
Fruit, to improve the efficiency of Log Interpretation and precision.
In order to solve the above technical problems, the embodiment of the present application provides a kind of AUTOMATIC SORTING IN LAYERS BY WELL-LOGGING CURVES value and evaluation side
Method, comprise the following steps:
First, AUTOMATIC ZONING
Step S1:Some logs are selected, the weights of each curve are set, are one further according to Curves compilation ruled synthesis
Bar synthesis layering curve;
Step S2:Smothing filtering is carried out with moving weighted average method to comprehensive layering curve;
Step S3:Resultant curve application bathmometry is layered;
Step S4:A crucial curve is selected, is layered using cutoff method, while utilizes its curve in interface
Monotonicity determines the top bottom attribute at interface;
Step S5:The rule limited according to minimum thickness, resultant curve slope threshold values, crucial curve cutoff, to step
3rd, layer is accepted or rejected caused by 4;
Step S6:Above-mentioned interface is built into non-reservoir units, reservoir units and small layer unit;
2nd, automatic value
Step S7:All flex points and extreme point are obtained with second dervative to every curve structure first derivative;
Step S8:By counting flex point, curve monotonicity judges song at extreme point quantity and flex point in small layer unit
Line morphology, tracing pattern are divided into four kinds, i.e. monotonous curve, sag vertical curve, convex curve and the curve of cyclical fluctuations;
Step S9:Minimum slope point value, average value, minimum, maximum are calculated in small layer unit;
Step S10:Value mode is determined according to tracing pattern, i.e. monotonous curve takes minimum slope point value or average value, ripple
Moving curve is averaged, and sag vertical curve minimalization, convex curve takes maximum, and substratum assignment is carried out to every curve;
3rd, automatic Evaluation
Step S11:Successively according to the standard set up, concluded judgement for substratum.
As the preferred embodiment of this programme, described step S1 selection curves are to user's formation characteristics of concern as
Layer lithology, physical property, oiliness parameter sensitivity, the response of formation characteristics can be reflected on comprehensive layering curve, various characteristics
Response it is separate, or with superposition amplification effect, the weights of all curves have positive and negative attribute.
Curves compilation formula is:
As the preferred embodiment of this programme, it is oblique that rule given described step S5 includes minimum thickness, resultant curve
Rate threshold values, the limitation of crucial curve cutoff.
As the preferred embodiment of this programme, described step S11 criterion includes shale content higher limit, hole
Spend lower limit, water saturation higher limit, and hydrocarbon types.
The one or more technical schemes provided in the embodiment of the present application, have at least the following technical effects or advantages:
AUTOMATIC ZONING, automatic value and automatic Evaluation are carried out to log by computer simulation manual type, namely
Before follow-up interactive layering value and evaluation is carried out, to be provided previously by a set of automatic explanation results;This method extremely accords with
Requirement and the feature of log manual zoning's value are closed, control parameter is flexible, as a result rationally reliable, can be follow-up interaction
Formula, which is explained, greatly reduces workload and raising efficiency.Auto-layering data with the result of evaluation compared with final explanation results,
Coincidence rate can reach more than 80%, and ideally, auto-layering data and the achievement evaluated are even directly available.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
There is the required accompanying drawing used in technology description to be briefly described, it should be apparent that, drawings in the following description are only this hairs
Some bright embodiments, for those of ordinary skill in the art, without having to pay creative labor, can be with
Other accompanying drawings are obtained according to these accompanying drawings.
Fig. 1 is the step flow chart of the embodiment of the present application.
Embodiment
In order to be better understood from above-mentioned technical proposal, below in conjunction with Figure of description and specific embodiment to upper
Technical scheme is stated to be described in detail.
As shown in figure 1, the application present invention carries out well log interpretation appraisal, comprise the following steps:
First, AUTOMATIC ZONING
(1)By achievement curve-porosity after well-log information digital processing(POR), water saturation(SW), shale content
(VSH)Curve weighting synthesizes a resultant curve;
S= SW*MWsw + VSH*MWvsh + POR*MWpor
(2)The resultant curve of synthesis is subjected to smothing filtering(5 moving weighted average methods, it can become according to weight coefficient and be divided into
Multiple ranks);
S=MW1*S[n-2]+MW2*S[n-1]+MW3*S[n]+MW4*S[n+1]+MW5*S[n+2]
(3)The first derivative and second dervative of resultant curve are calculated, and determines flex point depth as bed boundary.
Log is equal interval sampling data, therefore the first derivative of curve and second dervative are respectively:
S’[n]=S[n+1]-S[n]
S”[n]=S’[n+1]-S’[n]
When judging adjacent 2 points of second dervative product to bear, you can determine a flex point interface.
(4)Resultant curve slope at all flex points of statistic procedure (3), obtain reference of the G-bar as roughness
Value, compare the slope cutoff that the slope of curve at flex point calculates with roughness, determine the choice at flex point interface(Reservoir parameter reaches
Just it is layered to certain intensity of variation).
(5)Using shale content curve as crucial curve, cutoff layering is carried out(The friendship of shale content curve and cutoff
Point depth is cutoff interface).
(6)Give up all flex point interfaces that all shale contents are more than shale content cutoff(Mud stone is not finely divided
Layer).
(7)Give up the flex point interface for being unsatisfactory for minimum thickness requirement.
(8)Give up the cutoff interface for being unsatisfactory for minimum thickness requirement again.
(9)According to monotonicity of the shale content curve in interface(Judged by first derivative)Determine the top bottom category at interface
Property, it is reservoir top to subtract curve, and increasing curve is reservoir bottom.
(10)According to the top bottom attribute of bed boundary, build reservoir, reservoir formed between top and bottom, push up bottom and flex point interface it
Between form substratum, give up between bottom and top as non-reservoir(Well logging does not do interpretation and evaluation to non-reservoir).
2nd, automatic value
(11)All curves in substratum are scanned, by first derivative and second derivative method, extreme value points and flex point number are entered
Row counts, and in the form of this judgment curves, is specifically divided into convex curve, sag vertical curve, the curve of cyclical fluctuations and the type of monotonous curve 4.
(12)The maximum value minimum of all curves, average value, minimum slope curve values in statistics layer.
(13)According to tracing pattern, all curves are distinguished with assignment, specially convex curve assigns maximum, and sag vertical curve assigns
Minimum, the curve of cyclical fluctuations assign average value, and monotonous curve assigns minimum slope curve values(Or assign average value and also may be used).
3rd, automatic Evaluation
(14)According to geological layering, the shale content of subsection setup, porosity, water saturation effective thickness lower limit or
Upper limit standard, and hydrocarbonaceous type;
(15)According to above-mentioned standard, the differentiation of well log interpretation conclusion is successively made to substratum.
The above described is only a preferred embodiment of the present invention, any formal limitation not is made to the present invention, though
So the present invention is demonstrated as above with preferred embodiment, but is not limited to the present invention, any to be familiar with this professional technology people
Member, without departing from the scope of the present invention, when the technology contents using the disclosure above make a little change or modification
For the equivalent embodiment of equivalent variations, as long as being the content without departing from technical solution of the present invention, the technical spirit according to the present invention
Any simple modification, equivalent change and modification made to above example, in the range of still falling within technical solution of the present invention.
Claims (2)
1. a kind of AUTOMATIC SORTING IN LAYERS BY WELL-LOGGING CURVES value and evaluation method, it is characterised in that comprise the following steps:
First, AUTOMATIC ZONING
Step S1:Some logs are selected, the weights of each curve are set, are one comprehensive further according to Curves compilation ruled synthesis
Close layering curve;
Step S2:Smothing filtering is carried out with moving weighted average method to comprehensive layering curve;
Step S3:Resultant curve application bathmometry is layered;
Step S4:A crucial curve is selected, is layered using cutoff method, while is increased and decreased using its curve in interface
Property determines the top bottom attribute at interface;
Step S5:The rule limited according to minimum thickness, resultant curve slope threshold values, crucial curve cutoff, to step S3, S4
Caused layer is accepted or rejected;
Step S6:Above-mentioned interface is built into non-reservoir units, reservoir units and small layer unit;
2nd, automatic value
Step S7:All flex points and extreme point are obtained with second dervative to every curve structure first derivative;
Step S8:In small layer unit by count flex point, at extreme point quantity and flex point curve monotonicity come Detection curve shape
State, tracing pattern are divided into four kinds, i.e. monotonous curve, sag vertical curve, convex curve and the curve of cyclical fluctuations;
Step S9:Minimum slope point value, average value, minimum, maximum are calculated in small layer unit;
Step S10:Value mode is determined according to tracing pattern, i.e. monotonous curve takes minimum slope point value or average value, fluctuation song
Line is averaged, and sag vertical curve minimalization, convex curve takes maximum, and substratum assignment is carried out to every curve;
3rd, automatic Evaluation
Step S11:Successively according to the standard set up, concluded judgement for substratum.
2. a kind of AUTOMATIC SORTING IN LAYERS BY WELL-LOGGING CURVES value according to claim 1 and evaluation method, it is characterised in that:Described
Step S11 criterion includes shale content higher limit, PLL, water saturation higher limit, and hydrocarbon types.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510222975.3A CN104793263B (en) | 2015-05-05 | 2015-05-05 | A kind of AUTOMATIC SORTING IN LAYERS BY WELL-LOGGING CURVES value and evaluation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510222975.3A CN104793263B (en) | 2015-05-05 | 2015-05-05 | A kind of AUTOMATIC SORTING IN LAYERS BY WELL-LOGGING CURVES value and evaluation method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104793263A CN104793263A (en) | 2015-07-22 |
CN104793263B true CN104793263B (en) | 2018-03-02 |
Family
ID=53558224
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510222975.3A Active CN104793263B (en) | 2015-05-05 | 2015-05-05 | A kind of AUTOMATIC SORTING IN LAYERS BY WELL-LOGGING CURVES value and evaluation method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104793263B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105225018A (en) * | 2015-11-04 | 2016-01-06 | 上海斐讯数据通信技术有限公司 | Based on Forecasting Methodology and the prediction unit of the mobile phone flow of BP neural network |
CN109182463B (en) * | 2018-09-21 | 2022-04-29 | 博奥生物集团有限公司 | Fluorescence amplification curve inflection point determination method and device |
CN111173505B (en) * | 2018-10-23 | 2023-08-22 | 中国石油天然气股份有限公司 | Method and apparatus for determining a reservoir lower limit |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102454399B (en) * | 2010-10-26 | 2015-01-14 | 中国石油化工股份有限公司 | Method for correcting time difference signals of well logging acoustic waves |
CN102116871B (en) * | 2011-01-24 | 2013-05-01 | 中国石油天然气股份有限公司 | Method for improving capability of seismic inverted data prediction storage layer |
CN102830441A (en) * | 2012-05-16 | 2012-12-19 | 陶冶 | Comprehensive interpretation and evaluation method of hydrocarbon reservoir including data mining |
CN104047597B (en) * | 2014-01-26 | 2016-06-08 | 中国石油大学(华东) | Oily mud shale stratum log standardized method |
CN104295293B (en) * | 2014-10-23 | 2017-04-12 | 中国石油天然气股份有限公司 | Method for acquiring logging density curve |
-
2015
- 2015-05-05 CN CN201510222975.3A patent/CN104793263B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN104793263A (en) | 2015-07-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7873476B2 (en) | Well log correlation | |
US20220170366A1 (en) | 3d in-situ characterization method for heterogeneity in generating and reserving performances of shale | |
CN102116871B (en) | Method for improving capability of seismic inverted data prediction storage layer | |
CN102759748B (en) | The cross analysis Seismic Reservoir Prediction method that based target is analyzed | |
CN104040377A (en) | Integrated workflow or method for petrophysical rock typing in carbonates | |
CN106150477A (en) | A kind of method determining single well controlled reserves | |
CN104899928B (en) | Three-dimensional geological modeling method based on sparse boring point | |
CN103472484A (en) | Horizontal well track optimization method based on RS three-dimensional sensitivity seismic attribution analysis | |
CN103993871A (en) | Method and device for processing well logging information of thin interbed stratums in standardization mode | |
CN105549084A (en) | Method and system for building three-dimensional high-precision velocity model | |
CN104793263B (en) | A kind of AUTOMATIC SORTING IN LAYERS BY WELL-LOGGING CURVES value and evaluation method | |
CA2911107C (en) | Local layer geometry engine with work zone generated from buffer defined relative to a wellbore trajectory | |
CN108952699A (en) | A kind of complicated geological drilling process formation lithology intelligent identification Method | |
CN103376468A (en) | Reservoir parameter quantitative characterization method based on neural network function approximation algorithm | |
CN110988997A (en) | Hydrocarbon source rock three-dimensional space distribution quantitative prediction technology based on machine learning | |
CN112925017A (en) | Fracture type reservoir geological modeling method | |
CN114609675A (en) | Quantitative recovery method for carbonate rock stratum sedimentary micro-landform based on high-frequency cycle | |
CN105676288B (en) | A kind of sand body stacked relation method of discrimination based on reflectance factor analysis | |
US20140345946A1 (en) | Analysis of Geological Objects | |
Singleton | Geophysical data processing, rock property inversion, and geomechanical model building in a Midland Basin development project, Midland/Ector counties, Texas | |
CN108549101A (en) | Seismic data processing technique and device | |
CN113792499A (en) | Deposition characteristic-based dynamic determination method for water-rich property of loose confined aquifer | |
CN110389381A (en) | A kind of Sandbody Reservoirs prediction technique and device based on seismic properties | |
Coste et al. | Data mining techniques for optimizing fast track re-engineering of mature fields | |
RU2480795C2 (en) | Method of identifying structure of reservoirs of geologic strata and machine-readable medium |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP03 | Change of name, title or address | ||
CP03 | Change of name, title or address |
Address after: 100010 Chaoyangmen North Street, Dongcheng District, Dongcheng District, Beijing Co-patentee after: CNOOC (China) Limited Zhanjiang Branch Patentee after: China Offshore Oil Group Co., Ltd. Address before: 100010 China oil tower, 25 Chaoyangmen North Street, Dongcheng District, Beijing Co-patentee before: CNOOC (China) Limited Zhanjiang Branch Patentee before: China National Offshore Oil Corporation |