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 PDF

Info

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
Application number
CN201510222975.3A
Other languages
Chinese (zh)
Other versions
CN104793263A (en
Inventor
谢玉洪
林德明
钟华明
曾少军
何胜林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Offshore Oil Corp CNOOC
CNOOC China Ltd Zhanjiang Branch
Original Assignee
China National Offshore Oil Corp CNOOC
CNOOC China Ltd Zhanjiang Branch
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China National Offshore Oil Corp CNOOC, CNOOC China Ltd Zhanjiang Branch filed Critical China National Offshore Oil Corp CNOOC
Priority to CN201510222975.3A priority Critical patent/CN104793263B/en
Publication of CN104793263A publication Critical patent/CN104793263A/en
Application granted granted Critical
Publication of CN104793263B publication Critical patent/CN104793263B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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

A kind of AUTOMATIC SORTING IN LAYERS BY WELL-LOGGING CURVES value and evaluation method
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.
CN201510222975.3A 2015-05-05 2015-05-05 A kind of AUTOMATIC SORTING IN LAYERS BY WELL-LOGGING CURVES value and evaluation method Active CN104793263B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

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