CN110632665B - Sand body configuration quantitative analysis method - Google Patents

Sand body configuration quantitative analysis method Download PDF

Info

Publication number
CN110632665B
CN110632665B CN201910793343.0A CN201910793343A CN110632665B CN 110632665 B CN110632665 B CN 110632665B CN 201910793343 A CN201910793343 A CN 201910793343A CN 110632665 B CN110632665 B CN 110632665B
Authority
CN
China
Prior art keywords
sand body
configuration
relation
well
sand
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.)
Expired - Fee Related
Application number
CN201910793343.0A
Other languages
Chinese (zh)
Other versions
CN110632665A (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 University of Geosciences
Original Assignee
China University of Geosciences
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 University of Geosciences filed Critical China University of Geosciences
Priority to CN201910793343.0A priority Critical patent/CN110632665B/en
Publication of CN110632665A publication Critical patent/CN110632665A/en
Application granted granted Critical
Publication of CN110632665B publication Critical patent/CN110632665B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/306Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/61Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616Data from specific type of measurement
    • G01V2210/6169Data from specific type of measurement using well-logging

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention provides a sand body configuration quantitative analysis method, which comprises the following steps: identifying a study area gyrus interface; determining a sand body configuration interface according to the gyrating interface; analyzing the sand body configuration interface to obtain single-well sand body configuration elements and contact relations; combining the single-well sand body configuration elements and the contact relation to obtain a longitudinal sand body splicing relation; predicting the transverse sand body splicing relation according to the longitudinal sand body splicing relation; obtaining a sand body plane distribution diagram according to the longitudinal sand body splicing relation and the transverse sand body splicing relation; and (3) according to the sand body plane distribution diagram, utilizing Monte Carlo simulation to realize fine depiction of the sand body form and scale. The invention has the beneficial effects that: the problem of determining the sand body configuration parameters under the condition of less constraint can be solved.

Description

Sand body configuration quantitative analysis method
Technical Field
The invention relates to the technical field of mining, in particular to a sand body configuration quantitative analysis method.
Background
In the process of oil reservoir development, the research on the internal configuration of a reservoir sand body, particularly a single sand body, is of great significance to the research on the distribution of residual oil. Therefore, the development of fine research on the internal configuration of sand bodies has become an important direction for reservoir research in recent years. The main contents of the sand body configuration research generally comprise a configuration interface, configuration elements, sand body overlapping relation and extension range.
A plurality of scholars conduct detailed research on the aspects of single river channel and point dam identification, lateral volume quantitative characterization, river channel evolution analysis and the like, and provide basis for the close well network to finely depict reservoir configurations of different levels by utilizing geological data such as field outcrop, modern sediment and the like, but due to the uncertainty of inter-well sand body depicting, the research is conducted on the fine depicting of river channel characteristics of a thin well zone (large well spacing).
The seismic sedimentology has great advantages in the aspect of depicting the transverse spreading characteristics and the longitudinal evolution of sedimentary bodies, can be used for finely determining the boundaries of sand bodies, and is favorable for determining the transverse spreading characteristics and the configuration characteristics of the sand bodies among wells of a dilute well zone. In addition, methods such as dynamic combination and numerical simulation are also applied to reservoir configuration research. The research on the sand body configuration under the condition of no seismic data control is lacked at present.
Disclosure of Invention
In order to solve the problems, the invention provides a sand body configuration quantitative analysis method, which mainly comprises the following steps:
s101: identifying a gyrus interface of a research area according to a preliminary sedimentary facies research result by combining a rock core observation method and a logging interpretation method;
s102: dividing the research horizon according to the cyclic interface of the research area to determine a sand body configuration interface;
s103: carrying out microphase analysis on the sand body configuration interface, and determining the single-well sand body configuration elements and the contact relation thereof;
s104: combining the single-well sand body configuration elements and the contact relationship thereof to obtain a longitudinal sand body splicing relationship;
s105: predicting to obtain a transverse sand body splicing relation according to the longitudinal sand body splicing relation;
s106: finishing a sand body plane distribution map according to the transverse sand body splicing relation and the longitudinal sand body splicing relation;
s107: and substituting the actual data sample into Monte Carlo simulation in combination with the sand body plane distribution diagram to realize fine depiction of the sand body form and scale.
Further, step S103 specifically includes: and carrying out microphase analysis on the single well and well connecting sections of the sand body configuration interface, carrying out configuration element division on the single well according to well lithology combination and electrical logging curve characteristics, and determining the single well sand body configuration elements and the contact relation thereof.
Further, step S104 specifically includes: and (4) statistically calculating the thickness, the layer thickness ratio and the overlapping relation probability of each configuration element of the single-well sand body to obtain the splicing relation of the longitudinal sand bodies.
Further, step S105 specifically includes: and according to the longitudinal sand body splicing relation, combining well spacing, and predicting and simulating the extension width and splicing probability of each sand body on a transverse section by using a Walter phase law, an equiprobable principle and a Monte Carlo simulation method to obtain the transverse sand body splicing relation.
Further, step S106 specifically includes simulating an extension range and a splicing probability of each sand body between the cross sections by using a monte carlo simulation method according to the horizontal sand body splicing relationship and the longitudinal sand body splicing relationship, and completing the sand body plane distribution diagram.
Further, in step S107, the specific manner of acquiring the actual data sample is as follows: and obtaining an actual data sample of the sand body configuration by using the existing deposition data, the similar oil reservoir data and the actual configuration data of the earth surface outcrop.
The existing deposition data comprising a length and a width of an observable statistical configuration element; the similar reservoir data including length, width and thickness of a conformational element that can be interpreted by both seismic and well logging; the surface outcrop data comprises the width and thickness of observable statistical configuration elements; and the sand body configuration actual data samples comprise an aspect ratio and an aspect ratio.
The technical scheme provided by the invention has the beneficial effects that: the problem of determining the sand body configuration parameters under the condition of less constraint can be solved.
Drawings
The invention will be further described with reference to the accompanying drawings and examples, in which:
fig. 1 is a flow chart of a sand body configuration quantitative analysis method in the embodiment of the invention.
Detailed Description
For a more clear understanding of the technical features, objects and effects of the present invention, embodiments of the present invention will now be described in detail with reference to the accompanying drawings.
The embodiment of the invention provides a sand body configuration quantitative analysis method.
Referring to fig. 1, fig. 1 is a flowchart of a sand body configuration quantitative analysis method in an embodiment of the present invention, which specifically includes the following steps:
s101: identifying a gyrus interface of a research area according to a preliminary sedimentary facies research result by combining a rock core observation method and a logging interpretation method;
s102: dividing the research horizon according to the cyclic interface of the research area to determine a sand body configuration interface;
s103: carrying out microphase analysis on the sand body configuration interface, and determining the single-well sand body configuration elements and the contact relation thereof;
s104: combining the single-well sand body configuration elements and the contact relationship thereof to obtain a longitudinal sand body splicing relationship;
s105: predicting to obtain a transverse sand body splicing relation according to the longitudinal sand body splicing relation;
s106: finishing a sand body plane distribution map according to the transverse sand body splicing relation and the longitudinal sand body splicing relation;
s107: and substituting the actual data sample into Monte Carlo simulation in combination with the sand body plane distribution diagram to realize fine depiction of the sand body form and scale.
Step S103 specifically includes: and carrying out microphase analysis on the single well and well connecting sections of the sand body configuration interface, carrying out configuration element division on the single well according to well lithology combination and electrical logging curve characteristics, and determining the single well sand body configuration elements and the contact relation thereof.
Step S104 specifically includes: and (4) statistically calculating the thickness, the layer thickness ratio and the overlapping relation probability of each configuration element of the single-well sand body to obtain the splicing relation of the longitudinal sand bodies.
Step S105 specifically includes: and according to the longitudinal sand body splicing relation, combining well spacing, and predicting and simulating the extension width and splicing probability of each sand body on a transverse section by using a Walter phase law, an equiprobable principle and a Monte Carlo simulation method to obtain the transverse sand body splicing relation.
And S106, simulating the extension range and the splicing probability of each sand body between the sections by using a Monte Carlo simulation method according to the transverse sand body splicing relation and the longitudinal sand body splicing relation, and completing the plane distribution diagram of the sand bodies.
In step S107, the specific manner of acquiring the actual data sample is as follows: and obtaining an actual data sample of the sand body configuration by using the existing deposition data, the similar oil reservoir data and the actual configuration data of the earth surface outcrop.
The existing deposition data comprising a length and a width of an observable statistical configuration element; the similar reservoir data including length, width and thickness of a conformational element that can be interpreted by both seismic and well logging; the surface outcrop data comprises the width and thickness of observable statistical configuration elements; and the sand body configuration actual data samples comprise an aspect ratio and an aspect ratio.
The invention has the beneficial effects that: the problem of determining the sand body configuration parameters under the condition of less constraint can be solved.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.

Claims (5)

1. A sand body configuration quantitative analysis method is characterized in that: the method comprises the following steps:
s101: identifying a gyrus interface of a research area according to a preliminary sedimentary facies research result by combining a rock core observation method and a logging interpretation method;
s102: dividing the research horizon according to the cyclic interface of the research area to determine a sand body configuration interface;
s103: and (3) carrying out microphase analysis on the sand body configuration interface to determine the single-well sand body configuration elements and the contact relation thereof, specifically comprising the following steps:
carrying out microphase analysis on the single well and well connecting sections of the sand body configuration interface, carrying out configuration element division on the single well according to well lithology combination and electrical logging curve characteristics, and determining the single well sand body configuration elements and the contact relation thereof;
s104: combining the single-well sand body configuration elements and the contact relationship thereof to obtain a longitudinal sand body splicing relationship;
s105: predicting to obtain a transverse sand body splicing relation according to the longitudinal sand body splicing relation;
s106: finishing a sand body plane distribution map according to the transverse sand body splicing relation and the longitudinal sand body splicing relation;
s107: substituting the actual data sample into Monte Carlo simulation in combination with the sand body plane distribution diagram to realize fine depiction of the sand body shape and scale;
and according to the longitudinal sand body splicing relation, combining well spacing, and predicting and simulating the extension width and splicing probability of each sand body on a transverse section by using a Walter phase law, an equiprobable principle and a Monte Carlo simulation method to obtain the transverse sand body splicing relation.
2. The sand body configuration quantitative analysis method as claimed in claim 1, characterized in that: step S104 specifically includes: and (4) statistically calculating the thickness, the layer thickness ratio and the overlapping relation probability of each configuration element of the single-well sand body to obtain the splicing relation of the longitudinal sand bodies.
3. The sand body configuration quantitative analysis method as claimed in claim 1, characterized in that: and S106, simulating the extension range and the splicing probability of each sand body between the sections by using a Monte Carlo simulation method according to the transverse sand body splicing relation and the longitudinal sand body splicing relation, and completing the plane distribution diagram of the sand bodies.
4. The sand body configuration quantitative analysis method as claimed in claim 1, characterized in that: in step S107, the specific manner of acquiring the actual data sample is as follows: and obtaining an actual data sample of the sand body configuration by using the existing deposition data, the similar oil reservoir data and the actual configuration data of the earth surface outcrop.
5. The sand body configuration quantitative analysis method as claimed in claim 4, characterized in that: the existing deposition data comprising a length and a width of an observable statistical configuration element; the similar reservoir data including length, width and thickness of a conformational element that can be interpreted by both seismic and well logging; the surface outcrop data comprises the width and thickness of observable statistical configuration elements; and the sand body configuration actual data samples comprise an aspect ratio and an aspect ratio.
CN201910793343.0A 2019-08-27 2019-08-27 Sand body configuration quantitative analysis method Expired - Fee Related CN110632665B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910793343.0A CN110632665B (en) 2019-08-27 2019-08-27 Sand body configuration quantitative analysis method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910793343.0A CN110632665B (en) 2019-08-27 2019-08-27 Sand body configuration quantitative analysis method

Publications (2)

Publication Number Publication Date
CN110632665A CN110632665A (en) 2019-12-31
CN110632665B true CN110632665B (en) 2021-06-25

Family

ID=68969077

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910793343.0A Expired - Fee Related CN110632665B (en) 2019-08-27 2019-08-27 Sand body configuration quantitative analysis method

Country Status (1)

Country Link
CN (1) CN110632665B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112034528A (en) * 2020-09-07 2020-12-04 陕西延长石油(集团)有限责任公司 Method for depicting single sand body configuration of clastic rock reservoir

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106437690A (en) * 2015-08-07 2017-02-22 中国石油天然气股份有限公司 Method and device for determining three-dimensional spatial distribution of trigonal facies of braided river
CN105373648B (en) * 2015-10-09 2018-10-12 西南石油大学 A kind of modeling method of sand-mud interbed type lithosomic body reservoir architecture
CN108957549A (en) * 2018-06-26 2018-12-07 中国石油天然气股份有限公司 Braided river sediment heterogeneous compact sandstone gas reservoir geological modeling method

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10816684B2 (en) * 2017-02-02 2020-10-27 Exxonmobil Upstream Research Company Method for rapid calibration of seismic interval attributes to rock properties for assessment uncertainty

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106437690A (en) * 2015-08-07 2017-02-22 中国石油天然气股份有限公司 Method and device for determining three-dimensional spatial distribution of trigonal facies of braided river
CN105373648B (en) * 2015-10-09 2018-10-12 西南石油大学 A kind of modeling method of sand-mud interbed type lithosomic body reservoir architecture
CN108957549A (en) * 2018-06-26 2018-12-07 中国石油天然气股份有限公司 Braided river sediment heterogeneous compact sandstone gas reservoir geological modeling method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
砂质辫状河岩相与构型特征-以山西大同盆地中侏罗统云冈组露头为例;陈彬滔,等;《石油与天然气地质》;20150228;第111-117页 *
苏里格气田盒8气层组厚层辫状河道砂体构型分析;卢海娇,等;《东北石油大学学报》;20140228;第54-63页 *
辫状河三角洲露头构型-以塔里木盆地库车坳陷三叠系黄山街组为例;朱卫红,等;《石油勘探与开发》;20160630;第482-488页 *

Also Published As

Publication number Publication date
CN110632665A (en) 2019-12-31

Similar Documents

Publication Publication Date Title
CN108957549B (en) Braided river sediment heterogeneous compact sandstone gas reservoir geological modeling method
CN101487390B (en) Archie mode method for confirming initial oil saturation of oil layer
CN110318745B (en) Particle size lithology logging evaluation method under deposition microphase constraint
CN103645516B (en) Based on the method for rock physics phased oil-gas reactivation determination oil-production capacity
CN102434152A (en) Method for calculating oil saturation of reservoir
CN103256046A (en) Unconventional oil and gas reservoir horizontal well section full-fracture-length fracturing parameter analog method and device
CN106126936B (en) A kind of integrated evaluating method of densification low permeability reservoir fracture effectivity
CN103869052A (en) Method for quantitatively expressing diagenetic logging facies of sandstone reservoir
CN109375283B (en) Analysis method for 3D permeability evolution history of sandstone reservoir
CN105629308A (en) Phase control heterogeneous mechanical parameter crustal stress prediction method
US20190203593A1 (en) Method and System for Modeling in a Subsurface Region
CN104251137A (en) Interlayer identification prediction method in fluvial facies thick oil layer
CN110286416B (en) Fast two-dimensional density inversion method based on physical property function
CN105134185A (en) Reservoir fluid property identification method
CN105257285A (en) Evaluation method of oil gas filling capacity from source rock layer to karstic reservoir
CN105257284A (en) Method and device for determining content of coagulated ash by utilizing element capture energy spectrum logging
CN105205296A (en) Method for calculating porosity of shale gas reservoir
CN110632665B (en) Sand body configuration quantitative analysis method
CN108227033B (en) Method and device for predicting porosity of sandstone reservoir and computer storage medium
CN105804732A (en) Oil-gas distribution prediction method based on relative resistivity stochastic simulation inversion
CN103698494B (en) Method and device for determining saturation degree of hydrocarbon in lithologic trap
CN100349014C (en) Method for processing varying density terrain correction by heavy prospecting data
CN103197348A (en) Method using internal samples at reservoirs to carry out weighting and compile logging crossplot
Liu et al. Geostatistical modeling for fine reservoir description of Wei2 block of Weicheng oilfield, Dongpu depression, China
CN111608649B (en) Method for predicting beneficial area of exogenous supply type shale gas exploration

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
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210625

CF01 Termination of patent right due to non-payment of annual fee