CN104793263A - Well logging curve automatic tiered value extraction and evaluation method - Google Patents
Well logging curve automatic tiered value extraction and evaluation method Download PDFInfo
- Publication number
- CN104793263A CN104793263A CN201510222975.3A CN201510222975A CN104793263A CN 104793263 A CN104793263 A CN 104793263A CN 201510222975 A CN201510222975 A CN 201510222975A CN 104793263 A CN104793263 A CN 104793263A
- Authority
- CN
- China
- Prior art keywords
- curve
- value
- well logging
- curves
- evaluation
- 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
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 well logging curve automatic tiered value extraction and evaluation method and aims to perform tiering, value extraction and evaluation automatically according to well logging original curves and result curves. By the aid of the method, three specific targets are achieved, in other words, 1, the tiering is implemented automatically according to the change features of multiple well logging curves and boundary conditions comprehensively; 2, according to the curve shape, the representative value of each curve is acquired in the tiered manner; 3, the formation properties are evaluated automatically according to the value extraction result and criterion. According to the method, the initial result is provided for the subsequent well logging interactive interpretation and evaluation so as to improve the interpreting efficiency. The method can be applied to well logging curve square wave promotion and resistivity well logging 2D inversion initialization directly.
Description
Technical field
The present invention relates to Geophysical Logging field, specifically a kind of AUTOMATIC SORTING IN LAYERS BY WELL-LOGGING CURVES value and evaluation method.
Background technology
Along with the fast development of human civilization, day by day increase the demand of the energy, and modern society's consumer demand also becomes more diverse, demand of industrial production product can fast transformation, and this all increases the demand to oil and natural gas.In order to effectively detect hydrocarbon zone, in the process of oilfield prospecting developing, well-log information digital processing and integrated interpretation become one of main task of vast well logging interpretation personnel, in the process, explanation personnel are by the response according to original logging trace, and the achievement curve after well-log information digital processing, simultaneously with reference to drilling well, geology, well logging, other data such as rock core, carry out formation evaluation, be intended to determine reservoir quality (oil reservoir, gas-bearing formation, water layer or dried layer etc.), and acquisition reservoir lithology, physical property (factor of porosity, permeability) and the reservoir parameter such as oiliness (oil saturation).
Past, traditional logging Reservoir Evaluation method was generally divided into three steps, 1. utilized computing machine carry out digital processing to log data and export drafting pattern due to the reason such as restriction of computer software technology; 2. in the mode of desktop drawing operation, reservoir division is carried out to logging trace, curve reads value and reservoir quality evaluation; 3. manual entry the above results forms well logging interpretation achievement form and printout.Complete the logging evaluation work of a bite well, generally need the cost 2 people time of more than 3 days, there is the work period long, poor in timeliness, precision is not high, the problems such as easily to make mistakes.Nowadays logging evaluation work can complete alternately from data process completely to integrated interpretation under computer graphical interface, but though interactive interpretation mode of operation under graphical interfaces obtains great improvement, to be still difficult to avoid in manual explanation problem, the problem includes: most of problem.
Summary of the invention
Present applicant proposes a kind of AUTOMATIC SORTING IN LAYERS BY WELL-LOGGING CURVES value and evaluation method, be intended to the algorithm for well-log information digital processing and integrated interpretation software platform provide a set of well logging automatically to explain, for interactive interpretation provides a set of initial explanation results automatically, to improve efficiency and the precision of Log Interpretation.
For solving the problems of the technologies described above, the embodiment of the present application provides a kind of AUTOMATIC SORTING IN LAYERS BY WELL-LOGGING CURVES value and evaluation method, comprises the following steps:
One, AUTOMATIC ZONING
Step S1: select some logging traces, arranges the weights of each curve, then is a comprehensive layering curve according to Curves compilation ruled synthesis;
Step S2: to the smoothing filtering of comprehensive layering curve moving weighted average method;
Step S3: layering is carried out to resultant curve application bathmometry;
Step S4: select a crucial curve, application cutoff method carries out layering, utilizes it to determine attribute at the bottom of the top at interface at the curve monotonicity of interface simultaneously;
Step S5: according to the rule of minimum thickness, resultant curve slope threshold values, the restriction of crucial curve cutoff, accepts or rejects step 3,4 layers produced;
Step S6: above-mentioned interface is built into non-reservoir units, reservoir units and little layer unit;
Two, automatic value
Step S7: all flex points and extreme point are obtained to every bar curve structure first order derivative and second derivative;
Step S8: carry out Detection curve form by statistics flex point, extreme point quantity and flex point place curve monotonicity in little layer unit, tracing pattern is divided into four kinds, i.e. monotonous curve, concave curve, convex curve and the curve of cyclical fluctuations;
Step S9: calculate minimum slope point value, mean value, minimal value, maximum value in little layer unit;
Step S10: determine value mode according to tracing pattern, namely monotonous curve gets minimum slope point value or mean value, and the curve of cyclical fluctuations is averaged, concave curve minimalization, and convex curve gets maximum value, carries out substratum assignment to every bar curve;
Three, automatic Evaluation
Step S11: successively according to the standard set up, for substratum is concluded judgement.
As the preferred embodiment of this programme, described step S1 chooses formation characteristics that curve is concerned about user as formation lithology, physical property, oiliness parameter sensitivity, the response of formation characteristics can both be reflected on comprehensive layering curve, the response of various characteristic is separate, or there is superposition amplification effect, the weights of all curves have positive and negative attribute.
Curves compilation formula is:
As the preferred embodiment of this programme, the rule that described step S5 is given comprises minimum thickness, resultant curve slope threshold values, the restriction of crucial curve cutoff.
As the preferred embodiment of this programme, the criterion of described step S11 comprises shale index higher limit, PLL, water saturation higher limit, and hydrocarbon types.
The one or more technical schemes provided in the embodiment of the present application, at least have following technique effect or advantage:
By computer simulation manual type, AUTOMATIC ZONING, automatically value and automatic Evaluation are carried out to logging trace, be also before carrying out follow-up interactive layering value and evaluating, a set of automatic explanation results is provided in advance; This method extremely meets requirement and the feature of logging trace manual zoning value, and controling parameters is flexible, and result is rationally reliable, can explain greatly reduce workload and raise the efficiency for follow-up interactive mode.The result of auto-layering data and evaluation is compared with final explanation results, and coincidence rate can reach more than 80%, and ideally, the achievement of auto-layering data and evaluation is even directly available.
Accompanying drawing explanation
In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.
Fig. 1 is the flow chart of steps of the embodiment of the present application.
Embodiment
In order to better understand technique scheme, below in conjunction with Figure of description and concrete embodiment, technique scheme is described in detail.
As shown in Figure 1, application the present invention carries out well logging interpretation appraisal, comprises the following steps:
One, AUTOMATIC ZONING
(1) the achievement curve-factor of porosity (POR) after well-log information digital processing, water saturation (SW), the weighting of shale index (VSH) curve are synthesized a resultant curve;
S= SW*MWsw + VSH*MWvsh + POR*MWpor
(2) by the smoothing filtering of resultant curve (5 moving weighted average methods can be divided into multiple rank according to weight coefficient change) of synthesis;
S=MW1*S[n-2]+MW2*S[n-1]+MW3*S[n]+MW4*S[n+1]+MW5*S[n+2]
(3) calculate first order derivative and the second derivative of resultant curve, and determine that the flex point degree of depth is as bed interface.
Logging trace is equal interval sampling data, and therefore the first order derivative of curve and second derivative are respectively:
S’[n]=S[n+1]-S[n]
S”[n]=S’[n+1]-S’[n]
Judge that the second derivative product of adjacent 2 is as time negative, can determine a flex point interface.
(4) the resultant curve slope at all flex point places of statistic procedure (3), obtain the reference value of average gradient as roughness, the slope cutoff that the rate of curve at comparison flex point place and roughness calculate, determines the choice (reservoir parameter reaches certain intensity of variation and just carries out layering) at flex point interface.
(5) using shale index curve as crucial curve, cutoff layering (the intersection point degree of depth of shale index curve and cutoff is cutoff interface) is carried out.
(6) all flex point interfaces (subdivision of reservoir not being carried out to mud stone) that all shale indexs are greater than shale index cutoff are given up.
(7) the flex point interface not meeting the thick requirement of smallest tier is given up.
(8) the cutoff interface not meeting the thick requirement of smallest tier is given up again.
(9) determine attribute at the bottom of the top at interface according to shale index curve at the monotonicity (being judged by first order derivative) of interface, subtracting curve is reservoir top, and increasing curve is at the bottom of reservoir.
(10) according to attribute at the bottom of the top of bed interface, build reservoir, at top and the end, form reservoir, between the end, top and flex point interface, form substratum, give up (well logging does not do interpretation and evaluation to non-reservoir) as non-reservoir between the end and top.
Two, automatic value
(11) scan all curves in substratum, by first order derivative and second derivative method, extreme value to be counted and flex point number counts, with this judgment curves form, concrete be divided into convex curve, concave curve, the curve of cyclical fluctuations and monotonous curve 4 type.
(12) maximum value minimal value, mean value, the minimum slope curve values of all curves in statistics layer.
(13) according to tracing pattern, to all curves assignment respectively, be specially convex curve and give maximum value, concave curve gives minimal value, the curve of cyclical fluctuations gives mean value, and monotonous curve gives minimum slope curve values (or give mean value also can).
Three, automatic Evaluation
(14) according to geological layering, the net thickness lower limit of the shale index of subsection setup, factor of porosity, water saturation or upper limit standard, and hydrocarbonaceous type;
(15) according to above-mentioned standard, successively well logging interpretation conclusion is made to substratum and differentiate.
The above, it is only preferred embodiment of the present invention, not any pro forma restriction is done to the present invention, although the present invention with preferred embodiment demonstration as above, but and be not used to limit the present invention, any those skilled in the art, do not departing within the scope of technical solution of the present invention, make a little change when the technology contents of above-mentioned announcement can be utilized or be modified to the Equivalent embodiments of equivalent variations, in every case be the content not departing from technical solution of the present invention, according to any simple modification that technical spirit of the present invention is done above embodiment, equivalent variations and modification, all still belong in the scope of technical solution of the present invention.
Claims (3)
1. AUTOMATIC SORTING IN LAYERS BY WELL-LOGGING CURVES value and an evaluation method, is characterized in that, comprises the following steps:
(1) AUTOMATIC ZONING:
Step S1: select some logging traces, arranges the weights of each curve, then is a comprehensive layering curve according to Curves compilation ruled synthesis;
Step S2: to the smoothing filtering of comprehensive layering curve moving weighted average method;
Step S3: layering is carried out to resultant curve application bathmometry;
Step S4: select a crucial curve, application cutoff method carries out layering, utilizes it to determine attribute at the bottom of the top at interface at the curve monotonicity of interface simultaneously;
Step S5: according to the rule of minimum thickness, resultant curve slope threshold values, the restriction of crucial curve cutoff, accepts or rejects step 3,4 layers produced;
Step S6: above-mentioned interface is built into non-reservoir units, reservoir units and little layer unit;
(2) automatic value:
Step S7: all flex points and extreme point are obtained to every bar curve structure first order derivative and second derivative;
Step S8: carry out Detection curve form by statistics flex point, extreme point quantity and flex point place curve monotonicity in little layer unit, tracing pattern is divided into four kinds, i.e. monotonous curve, concave curve, convex curve and the curve of cyclical fluctuations;
Step S9: calculate minimum slope point value, mean value, minimal value, maximum value in little layer unit;
Step S10: determine value mode according to tracing pattern, namely monotonous curve gets minimum slope point value or mean value, and the curve of cyclical fluctuations is averaged, concave curve minimalization, and convex curve gets maximum value, carries out substratum assignment to every bar curve;
(3) automatic Evaluation:
Step S11: successively according to the standard set up, for substratum is concluded judgement.
2. a kind of AUTOMATIC SORTING IN LAYERS BY WELL-LOGGING CURVES value according to claim 1 and evaluation method, it is characterized in that: the formation characteristics that the curve that described step S1 chooses is concerned about user is as formation lithology, physical property, oiliness parameter sensitivity, the response of formation characteristics can both be reflected on comprehensive layering curve, the response of various characteristic is separate, or has superposition amplification effect; The weights of all curves have positive and negative attribute.
3. a kind of AUTOMATIC SORTING IN LAYERS BY WELL-LOGGING CURVES value according to claim 1 and evaluation method, is characterized in that: the criterion of described step S11 comprises shale index 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 true CN104793263A (en) | 2015-07-22 |
| CN104793263B 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) |
Cited By (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 |
| CN109182463A (en) * | 2018-09-21 | 2019-01-11 | 博奥生物集团有限公司 | A kind of fluorescent amplification curve inflection point determines method and device |
| CN111173505A (en) * | 2018-10-23 | 2020-05-19 | 中国石油天然气股份有限公司 | Method and apparatus for determining a reservoir lower bound |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102116871A (en) * | 2011-01-24 | 2011-07-06 | 中国石油天然气股份有限公司 | Method for improving capability of seismic inverted data prediction storage layer |
| CN102454399A (en) * | 2010-10-26 | 2012-05-16 | 中国石油化工股份有限公司 | Method for correcting time difference signals of well logging acoustic waves |
| CN102830441A (en) * | 2012-05-16 | 2012-12-19 | 陶冶 | Comprehensive interpretation and evaluation method of hydrocarbon reservoir including data mining |
| CN104047597A (en) * | 2014-01-26 | 2014-09-17 | 中国石油大学(华东) | Fat gas mud shale stratum well log standardizing method |
| CN104295293A (en) * | 2014-10-23 | 2015-01-21 | 中国石油天然气股份有限公司 | Method for acquiring logging density curve |
-
2015
- 2015-05-05 CN CN201510222975.3A patent/CN104793263B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102454399A (en) * | 2010-10-26 | 2012-05-16 | 中国石油化工股份有限公司 | Method for correcting time difference signals of well logging acoustic waves |
| CN102116871A (en) * | 2011-01-24 | 2011-07-06 | 中国石油天然气股份有限公司 | 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 |
| CN104047597A (en) * | 2014-01-26 | 2014-09-17 | 中国石油大学(华东) | Fat gas mud shale stratum well log standardizing method |
| CN104295293A (en) * | 2014-10-23 | 2015-01-21 | 中国石油天然气股份有限公司 | Method for acquiring logging density curve |
Cited By (5)
| 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 |
| CN109182463A (en) * | 2018-09-21 | 2019-01-11 | 博奥生物集团有限公司 | A kind of fluorescent amplification curve inflection point determines method and device |
| CN109182463B (en) * | 2018-09-21 | 2022-04-29 | 博奥生物集团有限公司 | Fluorescence amplification curve inflection point determination method and device |
| CN111173505A (en) * | 2018-10-23 | 2020-05-19 | 中国石油天然气股份有限公司 | Method and apparatus for determining a reservoir lower bound |
| CN111173505B (en) * | 2018-10-23 | 2023-08-22 | 中国石油天然气股份有限公司 | Method and apparatus for determining a reservoir lower limit |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104793263B (en) | 2018-03-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104950334B (en) | A kind of method and device of predicting reservoir distribution | |
| KR101148835B1 (en) | Prediction system and method for subsurface lithology in oil sands reservoir using statistical analysis of well logging data | |
| CN109884707B (en) | Static correction method for near-surface layered time-depth curve | |
| CN103513278B (en) | Seismic wave groups thickness is utilized to carry out the method for reservoir prediction | |
| CN103775075A (en) | Method for identifying lithology of whole well section | |
| CN103993871A (en) | Method and device for processing well logging information of thin interbed stratums in standardization mode | |
| Euzen et al. | Well log cluster analysis: an innovative tool for unconventional exploration | |
| CN105549084A (en) | Method and system for building three-dimensional high-precision velocity model | |
| CN103376468A (en) | Reservoir parameter quantitative characterization method based on neural network function approximation algorithm | |
| CN106226841A (en) | A kind of fluvial facies three-dimensional sedimentary facies model Decided modelling method | |
| CN108072915A (en) | A kind of method for identifying carbonate rock particle beach phase | |
| CN111257962A (en) | Method for positioning and predicting ore body by using three-dimensional geological modeling technology | |
| CN104793263A (en) | Well logging curve automatic tiered value extraction and evaluation method | |
| CN111487691B (en) | Method for quantitatively identifying lithology and sand body structure of tight sandstone reservoir | |
| CN114609675A (en) | Quantitative recovery method for carbonate rock stratum sedimentary micro-landform based on high-frequency cycle | |
| CN110988997A (en) | Hydrocarbon source rock three-dimensional space distribution quantitative prediction technology based on machine learning | |
| CN103376463A (en) | Inversion modeling method based on fault control | |
| CN103048689A (en) | Seismic multi-attribute HLSO (Hue, Lightness, Saturation, Opacity) fusion method | |
| Gupta et al. | Seismic attribute driven integrated characterization of the Woodford Shale in west-central Oklahoma | |
| CN105676288B (en) | A kind of sand body stacked relation method of discrimination based on reflectance factor analysis | |
| US20140345946A1 (en) | Analysis of Geological Objects | |
| CN104360386B (en) | A kind of panel method for granitic formation Division and contrast | |
| CN109991663B (en) | Work area seismic velocity body correction method and device | |
| Iltaf et al. | Facies and petrophysical modeling of Triassic Chang 6 tight sandstone reservoir, Heshui oil field, Ordos basin, China | |
| Singleton | Geophysical data processing, rock property inversion, and geomechanical model building in a Midland Basin development project, Midland/Ector counties, Texas |
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 |