CN107329181B - A method of seeking muddy ore formation water resistivity and litho-electric parameters - Google Patents

A method of seeking muddy ore formation water resistivity and litho-electric parameters Download PDF

Info

Publication number
CN107329181B
CN107329181B CN201710705083.8A CN201710705083A CN107329181B CN 107329181 B CN107329181 B CN 107329181B CN 201710705083 A CN201710705083 A CN 201710705083A CN 107329181 B CN107329181 B CN 107329181B
Authority
CN
China
Prior art keywords
water
resistivity
formation water
litho
electric parameters
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
CN201710705083.8A
Other languages
Chinese (zh)
Other versions
CN107329181A (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 CN201710705083.8A priority Critical patent/CN107329181B/en
Publication of CN107329181A publication Critical patent/CN107329181A/en
Application granted granted Critical
Publication of CN107329181B publication Critical patent/CN107329181B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)

Abstract

The present invention provides a kind of methods for seeking muddy ore formation water resistivity and litho-electric parameters, using this method can by litho-electric parameters a, m and formation water resistivity Rw by the method for iteration fast and accurately solve come, through the invention, it only needs to provide the water layer that a water saturation is 100%, so that it may accurately acquire litho-electric parameters and formation water resistivity by the method for iteration.Both solved the problems, such as that conventional Pickett plate was not used to the calculating of argillaceous sandstone formation water resistivity, it helps prospect pit, development well water saturation quickly calculate, and improve water saturation computational accuracy.

Description

A method of seeking muddy ore formation water resistivity and litho-electric parameters
Technical field
The present invention relates to geophysical well logging technology fields more particularly to one kind to seek muddy ore stratum water resistance The method of rate and litho-electric parameters.
Background technique
During oil gas field logging evaluation, litho-electric parameters and the accurate of formation water resistivity seek being water saturation The key point whether key of calculating is accurate.Litho-electric parameters acquisition methods mainly pass through coring rock-electric test and obtain at present, Formation water resistivity acquisition methods are mainly determined by water analysis of data.Pickett proposes classical Pickett within 1966 Plate need to only find the pure water layer that a water saturation is 100%, can be in resistivity-porosity log-log coordinate axis On calculate formation water resistivity and litho-electric parameters.However, technical staff's discovery passes through experimental data in actual production practice The method for determining litho-electric parameters and formation water resistivity takes a long time, and does not meet the demand of live quick look evaluation；And Although Pickett plate can quickly calculate formation water resistivity, it is not particularly suited for muddy ore.
Summary of the invention
The present invention provides one kind can quick and precisely seek muddy ore formation water resistivity and litho-electric parameters Method, using this method can by litho-electric parameters a, m and formation water resistivity Rw by the method for iteration fast and accurately It solves and, be conducive to the fine evaluation of prospect pit, development well.
In order to solve the above technical problems, the present invention provides one kind to seek muddy ore formation water resistivity and rock The method of electrical parameter, comprising the following steps:
S1: finding the water layer that a water saturation is 100% by log (lithology, electrical property, physical property), and reading should Layer resistivity Rt；
S2: calculating the shale content Vsh of the water layer by lithology curve (GR, SP curve), by physical property curve (DEN, CNL, AC curve) calculate water-yielding stratum porosity φ；
S3: given lithologic index and formation water resistivity product aRw and cementation factor m initial value are respectively as follows: aRw1, m1, Shale Correction item Ash, calculation formula are calculated with the parameter are as follows:
S4: being corrected resistivity using Shale Correction item Ash, calculates Rt/Ash and is intersected with porosity φ, Argillaceous sandstone Pickett plate is drawn out, water layer Rt/Ash and porosity φ are established into power function relationship Rt/Ash=(aRw)' φ-m', calculateIf Δ aRw< 1%, Δ m < 1%, then illustrate lithologic index It is true value with formation water resistivity product aRw and cementation factor m, otherwise by (aRw) ' brought into m' as step S3 initial value And step S3, step S4 are repeated until meeting the condition of convergence, solve lithologic index and formation water resistivity product aRw and glue Tie the true value of exponent m.
One or more technical solutions provided in the embodiments of the present application have at least the following technical effects or advantages:
Through the invention, it is only necessary to provide the water layer that a water saturation is 100%, so that it may by the method for iteration compared with Accurately to acquire litho-electric parameters and formation water resistivity.Both it had solved conventional Pickett plate and has been not used to argillaceous sandstone The problem of formation water resistivity calculates, it helps prospect pit, development well water saturation quickly calculate, and improve water saturation Computational accuracy.
Detailed description of the invention
Fig. 1 first time iteration effect picture；
Fig. 2 third time iteration effect picture；
Fig. 3 sixth iteration effect picture.
Specific embodiment
In order to better understand the above technical scheme, being carried out below in conjunction with specific embodiment to above-mentioned technical proposal Detailed description.
A kind of method for seeking muddy ore formation water resistivity and litho-electric parameters described in the present embodiment, including Following steps:
Step 1: the water layer that a water saturation is 100% being found by log (lithology, electrical property, physical property), is read Layer resistivity Rt out；
Step 2: calculating the shale content Vsh of the water layer by lithology curve (GR, SP curve), pass through physical property curve The porosity φ of (DEN, CNL, AC curve) calculating water-yielding stratum；
Step 3: given lithologic index is respectively as follows: 2,4 with formation water resistivity product aRw and cementation factor m initial value, uses The parameter calculates Shale Correction item Ash, calculation formula are as follows:
Step 4: resistivity being corrected using Shale Correction item Ash, Rt/Ash and porosity φ is calculated and is handed over Argillaceous sandstone Pickett plate is drawn out in meeting, and water layer Rt/Ash and porosity φ are established power function relationship Rt/Ash= (aRw)'φ-m', calculateIf Δ aRw< 1%, Δ m < 1%, then illustrate lithology Index and formation water resistivity product aRw and cementation factor m are true value, otherwise by (aRw) ' and m' as step 3 initial value Step 3, step 4 are brought and repeated into up to meeting the condition of convergence, iteratively solves out lithologic index and formation water resistivity by 6 times The true value of product aRw and cementation factor m is respectively 0.067,1.984.
Basic principle of the invention is to be incorporated Pickett plate thought based on Indonesia's formula, obtained stratum water power Resistance rate and litho-electric parameters relational expression, it may be assumed that
Have and only 1 group of lithologic index and formation water resistivity product aRw and cementation factor m can satisfy the formula, therefore Using the formula as iterative formula, initial value is inputted, passes through iteration, so that it may find out lithologic index and formation water resistivity product aRw And the true value of cementation factor m.
The above described is only a preferred embodiment of the present invention, be not intended to limit the present invention in any form, though So the present invention has been disclosed as a preferred embodiment, and however, it is not intended to limit the invention, any technology people for being familiar with this profession 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, but anything that does not depart from the technical scheme of the invention content, according to the technical essence of the invention Any simple modification, equivalent change and modification to the above embodiments, all of which are still within the scope of the technical scheme of the invention.

