CN105370273B - A kind of method of carbon/oxygen log environmental impact factor correction - Google Patents

A kind of method of carbon/oxygen log environmental impact factor correction Download PDF

Info

Publication number
CN105370273B
CN105370273B CN201510905675.5A CN201510905675A CN105370273B CN 105370273 B CN105370273 B CN 105370273B CN 201510905675 A CN201510905675 A CN 201510905675A CN 105370273 B CN105370273 B CN 105370273B
Authority
CN
China
Prior art keywords
values
gravel
annular fluid
layer
carbon
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
CN201510905675.5A
Other languages
Chinese (zh)
Other versions
CN105370273A (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 CN201510905675.5A priority Critical patent/CN105370273B/en
Publication of CN105370273A publication Critical patent/CN105370273A/en
Application granted granted Critical
Publication of CN105370273B publication Critical patent/CN105370273B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a kind of method of carbon/oxygen log environmental impact factor correction, the theory relation plate gone out using Monte Carlo simulation between influence factor and measured value, based on theoretical plate, the C/O values and the C/O values of uninfluenced layer position of affected layer position are chosen by reference lamina method, impacted C/O values and unaffected C/O values number are subjected to linear function fit using trend surface analysis approximating method, obtain the fit correlation formula of impacted C/O values and unaffected C/O values.The inventive method can fill the two influence factors to annular fluid and gravel and be corrected, the reliability of the log response value of impacted interval can be significantly improved using this method, and then the precision of C/O Log Interpretation evaluation remaining oil saturation is improved, there is stronger versatility.

Description

A kind of method of carbon/oxygen log environmental impact factor correction
Technical field
The present invention relates to a kind of logging evaluation technology for improving carbon-to-oxygen ratio and calculating remaining oil saturation Explanation Accuracy, it is specifically The bearing calibration that annular fluid and gravel filling influence during a kind of carbon/oxygen log, belongs to oil-gas exploration technical field.
Background technology
Carbon-to-oxygen ratio spectrum logging is also in carbon/oxygen log, is a kind of new pulsed neutron log method.It is main in oil Carbon elements and protium, mainly contain oxygen element and protium in water, carbon oxygen phase in stratum can be obtained by carbon/oxygen log To content ratio, the oil-gas Layer omitted can be found in the well for descended sleeve pipe, oil reservoir is determined in the oil well to have recovered the oil Residual saturation etc..
, it is necessary to which logging technique personnel enter to the C/O curves of measurement during carbon/oxygen log calculates remaining oil saturation Row influences correction, and so as to obtain the curve that can reflect stratum real information, the key of work is the school for determining influence factor Correction method.Main influence factor includes sleeve pipe, cement sheath, annular fluid, shale content and gravel during carbon/oxygen log Filling.Existing bearing calibration is primarily directed to these three factors of sleeve pipe, cement sheath and shale content, and these methods are main It is to be corrected according to experiment plate.And annular fluid and gravel fill the two influence factors and corrected there is presently no corresponding Method.
The content of the invention
In order to overcome above mentioned problem existing for prior art, the present invention provides a kind of carbon/oxygen log environmental impact factor school Positive method, single influence factor is carried out by the method for sharp Monte Carlo simulation and carries out theoretical modeling, the correction to theorize Plate, the relation between each influence factor and measured value is summarized, actual calibration model is then established according to reference lamina method.
Technical scheme is as follows:A kind of method of carbon/oxygen log environmental impact factor correction, including following step Suddenly:
Step 101, log obtains and layer choosing to be corrected takes:Using conventional carbon oxygen logging method, sent out using neutron source Project fast neutron and inelastic scattering occurs with stratum each element, stratum each element atomic nucleus discharges inelastic scattering gamma, profit Spectrum unscrambling is carried out with the feature spectral peak of each element, extracts log;Played back according to log, find and annular fluid influence be present The annular fluid factor correction reference lamina and gravel filling factor correction reference lamina influenceed with gravel filling, and by the two layers The C/O logs value of position is as a variable in data analysis;
Step 102, Monte Carlo theoretical modeling:The method of sharp Monte Carlo simulation carries out theoretical mould to single influence factor Intend, the correcting plate to theorize, analyze the relation between each influence factor and measured value;
Step 103, actual calibration model is established:Data fitting is carried out using reference lamina method, draws empirical equation COR=A × COR+B, COR is carbon-to-oxygen ratio in formula, and A, B are correction coefficient, is determined generally according to the actual geological condition in different oil fields;
Step 104:Influence factor corrects:On the basis of step 103, the empirical equation obtained using fitting is to carbon-to-oxygen ratio Carry out annular fluid and gravel filling is corrected, quantitative interpretation is carried out by the carbon-to-oxygen ratio value after correction.
In step 102, the influence factor includes:Sleeve pipe, cement sheath, annular fluid, shale content and gravel filling.
Further, it is the reference lamina based on step 102 selection to carry out data fitting using reference lamina method described in step 103 Carry out, wherein:
(1)The choosing method of annular fluid factor correction reference lamina:Select the water layer not influenceed by annular space, and the layer Shale content, porosity and whether have gravel filling etc. parameter try one's best it is roughly the same with the holding of layer to be corrected, take correspondingly The C/O values of interval are as annular space calibration mark layer characteristic;
(3)Gravel fills the choosing method of factor correction reference lamina:Select the minimum water layer of shale content, and the stratum Shale content, porosity and whether by annular space to be influenceed these three parameters roughly the same with layer holding to be corrected, take respective layer The C/O values of section fill calibration mark layer characteristic as gravel.
Further, the establishment step of calibration model described in step 103 specifically includes:Annular fluid and gravel are carried out to individual well Stone fills two influence factors and analyzed, according to the form of curve, as a certain interval C/O and Si/Ca envelope of curves areas Stepped increase is showed compared with the envelope size of adjacent layer position, then it is assumed that the layer has annular fluid influence.Select first Taking is influenceed the C/O value X of interval by annular fluid and gravel filling, then chooses the interval not influenceed by annular fluid and gravel filling C/O values Y.The data of selection are handled, are that ordinate establishes scatter diagram using X as abscissa, Y, by Y data and X parameter Data relationship is established, according to analysis result, it be Y=aX+b to establish correction fitting formula, wherein, Y is the C/O values for indicating interval, nothing Dimension;X be impacted interval C/O values, dimensionless;A is the slope of curve, and b is intercept.
Having the beneficial effect that relative to prior art of the invention:The present invention is to combine the Annular cutting to carbon/oxygen log Body and gravel fill the two influence factors, by sharp Monte Carlo theoretical modeling, the correcting plate that theorizes, according to reference lamina Method establishes calibration model;The inventive method can be accurately handled the well-log information influenceed by logging environment, by production It was verified that the true formation information in C/O curves can be obtained using the inventive method, and then improve stratum remaining oil saturation The precision of well log interpretation evaluation porosity is spent, there is stronger versatility;The present invention provides for RPM well logging influence factor corrections Method;The inventive method can integrally lift remaining oil saturation well log interpretation ability, pull operation, instruct the economy of offshore oilfield Efficient Development.
Brief description of the drawings
Fig. 1 is that annular fluid interval identifies schematic diagram.
Fig. 2 is influence of the annular fluid of Monte Carlo simulation to C/O values.
Fig. 3 is that the gravel of Monte Carlo simulation fills the influence to C/O values.
Fig. 4 is the annular fluid calibration model obtained in present example using data fit analysis.
Fig. 5 is that the gravel obtained in present example using data fit analysis fills calibration model.
Fig. 6 is that the annular fluid in the present invention corrects example, and wherein ALAVGCO is the C/O curves after correction, with measured value Compared to reduction, it is consistent with plate understanding.
Fig. 7 is the gravel filling correction example in the present invention, and wherein ALAVGCO is the C/O curves after correction, with measured value Compared to increase, it is consistent with plate understanding.
Embodiment
Carbon/oxygen log environmental impact factor bearing calibration of the present invention is made into one below in conjunction with the drawings and specific embodiments The detailed description of step, specifically comprises the following steps:
Step 101, operated using carbon oxygen logging technology conventional method and obtain log data, according to log data, picked out C/O and Si/Ca logs, according to the form of curve, the interval of annular fluid and gravel filling is selected, takes the two layers The C/O curve values of section;Annular fluid interval identification schematic diagram such as Fig. 1.
Step 102, Monte Carlo theoretical modeling:The actual geology and stratum water analysis data in combining target region, establish Model.Concrete model is as follows:Wellbore section is made up of wellbore fluid, sleeve pipe, cement sheath, and well is full of fresh water, sleeve pipe internal-and external diameter Respectively 224.5mm and 244.5mm, oil pipe internal-and external diameter are respectively 62.0mm and 73.0mm;The thickness of cement sheath is 33.3cm, into It is divided into CaSiO3, density 1.95g/cm3;The radial direction radius on stratum is 70cm, a height of 102cm, and fills different fluid respectively And lithology material;Formation water salinity is 35000ppm;The distance of long short space is respectively 159.766mm and 235.966mm.Cover Influence of the annular fluid of special Monte Carlo Simulation of Ions Inside to C/O values is as shown in Figure 2;Shadow of the gravel filling of Monte Carlo simulation to C/O values Ring as shown in Figure 3.
Step 103, the layer to be corrected influenceed by annular space and the reference lamina not influenceed by annular space are found.
Step 104, the C/O data picked out from step 102 are fitted using data fit analysis, be fitted The annular fluid calibration model arrived is as shown in Figure 4;It is as shown in Figure 5 to be fitted obtained gravel filling calibration model;Using fitting The model formation come is corrected to annular fluid and gravel filling respectively.
Analyzed specifically, carrying out two influence factors of annular fluid and gravel filling to individual well, according to the shape of curve State, when a certain interval C/O and Si/Ca envelopes of curves area show stepped increasing compared with the envelope size of adjacent layer position Greatly, then it is assumed that the layer has annular fluid influence(Such as the red block interval in Fig. 1).Choose and filled by annular fluid and gravel first Filling out influences the C/O value X of interval, then chooses the C/O values Y for the interval not influenceed by annular fluid and gravel filling.To the number of selection Be that ordinate establishes scatter diagram using X as abscissa, Y according to being handled, Y data and X parameter established into data relationship, according to point Result is analysed, annular fluid correction fitting formula is Y=- 0.069X+1.978, gravel filling correction fitting formula Y=- 2.688X+ 7.428, wherein, Y be mark interval C/O values, dimensionless;X be impacted interval C/O values, dimensionless.
The annular fluid calibration model obtained using step 104 fitting is corrected to annular fluid, as shown in fig. 6, its Middle ALAVGCO is the C/O curves after correction, is reduced compared with measured value, is consistent with plate understanding, illustrates established straightening die Type has good precision.
The gravel filling calibration model obtained using step 104 fitting is corrected to gravel filling, as shown in fig. 7, its Middle ALAVGCO is the C/O curves after correction, is increased compared with measured value, is consistent with plate understanding, illustrates established straightening die Type has good precision.
Above-described embodiment, the purpose of the present invention, technical scheme and beneficial effect are carried out further Describe in detail, should be understood that the embodiment that the foregoing is only the present invention, be not intended to limit the present invention Protection domain, within the spirit and principles of the invention, any modification, equivalent substitution and improvements done etc., all should include Within protection scope of the present invention.The content that this specification is not described in detail belongs to known to professional and technical personnel in the field Prior art.

Claims (4)

  1. A kind of 1. method of carbon/oxygen log environmental impact factor correction, it is characterised in that:Comprise the following steps:
    Step 101, log obtains and layer choosing to be corrected takes:Using conventional carbon oxygen logging method, launched using neutron source With stratum each element inelastic scattering occurs for fast neutron, and stratum each element atomic nucleus discharges inelastic scattering gamma, using each The feature spectral peak of element carries out spectrum unscrambling, extracts log;Played back according to log, find and annular fluid influence and gravel be present The annular fluid factor correction reference lamina and gravel filling factor correction reference lamina that stone filling influences, and by the two layer of position C/O logs value is as a variable in data analysis;
    Step 102, Monte Carlo theoretical modeling:Theoretical mould is carried out to single influence factor using the method for Monte Carlo simulation Intend, the correcting plate to theorize, analyze the relation between each influence factor and measured value;
    Step 103, actual calibration model is established:Data fitting is carried out using reference lamina method, draws empirical equation COR=A × COR + B, COR is carbon-to-oxygen ratio in formula, and A, B are correction coefficient, is determined according to the actual geological condition in different oil fields;
    Step 104:Influence factor corrects:On the basis of step 103, carbon-to-oxygen ratio is carried out using the empirical equation that fitting obtains Annular fluid and gravel filling are corrected, and quantitative interpretation is carried out by the carbon-to-oxygen ratio value after correction.
  2. 2. method according to claim 1, it is characterised in that:In step 102, the influence factor is:Sleeve pipe, cement sheath, Annular fluid, shale content and gravel filling.
  3. 3. method according to claim 1, it is characterised in that:Carrying out data fitting using reference lamina method described in step 103 is Carried out based on the reference lamina that step 102 is chosen, wherein:
    (1)The choosing method of annular fluid factor correction reference lamina:Select the water layer not influenceed by annular space, and the mud of this layer Matter content, porosity and whether there is gravel to fill these three parameters are roughly the same with layer holding to be corrected, take corresponding interval C/O values are as annular space calibration mark layer characteristic;
    (2)Gravel fills the choosing method of factor correction reference lamina:Select the minimum water layer of shale content, and the mud on the stratum Matter content, porosity and whether influenceed that these three parameters are roughly the same with layer holding to be corrected by annular space, take corresponding interval C/O values fill calibration mark layer characteristic as gravel.
  4. 4. method according to claim 1, it is characterised in that:The establishment step of calibration model described in step 103 specifically includes: Annular fluid is carried out to individual well and gravel fills two influence factors and analyzed, according to the form of curve, as a certain interval C/O Stepped increase is showed compared with the envelope size of adjacent layer position, then it is assumed that the layer is deposited with Si/Ca envelopes of curves area Influenceed in annular fluid;Choosing first is influenceed the C/O value X of interval by annular fluid and gravel filling, then is chosen not by Annular cutting The C/O values Y for the interval that body and gravel filling influence;The data of selection are handled, are that ordinate is established using X as abscissa, Y Scatter diagram, Y data and X parameter are established into data relationship, according to analysis result, it is Y=aX+b to establish correction fitting formula, wherein, Y be mark interval C/O values, dimensionless;X be impacted interval C/O values, dimensionless;A is the slope of curve, and b is intercept.
CN201510905675.5A 2015-12-10 2015-12-10 A kind of method of carbon/oxygen log environmental impact factor correction Active CN105370273B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510905675.5A CN105370273B (en) 2015-12-10 2015-12-10 A kind of method of carbon/oxygen log environmental impact factor correction

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510905675.5A CN105370273B (en) 2015-12-10 2015-12-10 A kind of method of carbon/oxygen log environmental impact factor correction

Publications (2)

Publication Number Publication Date
CN105370273A CN105370273A (en) 2016-03-02
CN105370273B true CN105370273B (en) 2018-01-26

Family

ID=55372831

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510905675.5A Active CN105370273B (en) 2015-12-10 2015-12-10 A kind of method of carbon/oxygen log environmental impact factor correction

Country Status (1)

Country Link
CN (1) CN105370273B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106958443B (en) * 2017-05-02 2020-10-09 长江大学 Method and system for correcting carbon-oxygen ratio logging drilling fluid invasion influence factors
CN107152277B (en) * 2017-06-07 2020-11-10 长江大学 Method and system for calculating residual oil saturation degree through carbon-oxygen ratio logging
CN108661632B (en) * 2018-05-07 2022-02-11 何巍巍 High-precision carbon-oxygen ratio logging method
CN109779618B (en) * 2019-01-10 2020-07-07 中国石油大学(北京) Method and device for determining formation residual oil saturation
CN109779611A (en) * 2019-02-21 2019-05-21 中陕核工业集团地质调查院有限公司 Application of the petroleum hole cementing concrete ring property parameters in the amendment of old well gamma repetition measurement coefficient
CN109915120B (en) * 2019-04-16 2022-06-03 北京恒泰万博石油技术股份有限公司 Correction method of resistivity logging while drilling system based on environmental factors
CN112031742B (en) * 2020-09-03 2023-07-04 中海油田服务股份有限公司 Carbon-oxygen ratio energy spectrum logging saturation interpretation method based on database

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4136278A (en) * 1977-07-14 1979-01-23 Dresser Industries, Inc. Method and apparatus for pulsed neutron spectral analysis using spectral stripping
CN101906963A (en) * 2010-07-23 2010-12-08 中国石油化工集团公司 Method for determining saturation by using C/O and formation macro capture cross-section intersecting technique
US8011238B2 (en) * 2008-10-09 2011-09-06 Chevron U.S.A. Inc. Method for correcting the measured concentrations of gas components in drilling mud
WO2014080107A1 (en) * 2012-11-20 2014-05-30 Total Sa Method for assessing a quantity of hydrocarbons in a deposit
CN104929627A (en) * 2015-06-25 2015-09-23 中国海洋石油总公司 Method for calculating wellbore water-holding capacity and formation water saturation in RPM well logging

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4136278A (en) * 1977-07-14 1979-01-23 Dresser Industries, Inc. Method and apparatus for pulsed neutron spectral analysis using spectral stripping
US8011238B2 (en) * 2008-10-09 2011-09-06 Chevron U.S.A. Inc. Method for correcting the measured concentrations of gas components in drilling mud
CN101906963A (en) * 2010-07-23 2010-12-08 中国石油化工集团公司 Method for determining saturation by using C/O and formation macro capture cross-section intersecting technique
WO2014080107A1 (en) * 2012-11-20 2014-05-30 Total Sa Method for assessing a quantity of hydrocarbons in a deposit
CN104929627A (en) * 2015-06-25 2015-09-23 中国海洋石油总公司 Method for calculating wellbore water-holding capacity and formation water saturation in RPM well logging

Also Published As

Publication number Publication date
CN105370273A (en) 2016-03-02

Similar Documents

Publication Publication Date Title
CN105370273B (en) A kind of method of carbon/oxygen log environmental impact factor correction
KR101782294B1 (en) Modeling method for gas production of cbm reservoir rocks
CN104047597B (en) Oily mud shale stratum log standardized method
Sullivan et al. Architectural analysis of deep-water outcrops: Implications for exploration and development of the Diana sub-basin, western Gulf of Mexico
França et al. Secondary porosity formed by deep meteoric leaching: Botucatu eolianite, southern South America
Budai et al. Fracture‐fill calcite as a record of microbial methanogenesis and fluid migration: a case study from the Devonian Antrim Shale, Michigan Basin
CN103670388A (en) Method for evaluating organic carbon content of shale
CN104806232B (en) A kind of method for determining porosity lower limit of fracture
Miller et al. Sequence boundaries are impedance contrasts: Core-seismic-log integration of Oligocene–Miocene sequences, New Jersey shallow shelf
CN105317431A (en) Method for explaining and evaluating logging parameters of horizontal well
CN101906963A (en) Method for determining saturation by using C/O and formation macro capture cross-section intersecting technique
CN103867198B (en) Method for distinguishing formation density of carbonatite natural gas reservoir
CN110056346B (en) Oil reservoir three-dimensional original water saturation simulation method based on trend change function
CN103744109B (en) Method for identifying a weathering crust structure of clastic rock in a area of covering no well
KR101708599B1 (en) Modeling method for gas content of cbm reservoir rocks
CN103993871A (en) Method and device for processing well logging information of thin interbed stratums in standardization mode
CN108073745A (en) Glutenite lithology physical property Logging Identification Method based on three air void curves
CN107288607A (en) A kind of method of the nearly well fracturing fracture of utilization Gd neutrons spike yield imaging in evaluation
CN112764122A (en) Delineation method of ancient river channel type sandstone uranium reservoir
KR101731221B1 (en) Method for estimating facies of cbm reservoir rocks
CN104251135B (en) Highly-deviated well space in-place method
CN109343099A (en) A method of with GR and SP joint account stratum containing radioactive substance shale content
Jaeger et al. Site U1417
CN109444980B (en) Efficient deep geothermal resource exploration method
CN105865999A (en) Analysis method for diagenetic stages of sedimentary rock

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No.

Applicant after: China Offshore Oil Group Co., Ltd.

Applicant after: CNOOC (China) Limited Zhanjiang Branch

Address before: China CNOOC Building No. 25 Beijing 100000 Chaoyang District Chaoyangmen North Street

Applicant before: China National Offshore Oil Corporation

Applicant before: CNOOC (China) Limited Zhanjiang Branch

GR01 Patent grant
GR01 Patent grant
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Li Yuelin

Inventor after: Guo Haimin

Inventor after: He Shenglin

Inventor after: Zheng Yongjian

Inventor after: Wang Lijuan

Inventor before: Li Yuelin

Inventor before: Guo Haimin

Inventor before: He Shenglin

Inventor before: Zheng Yongjian

Inventor before: Wang Lijuan