Claims (1)

1. a kind of method for seeking muddy ore formation water resistivity and litho-electric parameters, which is characterized in that including as follows Step:
S1: pass through log: lithology, electrical property, physical property find the water layer that a water saturation is 100%, read this layer of electricity Resistance rate Rt；
S2: passing through lithology curve: GR, SP curve, calculates the shale content Vsh of the water layer, passes through physical property curve: DEN, CNL, AC Curve calculates the porosity φ of water-yielding stratum；
S3: given lithologic index and formation water resistivity product aRw and cementation factor m initial value are respectively as follows: aRw1, m1, with this Parameter calculates Shale Correction item Ash, calculation formula are as follows:
S4: being corrected resistivity using Shale Correction item Ash, calculates Rt/Ash and is intersected with porosity φ, draws Water layer Rt/Ash and porosity φ are established power function relationship R by argillaceous sandstone Pickett plate outt/Ash=(aRw)'φ-m', It calculatesIf Δ aRw< 1%, Δ m < 1%, then illustrate lithologic index and stratum Water resistance rate product aRw and cementation factor m are true value, otherwise by (aRw) ' brought into as S3 initial value with m' and repeat step S3, S4 solve the true of lithologic index and formation water resistivity product aRw and cementation factor m until meet the condition of convergence Value.
CN201710705083.8A 2017-08-17 2017-08-17 A method of seeking muddy ore formation water resistivity and litho-electric parameters Active CN107329181B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710705083.8A CN107329181B (en) 2017-08-17 2017-08-17 A method of seeking muddy ore formation water resistivity and litho-electric parameters

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710705083.8A CN107329181B (en) 2017-08-17 2017-08-17 A method of seeking muddy ore formation water resistivity and litho-electric parameters

Publications (2)

Publication Number Publication Date
CN107329181A CN107329181A (en) 2017-11-07
CN107329181B true CN107329181B (en) 2019-03-01

Family

ID=60228107

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710705083.8A Active CN107329181B (en) 2017-08-17 2017-08-17 A method of seeking muddy ore formation water resistivity and litho-electric parameters

Country Status (1)

Country Link
CN (1) CN107329181B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111720114B (en) * 2019-03-22 2023-06-02 中海油能源发展股份有限公司 Well logging saturation calculation method for sandstone hydrocarbon reservoir

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101725344A (en) * 2008-10-29 2010-06-09 中国石油天然气集团公司 Method for determining litho-electric parameters
CN103225500A (en) * 2013-05-02 2013-07-31 中国石油大学(华东) Novel water flooding layer logging evaluation method applying three parameters self-consistent iterative algorithm
CN104181090A (en) * 2013-05-27 2014-12-03 中国石油化工股份有限公司 Hydrocarbon saturation evaluation method based on rock electrical structure coefficient

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101358943A (en) * 2007-07-31 2009-02-04 中国石油天然气集团公司 Method for correcting electrical experiment parameter m and n of undersaturated core
US20120143508A1 (en) * 2010-12-01 2012-06-07 Conocophillips Company Automatic estimation of source rock petrophysical properties
CN102156297B (en) * 2011-05-16 2012-10-10 中国石油大学(北京) Fluid substitution method based on sandstone reservoir post-stack seismic data

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101725344A (en) * 2008-10-29 2010-06-09 中国石油天然气集团公司 Method for determining litho-electric parameters
CN103225500A (en) * 2013-05-02 2013-07-31 中国石油大学(华东) Novel water flooding layer logging evaluation method applying three parameters self-consistent iterative algorithm
CN104181090A (en) * 2013-05-27 2014-12-03 中国石油化工股份有限公司 Hydrocarbon saturation evaluation method based on rock electrical structure coefficient

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
高温高压条件下不同气体组分储层岩电实验及应用;何胜林，等;《天然气地球科学》;20170430;第28卷(第4期);575-581

Also Published As

Publication number Publication date
CN107329181A (en) 2017-11-07

Similar Documents

Publication Publication Date Title
CN109838230B (en) Quantitative evaluation method for oil reservoir water flooded layer
CN103775057B (en) A kind of fine and close oil-gas reservoir Effective Reservoirs sentence knowledge method and device
CN104500055B (en) A kind of extra-low permeability oil reservoirs Water Flooding Layer water saturation computational methods
CN106951660A (en) A kind of marine clastics horizontal well reservoir log interpretation method and device
CN103437760B (en) A kind of method using array induction data Fast Evaluation oil-water-layer
CN105447762B (en) A kind of calculation method of the low-permeability oil deposit water logging information of fluid replacement
CN104375172B (en) Volcanic underlayer structural configuration correct imaging method
CN109100793A (en) The method that a kind of quantitative analysis crack factor influences reservoir
CN109298464B (en) Tight sandstone reservoir diagenetic facies logging identification method and device
Wimmers et al. Integration of sedimentology, petrophysics and rock typing as key to understanding a tight gas reservoir
CN109611073A (en) A kind of densification grease horizontal well mud stone can the preferred method and device in fractured layer position
CN110532507B (en) Method for improving drilling rate of compact oil reservoir of horizontal well
CN107329181B (en) A method of seeking muddy ore formation water resistivity and litho-electric parameters
CN112861890A (en) Reservoir evaluation model construction method and reservoir identification method
Li et al. Discrimination of effective source rocks and evaluation of the hydrocarbon resource potential in Marsel, Kazakhstan
Shchipanov et al. A new approach to deformable fractured reservoir characterization: case study of the Ekofisk field
Kelkar et al. Upgridding method for tight gas reservoirs
CN109343099B (en) Method for calculating mud content of stratum containing radioactive substances by using GR and SP in combination mode
CN109973068B (en) Method and device for identifying oil reservoir water injection induced cracks
CN109298462B (en) Method and system for estimating shale content of stratum
Yang et al. Prediction and research of fracture distribution characteristics of the Chang 8 reservoir in Zhenjing area, Ordos Basin
CN107808229A (en) Array induction, the sandstone reservoir mobile water saturation quantitative evaluation method of laterolog
WANG et al. Research on predicting S-wave velocity in tight reservoir: Sulige gas field for example
Deng et al. Estimation of porosity and permeability from conventional logs in tight sandstone reservoirs of north Ordos basin
ZHANG et al. A saturation evaluating method based on pore structure

Legal Events

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

Address after: 100010 Chaoyangmen North Street, Dongcheng District, Dongcheng District, Beijing

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

Applicant after: CNOOC (China) Limited Zhanjiang Branch

Address before: 100010 China oil tower, 25 Chaoyangmen North Street, Dongcheng District, Beijing

Applicant before: China National Offshore Oil Corporation

Applicant before: CNOOC (China) Limited Zhanjiang Branch

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant